Seite 1 von 1
Kaplan Meier Schätzer
Verfasst: Mi Aug 09, 2017 6:51 pm
von Hendrik711
Hallo Leute,
ich bin Neu hier und weis nicht so genau ob mein Thema hier richtig ist. Ich möchte den Kaplan-Meier-Schätzer anwenden.
Ganz einfaches Beispiel:
Code: Alles auswählen
rm(list=ls(all=TRUE))
library(survival)
zeiten <- c(7, 3, 24, 19, 8, 7, 11, 19)
delta <- c(1, 1, 1, 1, 1, 1, 0, 0)
survivalobjekt <- Surv(zeiten, delta)
km1 <- survfit(survivalobjekt ~ 1)
plot(x=km1, col=c("black", "red", "red"), xlab="x", ylab="S(x)", lty=1)
Jetzt habe ich zwei Probleme.
1) Die zensierten Eintrage werden nicht angezeigt.
2) Konfidenzintervall bricht ab
Das habe ich:
Das sollte raus kommen:
Wäre echt super wenn einer sagen kann, woran das liegt.
Mfg Hendrik
Re: Kaplan Meier Schätzer
Verfasst: Fr Aug 11, 2017 6:40 pm
von Hendrik711
So das mit dem Konfidenzintervallen:
Code: Alles auswählen
Call: survfit(formula = survivalobjekt ~ 1)
time n.risk n.event survival std.err lower 95% CI upper 95% CI
3 8 1 0.875 0.117 0.673 1.000
7 7 2 0.625 0.171 0.365 1.000
8 5 1 0.500 0.177 0.250 1.000
19 3 1 0.333 0.180 0.116 0.961
24 1 1 0.000 NaN NA NA
Die Konfidenzintervalle klappen nicht, wegen diesen NA. Wie ignoriere ich die? ich würde gerne gerade
Linien haben, wie bei dem Bild im ersten Beitrag.
Liebe Grüße, Hendrik
Re: Kaplan Meier Schätzer
Verfasst: Sa Aug 12, 2017 11:24 am
von jogo
Hallo Hendrik,
willkommen im Forum!
Hendrik711 hat geschrieben: ↑Fr Aug 11, 2017 6:40 pm
So das mit dem Konfidenzintervallen:
Code: Alles auswählen
Call: survfit(formula = survivalobjekt ~ 1)
time n.risk n.event survival std.err lower 95% CI upper 95% CI
3 8 1 0.875 0.117 0.673 1.000
7 7 2 0.625 0.171 0.365 1.000
8 5 1 0.500 0.177 0.250 1.000
19 3 1 0.333 0.180 0.116 0.961
24 1 1 0.000 NaN NA NA
Die Konfidenzintervalle klappen nicht, wegen diesen NA. Wie ignoriere ich die?
Vielleicht hilft
na.omit() oder Du könntest in diesem Fall versuchen, die letzte Beobachtung explizit zu löschen.
Gruß, Jörg
Re: Kaplan Meier Schätzer
Verfasst: Sa Aug 12, 2017 6:40 pm
von Hendrik711
Danke für deine Antwort Jörg!
Sorry ich bin echt nicht geübt mit R. Kannst du vielleicht genauer sagen, wie das funktionieren soll?
Ich habe jetzt einfach na.omit(km1) probiert, aber die NA´s bleiben trotzdem da.
Liebe grüße
Hendrik
Re: Kaplan Meier Schätzer
Verfasst: So Aug 13, 2017 7:55 am
von jogo
Hallo Hendrik,
zur Zeit kämpfe ich noch damit, die gleichen Ergebnisse herauszubekommen, wie Du sie hier präsentierst. Bei mir gibt es immer nur:
Code: Alles auswählen
> survfit(formula = survivalobjekt ~ 1)
Call: survfit(formula = survivalobjekt ~ 1)
n events median 0.95LCL 0.95UCL
8.0 6.0 13.5 7.0 NA
Gruß, Jörg
Re: Kaplan Meier Schätzer
Verfasst: Mi Aug 16, 2017 7:01 pm
von Hendrik711
Da musst du einfach nur summary(km1) eingeben. Aber Danke echt für dein Bemühen! Ich habe schon die Lösungen gefunden. Package ggfortify löst das Problem!