分佈擬閤中的AIC的R代碼

Question

我在R中使用fitdistr對我的索賠量數據擬合了Weibull，對數正態和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。這是最合適的。