"random effects for country code" is not a very precise description, but my
guess is that in Stata he did:
xtlogit y x1 x2 , fe i(ccode)
this should mirror:
a<-zelig(y~x1+x2+tag(1|ccode), data=dat, model="logit.mixed")
this gives you a random country intercept (as opposed to a country fixed
effect). Best way to grasp RE models is read up on them in the Bates books
listed in the zelig references on "logit.mixed" or lme4 package.
hth,
j.
-----Original Message-----
From: gov2001-l-bounces at
lists.fas.harvard.edu [mailto:gov2001-l-
bounces at
lists.fas.harvard.edu] On Behalf Of Didi Kuo
Sent: Thursday, March 27, 2008 9:13 AM
To: gov2001-l at
lists.fas.harvard.edu
Subject: [gov2001-l] Random effects logit
Hey yall,
The author of the paper we're replicating did a logit model in stata
using xtlogit with, in his words, "random effects for country code." I
haven't been able to replicate his numbers in either stata or R...
In R, I've been doing:
a<-zelig(y~x1+x2...+tag(1|ccode), data=dat, model="logit.mixed")
But I'm not sure if finding a different intercept for each country
(which is what the above code does) is the same as what the author did.
Can anybody clarify what the author meant? (we don't have his exact
stata code - only what he provided above.)
thanks,
didi & shahrzad
_______________________________________________
gov2001-l mailing list
gov2001-l at
lists.fas.harvard.edu
http://lists.fas.harvard.edu/mailman/listinfo/gov2001-l