4

我正在尝试为 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被认为是图片中的变量。

4

0 回答 0