[R] Competing risks Kalbfleisch & Prentice method
er339 at medschl.cam.ac.uk
Fri Mar 27 15:44:19 CET 2009
Dear Prof. Therneau,
Thank you for your views on this subject. I think all R users who play
with survival analysis are most grateful for the functions you have
already supplied us with.
I'm guessing Ravi is wondering why you have not implemented the
smoothing of the baseline hazard from the Cox model.
I actually tried to do this originally, inspired from this thread (i.e
use sm.spline to smooth the hazard):
but it overestimated the CI (perhaps I implemented it wrong). I was then
advised to treat CI as a step function, rather than continuous, which
means that F(t+1, cause k)-F(t, cause k) will be 0 unless an event of
cause k has occurred in that interval (see also "Competing Risks, by
Melanie Pintilie, page 62). This is obviously problematic if one wants
to estimate the CI at times that are not close to observed events for
either cause (perhaps a parametric model could be used in this case).
But then again, this was not an issue wtih my data.
Strangeways Research Laboratory
Department of Public Health and Primary Care
University of Cambridge
From: Terry Therneau [mailto:therneau at mayo.edu]
Sent: 27 March 2009 13:53
To: Eleni Rapsomaniki; tuechler at gmx.at; Ravi Varadhan
Cc: r-help at r-project.org
Subject: RE: Competing risks Kalbfleisch & Prentice method
Ravi's last note finished with
> I am wondering why Terry Therneau's "survival" package doesn't
> have this option.
The short answer is that there are only so many hours in a day.
I've recently moved the code base from an internal Mayo repository to
one long term goal with this is to broaden the developer base to n>2 (me
A longer statistical answer:
I'm not sure if the "this" of Ravi's question is a. smoothed hazards,
K&P cumulative incidence or c. the Fine & Gray model.
b. I like the CI model and am using it more. We also have local code.
latest version of survival (on rforge, likely in the next default R
added simple CI curves to the survfit function. Adding code for survfit
models is on the todo list. But -- this release also fixes up
handle weighted Cox models and that was on my list for approx 10 years,
don't hold your breath. I don't release something until it also has a
worked out test cases to add to the 'tests' directory.
a. smoothed hazards. For the case at hand I don't see any particular
advantage of this. On the other hand, I often would like to display
functions instead of CI functions for Cox models; with time dependent
I don't think a survival curve makes sense. But I haven't had the time
through exactly which methods should be added.
c. Fine & Gray model, i.e., where covariates have a direct influence
competing risk. I find the model completely untenable from a biologic
view, so have no interest in adding it. (Due to finite time, everything
survival package is code that I needed for an analysis; medical research
pays my salary.) Assume that I have competing processes/risks, say
of a tumor and heart disease; I expect that the tumor process pays no
whatsoever to what is going on in the heart. But this is necessary if
"type=squamous" is modeled as an absolute beta=__ increase in the CI for
The squamous cells need to "step up the pace" of invasion if heart
threatens, like jockeys in a horse race.
More information about the R-help