> #4.6.1 Ice Cream (Off Flavour) Example 
> x<-matrix(c(
+ 6,1,1,2,NA,NA,
+ 6,NA,NA,1,3,3,
+ NA,4,2,NA,5,2,
+ 7,2,3,NA,2,NA,
+ 3,5,NA,1,NA,1,
+ NA,NA,1,1,3,2,
+ 7,4,4,NA,NA,3,
+ 2,NA,1,1,1,NA,
+ NA,2,NA,2,2,3,
+ 4,2,NA,2,5,NA,
+ 5,NA,3,NA,1,1,
+ NA,3,2,1,NA,2,
+ 4,2,NA,NA,1,1,
+ 5,NA,2,2,NA,1,
+ NA,2,4,5,3,NA),ncol=6,byrow=T)
> colnames(x)<-c("A","B","C","D","E","F")
> 
> xtotab(x)
         [,1]     [,2]     [,3] [,4]
[1,] 0.000000 0.000000 2.000000    8
[2,] 2.833333 2.333333 2.833333    2
[3,] 1.833333 3.333333 4.833333    0
[4,] 4.166667 2.666667 2.166667    1
[5,] 2.166667 3.666667 1.166667    3
[6,] 4.000000 3.000000 2.000000    1
> U(x)
          [,1]     [,2]       [,3] [,4]
[1,] 9.6666667 4.666667  0.6666667    0
[2,] 4.6666667 7.666667  2.6666667    0
[3,] 0.6666667 2.666667 11.6666667    0
[4,] 0.0000000 0.000000  0.0000000   15
> BlockAnalysis(x)
$Cri
       Location  Dispersion
A  3.605551e+00  1.33330989
B -2.773501e-01 -0.05507975
C -5.547002e-01 -2.20888803
D -1.386750e+00  0.46532895
E  8.942271e-17  0.12345462
F -1.386750e+00  0.34187433

$partition
           df        SS      pvalue
Location    5 17.230769 0.004082306
Dispersion  5  7.008585 0.220002577
Residual    5  4.295180 0.507745895
Total      15 28.534535 0.018452738

>