Hi,
I have a problem with using streg for weighted data. Unlike with
unweighted data, streg with unshared gamma frailty for weighted data
does not report a likelihood ratio test for unobserved
heterogeneity. I used the code [pw=sampleweight] in stset to
incorporate sample weights. I tried to do the likelihood ratio test
manually, by running models with and without the frailty
parameter. But this seemed not to work: The loglikelihoods from streg
are substantially different between when frailty(gammma) is included
(ll=-464), and when frailty is not included (ll=448). With only 1
degree of freedom, the difference in ll seems to be too big,
considering that the theta estimate is .45 with a standard error of
.48.
Do you have any idea about what is going on? Could you suggest some
appropriate test for heterogeneity?
Thank you,
Yan
--
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/