R coxph返回时变协变量的NA

我使用

coxph

来估计时变变量（

tvc

）对生存事件（

事件

）的影响。 按照互联网上的几个说明（考虑到

tmerge

），我准备数据（

df

），如下所示：

然而，当我运行

fit=coxph（Surv（开始、停止、事件）~country+tvc，data=df）

时，它只返回每个国家的适当系数，tvc结果的系数为NA。 因此，我不知道我的公式是否有问题，或者它在某种程度上是错误的，将tvc作为时变协变量考虑（因为它在精确的时间步长上变化（0-3=0.1，3-6=0.2，6-9=0.3），并且在这些时间步长中，无论id如何，它都具有相同的值）

任何帮助都将不胜感激

谢谢，

---

阿尔贝托

你有没有偶然发现这个问题？我还是不确定。问题是时变协变量在IDs中是静态的，因此（我假设）代码会为此返回NAN。当我尝试一个tvc，它不仅在不同的时间，而且在不同的ID上呈现不同的值时，这个公式起了作用。

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