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)

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

1 个答案: