为什么我的生存曲线没有显示为分层类别?

Why are my survival curves not displaying as stratified categories?

第一次发帖,希望这里有足够的信息。

我试图在 4 个类别中显示我的生存曲线。分析是根据我在生存表中的 4 个类别进行分层的,但是生存图并没有描绘这 4 个类别,而是显示了许多不同的生存曲线。我在这里做错了什么?

Survival curve

# categorise ADAMTS13 levels
TMAdata$ADAMTS13level.f<-cut(TMAdata$ADAMTS13level, 
                            breaks=c(0.0,10.0,40.0, 60.0,160.0),
                             labels=c('0-10.0',
                                      '10.1-40.0',
                                      '40.1-60.0',
                                      '60.1-160.0'))
summary(TMAdata$ADAMTS13level.f)

# use 10-40% ADAMTS13 level as reference point
TMAdata$ADAMTS13level.f = relevel(TMAdata$ADAMTS13level.f, ref="10.1-40.0")

# platelet recovery according to ADAMTS13 level (reference point is 10.1-40.0)

    pltrecovery_ADAMTS13_table <- survfit(Surv(TMAdata$Daysplateletrecovery, TMAdata$Recoveredplatelets)~TMAdata$ADAMTS13level.f)
    summary(pltrecovery_ADAMTS13_table)

    plot(pltrecovery_ADAMTS13_table, conf.int=0, 
     xlab = "Days", 
     ylab = "Probability of not achieving platelet count =>150")
     legend("topright", inset=0.03,
            c("0-10.0",
              "10.1-40.0",
              "40.1-60.0",
              "60.1-160.0"),
            lty=1:2,
            lwd=2,
            cex=1)

额外的线是置信边界。指定 conf.int=0 不会抑制置信区间绘图。这可以说是不正确的,因为使用 ?survfit.formula 中的第一个示例很容易证明这一点。如果您不想要置信边界,则将 conf.int 参数一起省略。

图例将只有两种类型的线条,它们可能与绘制的生存图类型不匹配。