Can you try to roll back to 3.5?
Best,
Kosuke
Department of Politics
Princeton University
http://imai.princeton.edu
On Sep 5, 2013, at 10:11 AM, Eva Bernauer <eva.bernauer(a)uni-mannheim.de<mailto:eva.bernauer@uni-mannheim.de>> wrote:
Dear Kosuke Imai,
may I ask you a question on Zelig? (My mail to Matt Owen somehow could not be delivered).
I found that the logit and the relogit model in zelig do not appropriately run
with robust standard errors. In fact, the results do not change when robust=TRUE. As you might know, this issue has already been discussed in the Zelig forum.
Is there any solution to the problem?
I would be very glad, if you could tell me, as such a fix would faciliate everything so much!
Best,
Eva
Hi Zelig list,
I've been having issues using plot.ci for a Poisson GEE. I originally thought it was due to a unique issue with my dataset, but I've tried doing the same with the poisson GEE example from the Zelig documentation, and get the same error:
Warning messages:
1: In simulation.matrix.sim(obj[[label]], which = which, exact.match = FALSE) :
The "which" parameter does not exist. Valid titles are:
"Expected Values (for x): E(Y|X)", "Expected Values (for x1): E(Y|X1)", "First Differences: E(Y|X1) - E(Y|X)"
The code that I'm using is:
data(sanction)
sanction$cluster <- c(rep(c(1:15),5),rep(c(16),3))
sorted.sanction <- sanction[order(sanction$cluster),]
z.out <- zelig(num ~ target + coop, model = "poisson.gee", id = "cluster", data = sorted.sanction, robust=TRUE, corstr="exchangeable")
summary(z.out)
coop.range<-1:4
x.low <- setx(z.out, target=2, coop=coop.range)
x.high <- setx(z.out, target=2, coop=coop.range)
s.out <- sim(z.out, x = x.low, x1 = x.high)
plot.ci(s.out, var="coop", ylim=c(1, 150))
Any insights into what might be going wrong?
Cheers,
Laura
Laura R. Crothers
University of Texas at Austin
Section of Integrative Biology
1 University Station, C0990
Austin, TX 78712