[R] A question about external time-dependent covariates in cox model
Rui Song
rsong at stat.wisc.edu
Thu Aug 19 05:20:30 CEST 2004
Dear Sir or Madam:
I am a graduate student in UW-Madison statistics department. I have a
question about fitting a cox model with external time-dependent
covariates.
Say the original data is in the following format:
Obs Eventtime Status Cov(time=5) Cov(time=8) Cov(time=10) Cov(time=12)
1 5 1 2
2 8 0(censored) 2 4
3 10 1 2 4 6
4 12 1 2 4 6 8
....
Notice that the time-dependent covariates are identical at the same
time points for all obs since they are external to the failure process.
process.
Then I organized the data as the following:
obs start end eventtime status cov
1 0 5 5 1 2
2 0 5 8 0 2
2 5 8 8 0 4
3 0 5 10 1 2
3 5 8 10 1 4
3 8 10 10 1 6
4 0 5 12 1 2
4 5 8 12 1 4
4 8 10 12 1 6
4 10 12 12 1 8
And fit the model using:
fit<-coxph(Surv(start, end, status)~cov);
When I fit the model to my data set (Which has 89 observations and 81
distinct time points, sort of large.), I always got a message that
"Process R segmentation fault (core dumped)". Would you let me know if it
is due to the matrix sigularity in the computation of the partial
likelihood or something else? And how should I fit a cox model with
external time-dependent covariates?
Thanks a lot for your time and help!
Sincerely,
Rui Song
More information about the R-help
mailing list