Блог

Я использую coxph для оценки влияния изменяющейся во времени переменной ( tvc ) на событие выживания ( event ). Следуя нескольким инструкциям в Интернете (учитывая tmerge ), я подготовил data ( df ) как:

 id  country   start stop    event   endpoint tvc    time
1   a         0     3       0       0        0.1    9
1   a         3     6       0       0        0.2    9
1   a         6     9       0       0        0.3    9
2   a         0     3       1       0        0.1    8
2   a         3     6       1       0        0.2    8
2   a         6     8       1       1        0.3    8
3   b         0     3       1       0        0.1    5
3   b         3     5       1       1        0.2    5
  

Тем не менее, когда я запускаю fit = coxph(Surv(start, stop, event) ~ country tvc, data=df) , он возвращает правильный коэффициент только для каждой страны, а коэффициент для результатов tvc NA.
Поэтому я не понимаю, есть ли что-то не так с моей формулой или это каким-то образом неправильно рассматривать tvc как изменяющуюся во времени ковариацию (поскольку она изменяется с точными временными шагами (0-3 = 0.1, 3-6 = 0.2, 6-9 = 0.3) и на этих временных этапах он имеет одно и то же значение независимо отидентификатор).

Любая помощь будет высоко оценена.

Спасибо, Альберто