> #6.8 Page Analysis of Ordinal Trinary Data
> 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")
> x
       X Y   Z
[1,] 2.5 1 2.5
[2,] 1.5 3 1.5
[3,] 2.5 1 2.5
[4,] 1.5 3 1.5
[5,] 1.5 3 1.5
[6,] 1.0 3 2.0
[7,] 1.0 3 2.0
[8,] 1.5 3 1.5
> 
> xtotab3(x)
  1 1.5 2 2.5 3
X 2   4 0   2 0
Y 2   0 0   0 6
Z 0   4 2   2 0
> OrderedBlockAnalysis(x)
$Z
           [,1]      [,2]       [,3]       [,4]
[1,]  0.5547002 -1.186436  0.2139399  1.1130624
[2,] -1.9215378 -2.952539 -0.5601411 -0.5257856

$partition
         df      Value      pvalue
Page      1  0.3076923 0.579099742
Umbrella  1  3.6923077 0.054663936
Residual  2 12.0000000 0.002478752
Total     4 16.0000000 0.003019164

> 
> #Note see example 6.6 for CMH "Correlation Statistic" output
>