使用SPEI的标准化降水指数

Question

我试图使用R中提供的Wichita数据开发标准化降水指数(SPI)，为此我使用了 "SPEI "包。我可以生成SPI数据，但质量不好，不适合发布。我在这里发现了一个非常类似的问题。如何在R中格式化SPEI包的硬编码绘图函数的X轴？.

当我试着用上面链接中提供的代码时，我可以为SPEI生成质量很好的图，但不能为SPI生成。我在下面一行出现了问题。dplyr::mutate(sign = ifelse(ET0_har >= 0, "pos", "neg")) 代码的下面一行出现了问题。我的具体问题是我应该用什么代替 "ET0_har "来生成SPI？

如果能像上面的链接那样生成SPI数据，我将感激您的意见。下面是我的代码。

install.packages("SPEI")
library(SPEI)
data("wichita")
SPI1<-spi(wichita$PRCP,scale = 1, distribution = 'Gamma')
plot(SPI1)

Answer 1

把这一行改成。

dplyr::mutate(sign = ifelse(`Series 1` >= 0, "pos", "neg")) 

这样就可以了。