1
 2
 3
 4
 5
 6
 7
 8
 9
10
\message{ !name(Introduction.Rnw.tex)}
\message{ !name(Introduction.Rnw) !offset(291) }
fit <- stpm2(Surv(rectime,censrec==1)~hormon,
             data=brcancer, link.type="AH")
summary(fit)
plot(fit, newdata=data.frame(hormon=0))
lines(fit, newdata=data.frame(hormon=1), lty=2)
lines(survfit(Surv(rectime,censrec==1)~hormon, data=brcancer), col="blue", lty=1:2)

\message{ !name(Introduction.Rnw.tex) !offset(-10) }