[R] Survival analysis

David Winsemius dwinsemius at comcast.net
Wed Jan 29 01:30:37 CET 2014


On Jan 28, 2014, at 4:19 PM, Dennis Murphy wrote:

> Hi David:
> 
>>> Do you meant [sic] the (cumulative) hazard function that is 1 minus the survival function?
> 
> Um, 1 - S(t) = F(t), the cumulative *distribution function* of the
> random variable T. The cumulative hazard function is -ln S(t) and the
> hazard function is
>       h(t) = d(-ln S(t))/dt = f(t)/S(t)  [Klein and Moeschberger, p.27].
> 
> I know you know this, but we all need a reminder occasionally...

Actually I had completely forgotten the distinction, so I hope you don't mind my forwarding it to rhelp. I do appreciate being corrected when I make errors.

-- 
david

> 
> Dennis
> 
> On Tue, Jan 28, 2014 at 4:01 PM, David Winsemius <dwinsemius at comcast.net> wrote:
>> 
>> On Jan 28, 2014, at 7:58 AM, Endy BlackEndy wrote:
>> 
>>> Dear R users, how I can implement Breslow and Tarone-Ware tests for
>>> survival analysis?
>> 
>> Tests ... of what?
>> 
>> Perhaps you should learn to search:
>> 
>> install.packages('sos')
>> library(sos)
>> findFn("Tarone-Ware")
>> 
>> The 'comp' function in {survMisc} package promises several log-rank tests, one of them named Tarone-Ware, another named Gehan-Breslow. It would probably help to refer to the section of Klein J, Moeschberger M 2003 Survival Analysis, 2nd edition, that is cited on the help page. These test are described in terms of Fleming-Harrington weights and Klein and Moeschberger has several pages describing this family of tests.
>> 
>>> Also, is there any way I can estimate and plot the hazard function,along
>>> the lines of the survival function in Survival Package?
>> 
>> Do you meant the (cumulative) hazard function that is 1 minus the survival function? (And you should not capitalize the name of the 'survival' package.)
>> 
>>> (The package muhaz
>>> seems that does not do what it promises, except if I could not understand
>>> how it works.)
>>> 
>> --
>> 
>> David Winsemius
>> Alameda, CA, USA
>> 
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.

David Winsemius
Alameda, CA, USA




More information about the R-help mailing list