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/
