[R] understanding loops for "loop-plotting"

[Christian Schulz]

yes, yesterday when i go to
sleep i recognize my stupid error!

Thanks,Christian

Am Dienstag, 17. Februar 2004 03:42 schrieb Paul Murrell:
> Hi
>
> If you make the pdf() call BEFORE entering the loop does it do what you
> want?
>
> Paul
>
> Christian Schulz wrote:
> > Hi,
> >
> > i make some practice on "loops"  to understand this important
> > issue in "difference" to the wrappers like  apply  finally!
> > The point is that i want plot's for different different const  and wish
> > in one pdf the 5 plots . So i want  plots for const=1:5 which change the
> > undelying data expressed in t1 and t2.
> >
> > I attempt  some trials , while, for  and now
> > repeat loop now , but didn't getting success . I get only the first plot 
> > and the loop hang?
> >
> > Many thanks for help ,Christian.
> >
> > #modeldat  and YS are data.frames.
> > special <-  function(const,modeldat,YS) {
> > const=const
> > repeat {
> >
> > t1  <-  apply(YS,1, function(x) {  ifelse(all(is.na(x)) | all(na.omit(x)
> > < 0 ,NA, which( x > const ))})
> >
> > t1[is.na(t1)] <-  13
> >
> > t2  <-  sapply(t1,function(x) { ifelse(x ==13,0,1)})
> >
> > modeldat$MONTH  <- t1
> > modeldat$ACTIVE <- t2
> > modeldats <- na.omit(modeldat)
> > mod1 <-  coxph(Surv(MONTH,ACTIVE) ~  ALTER+RISIKO,data=modeldats)
> > pdf(file = "c:/Survivalx.pdf",    width = 6, height = 6, onefile = TRUE,
> > family = "Helvetica",title = "R Graphics Output")
> > plot(survfit(mod1),ylim=c(.7,1),xlab='Month',ylab='Proportion not
> > Active') const=const+1
> > if(const > 5) break
> > }
> > dev.off()
> > }
> >
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide!
> > http://www.R-project.org/posting-guide.html