Can somebody give a hint on how to speed-up the following loop:


for(j in 0:KM1)
  {
   k=j*60
   for(i in 1:60)
   {
  dat$yvac[k+i]= rbinom(1,dat$nvac[k+i],dat$p.trt[j+i])
  }
 }

K1=999



-- 
-Tony

	[[alternative HTML version deleted]]


