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: predicting survival with a semiparameteric model


From   "Feiveson, Alan H. (JSC-SK311)" <alan.h.feiveson@nasa.gov>
To   "statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu>
Subject   RE: st: predicting survival with a semiparameteric model
Date   Wed, 25 Aug 2010 08:06:16 -0500

Maarten - 

Thanks for your suggestion and the example. Since then, John Moran has suggested another possibility  - Yulia Marchenkp's -survci-, which uses an asymptotic variance for the covariate-adjusted survivor function estimate after Cox regression.

Al F.

-----Original Message-----
From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Maarten buis
Sent: Tuesday, August 24, 2010 2:03 PM
To: statalist@hsphsun2.harvard.edu
Subject: Re: st: predicting survival with a semiparameteric model

--- On Tue, 24/8/10, Feiveson, Alan H. (JSC-SK311)  wrote:
> Of course, I can always get a standard error of S(t) with a
> fully parametric model using -streg-, but if possible, I'd
> like to use a PH model without having to specify a
> distribution.

The key issue with specifying the distribution that -stcox- 
tries to solve is that you don't have to specify the baseline
hazard rate. This is a strength and a weakness, you make the
model more robust, but you limit what you can say about the
distribution, as you find out.

A solution is to choose a parametric model, and add a 
parametric form for the baseline, but keep it flexible:

A common method is the piecewise constant model, but you can
go further and model it with (restricted) cubic splines, like
this:

*------------- begin example --------------
sysuse cancer, clear
gen int id = _n
stset studytime, failure(died) id(id)
stsplit time, every(1)
mkspline sp = time, cubic knots(5 10 20 30)
streg i.drug age sp*, dist(exp) hr
*-------------- end example ---------------

I believe Paul Lambert and Patrick Royston has done and published 
stuff like this, see:


P. C. Lambert and P. Royston (2009) Further development of flexible 
parametric models for survival. The Stata Journal, 9(2):265--290.
http://www.stata-journal.com/article.html?article=st0165

Hope this helps,
Maarten

--------------------------
Maarten L. Buis
Institut fuer Soziologie
Universitaet Tuebingen
Wilhelmstrasse 36
72074 Tuebingen
Germany

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


      

*
*   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/

*
*   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/


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