Hi;
Is there any faster way to run this code :
## b0, beta.r  and stud  are already defined.
 
  Ab<- matrix(numeric(1),nrow=2,ncol=1000)
    Bb<- matrix(numeric(1),nrow=2,ncol=1000)
  for( j in 1:1000) { 
 
   for(i in 1:1000) {
  
      eq<- sample(stud,replace=T)
      ystar<-  beta.r*x+eq+b0    
      
      Ab[,i]<- wwfit(x,ystar)$coef 
                     }
      Bb[1,j]<-sd(Ab[1,])    
      Bb[2,j]<-sd(Ab[2,])
        
                          }               
      se1<- mean(Bb[1,])      
      se2<- mean(Bb[2,])
      se1
      se2    

thanks alot for your help.
Rebhi Bsharat.
dept of stat
KSU.


__________________________________________________



	[[alternative HTML version deleted]]


