News

 

 

 

 

    perc.fun.all <- function(data,i) {
         m <- apply(data[i,],2,mean)
         if((m[1]<=m[2])&&(m[2]<=m[3])&&(m[3]<=m[4])) return(1)
         if((m[1]<=m[2])&&(m[2]<=m[4])&&(m[4]<=m[3])) return(2)
         if((m[1]<=m[3])&&(m[3]<=m[2])&&(m[2]<=m[4])) return(3)
         if((m[1]<=m[3])&&(m[3]<=m[4])&&(m[4]<=m[2])) return(4)
         if((m[1]<=m[4])&&(m[4]<=m[2])&&(m[2]<=m[3])) return(5)
         if((m[1]<=m[4])&&(m[4]<=m[2])&&(m[2]<=m[3])) return(6)
         if((m[2]<=m[1])&&(m[1]<=m[3])&&(m[3]<=m[4])) return(7)
         if((m[2]<=m[1])&&(m[1]<=m[4])&&(m[4]<=m[3])) return(8)
         if((m[2]<=m[3])&&(m[3]<=m[1])&&(m[1]<=m[4])) return(9)
         if((m[2]<=m[3])&&(m[3]<=m[4])&&(m[4]<=m[1])) return(10)
         if((m[2]<=m[4])&&(m[4]<=m[1])&&(m[1]<=m[3])) return(11)
         if((m[2]<=m[4])&&(m[4]<=m[3])&&(m[3]<=m[1])) return(12)
         if((m[3]<=m[1])&&(m[1]<=m[2])&&(m[2]<=m[4])) return(13)
         if((m[3]<=m[1])&&(m[1]<=m[4])&&(m[4]<=m[2])) return(14)
         if((m[3]<=m[2])&&(m[2]<=m[1])&&(m[1]<=m[4])) return(15)
         if((m[3]<=m[2])&&(m[2]<=m[4])&&(m[4]<=m[1])) return(16)
         if((m[3]<=m[4])&&(m[4]<=m[1])&&(m[1]<=m[2])) return(17)
         if((m[3]<=m[4])&&(m[4]<=m[2])&&(m[2]<=m[1])) return(18)
         if((m[4]<=m[1])&&(m[1]<=m[2])&&(m[2]<=m[3])) return(19)
         if((m[4]<=m[1])&&(m[1]<=m[3])&&(m[3]<=m[2])) return(20)
         if((m[4]<=m[2])&&(m[2]<=m[1])&&(m[1]<=m[3])) return(21)
         if((m[4]<=m[2])&&(m[2]<=m[3])&&(m[3]<=m[1])) return(22)
         if((m[4]<=m[3])&&(m[3]<=m[1])&&(m[1]<=m[2])) return(23)
         if((m[4]<=m[3])&&(m[3]<=m[2])&&(m[2]<=m[1])) return(24) }

           can be replaced by 

    perc.fun.all <- function(data,i) {
         m <- apply(data[i,],2,mean)
         b <- c(1000,100,10,1)
         return(t(b)%*%rank(m)) }

     The returned value then reflects the ranks of the elements of the mean vector.


Up Ordering the Book News References R Code Data Corrections About the Author