01-05-Stat-Methods-Survival.Rmd

## Survival Analysis {-}

```{r a12, comment=NA, fig.align='center'}

library(survival)

summary(lung)

mdl = survfit(Surv(time, status) ~ 1, data = lung)

plot(mdl, conf.int = FALSE,
     main="Kaplan-Meier Estimator of Survival Function",
     xlab="Survival Time",
     ylab="Survival Function")

summary(mdl)

mdl = survfit(Surv(time, status) ~ strata(sex), data = lung)

summary(mdl)

plot(mdl, conf.int = FALSE, col = c(2,4),
     main="Survival Function",
     xlab="Survival Time",
     ylab="Survival Function")

legend(800, 1, c("men", "women"), lty = c(1, 1), col = c(2, 4))


summary(survreg(Surv(time, status) ~ ph.ecog + age + strata(sex), lung))
```