r - 使用 R 包 "Survival"的多态生存分析

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

谁能解释一下我哪里出错了？

问候，

J

Answer 1

我遇到了同样的问题，但最终意识到 frailty 术语是累加的!也就是说，它不应该出现在逗号之后，而是出现在加号之后。

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

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