Unterscheiden Sie jede Zeile mit verschiedenen Typen in `ggsurv` Plots (oder in` Plot`)

Question

Ich verwende Rstudio. Ich verwende ggsurv Funktion von GGally Paket zum Zeichnen von Kaplan-Meier-Kurven für meine Daten (für Überlebensanalyse), aus Tutorial here. Ich benutze es anstelle von plot weil ggsurv kümmert sich um Legenden von selbst.

Wie auf dem Link gezeigt, werden mehrere Kurven farblich unterschieden. Ich möchte anhand des Linientyps unterscheiden. Das Tutorial scheint keine Option dafür zu haben. mein Befehl folgend:

surv1 <- survfit(Surv(DaysOfTreatment,Survived)~AgeOnFirstContactGroup) 
print(ggsurv(surv1, lty.est = 3)+ ylim(0, 1)) 

lty.est=3 (oder 2) ergibt gleiche gestrichelte Linien für alle Linien. Ich möchte für jede Zeile eine andere gestrichelte Linie haben. Mit lty=type gibt Fehler: object 'type' not found. Und lty=type würde in ggplot arbeiten, aber ggplot nicht direkt mit survfit Plots befassen.

Bitte zeigen Sie mir, wie Kurven in entweder ggsurv oder einfach plot von Linientyp unterscheiden (obwohl ich ggsurv bevorzugen würde, weil sie kümmert sich der Legenden nimmt)

Answer 1

Aus der Dokumentation für ggsurv

lty.est: linetype of the survival curve(s). Vector length should be either 1 or equal to the number of strata.

So, Um für jede Schicht einen anderen Linientyp zu erhalten, setzen Sie lty.est gleich einem Vektor gleicher Länge wie die Anzahl der Linien, die Sie zeichnen, wobei jeder Wert einem anderen Linientyp entspricht.

Zum Beispiel unter Verwendung der Lungen Daten aus dem survival Paket

library(GGally) 
library(survival) 
data(lung) 
surv1 <- survfit(Surv(time,status) ~ sex, data = lung) 
ggsurv(surv1, lty.est=c(1,2), surv.col = 1) 

gibt die folgende Handlung

Sie können auch ggplot Themen oder andere ggplot Elemente der Handlung hinzufügen . Zum Beispiel können wir das Aussehen mit dem cowplot Thema verbessern

library(ggplot2) 
library(cowplot) 
ggsurv(surv1, lty.est=c(1,2), surv.col = 1) + theme_cowplot() 

Wenn Sie die Legende folgt Etiketten nach Differenzierung durch Linientyp ändern müssen, dann können Sie es auf diese Weise

ggsurv(surv1, lty.est=c(1,2), surv.col = 1) + 
    guides(colour = FALSE) + 
    scale_linetype_discrete(name = 'Sex', breaks = c(1,2), labels = c('Male', 'Female')) 

tun