Hello Zelig people,
I am having difficulty simulating from ls.mixed. My code and errors
are below:
model <- zelig(y.sn ~
violator.democracy7*duration*violator.ratified +
violator.ratified*victim.ratified + Aerial + Armistice + CBW + POWs +
Civilians + Cultural + HighSeas + powerbal + tag(1 |violatorccode/
victimccode), model="ls.mixed", data=morrow)
x <- setx(model, fn=list(ordered = median, others = mode,
numeric=mean),
duration=log(46):log(1182), violator.democracy7=1,
violator.ratified=0)
Warning messages:
1: In FUN(X[[1L]], ...) :
There is more than one mode. The first level is selected.
2: In FUN(X[[1L]], ...) :
There is more than one mode. The first level is selected.
x1 <- setx(model, fn=list(ordered = median, others
= mode,
numeric=mean), duration=log(46):log(1182), violator.democracy7=1,
violator.ratified=1)
Warning messages:
1: In FUN(X[[1L]], ...) :
There is more than one mode. The first level is selected.
2: In FUN(X[[1L]], ...) :
There is more than one mode. The first level is selected.
s.out1 <- sim(model, x)
(Intercept)
violator.democracy7 duration violator.ratified
victim.ratified Aerial Armistice CBW POWs Civilians Cultural HighSeas
powerbal violator.democracy7:duration
violator.democracy7:violator.ratified duration:violator.ratified
violator.ratified:victim.ratified
violator.democracy7:duration:violator.ratified
Error in gammas[[names(rTerms[i])]] %*% t(as.matrix(rTerms[[i]])) :
requires numeric matrix/vector arguments
s.out2 <- sim(model, x1)
(Intercept)
violator.democracy7 duration violator.ratified
victim.ratified Aerial Armistice CBW POWs Civilians Cultural HighSeas
powerbal violator.democracy7:duration
violator.democracy7:violator.ratified duration:violator.ratified
violator.ratified:victim.ratified
violator.democracy7:duration:violator.ratified
Error in gammas[[names(rTerms[i])]] %*% t(as.matrix(rTerms[[i]])) :
requires numeric matrix/vector arguments
When I use a nearly identical specification with plain old least
squares, I do not get the error with sim or the warnings with setx.
It is clearly not a bug with ls.mixed, as I get no errors when I work
through the voteincome example from the manual.
Ideas?
Thanks,
Keith Schnakenberg
-
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/