[R] cumulative incidence plot vs survival plot

David Winsemius dwinsemius at comcast.net
Tue Jun 28 00:38:36 CEST 2011


On Jun 27, 2011, at 5:04 PM, array chip wrote:

> Hi David,
>
> Thanks for responding, and plain text ...(didn't realized I was in  
> rich text).
>
> The endpoint is disease recurrence, I was producing a regular KM  
> plot of
> recurrence-free probability. Then someone recommend using cumulative  
> incidence
> is preferred because death was censored in the dataset. I did a  
> little googling,
> I found CI was used often in the context of competing risk. I am  
> totally new to
> competing risk and trying to understand what competing risk means  
> and why CI is
> preferred than KM survival in this context. If you could share your  
> thoughts
> helping me to understand, greatly appreciated.
>
> Searched archive, found people talking about cmprsk package for  
> estimating and
> plotting CI. would that be the same as the code you suggested:  
> plot(time,
> cumsum(dead))
>

I do not think it would. I think this morning's post from Dr Geskus is  
on point here. You should read his 2011 'Biometrics' paper.

-- 
David.

> Thanks very much!
>
> John
>
>
>
>
> ________________________________
> From: David Winsemius <dwinsemius at comcast.net>
> To: array chip <arrayprofile at yahoo.com>
> Cc: r-help at r-project.org
> Sent: Mon, June 27, 2011 1:45:35 PM
> Subject: Re: [R] cumulative incidence plot vs survival plot
>
>
> On Jun 27, 2011, at 4:31 PM, array chip wrote:
>
>> Hi, I am wondering if anyone can explain to me if cumulative  
>> incidence (CI) is
>> just "1 minus kaplan-Meier survival"?
>
> First tell us what you think CI is defined as. I suspect it is not  
> the same. The
>
> KM estimator is cumulative product of (alive-n(dead))/alive so is  
> the product of
>
> interval survival probabilities. I doubt that your definition of CI  
> has a
> similar denominator.
>
>
>> Under what circumstance, you should use
>> cumulative incidence vs KM survival? If the relationship is just CI =
>> 1-survival, then what difference it makes to use one vs. the other?
>>
>> And in R how I can draw a cumulative incidence plot.
>
> plot(time, cumsum(dead)) ...?
>
>> I know I can make a
>> Kaplan-Meier survival plot using plot(survfit()), for example:
>>
>> fit<-survfit(Surv(time,status)~group,data=data)
>> plot(fit, col=1:2)
>>
>> How to draw CI plot then?
>
> As above. Specify what you are seeking.
>
> There is a well-defined relationship between S(t) and the cumulative  
> hazard.
> Maybe you should do a little study of those terms in texts regarding  
> survival
> analysis.
>
>> Thanks very much!
>>
>> John
>>    [[alternative HTML version deleted]]
>
> Isn't it time you learned to post in plain text?
>
> --
> David Winsemius, MD
> West Hartford, CT

David Winsemius, MD
West Hartford, CT



More information about the R-help mailing list