小编Dar*_*n J的帖子

我有三个数据框如下.我希望根据Lon和Lat将它们组合成一个数据帧,并平均每个"单元"的3个值.我已经读过这个(计算多个数据帧的平均值)并试图利用聚合但无济于事....任何指针都赞赏.

> head(CSR.GRACE[,c(1:14)],10)
    Lon  Lat   January  February     March     April       May     June        July     August  September   October    November  December
1  28.5 -4.5 17.710425 13.855327 12.385712 13.558101 12.789865 6.913783  1.03770075 -5.3901741 -6.6351015 -7.661375 -3.09337944 6.0659410
2  29.5 -4.5 14.010154 10.257435  9.009641 10.275778  9.598241 5.166972  0.73570247 -4.2733162 -5.0861417 -5.850192 -2.93521806 4.1240150
3  30.5 -4.5 16.288443 10.467614  9.275714 10.904162 10.228808 5.364853  0.50089883 -4.7478741 -5.4320069 -6.316568 -3.80160315 3.8494745
4  31.5 -4.5 18.560677  9.932461  9.239592 11.037748 10.551886 5.281853  0.01181973 -4.9034324 -5.3504391 -6.438050 -4.41695714 3.3432301
5 …

我试图在Taylor图中绘制多个模型,并且在代码中稍微挣扎.我已设法生成图表(见图)但无法弄清楚如何减小轴,因为它们太大,标准化标记为1,2,3,4的轴并在相关性上添加刻度线 - 用勾号标记我希望每0.1一次有一个主要的刻度,每0.05到0.9就有一个小刻度,之后我试图在0.95处有一个主要的刻度,此时每0.01都有一个小刻度(如果这是有道理的).任何有关上述的帮助/建议都会有所帮助.我在'plotrix'包中使用了'taylor.diagram'(并阅读了它的指南 - 但我对R相对缺乏经验)并且附加了我的(有点基本的)代码到目前为止,但我的情节看起来相当混乱.谢谢

all.models <- as.data.frame(cbind(Sy.One, Sy.Two, Sy.Three, Sy.Four, Sy.Five, Sy.Six, Sy.Seven, Sy.Eight, Sy.Nine, Sy.Ten))

taylor.diagram(CSR, Sy.One, sd.arcs=T, ref.sd=T, pcex=1.5, main=NULL, pos.cor=F,
              xlab="Standard Deviation (cm)", ylab="Standard Deviation (cm)")

for (i in 1:dim(all.models)[2]) {
  model.wanted <- all.models[,i]
  taylor.diagram(CSR, model.wanted, sd.arcs=T, ref.sd=T, pcex=1.5, col=i, add=T, pos.cor=F)}

# Add legend
model.names <- c("Sy=1%","Sy=2%","Sy=3%","Sy=4%","Sy=5%","Sy=6%","Sy=7%","Sy=8%","Sy=9%","Sy=10%")
legend("top", model.names, pch=19, col=i, cex=1.0, bty="n", ncol=5)

我在下面为此做了一个快速示例数据框.基本上我想创建一个线图,其平均值为一条线,并且该线周围的阴影代表值的范围.我意识到我可能必须找到行最小值/最大值,但我不确定如何为行执行此操作,也不知道我将如何绘制此行

TEST <- data.frame(a=c(1,5,7,2), b=c(3,8,2,5), c=c(6,10,2,1))
TEST$mean <- rowMeans(TEST)

任何帮助表示感谢 - 谢谢

我有一个 0.25 * 0.25 度网格分辨率的数据集。我有另一个 1 * 1 度的分辨率，并希望使它们具有可比性（两者都更改为 1 * 1 分辨率）。该地区是 28.5 至 36.5 经度和 -4.5 至 4.5 纬度与 2003 年至 2012 年的月度数据。我已附加前几行作为例如，如果以下内容没有帮助，我已附加完整数据的 excel 文件！https://onedrive.live.com/redir?resid=9E74848E574367C6!2625&authkey=!AOcpftTyHVqU8CA&ithint=file%2ccsv我一直在研究“sp”（聚合）和“光栅”包，但无法弄清楚。一如既往地感谢任何帮助！

> head(cur.data, 40)
    Longitude Latitude Year Month DecimDate  Rainfall  RainAnom
1     28.625   -4.375 2012    12  2012.917 173.69343  73.74917
2     28.875   -4.375 2012    12  2012.917 186.91148  86.96723
3     29.125   -4.375 2012    12  2012.917 158.17921  58.23495
4     29.375   -4.375 2012    12  2012.917 194.71006  94.76581
5     29.625   -4.375 2012    12  2012.917 160.84774  60.90349
6 …