Hi everyone
I am running a mixed effects logit on a multiple imputed dataset, as follows:
G1m <- zelig(fm.mixed, data=a.out$imputations, model="logit.mixed")
this runs, but summary(G1m) gives:
summary(G1m)
Error in x$coef : $ operator is invalid for atomic vectors
Poking around, I discover that the model "object[i]" in the method
"summary.MI" is not responding nicely to summary() when that is called
to give res[i]:
summary(object[i])
Length Class Mode
[1,] 1 mer S4
class(object[i])
[1] "list"
by contrast, running the same thing on a normal dataset gives a "mer"
object which responds nicely to summary():
G1mtmp <- zelig(fm.mixed, data=ds,
model="logit.mixed")
summary(G1mtmp)
Generalized linear mixed model fit by the Laplace approximation
Formula: fm.mixed
Data: ds
AIC BIC logLik deviance
227 275 -103 205
[et cetera..........]
Is this just my setup, or can others reproduce this bug?
cheers,
David Hugh-Jones
Post-doctoral Researcher
Max Planck Institute of Economics, Jena
http://davidhughjones.googlepages.com
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)lists.gking.harvard.edu
[un]subscribe Options:
http://lists.gking.harvard.edu/?info=zelig
Zelig program information:
http://gking.harvard.edu/zelig/