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