我试图绘制 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 模型中添加了随机效应。任何人都知道如何解决这个问题?谢谢!