> #4.2.2 Cleaning materials Example
> x<-matrix(c(
+ 1,2,3,4,NA,   NA,   NA,   NA,NA,NA,
+ 1,3,NA,NA,   2,4,   NA,   NA,NA,NA,
+ 3,NA,4,NA,NA,NA,2,1,NA,NA,
+ 4,NA,NA,3,NA,NA,NA,NA,1,2,
+ 2,NA,NA,NA,3,NA,1,NA,4,NA,
+ 3,NA,NA,NA,NA,4,NA,1,NA,2,
+ NA,2,1,NA,NA,4,NA,NA,3,NA,
+ NA,2,NA,4,NA,NA,3,NA,NA,1,
+ NA,1,NA,NA,3,NA,NA,4,NA,2,
+ NA,3,NA,NA,NA,NA,2,4,1,NA,
+ NA,NA,1,NA,4,NA,NA,NA,3,2,
+ NA,NA,1,NA,NA,2,4,NA,NA,3,
+ NA,NA,1,3,2,NA,NA,4,NA,NA,
+ NA,NA,NA,4,2,3,1,NA,NA,NA,
+ NA,NA,NA,4,NA,2,NA,3,1,NA,
+ 3,1,2,4,NA,NA,NA,NA,NA,NA,
+ 1,2,NA,NA,3,4,NA,NA,NA,NA,
+ 2,NA,4,NA,NA,NA,3,1,NA,NA,
+ 3,NA,NA,4,NA,NA,NA,NA,1,2,
+ 1,NA,NA,NA,3,NA,2,NA,4,NA,
+ 3,NA,NA,NA,NA,4,NA,2,NA,1,
+ NA,1,3,NA,NA,4,NA,NA,2,NA,
+ NA,1,NA,4,NA,NA,3,NA,NA,2,
+ NA,1,NA,NA,2,NA,NA,4,NA,3,
+ NA,1,NA,NA,NA,NA,4,3,2,NA,
+ NA,NA,1,NA,3,NA,NA,NA,2,4,
+ NA,NA,1,NA,NA,2,4,NA,NA,3,
+ NA,NA,1,3,2,NA,NA,4,NA,NA,
+ NA,NA,NA,4,2,3,1,NA,NA,NA,
+ NA,NA,NA,4,NA,1,NA,2,3,NA),ncol=10,byrow=T)
> colnames(x)<-c("A","B","C","D","E","F","G","H","I","J")
> 
> xtotab(x)
      [,1] [,2] [,3] [,4]
 [1,]    4    2    5    1
 [2,]    6    4    2    0
 [3,]    7    1    2    2
 [4,]    0    0    3    9
 [5,]    0    6    5    1
 [6,]    1    3    2    6
 [7,]    3    3    3    3
 [8,]    3    2    2    5
 [9,]    4    3    3    2
[10,]    2    6    3    1
> U(x)
     [,1] [,2] [,3] [,4]
[1,]   30    0    0    0
[2,]    0   30    0    0
[3,]    0    0   30    0
[4,]    0    0    0   30
> BlockAnalysis(x)
$Cri
    Location    Dispersion
A -0.7348469 -5.477226e-01
B -2.4494897  2.723330e-15
C -1.7146428  1.643168e+00
D  3.6742346  1.643168e+00
E  0.2449490 -2.738613e+00
F  1.7146428  5.477226e-01
G  0.0000000  0.000000e+00
H  0.7348469  1.095445e+00
I -0.7348469  7.497978e-16
J -0.7348469 -1.643168e+00

$partition
           df   SS      pvalue
Location    9 27.6 0.001112224
Dispersion  9 17.4 0.042808282
Residual    9  5.4 0.798139062
Total      27 50.4 0.004090349

>