> #4.2.3 Ice cream Flavour Example
> x<-matrix(c(
+ 2,1,NA,3,NA,NA,NA,
+ NA,2,3,NA,1,NA,NA,
+ NA,NA,1,2,NA,3,NA,
+ NA,NA,NA,1,3,NA,2,
+ 2,NA,NA,NA,1,3,NA,
+ NA,1,NA,NA,NA,2,3,
+ 1,NA,3,NA,NA,NA,2,
+ NA,NA,3,1,NA,2,NA,
+ NA,NA,NA,2,1,NA,3,
+ 2,NA,NA,NA,1,3,NA,
+ NA,1,NA,NA,NA,3,2,
+ 1,NA,2,NA,NA,NA,3,
+ 1,3,NA,2,NA,NA,NA,
+ NA,2,3,NA,1,NA,NA),ncol=7,byrow=T)
> colnames(x)<-c("A","B","C","D","E","F","G")
> 
> xtotab(x)
     [,1] [,2] [,3]
[1,]    3    3    0
[2,]    3    2    1
[3,]    1    1    4
[4,]    2    3    1
[5,]    5    0    1
[6,]    0    2    4
[7,]    0    3    3
> U(x)
     [,1] [,2] [,3]
[1,]   14    0    0
[2,]    0   14    0
[3,]    0    0   14
> BlockAnalysis(x)
$Cri
    Location    Dispersion
A -1.3887301 -8.017837e-01
B -0.9258201 -8.533101e-16
C  1.3887301  8.017837e-01
D -0.4629100 -8.017837e-01
E -1.8516402  1.603567e+00
F  1.8516402  1.706620e-15
G  1.3887301 -8.017837e-01

$partition
           df            SS     pvalue
Location    6  1.371429e+01 0.03299579
Dispersion  6  5.142857e+00 0.52562569
Residual    0 -1.554312e-13         NA
Total      12  1.885714e+01 0.09203482

>