> #3.5 Tomato Example
> 
> Nij<-rbind(
+ c( 4, 8, 9, 3),
+ c( 5, 8, 6, 5),
+ c( 4, 5, 7, 8),
+ c(11, 3, 2, 8))
> rownames(Nij)<-c("Floradade","Momotaro","Summit","Rutgers")
> colnames(Nij)<-1:4
> Nij
           1 2 3 4
Floradade  4 8 9 3
Momotaro   5 8 6 5
Summit     4 5 7 8
Rutgers   11 3 2 8
> RBD(Nij)
$Cri
           
              Location    Dispersion      Skewness
  Floradade -0.1581139 -1.767767e+00 -3.162278e-01
  Momotaro  -0.1581139 -7.071068e-01  4.743416e-01
  Summit     1.1067972 -1.570092e-16 -1.581139e-01
  Rutgers   -0.7905694  2.474874e+00  1.413083e-15

$partition
           df    SS     pvalue
Location    3  1.90 0.59341918
Dispersion  3  9.75 0.02081525
Residual    3  0.35 0.95036612
Total       9 12.00 0.21330931

>