Hello,
I am having problems with the "ls.mixed" model in zelig() and the weights
argument. This is similar to the problem documented here:
https://lists.gking.harvard.edu/pipermail/zelig/2011-July/000183.html. However, while it
seems like the issue in that thread was missing values on the dependent variable, I am
having the issue with no missing data. I can estimate the model with weights using lmer(),
but not using zelig(). I have an example that reproduces this error using fake data posted
here:
http://pastebin.com/29KbmRZD and listed below.
##
library(Zelig)
set.seed(1000)
x <- runif(1000, -1, 1)
y <- x + rnorm(1000)
wt <- rexp(1000)
id <- rep(1:20, times = 50)
d <- as.data.frame(cbind(y, x, wt, id))
m1 <- zelig(y ~ x + tag(1|id), model = "ls.mixed", weights = "wt",
data = d)
m2 <- lmer(y ~ x + (1|id), weights = wt, data = d)
##
Thanks for any insight with this problem.
Jeff Harden
---
Jeffrey J. Harden
Ph.D Candidate, Department of Political Science
Program Coordinator, American Politics Research Group (APRG)
University of North Carolina at Chapel Hill
jjharden@unc.edu<mailto:jjharden@unc.edu>
http://jjharden.web.unc.edu/