第一次海報，所以我希望我在這里有足夠的資訊。
我試圖在 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)
uj5u.com熱心網友回復：
額外的線是置信邊界。指定 conf.int=0 不會抑制置信區間繪圖。這可以說是不正確的，因為使用?survfit.formula. 如果您不想要置信邊界，則將 conf.int 引數一起省略。
圖例只有兩種型別的線，它們可能與繪制的生存型別不匹配。
轉載請註明出處，本文鏈接：https://www.uj5u.com/yidong/343040.html
上一篇：將資料幀中每一行中的向量分成兩半