R中的“极值”包出错:无法绘制重现期

您得到的错误来自 aggregate 电话:

annmax <- aggregate(dat$Rain, by = list(dat$Year),max, na.rm=TRUE,na.action=NULL)
# Warning messages:
# 1: In FUN(X[[i]], ...) : no non-missing arguments to max; returning -Inf
# 2: In FUN(X[[i]], ...) : no non-missing arguments to max; returning -Inf
# [...]

在这个调用中,您正在计算 Rain 对于每个 Year

annmax
#    Group.1     x
# 1     1979   5.4
# 2     1980  27.1
# 3     1981  62.5
# [...]
# 33    2011  58.2
# 34    2012   5.7
# 35    2013  74.9
# 36    2014  -Inf
# 37    2015  -Inf
# 38    2016  -Inf
# 39    2017  -Inf

这些年,当 骨料 返回 -Inf

dat[dat$Year %in% 2014:2017,]
#      Year Month Day Rain
# 1086 2014     1   1   NA
# 1087 2014     1   2   NA
# 1088 2014     1   3   NA
# 1089 2014     1   4   NA
# 1090 2014     1   5   NA
# 1091 2014     1   6   NA
# 1092 2014     1   7   NA
# 1093 2014     1   8   NA
# 1094 2014     1   9   NA
# 1095 2014     1  10   NA
# 1096 2014     1  11   NA
# [...]

所以,你要决定如何处理那些失踪的岁月。这取决于分析。分析是否需要缺失年份的数据?

fevd )需要一些缺失年份的数据,你需要一些方法从其他年份中插值。

2) 如果分析不需要缺失的年份,就删除它们。这是最简单的解决方案:

annmax2 <- annmax[is.finite(annmax[,2]),]
gevfit1 <- fevd(annmax2[ ,2])