[R] A stats question -- about survival analysis and censoring
David Winsemius
dwinsemius at comcast.net
Tue Mar 11 01:29:28 CET 2008
"Geoff Russell" <geoffrey.russell at gmail.com> wrote in
news:93c3eada0803100546l2a604d53jead5b22bf39bf26c at mail.gmail.com:
> Dear UseRs,
>
> Suppose I have data regarding smoking habits of a prospective cohort
> and wish to determine the risk ratio of colorectal cancer in the
> smokers compared to the non-smokers. What do I do at the end of the
> study with people who die of heart disease? Can I just censor them
> exactly the same as people who become uncontactable or who die in a
> plane crash? If not, why not?
>
> I'm thinking that heart disease isn't independent of smoking even
> though a death from heart disease is probably uninformative about
> colorectal cancer risk. Hence,I suspect simply censoring these deaths
> will introduce a bias, but I don't know how to correct for it.
>
This isn't really an R question, of course. Quick answer; don't censor if
the subject survives a heart disease event. One is still at risk for
colorectal cancer.
I am not sure there is a completely correct statitical method for
correcting this sort of informative censoring induced by competing risks
if your focus is colonic neoplasms. There would also be other cancers
(oropharyngeal, laryngeal, lung, kidney, bladder) that would also pose
similar potentially informative censoring events, especially if the
subject dies, eh. I think you can take solace in that fact that any
positive association you find will be an underestimate regarding
smoking's induction of colorectal neoplasms.
A full treatment might be a Markov model that would allow transitions to
"heart disease", "colorectal cancer", "other smoking related cancers" and
"death". You can still transition to colorectal cancer after heart
disease and I doubt that the risk of colorectal cancer is much changed
after making the transition into "heart disease" except for the fact that
such subjects will be bombarded by messages to stop smoking. I would
predict that one's risk of staying in the surviving risk set would be
very much loweer after transitioning to lung or oropharyngeal cancer,
however.
--
David Winsemius
More information about the R-help
mailing list