Question

我正在尝试使用L-Moment方法拟合GEV df用于年度最大值，但是当使用不同的包时，拟合GEV分布的参数（位置，比例和形状）不相似。

我使用了不同的软件包而奇怪的软件包是lmomRFA，是否有人可以帮我识别问题？

这里是我正在做的样本：

install.packages("extRemes")
install.packages("lmom")
install.packages("lmomRFA")
install.packages("nsRFA")
library (nsRFA)
library (extRemes)
library (lmomRFA)
library (lmom)

amax_standard <- c (0.6274510, 0.7545455, 0.6521739, 1.5102041, 2.0937500, 1.0000000, 0.7094017, 1.0315789, 1.3207547)
amax_standard_matrix <- as.matrix(amax_standard)## The AMAX as a matrix

amax_standard_vector <-  as.vector(amax_standard_matrix) ## The AMAX as a vector
amax_standard_vector_order <- amax_standard_vector[order (amax_standard_vector)] ## The AMAX as an ordered vector

fit1 <- fevd(amax_standard_vector_order,type =("GEV"),method=("Lmoments")) #fittig EVD using the "extRemes" package

lmom <- samlmu(amax_standard_vector_order)  # calculating the L-Moment ratios for the region. This function deals with all the observation as one vector
fit2 <- pelgev(lmom) #fittig EVD using the "lmom" package

regdata <- regsamlmu(amax_standard_matrix, nmom = 4, sort.data = TRUE, lcv = FALSE)
colnames (regdata) [4]<- c("t") # to make the column name acceptable for the "regfit" function
fit3 <- regfit(regdata, "gev") #fittig EVD using the "lmomRFA" package

extRemes和lmom包中的前两个函数给出了GEV df估计参数的类似结果，而来自lmomRFA包的第三个函数提供了不同的答案。

此外，我尝试使用不同的包，nsRFA只是为了估计参数并使用函数par.GEV而我提供的函数具有估计的l-力矩比率，它们是相同的所有包都可以在变量lmom和regdata中找到，以估算参数，结果与extRemes和lmom包相似，那么问题是什么使用lmomRFA

中的函数

check <- par.GEV (1.0777622, 0.2760841, 0.3553799)

适用于不同包装的GEV df参数/不同结果（lmomRFA，lmom，extRemes，nsRFA）

0 个答案: