\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) }