时间戳

首页 > 解决方案 > 在 R 中使用 mstate

r - 在 R 中使用 mstate

问题描述

我正在使用该mstate程序来估计疾病-死亡模型，并有一些问题希望经验丰富的 mstate 用户能够对此有所了解。

我的第一个问题很基本：

如果开始时间并不总是大于停止时间，这会是一个问题吗？例如，以下是数据中 Tstart>=Tstop 的一些案例：

       ID   Tstart    Tstop trans status
480   228 5.111874 5.111874     3      0
587   277 4.996568 4.996568     3      0
1184  567 5.015786 5.015786     3      0
1258  603 0.000000 0.000000     1      1
1259  603 0.000000 0.000000     2      0
1397  667 0.000000 0.000000     1      1
1398  667 0.000000 0.000000     2      0

如您所见，问题实际上是我的 Tstart=Tstop。

我像这样设置我的模型： model <- coxph(Surv(Tstart, Tstop, status) ~ x.1 + x.2 + x.3 .... + strata(strata), data = example, method = "breslow")，在运行我的模型（确实运行并提供估计）之后，这是我收到的错误消息：

Stop time must be > start time, NA created

我的第二个问题有点复杂。也就是说，有没有办法包含一个不是新状态的时间更新变量？例如，如果我对疾病 X 的风险感兴趣，但感兴趣的暴露可能会在研究进入时间和疾病 X 诊断之间发生变化。

标签： rsurvival-analysis

解决方案

推荐阅读