我正在尝试为 cox 模型的结果绘制森林图,并考虑受试者 ID 进行重复测量;这是我的大数据的一部分;
data <- read.table(header=TRUE, text="
Id start end QS Age Event
01 0 70 1 25 1
01 70 78 2 25 1
01 78 85 3 25 1
02 0 92 4 23 1
02 92 98 5 23 1
02 98 105 6 23 1
02 105 106 7 23 0
")
考克斯模型为:
set.seed(1234)
model <- coxph(Surv(start, end, Event) ~ QS + Age
+ cluster(Id), data = data, id=Id)
结果:
> model
Call:
coxph(formula = Surv(start, end, Event) ~ QS, data = data, id = Id,
cluster = Id)
coef exp(coef) se(coef) robust se z p
QS -1.896e+01 5.829e-09 1.310e+04 1.000e+00 -18.96 <2e-16
Likelihood ratio test=4.16 on 1 df, p=0.04142
n= 7, number of events= 6
森林图的功能
ggforest(model, data = data)
提供如下图:
这不是我的理想,我不确定这是否是正确的图片。例如为什么id被认为是图片中的变量。