使用R包的多状态生存分析；“生存”；

我试图使用R

survival

package估算从状态1过渡到状态2的Cox比例风险模型，如下所示：

Altman <- coxph(Surv(Tstart, Tstop, to == 2) ~ wWCTA + wRETA + wEBITTA + wMETL + 
            wSTA, data=Multistate, subset = from == 1, 
          frailty(x=id, distribution= "gamma"))

谁能解释一下我哪里出了问题

问候,

---

我也有同样的问题，但最终意识到脆弱性这个术语是加性的！也就是说，它不应该在逗号之后，而应该在加号之后

因此，将代码更改为以下应该可以解决问题：

Altman <- coxph(Surv(Tstart, Tstop, to == 2) ~ wWCTA + wRETA + wEBITTA + wMETL + wSTA + frailty(x=id, distribution= "gamma"), data=Multistate, subset = from == 1)

---

Altman通常Therneau建议现在对脆弱模型使用coxme。
Altman <- coxph(Surv(Tstart, Tstop, to == 2) ~ wWCTA + wRETA + wEBITTA + wMETL + wSTA + frailty(x=id, distribution= "gamma"), data=Multistate, subset = from == 1)