Dear Zelig users,
I just started out with Zelig and have a *very* newbie-ish question: how
do I extract the standard errors of fixed effects coefficients from the
mer object returned by a mixed-effects zelig()? For instance:
set.seed(2009)
foo <-
data.frame(dv=runif(20)<.3,iv=rnorm(20),group=rep(c("A","B","C","D"),times=5))
z.out <- zelig(dv~iv+tag(1|group), data=foo, model="logit.mixed")
summary(z.out)
I would now like to extract the estimated fixed effects coefficients and
their standard errors (which are printed by summary() under "Fixed
effects") - but the @fixef slot only contains the estimates, not the
standard errors:
summary(z.out)@fixef
Did I overlook some slot? Or do I need to calculate the standard errors
myself (if so, any pointers on how? At this point, I would probably just
bootstrap them out of the original sample...)?
I'm using Zelig 3.4-5 in R 2.9.2 (yes, these are not the newest
versions, long story) and have read the FAQ and the "Statistical
Commands" part of the documentation and have searched both on this
list's and R-help's archive, to no success...
Best,
Stephan
-
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/