plot( survfit (Surv(time, cens)  ~ horTh, data = GBSG2), lty = 1:2, mark.time = FALSE, ylab = "Probabilità",  xlab = "tempo di sopravvivenza in giorni")
legend (250, 0.2, legend = c( "sì", "no") , lty = c(2,1),  title = "Terapia ormonale", bty = "n")


####  interessante scoprire qual è l'output del comando


survfit (Surv(time, cens)  ~ horTh, data = GBSG2)