如何在zelig中獲得多重估算數據的模型擬合（AIC，F-統計）測量？

Question

跟進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）有沒有人知道如何獲得AIC，F-統計和度MI數據自由的值是多少？ （2）我找到了一個similar question，Koskuke Imai回答說，可以簡單地將單個數據集的摘要產生的度量值取平均值。 「可能」讓我有點懷疑。對此有何想法？

非常感謝！

Answer 1

我跟同事討論了這個問題，我們都同意去找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