it might help to look at the Zelig manual, particularly the setx and sim
functions:
Hi,
Maybe I?m missing something in my understanding, but how do you evaluated
the difference in the expected probability in the data when checking for
Zelig? Do you mind describing how you specified the inputs y and X?
Thanks!
Clarence
*From:* gov2001-l-bounces at
lists.fas.harvard.edu [mailto:
gov2001-l-bounces at
lists.fas.harvard.edu] *On Behalf Of *Z. A. Townsend
*Sent:* Thursday, March 12, 2009 2:14 AM
*To:* gov2001-l at
lists.fas.harvard.edu
*Subject:* Re: [gov2001-l] Using Zelig to check our results in 2.3-2.5
for 2.3, for example, something like
data <- as.data.frame(cbind(y,X)) %which I defined earlier to be the vector
and matrix of interest
names(data) <- c("incident", "intercept", "temperature",
"pressure")
output <- zelig(incident ~ temperature +
pressure, data = data, model = "logit")
Should confirm your results.
On Thu, Mar 12, 2009 at 12:33 AM, Lee, Clarence <clee at hbs.edu> wrote:
Does anyone know the commands in Zelig to quickly check our results for
2.3-2.5?
Thanks,
Clarence
_______________________________________________
gov2001-l mailing list
gov2001-l at
lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l
_______________________________________________
gov2001-l mailing list
gov2001-l at
lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l
--
Patrick Lam
Department of Government and Institute for Quantitative Social Science,
Harvard University