通过AIC选择型号

Question

我已经使用R中的fitdistr将正态分布拟合到索赔金额数据中。如何拟合多元正态分布（二维正态分布）？ 我想使用R中的AIC选择最适合我的样本数据的一个。如何进行？ 我努力了

IC<-Mclust(data,G=1) 
IC1<-Mclust(data,G=2)

BIC越小越好。 但是如何根据此Mclust结果计算AIC。

Answer 1

赤池信息准则（AIC）定义为2 * k-2 * ln（L），其中k是模型中参数的数量，L是模型的似然性（通过拟合最佳参数来最大化数据） 。

Mclust()将这两个都返回给您，因此您可以从返回的Mclust对象中的字段计算AIC。

require(mclust)
data(iris)
IC <- Mclust(data=iris, G=1)
IC1 <- Mclust(data=iris, G=2)
aic <- 2*IC$df - 2*IC$loglik
aic1 <- 2*IC1$df - 2*IC1$loglik

请执行?Mclust以获取更多信息。