r语言 - 为什么我的生存曲线没有显示为分层分类?

Question

第一次贴海报，希望这里有足够的信息。

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

存活曲线

# 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)

Answer 1

额外的行是置信界限。指定conf.int=0不会抑制置信区间绘图。这可能是不正确的，因为使用?survfit.formula中的第一个示例很容易演示。如果您不想要置信边界，那么请将conf.int参数全部省略。

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