使用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)