跟进an earlier post ,我有兴趣学习如何使用多重插补数据(使用 Amelia 创建)来获取 zelig 中统计模型相对质量的常用度量,以进行回归。
require(Zelig)
require(Amelia)
data(freetrade)
#Imputation of missing data
a.out <- amelia(freetrade, m=5, ts="year", cs="country")
# Regression model
z.out <- zelig(polity~tariff+gdp.pc, model="ls", data=a.out$imputations)
summary(z.out)
Model: ls
Number of multiply imputed data sets: 5
Combined results:
Call:
lm(formula = formula, weights = weights, model = F, data = data)
Coefficients:
Value Std. Error t-stat p-value
(Intercept) 1.6740501340 1.0270535468 1.6299541 0.10342186
tariff 0.0196015092 0.0233789523 0.8384255 0.40234214
gdp.pc 0.0003296261 0.0001844909 1.7866798 0.07409327
For combined results from datasets i to j, use summary(x, subset = i:j).
For separate results, use print(summary(x), subset = i:j).
问题
(1) 有谁知道如何获取 MI 数据的 AIC、F-statistics 和 自由度 值?
(2) 我找到了 similar question Koskuke Imai 回答说,人们可以对各个数据集摘要产生的测量值进行简单平均。这个“可能”让我有点怀疑。对此有什么想法吗?
非常感谢!!
最佳答案
我与同事讨论了这个问题,我们都同意采用 Koskuke 建议的解决方案。
x1 <- summary(z.out[[1]])
x2 <- summary(z.out[[2]])
x3 <- summary(z.out[[3]])
x4 <- summary(z.out[[4]])
x5 <- summary(z.out[[5]])
#F-statistics
(x1$fstatistic[1]+x1$fstatistic[1]+x3$fstatistic[1]+x4$fstatistic[1]+x5$fstatistic[1])/5
#AIC
(AIC(z.out[[1]])+AIC(z.out[[2]])+AIC(z.out[[3]])+AIC(z.out[[4]])+AIC(z.out[[5]]))/5
关于r - 如何在 zelig 中获取多重插补数据的模型拟合度(AIC、F 统计量)度量?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/16694587/