Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down at the end of May, and its replacement, statalist.org is already up and running.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: Cox model with unobserved heterogeneity


From   Steve Samuels <sjsamuels@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: Cox model with unobserved heterogeneity
Date   Fri, 28 Sep 2012 17:58:17 -0400

I agree with Maarten. I would only add that, as a matter of showing professional respect, Statalist practice is to use full names when posting (FAQ 2.1.3).

Steve

On Sep 28, 2012, at 5:27 AM, Maarten Buis wrote:

On Fri, Sep 28, 2012 at 11:12 AM, Nynke  wrote:
> I would like to estimate a mixed proportional hazard model with a flexibel specification of the hazard in Stata, comparable to a Cox model with unobserved heterogeneity. I can estimate a MPH with a parametric form of the hazard and unobserved heterogeneity (streg with frailty) or a Cox model without accounting for unobserved heterogeneity (stcox), but how do I sort of combine these two? And when I have estimated this model, is there a way to check the PH assumption?

I suspect that the model you are looking for is just not identified.
You can estimate a Cox model with group level unobserved
heterogeneity, with the -shared()- option in -stcox-. In that cases it
uses information from how similar observations are within groups. For
individual level unobserved heterogeneity the parametric models use
deviations from theoretical form of the baseline hazard function to
identify the unobserved heterogeneity component. Since Cox regression
has no theoretical form for the baseline hazard function, there is no
information that can be used to estimate the individual level
unobserved heterogeneity component.

Hope this helps,
Maarten

---------------------------------
Maarten L. Buis
WZB
Reichpietschufer 50
10785 Berlin
Germany

http://www.maartenbuis.nl
---------------------------------

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/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/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index