When running summary() on a zelig() result computed from a list of imputed data sets, I get the individual results, not the combined results. This used to be different.
Steps to reproduce:
## Fake "imputed" data
library(MASS)
dat <- mvrnorm(100, c(0, 0), matrix(c(1, .4, .4, 1), nrow = 2))
dat <- data.frame(dat)
names(dat) <- c("x", "y")
naidx <- sample(c(T, F), 100, replace = T, prob = c(.2, .8))
imp.dat <- lapply(1:5, function(i) {
dat[naidx,2] <- dat[naidx,2] + rnorm(sum(naidx), 0, .2)
dat
})
## Run Zelig
library(Zelig)
z.out <- zelig(x ~ y, model = "normal", data = do.call(mi, imp.dat))
summary(z.out)
whatever
commented
11 years ago
When running summary() on a zelig() result computed from a list of imputed data sets, I get the individual results, not the combined results. This used to be different.
Steps to reproduce:
Is this intended?