20 Aug 08:30
cmprsk and a time dependent covariate in the model
From: <phguardiol <at> aol.com>
Subject: cmprsk and a time dependent covariate in the model
Newsgroups: gmane.comp.lang.r.general
Date: 2008-08-20 06:31:49 GMT
Subject: cmprsk and a time dependent covariate in the model
Newsgroups: gmane.comp.lang.r.general
Date: 2008-08-20 06:31:49 GMT
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)
RSS Feed