<>
The equivalance should be obvious from this:
***
webuse bangladesh, clear
// Random-intercept model,
// analogous to xtlogit
xtmelogit c_use urban age /*
*/child* || district:, nolog
xtlogit c_use urban /*
*/age child*, i(district)/*
*/ nolog
***
HTH
Martin
-----Original Message-----
From: [email protected]
[mailto:[email protected]] On Behalf Of Thomas Klausch
Sent: Montag, 10. August 2009 18:56
To: [email protected]
Subject: st: Reproducing xtlogit with xtmelogit
Hello stata list member,
As I understand it should be possible to reproduce a random effects
logistic regression model that was fitted with xtlogit by using the
multilevel syntax of xtmelogit. But for my case it just does not seem
to work.
I have longitudinal data, i.e. respondents nested in days. First step
is that I want to fit an empty random intercept model, that is:
xtset respondents days
xtlogit dummy, re
Which gives me a fit using Gauss-Hermite Quadrature in reasonable
time, say 3 minutes for NxT=60,000.
Then giving
xtmelogit dummy || respondents:
Does give no fit. I was not sure what would be the correct syntax, maybe it
is
xtmelogit dummy || respondents: ||days:
This takes forever to fit (I am not sure if the algorithm is ever
going to stop, It's still running after 1h waiting).
a. Am I specifying xtmelogit correctly?
b. Does xtmelogit use the same optimization algorithm as xtlogit?
c. What should be changed to yield the same results as with xtlogit?
Many thanks in advance.
Kind regards
Thomas
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
