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

Question

我正在尝试使用 Rsurvival包估计从状态 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)