You are right to be dissatisfied with your code.

I suspect you will get more response if you say
what the code is supposed to do, and give a
smaller example of what the answer should be.

Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")

Sofia Martinez wrote:
HI:
I need ideas on how to make this code shorter (maybe with a second loop?).
The code as it is works, but in this case I only have 14 samples, but it
will become insane with more, so I need a way to make it more automatic. The
problem is that the output from ts1, ts2, and so on is a vector with more
than one value, so I do not know how to solve this.
   Thanks
      Prenewbie

The code is the following:

duration<-c(750, 8940,17180, 8693,10100, 7990,24820, 6770,
7390,8450,18400,16252,6080,11030)

tmax<-0
dt<-0
for (m in 1:14)
 {
  tmax[m]<-duration[m]
   dtmin<-969
   dtmax<-9884
   i<-1
   dt[m]<-round(runif(1,dtmin,(tmax[m]-1)))}
      ts<-c(1)
      ts1<-ts
      ts2<-ts1
      ts3<-ts1
      ts4<-ts1
      ts5<-ts1
      ts6<-ts1
      ts7<-ts1
      ts8<-ts1
      ts9<-ts1
      ts10<-ts1
      ts11<-ts1
      ts12<-ts1
      ts13<-ts1
      ts14<-ts1
      for (m in 1:2){
      while(i+dt[m]<tmax[m])
      {
      ts1<-append(ts,i+dt[1])
      i<-i+dt[1]
      ifelse(dtmin<=(tmax[1]-i), dt[1]<-round(runif(1,dtmin,(tmax[1]-i))),
dt[1]<-(tmax[1]-i))
      ts2<-append(ts,i+dt[2])
      i<-i+dt[2]
      ifelse(dtmin<=(tmax[2]-i), dt[2]<-round(runif(1,dtmin,(tmax[2]-i))),
dt[2]<-(tmax[2]-i))
       ts3<-append(ts,i+dt[3])
      i<-i+dt[3]
      ifelse(dtmin<=(tmax[3]-i), dt[3]<-round(runif(1,dtmin,(tmax[3]-i))),
dt[3]<-(tmax[3]-i))
       ts4<-append(ts,i+dt[4])
      i<-i+dt[4]
      ifelse(dtmin<=(tmax[4]-i), dt[4]<-round(runif(1,dtmin,(tmax[4]-i))),
dt[4]<-(tmax[4]-i))
       ts5<-append(ts,i+dt[5])
      i<-i+dt[5]
      ifelse(dtmin<=(tmax[5]-i), dt[5]<-round(runif(1,dtmin,(tmax[5]-i))),
dt[5]<-(tmax[5]-i))
       ts6<-append(ts,i+dt[6])
      i<-i+dt[6]
      ifelse(dtmin<=(tmax[6]-i), dt[6]<-round(runif(1,dtmin,(tmax[6]-i))),
dt[6]<-(tmax[6]-i))
       ts7<-append(ts,i+dt[7])
      i<-i+dt[7]
      ifelse(dtmin<=(tmax[7]-i), dt[7]<-round(runif(1,dtmin,(tmax[7]-i))),
dt[7]<-(tmax[7]-i))
       ts8<-append(ts,i+dt[8])
      i<-i+dt[8]
      ifelse(dtmin<=(tmax[8]-i), dt[8]<-round(runif(1,dtmin,(tmax[8]-i))),
dt[8]<-(tmax[8]-i))
       ts9<-append(ts,i+dt[9])
      i<-i+dt[9]
      ifelse(dtmin<=(tmax[9]-i), dt[9]<-round(runif(1,dtmin,(tmax[9]-i))),
dt[9]<-(tmax[9]-i))
       ts10<-append(ts,i+dt[10])
      i<-i+dt[10]
      ifelse(dtmin<=(tmax[10]-i),
dt[10]<-round(runif(1,dtmin,(tmax[10]-i))), dt[10]<-(tmax[10]-i))
       ts11<-append(ts,i+dt[11])
      i<-i+dt[11]
      ifelse(dtmin<=(tmax[11]-i),
dt[11]<-round(runif(1,dtmin,(tmax[11]-i))), dt[11]<-(tmax[11]-i))
       ts12<-append(ts,i+dt[12])
      i<-i+dt[12]
      ifelse(dtmin<=(tmax[12]-i),
dt[12]<-round(runif(1,dtmin,(tmax[12]-i))), dt[12]<-(tmax[12]-i))
       ts13<-append(ts,i+dt[13])
      i<-i+dt[13]
      ifelse(dtmin<=(tmax[13]-i),
dt[13]<-round(runif(1,dtmin,(tmax[13]-i))), dt[13]<-(tmax[13]-i))
       ts14<-append(ts,i+dt[14])
      i<-i+dt[14]
      ifelse(dtmin<=(tmax[14]-i),
dt[14]<-round(runif(1,dtmin,(tmax[14]-i))), dt[14]<-(tmax[14]-i))
      } }
     ts1<-append(ts,i+dt[1])
     ts2<-append(ts,i+dt[2])
     ts3<-append(ts,i+dt[3])
     ts4<-append(ts,i+dt[4])
     ts5<-append(ts,i+dt[5])
     ts6<-append(ts,i+dt[6])
     ts7<-append(ts,i+dt[7])
     ts8<-append(ts,i+dt[8])
     ts9<-append(ts,i+dt[9])
     ts10<-append(ts,i+dt[10])
     ts11<-append(ts,i+dt[11])
     ts12<-append(ts,i+dt[12])
     ts13<-append(ts,i+dt[13])
     ts14<-append(ts,i+dt[14])

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.




______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to