пловец выживания
Есть ли простой способ создать пловец в R? Те же данные, что и на кривой КМ, но каждая выживаемость представлена в виде линии. Пример:
Я искал stackoverflow, список рассылки R-help и консультировался с доктором Google без очевидного ответа, хотя моя техника поиска может быть неоптимальной. Спасибо!
**** ДОБАВЛЕНО **** Прошу прощения за неправильное задание вопроса - это мой первый раз! Играя вокруг, я смог сделать следующее:
OS DeathYN TreatmentGroup
4 444 days 1 0
5 553 days 1 0
8 812 days 0 0
1 844 days 0 0
10 1071 days 0 0
9 1147 days 0 0
6 1349 days 0 0
3 1375 days 0 0
2 1384 days 0 1
7 1687 days 0 0
orderedData$GroupColor[orderedData$TreatmentGroup==0] <- "yellow"
orderedData$GroupColor[orderedData$TreatmentGroup==1] <- "red"
orderedData$YCoord <- barplot(as.numeric(orderedData$OS), horiz=TRUE, col=orderedData$GroupColor, xlim=c(0,max(orderedData$OS) + 50), xlab="Overall Survival")
points(x=20+as.numeric(orderedData$OS), y=orderedData$YCoord,pch=62, col="green")
legend(1000,2, c("Control", "Treatment", "still living"), col=c("yellow","red", "green"), lty=1, lwd=c(10,10,0),pch=62)
На данный момент это достаточно близко, но эстетика не идеальна. Если есть пакет или лучшее решение, кто-то может подсказать, я хотел бы увидеть его!