Как получить показатели соответствия модели (AIC, F-статистика) в zelig для многократных вмененных данных?

Question

В продолжение предыдущего поста мне интересно узнать, как получить обычные показатели относительного качества статистической модели в zelig для регрессии, используя многократные вмененные данные (созданные с помощью Amelia).

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) Я нашел похожий вопрос, на который Koskuke Imai ответил, что, вероятно, можно взять простое среднее из мер, полученных из сводок отдельных наборов данных. "Вероятно" делает меня немного подозрительным. Есть мысли по этому поводу?

Большое спасибо!!

Answer 1

Я обсудил этот вопрос с коллегами, и мы все согласились пойти на решение, предложенное Коскюке.

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