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

Question

我正在尝试使用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

任何人都可以解释我哪里出错了吗？

问候，

Ĵ

Answer 1

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

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

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