r - 使用R包“生存”进行多状态生存分析

我正在尝试使用R survival包估算从状态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"))

当我估计不带frailty参数的模型时，它可以正常工作。但是当我包含它时，出现以下错误:

Error in if (any(ord > 1)) stop("Penalty terms cannot be in an interaction") :
missing value where TRUE/FALSE needed

谁能解释我要去哪里错了？

问候，

Ĵ

最佳答案

我遇到了同样的问题，但最终意识到脆弱的术语是可加的!也就是说，它不应该在逗号之后，但要在加号之后。

将您的代码更改为以下代码应该可以解决问题:

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

关于r - 使用R包“生存”进行多状态生存分析，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/23688701/