[R] multi-level cox ph with time-dependent covariates

Terry Therneau therneau at mayo.edu
Thu Nov 4 14:26:40 CET 2010


Your question has two levels:
   1. What is the right model for this data
   2. Can model __ be fit

Wrt 2 and coxme: For a reliable fit you need to have more events than
random effects.  Thus for patient/tissue I would want to see multiple
events per patient/tissue pair.  This is statistical issue -- when there
are too few events the confidence intervals for the random effects end
up being a mile wide.  (Exception, if the number of events is very
large, >10^5 say as sometimes occurs in economics studies, the estimates
can work.)  
   coxme works fine with start,stop data.

Wrt question 1.  Your models assume that marker1, marker2, ... each have
the same effect across tissue types.  Adding a random effect gave per
subject or per subject/tissue intercepts.  Do you instead want to do
shrinkage of the marker1, .. coefficients?

Terry Therneau



More information about the R-help mailing list