> #3.2 Peach Example
> x<-rbind(
+ c(1,2,3,4,5),
+ c(1,2,3,5,4),
+ c(1,4,5,2,3),
+ c(1,3,4,5,2),
+ c(1,2,3,5,4),
+ c(1,3,5,2,4),
+ c(1,5,4,3,2),
+ c(1,2,3,5,4),
+ c(1,2,4,3,5),
+ c(1,5,2,4,3))
> colnames(x)<-c("A","B","C","D","E")
> 
> xtotab(x)
     [,1] [,2] [,3] [,4] [,5]
[1,]   10    0    0    0    0
[2,]    0    5    2    1    2
[3,]    0    1    4    3    2
[4,]    0    2    2    2    4
[5,]    0    2    2    4    2
> U(x)
     [,1] [,2] [,3] [,4] [,5]
[1,]   10    0    0    0    0
[2,]    0   10    0    0    0
[3,]    0    0   10    0    0
[4,]    0    0    0   10    0
[5,]    0    0    0    0   10
> BlockAnalysis(x)
$Cri
       Location    Dispersion
A -4.000000e+00  3.380617e+00
B -3.185596e-17 -1.014185e+00
C  1.200000e+00 -1.352247e+00
D  1.600000e+00  2.252948e-15
E  1.200000e+00 -1.014185e+00

$partition
           df       SS       pvalue
Location    4 21.44000 2.589946e-04
Dispersion  4 15.31429 4.091850e-03
Residual    8 11.24571 1.881616e-01
Total      16 48.00000 4.749992e-05

>