Re: [R] Cumulative dose in survival models

From: peter dalgaard <>
Date: Tue, 10 May 2011 15:51:24 +0200

On May 10, 2011, at 10:52 , wrote:

> Dear List,
> I want to fit a CPH survival model with a time-dependent cumulative exposure variable. More specifically, I have exposure measurements for all cohort members at discrete times (yearly measurements) during follow-up.
> I suppose that when estimating the model parameters the most correct would be to have the cumulative exposure for everybody at risk at each death time, as I cannot assume the cumulative exposure to be constant between the yearly measurements . Thus, I need to calculate (interpolate) cumulative exposure for everybody at risk at each death time and pass this to the estimation routine.
> Does anyone know how to do this in R? Any help would be much appreciated.

You didn't say what your data layout is. One record per yearly measurement?

Generic time splitting is discussed in the 2nd ed. of ISwR (or peek into the scripts in the ISwR package). That is for a simplified situation (5 year age groups) but should generalize to arbitrary cutpoints with some elaboration. Then you "just" need to figure out the exposure. If you can things narrowed down to

then it would be obvious how to proceed. One worry would be whether you have the current exposure for those who die. I.e., you do not want to use interpolations that depend on _next_ year's measurements.

Peter Dalgaard
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email:  Priv:

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Tue 10 May 2011 - 13:55:21 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 10 May 2011 - 14:20:06 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive