> #6.3 Ordinal Data Example
> x<-rbind(
+ c(2.5,1,2.5),
+ c(1.5,3,1.5),
+ c(2.5,1,2.5),
+ c(1.5,3,1.5),
+ c(1.5,3,1.5),
+ c(1,3,2),
+ c(1,3,2),
+ c(1.5,3,1.5))
> colnames(x)<-c("X","Y","Z")
> 
> #Have to get ranks from x*10 as xtotab assumes integer
> xtotab(x*10)
  1 2 3
X 4 3 1
Y 2 0 6
Z 2 5 1
> U(x*10)
  1 2 3
1 6 2 0
2 2 5 1
3 0 1 7
> BlockAnalysis(x*10)
$Cri
    Location  Dispersion
X -1.1766968 -0.09656091
Y  1.5689291  2.22090093
Z -0.3922323 -2.12434002

$partition
           df           SS      pvalue
Location    2 4.000000e+00 0.135335283
Dispersion  2 9.454545e+00 0.008850576
Residual    0 5.861978e-13          NA
Total       4 1.345455e+01 0.009255690

>