R коды для AIC в распределительной арматуре

Question

Я установил распределения Weibull, lognormal и Gamma для моих данных о сумме претензии, используя fitdistr в R. Я хочу выбрать, какой из них лучше всего подходит для моих выборочных данных, используя AIC в R. Как мне поступить?

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. Это как нельзя лучше.