[R] Question about 'logit' and 'mlogit' in Zelig
Mathew, Abraham T
amathew at ku.edu
Mon Mar 29 17:10:07 CEST 2010
I'm running a multinomial logit in R using the Zelig packages. According to str(trade962a), my dependent variable is a factor with three levels. When I run the multinomial logit I get an error message. However, when I run 'model=logit' it works fine. any ideas on whats wrong?
## MULTINOMIAL LOGIT
anes96two <- zelig(trade962a ~ age962 + education962 + personal962 + economy962 + partisan962 + employment962 + union962 + home962 + market962 + race962 + income962, model="logit", data=data96)
summary(anes96two)
#Error in attr(tt, "depFactors")$depFactorVar :
# $ operator is invalid for atomic vectors
## LOGIT
Call:
zelig(formula = trade962a ~ age962 + education962 + personal962 +
economy962 + partisan962 + employment962 + union962 + home962 +
market962 + race962 + income962, model = "logit", data = data96)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.021 -1.179 0.764 1.032 1.648
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.697675 0.600991 -1.161 0.2457
age962 0.003235 0.004126 0.784 0.4330
education962 -0.065198 0.038002 -1.716 0.0862 .
personal962 0.006827 0.072421 0.094 0.9249
economy962 -0.200535 0.084554 -2.372 0.0177 *
partisan962 0.092361 0.079005 1.169 0.2424
employment962 -0.009346 0.044106 -0.212 0.8322
union962 -0.016293 0.149887 -0.109 0.9134
home962 -0.150221 0.133685 -1.124 0.2611
market962 0.292320 0.128636 2.272 0.0231 *
race962 0.205828 0.094890 2.169 0.0301 *
income962 0.263363 0.048275 5.455 4.89e-08 ***
---
Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 1841.2 on 1348 degrees of freedom
Residual deviance: 1746.3 on 1337 degrees of freedom
(365 observations deleted due to missingness)
AIC: 1770.3
Number of Fisher Scoring iterations: 4
Thanks
Abraham
More information about the R-help
mailing list