我coxph用来估计时变变量 ( tvc) 对生存事件 ( event) 的影响。按照互联网上的几个说明（考虑tmerge），我将数据（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），并且在这些时间步长中，无论 id 是什么，它都具有相同的值）。
任何帮助将非常感激。
谢谢，阿尔贝托