## Lab 15 R Code

anemia <- read.delim("http://myweb.uiowa.edu/pbreheny/data/anemia.txt")
attach(anemia)

require(survival)
S <- Surv(Time,Status!=0)
fit <- survfit(S~1)
plot(fit, ylab = "Probability", xlab = "Time")

survdiff(S~Trt)
fit <- survfit(S~Trt)
plot(fit, ylab = "Probability", xlab = "Time", col = c("red","blue"))
legend("bottomleft", c("MTX","MTX + CSP"), text.col = c("red","blue"), bty = "n")