probability - 如何从 extRemes 输出中提取单个参数

Question

我已经为我的数据拟合了不同的概率分布。为了比较拟合优度 (GOF) 并确定最佳分布，我使用 AIC 作为标准。虽然我可以轻松地为“fitdist”对象提取 GOF，例如

fit.gamma <- fitdist(x, "gamma")
fit.norm <- fitdist(x, "norm")

aic.gamma <- fit.gamma$aic
shape.gamma <- fit.gamma$estimate[1]

等等 ...

但对于

fit.gev = fevd(x ,type="GEV") # 这使用 fevd（来自 'extRemes'）

虽然我们可以使用“summary(fit.gev)”显示 fevd 的对象，但似乎没有办法提取这些单独的参数（AIC 或比例等）以供后续计算使用。

任何帮助将不胜感激。谢谢你。

标签： probabilityprobability-densityprobability-distribution

Answer 1

答案在 Eric Gilleland 的 extRemes 手册第 55 页
https://cran.r-project.org/web/packages/extRemes/extRemes.pdf

1. 将 fevd 的摘要分配给一个对象，例如

look1 <- 摘要（fit.gev.12，静默 = TRUE）

2. 将对象提取为：

AIC <- c(look1$AIC)
BIC <- c(look1$BIC)

和参数

位置<-look1$par[1]
比例<-look1$par[2]
形状<-look1$par[3]

完毕！