我使用 R 中的以下编码绘制了回报水平图。
data(ftcanmax)
fit <- gev.fit(ftcanmax[,"Prec"])
return.level(fit)
如何显示对应回报期的回报水平值?比如return period为10,return value的值是多少?
我使用 R 中的以下编码绘制了回报水平图。
data(ftcanmax)
fit <- gev.fit(ftcanmax[,"Prec"])
return.level(fit)
如何显示对应回报期的回报水平值?比如return period为10,return value的值是多少?
该代码是 extRemes ?return.level 包中的示例:
library(extRemes)
#...code
rets<-return.level(fit)
rets$return.level[rets$return.period==10]
#[1] 281.459
仅当 'rperiods' 参数return.level
包含感兴趣的期间......在默认情况下它为 10 时才可能这样做。因此,对于非默认期间,这应该足够了:
> return.level(fit ,rperiod=12)$return.level
[1] 296.796