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

Продължение на по-ранна публикация, интересувам се да науча как да получа обичайните мерки за относителното качество на статистически модел в 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 отговори, че вероятно може да се вземе проста средна стойност на мерките, получени от обобщенията на отделните набори от данни. „Вероятно“ ме прави малко подозрителен. Някакви мисли по този въпрос?

Много благодаря!!


r missing-data regression model-fitting r-zelig
person TiF    schedule 22.05.2013    source източник

Отговори (1)


arrow_upward
1
arrow_downward

Обсъдих въпроса с колеги и всички се съгласихме да приемем решението, предложено от 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
person TiF    schedule 23.05.2013