R中的生存曲线与幸存者

Question

我想使用以下数据绘制生存曲线。我将数据文件称为 A.txt 和对象A

A <- read.table(text = "
Time Status Group
8 1 A
8 1 A
8 1 A
9 1 A
9 1 A
9 1 A
15 0 A
15 0 A
7 1 B
7 1 B
8 1 B
9 1 B
10 1 B
10 1 B
15 0 B
15 0 B", header = TRUE)

我尝试使用此代码绘制生存曲线：

title(main="Trial for Survival Curve")
fit <- survfit(Surv(Time, Status) ~ Group, data = A)
par(col.lab="red")
legend(10, .9, c("A", "B"), pch=c(2,3),col=2:3)
plot(fit, lty=2:3, col=2:3,lwd=5:5, xlab='Time(Days)',
  ylab='% Survival',mark.time=TRUE,mark=2:3)

每当生存率降低时（例如在第7天和第8天），我想为A和"+"添加标记（我需要B和{{1}}。我想在整个过程中使用此标记图表，但它仅在实验结束时添加标签。

Answer 1

首先，我建议重新安排绘图调用：

par(col.lab="red")
plot(fit, lty=2:3, col=2:3,lwd=5:5, xlab='Time(Days)',
  ylab='% Survival',mark.time=TRUE,mark=2:3)
title(main="Trial for Survival Curve")
legend(10, .9, c("A", "B"), pch=c(2,3),col=2:3)

您可以使用points功能为生存图添加点。但是，看起来有一个小错误，你可以很容易地解决这个问题：

firsty <- 1 ## Gets around bug
points(fit[1], col = 2, pch = 2) # Plots first group in fit
points(fit[2], col = 3, pch = 3) # Plots second group in fit

这些点被绘制在生存图中“悬崖”的底部。