我正在使用 R 中的 coxme 函数 {coxme} 运行混合效应 Cox 模型，并且我想检查比例风险的假设。
我知道 PH 假设可以通过 cox.ph 模型上的 cox.zph 函数 {survival} 进行验证。
但是，我找不到 coxme 模型的等价物。
2015 年在这里发布了一个类似的问题，但没有答案。
我的问题是：1）如何测试混合效应 cox 模型 coxme 的 PH 假设？2) 如果 coxme 模型没有等效的 cox.zph，在科学文章中发表运行混合效应 coxme 模型但在与 coxme 模型相同但没有随机效应的 cox.ph 模型上测试 PH 假设是否有效?
提前感谢您的回答。问候
您可以在函数中使用frailty选项。coxph假设您的随机效应变量是B，您的固定效应变量是A。然后你适合你的模型如下
myfit <- coxph( Surv(Time, Censor) ~ A + frailty(B) , data = mydata )
现在，您可以使用cox.zph(myfit)比例风险假设来检验。
我没有足够的声誉发表评论，但我认为frailty在函数中使用该选项coxph不会起作用。在cox.zph文档中，它说：
随机效应术语（如coxme 模型中的脆弱或随机效应）不会检查比例风险，而是将它们视为模型中的固定偏移量。
因此，在测试比例风险假设时，它没有考虑随机效应。