Hello All,
I am trying to use the unfold function in RcmdrPlugin.survival library,
which converts the survival data with time varying covariates to the
counting process notation. The problem is somehow, the event indicator
created is not correct.
Below is the data, I am trying to convert:
CASE TRT FAILTIME FAILCENS SEX AGE IGG0 IGG28 IGG42 IGG84 IGG364
26003 A 11.2033 2 1 43.9535 0 80 80
320 NA
IGG0 to IGG364 are my time varying covariates and FAILCENS = 2 indicates
observed times. So, when I tried to use,
igg2.long<-unfold(igg.2[3,],time="FAILTIME",event="FAILCENS",cov=9:13,cov.na
mes="iggtitre"), I obtained the following output:
start stop FAILCENS.time CASE TRT FAILTIME FAILCENS SURVTIME SURVCENS
SEX
3.1 0 1 0 26003 A 11.2033 1 17.5441 1
1
3.2 1 2 0 26003 A 11.2033 1 17.5441 1
1
3.3 2 3 0 26003 A 11.2033 1 17.5441 1
1
3.4 3 4 0 26003 A 11.2033 1 17.5441 1
1
AGE IGGrel iggtitre
3.1 43.9535 NA 0
3.2 43.9535 NA 80
3.3 43.9535 NA 80
3.4 43.9535 NA 320
>
As it can be seen, this data was observed, but my censoring indicator
(FAILCENS.time) for the last start stop line should have been one, which is
not.
Could somebody help me with this problem? I did try using the Rossi data
also, which is an example in the unfold function. Even there, the
arrest.time are all zero, even when the event happened.
Thanks a lot for your time.
Mathangi
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
