> x<-rbind( + c(2.5,1,2.5), + c(1.5,3,1.5), + c(2.5,1,2.5), + c(1,3,2), + c(1.5,3,1.5), + c(1,2.5,2.5), + c(1,3,2), + c(1.5,3,1.5)) > colnames(x)<-c("A","B","C") > > #Have to get ranks from x*10 as xtotab assumes integer > xtotab(x*10) [,1] [,2] [,3] [1,] 4.5 2.5 1.0 [2,] 2.0 0.5 5.5 [3,] 1.5 5.0 1.5 > U(x*10) [,1] [,2] [,3] [1,] 6.5 1.5 0.0 [2,] 1.5 5.0 1.5 [3,] 0.0 1.5 6.5 > BlockAnalysis(x*10) $Cri Location Dispersion A -1.372813e+00 0.1543033 B 1.372813e+00 2.0059435 C 4.050007e-17 -2.1602469 $partition df SS pvalue Location 2 3.769231e+00 0.15188747 Dispersion 2 8.714286e+00 0.01281495 Residual 0 -2.700062e-13 NA Total 4 1.248352e+01 0.01409558 >