我想使用 R 软件在 3 个月的时间尺度上计算网格级的 SPI 和 SPEI 。
图书馆(SPEI) 小蜘蛛<-read.csv("GriddedMonthlyRF.csv")
我尝试了以下两种选择
sapply(spigrid,函数(x){ spii<-spi(x, 3, na.rm = TRUE)}
和
for (x in names(spigrid)){ spi[x] <- spi(x, 3, na.rm = TRUE)
我期待网格级别的 SPI 在网格的每一列的 3 个时间尺度但失败了。