Dear List,
I am running a logit model on the following data, and I have truncated my
data set (a panel data set with country year as the unit of analysis) from
the original 418 country-year observations (rows) to 275 observations
(rows). This is the data.4 that I specify in model.4 below:
model.4 <- svydesign(ids=idcode, strata=NULL, weights=NULL, data=data.4)
However, when I attempt to run the logit model in the survey package, I
enter this
#Running Model 1 with three individual year lags in DV included
multilag.logit1 <-
svyglm(dintensity~lag1yr1+lag2yr1+lag3yr1+polity2l+gdpenl+oil+lpop+epop+p9010.prime+provcont+eprovcont+lsbextex+lsgrtex+lsbextex*
p9010.prime+epop*p9010.prime+epop*lsgrtex+provcont*eprovcont, design=model.4,
family=quasibinomial(link="logit"))
and get the following error:
Error in `$<-.data.frame`(`*tmp*`, ".survey.prob.weights", value = c(
0.00239234449760766, :
replacement has 418 rows, data has 275
I am confused because nothing in my data is 418 rows long any more. Any
clue as to what R is trying to tell me?
Many thanks,
Cat
On Sat, Mar 29, 2008 at 1:12 PM, <johns18 at fas.harvard.edu> wrote:
Does anyone know how to make R include missing data as
its own row in
contingency tables? thanks!
_______________________________________________
gov2001-l mailing list
gov2001-l at
lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l
--
Catherine Lena Kelly
Ph.D. Student
Harvard University Department of Government