分配拟合中的AIC的R代码

Question

我已经使用R中的fitdistr将Weibull，lognormal和Gamma分布拟合到我的索赔金额数据。我想使用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。这是最合适的。