[R] R question......

From: rebhi bsharat <rebhi2001_at_yahoo.com>
Date: Sat 31 Jul 2004 - 06:29:51 EST

 

Hi;
Is there a faster way to run this code?  

    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 A LOT.                 


        [[alternative HTML version deleted]]



R-help@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 Received on Sat Jul 31 06:35:04 2004

This archive was generated by hypermail 2.1.8 : Wed 03 Nov 2004 - 22:55:23 EST