Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at

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

Re: st: RE: Survival analysis - individual survival functions

From   Steven Samuels <>
Subject   Re: st: RE: Survival analysis - individual survival functions
Date   Thu, 16 Dec 2010 16:27:18 -0500


You are very welcome, Paul, but I think that using the post-estimation facilities of -stpm2-, as outlined by Paul Lambert, might be better, as well as easier. -stpm2- has the advantage of being parametric, but doesn't require that you pick (and defend!) a parametric distribution.

In any case, with -streg- you'll need -nlcom- to produce estimates and CIs on a suitable scale and then transform back to the probability scale. With -stpm2- it would be simpler to bootstrap, given the ease of producing S(t)/S(t_0).


On Dec 16, 2010, at 3:31 PM, Seed, Paul wrote:

Thanks Steve.
I had been hoping to avoid doing all the algebra (or at
least get some pre-calculated values to check my algebra

It6 seems I'm out of luck.

Date: Wed, 15 Dec 2010 20:09:30 -0500
From: Steven Samuels <>
Subject: Re: st: Survival analysis - individual survival functions


 -predict- gives values at the observations after every regression
command. To get it for arbitrary values of the inputs, you have to do
it yourself.   S(t| t_0) = S(t)/S(t_0), and the formulas for S(.) are
on page 358 of the Stata Manual (V 11).  You can get the linear
predictor for any values of the X's with  -margins-  or -adjust-
( Stata<11).  Add the appropriate auxiliary parameters (returned
results e(theta), e(aux_p), e(gamma), e(sigma), e(kappa)); your t and
t_0; then plug  everything into the formulas.


Steven J. Samuels
18 Cantine's Island
Saugerties NY 12477
Voice: 845-246-0774
Fax:    206-202-4783

*   For searches and help try:

*   For searches and help try:

© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index