小编CFo*_*ate的帖子

我正在使用 SPEI 包及其 32 年的每月样本数据。我想修改 x 轴标签以反映年份而不是数字。但是，硬编码绘图功能不允许我这样做。我厌倦了提取SPEI$fitted数据并尝试使用 ggplot 复制相同的图，但没有成功。这是示例代码

install.packages("SPEI")
library(SPEI)
data("wichita")
wichita$PET=hargreaves(Tmin=wichita$TMIN, Tmax = wichita$TMAX, lat = 37.64)
wichita$BAL=wichita$PRCP - wichita$PET
SPEI_12=spei(wichita[,"BAL"],12)
plot.spei(SPEI_12, main = 12-Month SPEI)

任何帮助，将不胜感激。我想生成一个如图所示的图表。

我想将数据分组到多个季节，这样我的季节就是冬季：12 月 - 2 月；春季：三月至五月；夏季：6 月 - 8 月，秋季：9 月 - 11 月。然后，我想绘制冬季和春季季节性数据的箱线图，比较 A 到 B，然后比较 A 到 C。这是迄今为止我的费力代码。我希望有一种有效的数据分组和绘图方法。

library(tidyverse)
library(reshape2)

Dates30s = data.frame(seq(as.Date("2011-01-01"), to= as.Date("2040-12-31"),by="day"))
colnames(Dates30s) = "date"
FakeData = data.frame(A = runif(10958, min = 0.5, max = 1.5), B = runif(10958, min = 1.6, max = 2), C = runif(10958, min = 0.8, max = 1.8))
myData = data.frame(Dates30s, FakeData)
myData = separate(myData, date, sep = "-", into = c("Year", "Month", "Day"))
myData$Year = as.numeric(myData$Year)
myData$Month = as.numeric(myData$Month) …

我想生成满足以下条件的 500 种不同的 a、b 和 c 组合

1. a+ b+ c = 1 和
2. a < b < c

这是生成随机数的基本示例，但是，我需要根据上述条件生成它。

Coeff = data.frame(a=runif(500, min = 0, max = 1),
b=runif(500, min = 0, max = 1),
c=runif(500, min = 0, max = 1))