Hi -
I don't know if have an answer to this already, but this might help. If
you are trying to replicate something like random effects logit in
Stata, you can use the glmmML package (available at CRAN) to run the
model. You need to get the NAs out of the dataset, but other than that it
is pretty simple. Here is the code that I used to replicate a table from
our paper - it was right to at least two decimal places.
library(glmmML)
data7.1 <- na.omit(data.frame(data$warsa, data$sxp, data$sxp2, data$secm,
data$gy1, data$peace, data$geogia, data$lnpop, data$frac, data$etdo4590,
data$country))
table7.1 <- glmmML(data.warsa ~ data.sxp + data.sxp2 + data.secm +
data.gy1 + data.peace + data.geogia + data.lnpop + data.frac +
data.etdo4590, data=data7.1, cluster=as.factor(data7.1$ data.country) )
table7.1
Hope this helps.
Cheers,
Mike
On Fri, 19 Nov 2004, Yanbo Wang wrote:
Hello guys,
After spending several hours at a failing course, I guess it is time to ask for
help! If you happen to have any clue on how to ran random effect models on R, DO
please share it with me!
Best, Yanbo
ÒýÓÃ Olivia Lau <olau(a)fas.harvard.edu>du>:
Hi, guys.
FAS has just upgraded R to 2.0.0 on the servers. Those of you who are
using non-standard packages should run update.packages() to get the clean
version.
If you have problems come to OH tonight.
Yours,
Olivia.
_______________________________________________
gov2001-l mailing list
gov2001-l(a)lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l
--
Yanbo Wang
Ph.D Student
International Management and Strategy
Sloan School, MIT
_______________________________________________
gov2001-l mailing list
gov2001-l(a)lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l