r - 分布拟合中 AIC 的 R 代码

Question

我已经使用 in R 将 Weibull、对数正态和 Gamma 分布拟合到我的索赔金额数据fitdistr中。我想使用 in R 选择最适合我的样本数据的分布。我该AIC如何进行？

Answer 1

您可以使用以下AIC功能：

set.seed(1)
x <- rlnorm(100) # random values from a log-normal distribution

# fit distributions
library(MASS)
weibull <- fitdistr(x, "weibull")
lognormal <- fitdistr(x, "lognormal")
gamma <- fitdistr(x, "gamma")

# compare AICs
AIC(weibull)
# [1] 300.9519
AIC(lognormal)
# [1] 287.0875
AIC(gamma)
# [1] 297.1818

毫不奇怪，对数正态拟合的AIC最低。这是最合适的。