具有survfit的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)

我想放標記（三角形為A和"+"為B ）每當存活％在7天和第8天我想整個圖表這個標記例如減小，但它僅在末端添加標簽本實驗。

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

這些點繪制在生存圖中“懸崖”的底部。