> #1.3 Cigarette Example
> 
> Nij<-rbind(
+ c( 6, 1, 6, 4, 0, 1, 0, 2),
+ c( 8, 2, 0, 1, 4, 1, 3, 1),
+ c( 2, 7, 0, 1, 2, 3, 3, 2),
+ c( 0, 3, 4, 5, 3, 4, 0, 1),
+ c( 2, 3, 2, 5, 3, 5, 0, 0),
+ c( 0, 1, 4, 3, 7, 1, 4, 0),
+ c( 0, 0, 0, 1, 0, 4, 3,12),
+ c( 2, 3, 4, 0, 1, 1, 7, 2))
> rownames(Nij)<-c("A","B","C","D","E","F","G","H")
> colnames(Nij)<-1:8
> Nij
  1 2 3 4 5 6 7  8
A 6 1 6 4 0 1 0  2
B 8 2 0 1 4 1 3  1
C 2 7 0 1 2 3 3  2
D 0 3 4 5 3 4 0  1
E 2 3 2 5 3 5 0  0
F 0 1 4 3 7 1 4  0
G 0 0 0 1 0 4 3 12
H 2 3 4 0 1 1 7  2
> RBD(Nij) #gives default level of detail
$Cri
   
      Location Dispersion    Skewness   Kurtosis 5-th order 6-th order
  A -2.3734644  0.7302967  0.87386290  1.3825208 -1.8988863  1.6749039
  B -1.7344548  1.8257419 -2.54876678  0.9534626 -0.1265924 -1.1651505
  C -0.4564355  0.6390097 -0.14564382 -2.0022714  2.0001602 -1.2379724
  D -0.4564355 -2.4647515  1.01950671  0.3813850  0.5823251  0.6553972
  E -1.0041580 -2.0083160 -0.50975336  0.6197507  0.9621024  0.3641095
  F  0.4564355 -2.7386128 -0.21846572  0.2383656 -0.7595545 -1.0923286
  G  5.0207901  3.1950483  1.60208198  1.0011357  0.8608285  1.0195067
  H  0.5477226  0.8215838 -0.07282191 -2.5743490 -1.6203830 -0.2184657
   
    7-th order
  A  0.3837461
  B -0.9896610
  C  0.5655206
  D  0.9290695
  E  1.5349844
  F -2.2620823
  G  1.1108440
  H -1.2724213

$partition
           df        SS       pvalue
Location    7  35.78333 7.964563e-06
Dispersion  7  32.76667 2.926183e-05
Skewness    7  11.20000 1.301299e-01
Kurtosis    7  15.04545 3.542114e-02
Residual   21  34.00455 3.619873e-02
Total      49 128.80000 4.279834e-09

> RBD(Nij,SSeffects.max=2) #to give same level of detail as in book specify the maximum number of effects to display
$Cri
   
      Location Dispersion    Skewness   Kurtosis 5-th order 6-th order
  A -2.3734644  0.7302967  0.87386290  1.3825208 -1.8988863  1.6749039
  B -1.7344548  1.8257419 -2.54876678  0.9534626 -0.1265924 -1.1651505
  C -0.4564355  0.6390097 -0.14564382 -2.0022714  2.0001602 -1.2379724
  D -0.4564355 -2.4647515  1.01950671  0.3813850  0.5823251  0.6553972
  E -1.0041580 -2.0083160 -0.50975336  0.6197507  0.9621024  0.3641095
  F  0.4564355 -2.7386128 -0.21846572  0.2383656 -0.7595545 -1.0923286
  G  5.0207901  3.1950483  1.60208198  1.0011357  0.8608285  1.0195067
  H  0.5477226  0.8215838 -0.07282191 -2.5743490 -1.6203830 -0.2184657
   
    7-th order
  A  0.3837461
  B -0.9896610
  C  0.5655206
  D  0.9290695
  E  1.5349844
  F -2.2620823
  G  1.1108440
  H -1.2724213

$partition
           df        SS       pvalue
Location    7  35.78333 7.964563e-06
Dispersion  7  32.76667 2.926183e-05
Residual   35  60.25000 5.030076e-03
Total      49 128.80000 4.279834e-09

>