Hallo EDi,
mega, das ist genau das was ich brauchte.
Dass es etwas unübersichtlich aussieht ist erstmal nicht schlimm. Ich möchte letztendlich nur eines der Events vergleichen und die Kurven trennen sich etwas besser. Sieht erstmal ganz gut aus
Könntest du mir noch helfen das gleiche nochmal mit dem normalen "survfit()" aus dem survival Package zu schaffen?
Hier wieder die Testdaten:
Code: Alles auswählen
set.seed(2)
ss1 <- rexp(100)
gg1 <- factor(sample(1:2,100,replace=TRUE),1:2,c('a','b'))
cc1 <- sample(0:1,100,replace=TRUE)
surv_fit1 <- survfit(Surv(ss1,cc1)~gg1)
set.seed(3)
ss2 <- rexp(100)
gg2 <- factor(sample(1:2,100,replace=TRUE),1:2,c('a','b'))
cc2 <- sample(0:1,100,replace=TRUE)
surv_fit2 <- survfit(Surv(ss2,cc2)~gg2)
plot(surv_fit1)
#mein Versuch
lines(surv_fit2$time[1:unname(surv_fit2$strata)[1]],surv_fit2$surv[1:unname(surv_fit2$strata)[1]])
lines(surv_fit2$time[(unname(surv_fit2$strata)[1]+1):length(surv_fit2$time)],surv_fit2$surv[(unname(surv_fit2$strata)[1]+1):length(surv_fit2$surv)])
Das Problem liegt wieder darin die eigentlichen Koordinaten für die KM-Kurve aus dem survfit-Objekt herauszubekommen. Ich habe auch mal mit seq_along() durchgeguckt, aber da sind diesmal deutlich mehr Listen enthalten und keine davon sieht wirklich brauchbar aus. Mit $time und $surv kriegt man grob den Verlauf der Kurve, aber eben nicht in der typischen KM-Form mit nur Bewegung in entweder x- oder y-Richtung. Hoffe es gibt hierfür auch so eine elegante Lösung.
Danke