我试图绘制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模型中添加了一个随机效果。谁知道如何解决这个问题?谢谢!
这表明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
的测量治疗效果一致,并且模型中不显着的“年龄”效应是非常浅的线性上升。