Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: Re: Interpreting Theta with gompertzhet


From   rgutierrez@stata.com (Roberto G. Gutierrez, StataCorp.)
To   statalist@hsphsun2.harvard.edu
Subject   st: Re: Interpreting Theta with gompertzhet
Date   Mon, 07 Oct 2002 11:07:05 -0500

A. A. Siyam <A.A.Siyam@lse.ac.uk> asks:

> A colleague of mine is asking about how to best interpret the "theta" value
> following the estimation of a "gompertzhet" model.

> As an example, we have the following model:

> xi:gompertzhet time i.sex i.tenureg i.newcohort i.marstat i.nmm i.socg
> i.educ i.bathg i.carg i.urind i.pob, dead(outcome) frailty(gamma) robust hr

[output omitted]

> We found a great deal of help in understanding the gompertz model from the
> book "Introduction to Survival Analysis Using Stata" by Cleves,Gould and
> Gutierrez. But, there seems to be very little explanation of how to use and
> interpret the "theta" coeffi cient estimated by the model.

First of all, the equivalent way to fit your model using the more modern 
-streg- is

    . stset time, failure(outcome)
    . xi: streg i.sex i.tenureg i.newcohort i.marstat i.nmm i.socg /*
    	  */ i.educ i.bathg i.carg i.urind i.pob, dist(gompertz) /*
	  */ frailty(gamma) robust

That being said, you'll find the discussion of the frailty variance, theta,
estimated by -streg, frailty()- in Section 15.3 of the aforementioned text
applicable.

Namely, theta is the estimated frailty variance.  If significantly different
from zero (in the case of robust no formal likelihood-ratio test is available,
so merely compare the parameter estimate to its standard error), this
indicates significant heterogeneity.  This heterogeneity can be thought of as
arising from some unobserved observation-specific effect, or frailty.  Theta
measures the variability of these frailties.  If close to zero, then the
frailties are virtually non-existent, and you can get by with the standard
Gompertz model.

--Bobby
rgutierrez@stata.com
*
*   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/



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