Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: st: equivalence of log-logistic survival estimation with gllamm
From
"JVerkuilen (Gmail)" <[email protected]>
To
[email protected]
Subject
Re: st: equivalence of log-logistic survival estimation with gllamm
Date
Tue, 26 Mar 2013 15:42:30 -0400
On Tue, Mar 26, 2013 at 2:43 PM, Nick Cox <[email protected]> wrote:
>
> The loglog and cloglog link functions have no application to survival
> times whatsoever. They are relevant _only_ to mean responses bounded
> by 0 and 1.
I'm with Nick. It's pretty clear there's some confusion going on.
However, there are discrete time proportional hazards survival models
that involve the cloglog link, and maybe that's what the original
poster wanted? I just checked Multilevel and Logitudinal Modeling
Using Stata, Volume II: Categorical Responses, Counts, and Survival,
Third Edition, S. Rabe-Hesketh and A. Skrondal, 2012, Stata Press.
They give an example using -xtcloglog- on p. 783, and discuss how this
could be fit using -gllamm-.
For a continuous time parametric survival model, I'm guessing that
some kind of censored normal model on the log-transformed time would
be necessary. That would be the lognormal model, not the log-logistic.
Out of my area, but I wonder if it would be possible to trick -gllamm-
to use the same basic "Poisson" likelihood discussed here, but with
censoring?
http://blog.stata.com/2011/08/22/use-poisson-rather-than-regress-tell-a-friend/
*
* 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/