postscript("l04.ps",paper="letter") larynx<-as.data.frame( scan("larynx.dat", what=list(stage=0,time=0,y1=0,y2=0,status=0))) stages14<-larynx[(larynx$stage==1)|(larynx$stage==4),] stages14$fs<-(stages14$status==0)*(abs(stages14$time-5.8)>1.4)+ (stages14$status==1)*(abs(stages14$time-2.5)>1.5) library(survival) source("l04f.R") fulltest(stages14,stages14$time,stages14$status,stages14$fs)