Did you try the same model with a small number of covariates? Also, try
it without the cluster variable to see if you can identify the exact
problem.
Kosuke
--
Department of Politics
Princeton University
http://imai.princeton.edu
On Sun, 23 May 2010, Prashant wrote:
Hi,
I apologize in advance if this is a basic question. I am using zelig and
logit.survey and trying to simulate simple quantities of interest using setx
and sim. After a standard sim command, I get the error
Argument eta must be a nonempty numeric vector
I have checked my syntax against standard examples but am still not sure
what the problem could be. Any help would be most appreciated.
An example of some of the code is:
impa.model <- zelig(impa ~ age + age_sqr + hh_income + urban + female, model
= "logit.survey", weights = ~weight, ids = ~county, data = my.data)
impa.urb20pf.age <- setx(impa.model, urban = 1, female = 1, age = 0:100,
hh_income = quantile(hh_income, .2))
impa.urb50pf.age <- setx(impa.model, urban = 1, female = 1, age = 0:100,
hh_income = quantile(hh_income, .5))
#(everything runs fine up to here)
urbf.age.out <- sim(impa.model, x = impa.urb20pf.age, x1 = impa.urb50pf.age)
Additional background: I have a large dataset (2 million plus) observations
but with no missing data. All of the variables in the logit equation are
numeric (binary or continuous), while the cluster identifier is a
string/categorical.
Thank you so much in advance,
Prashant
-
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/