ich möchte für meine Doktorarbeit eine Grafik der kumulierten Inzidenz für Competing Risks erstellen. Dafür nutze ich die Funktion cuminc() aus dem Paket 'cmprsk' und rufe danach plot() auf. Das funktioniert soweit auch schon ganz gut.
Ich würde gerne zusätzlich unter der Grafik eine "Number at risk" Tabelle darstellen, ganz einfach die Menge der verbliebenen Fälle in der Interventions- und Kontrollgruppe. Leider kann ich dazu keine Möglichkeit finden. Wie könnte ich das umsetzen?
Ich weiß, dass es auch noch weitere Pakete wie 'ggcompetingrisk()' etc. gibt, aber mir gefällt bisher die ganz einfache plot() Funktion graphisch am besten und deshalb würde ich am liebsten dabei bleiben.
Der Plot entsteht mit:
Code: Alles auswählen
cuminc_fit <- cuminc(ftime=time,fstatus=status,group=inter,cencode=0)
plot(cuminc_fit)
Code: Alles auswählen
surv_fit <- survfit(Surv(time,status) ~ inter)
risk_table <- get_nrisk_tbl(surv_fit,2)
Das ganze müsste ich jetzt nur noch unter die jeweiligen Zeitpunkte in dem Plot bekommen.strata time n.risk
1 inter=0 0 115
2 inter=0 2 81
3 inter=0 4 55
4 inter=0 6 33
5 inter=0 8 21
6 inter=0 10 11
7 inter=1 0 18
8 inter=1 2 8
9 inter=1 4 4
10 inter=1 6 2
11 inter=1 8 1
12 inter=1 10 1
Vielen Dank
So würde ich mir das in etwa vorstellen: