随机效应Cox模型

Question

我正在尝试绘制pm2.5暴露与高血压发病率之间的剂量反应关系。我拟合了一个随机效应cox模型，其中为研究中心添加了随机效应。然后，我使用plotHR函数绘制剂量反应关系。但是我遇到了一个错误。下面是我使用的示例R代码。

library(survival)
library(coxme) 
library(splines)
library(Greg)

data("eortc")
eortc$age<-rnorm(2323,40,10)

efit1 <- coxph(Surv(y, uncens) ~ ns(age) + trt , eortc)
plotHR(efit1,term = 1,xlim = c(30,70))

efit2 <- coxme(Surv(y, uncens) ~ ns(age) + trt + (1|center), eortc)
plotHR(efit2,term = 1,xlim = c(30,70))

我可以使用plotHR绘制efit1，但是在绘制efit2时遇到错误，其中我在cox模型中添加了随机效果。有人知道如何解决这个问题吗？谢谢！

Answer 1

这显示predict.coxme对于这种模型可以返回什么。我猜它不会完全令人满意，但是我声称它很有用，因为它演示了混合模型在“引擎盖下”的样子。每个中心都有单独的预测，但是由于没有人声称它们具有任何特定的分布，因此从统计上考虑将其汇总的努力被认为是可疑的。这将绘制模型的指数线性预测变量，其中涉及不带结的“样条曲线”拟合（这会给您带来一堆线条）。我正在按trt状态着色：黑色代表“ 0”，红色代表“ 1”

plot( eortc$age, exp( predict(efit2)), col=1+(eortc$trt==1) )

trt == 0和trt == 1组之间的“平均”差异确实出现了，并且与exp(0.705) -> [1] 2.023847的实测治疗效果一致，而模型中的“年龄”效果则没有重要的是非常浅的线性上升。