Good morning
If I run
<<<
susan.lsmixed.out <- zelig(formula = unprot_vag_sex ~ married + age + TREATMENT.ARM*time + highest_grade + income + tag(1|id),
data = susanMI.out$imputations, model = "ls.mixed")
summary(susan.lsmixed.out)
>>>>
I get an error
Error in x$coef : $ operator is invalid for atomic vectors
Searching the archives, I see that others have had similar problems. Is there a workaround?
summary(susan.lsmixed.out[[1]])
works fine; should I then average across the five imputed data sets?
thanks!
Peter
Peter L. Flom, PhD
Statistical Consultant
Website: http://www DOT statisticalanalysisconsulting DOT com/
Writing; http://www.associatedcontent.com/user/582880/peter_flom.html
Twitter: @peterflom
-
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/
Using the example in:
http://cran.r-project.org/web/packages/Zelig/vignettes/probit.pdf
If I add robust=TRUE, I get the following error:
> z.out <- zelig(vote ~ race +
educate,mode="probit",data=turnout,robust=TRUE)
Error in glm.control(robust = TRUE) : unused argument(s) (robust = TRUE)
Thoughts?
Thank you,
bp
I'm trying to run a relogit model thru zelig with 1.4 million observations,
using weights, with about 15 covariates. Has been running for quite some
time. Is this just way to big to run, or is there some way to run this with
large data? Memory should not be a problem, as I'm running it on a grid
with about 50 gb of allocated RAM.
Any guidance would be much appreciated. thanks!
EXL
--
------------------------------
ERIC LIN
elin(a)iq.harvard.edu
------------------------------
After looking around at the documentation and some examples, it isn't clear
to me how to use plot.ci with model="sur".
Using the example in Zelig,
>data(grunfeld)
>formula <- list(mu1 = Ige ~ Fge + Cge, mu2 = Iw ~ Fw + Cw)
Estimating the model using sur:
>z.out <- zelig(formula = formula, model = "sur", data = grunfeld)
> summary(z.out)
How would I set Fw=10:91 (one of the independent variables above) in the
setx command? Just putting it as in a single equation model leads to
errors.
>f.range=10:90
> x.out <- setx(z.out,Fw=f.range)
Error in X[, xtmp] <- Xs[[i]][, xtmp] :
It seems like I need to call 'mu2' with Fw in it, but I am not sure how to
do that.
Thanks,
Greg
--
Gregory E. McAvoy
Associate Professor
Department of Political Science
University of North Carolina at Greensboro
PO Box 26170
Greensboro, NC 27402
