# NOT RUN {
require(survival)
cancer$sex <- factor(ifelse(cancer$sex == 1, "male", "female"))
mod.1 <- coxph(Surv(time, status) ~ age + wt.loss, data=cancer)
plot(mod.1)
plot(mod.1, typical=function(x) quantile(x, c(.25, .75)))
mod.2 <- coxph(Surv(time, status) ~ age + wt.loss + sex, data=cancer)
plot(mod.2)
mod.3 <- coxph(Surv(time, status) ~ (age + wt.loss)*sex, data=cancer)
plot(mod.3)
mod.4 <- coxph(Surv(time, status) ~ age + wt.loss + strata(sex), data=cancer)
plot(mod.4)
mods.1 <- survreg(Surv(time, status) ~ age + wt.loss, data=cancer)
# }
Run the code above in your browser using DataLab