存活曲線在R與survfit

我想繪製一條生存曲線使用以下數據。我叫A.TXT數據文件和對象A存活曲線在R與survfit

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天減少時。我希望在整個圖中標記該標籤，但它只在實驗結束時添加標籤。

來源

2015-04-06 Lothlorian

沒有數據，因此沒有數據可以用於例子。看看'plot'和它的'type =「s」'以及'？stepfun'。並且學會使用代碼製作示例，以便您的答案更有幫助。 –

我很抱歉，因爲我是新用戶，系統不允許我添加一個數字來澄清我想要做什麼。此外，我的生存數據由標爲A和B的兩組組成.A中，第8天和第B天有3名受試者死亡，第8天2名受試者死亡，感染後第10天死亡2名。在第15天的實驗結束時，每組有兩名受試者存活。我試圖爲這些數據做一個圖表。試圖學習基本仍然。非常感謝！ – Lothlorian

仍然沒有數據。使用dput（）來產生R代碼，這是你添加到問題中的東西。（請仔細閱讀說明，不要對評論的內容使用評論。） –

首先，我建議你重新排列繪製電話：

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

的點在生存情節「懸崖」的底部繪製。 enter image description here

來源

2015-04-07 09:11:48 BenBarnes

非常感謝。這對克服這個錯誤非常有幫助。祝一切順利。 – Lothlorian

存活曲線在R與survfit

回答

相關問題