phguardiol | 20 Aug 08:30
Favicon

cmprsk and a time dependent covariate in the model


 Dear R users,

I d like to assess the effect of "treatment" covariate on a disease relapse risk with the package cmprsk. 

However, the effect of this covariate on survival is time-dependent
(assessed with cox.zph): no significant effect during the first year of follow-up,
then after 1 year a favorable effect is observed on survival (step
function might be the correct way to say that ?). 
For overall survival analysis
I have used a time dependent Cox model which has confirmed this positive effect after
1 year.

Now I m moving to disease relapse incidence and a similar time dependency seems to be present. 

what I d like to have is that: for
patients without "treatment" the code for "treatment" covariate is
always 0, and for patients who received "treatment" covariate I d like
to have it = 0 during time interval 0 to 1 year, and equal to 1 after 1
year. Correct me if I m wrong in trying to do so.

First, I have run the following script (R2.7.1 under XPpro) according to previous advices:

library(cmprsk)

attach(LAMrelapse)

fit1<- crr(rel.t, rel.s, treatment, treatment, function(uft)
cbind(ifelse(uft<=1,1,0),ifelse(uft>1,1,0)), failcode=1,
cencode=0, na.action=na.omit, gtol-06, maxiter)
(Continue reading)


Gmane