hi,
I need to extract qi values from s.out to bring the results to another app.
When I do this s.out$qi, I get all null values. Any idea what might be
wrong?
Regards,
Sean
Zelig package is what I have been looking for a long time. I slight problem
is that I dont understand how sim() is run. I read through the docs and I
think it is written for graduate students. We can run the numbers but at
the end we need to be able to explain the results to the audience.
Can somebody explain with an example, how does sim() work?
"Zelig simulates parameters from classical *maximum likelihood* models
using asymptotic normal approximation to the log-likelihood. This is the
same assumption as used for frequentist hypothesis testing (which is of
course equivalent to the asymptotic approximation of a Bayesian posterior
with improper uniform priors). See King, Tomz, and Wittenberg
(2000)<http://gking.harvard.edu/files/abs/making-abs.shtml>.
For *Bayesian models*, Zelig simulates quantities of interest from the
posterior density, whenever possible. For *robust Bayesian models*,
simulations are drawn from the identified class of Bayesian posteriors."
Regards,
Mike
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