#3.2 Peach Example

Nij<-rbind(
c(10, 0, 0, 0, 0),
c( 0, 5, 2, 1, 2),
c( 0, 1, 4, 3, 2),
c( 0, 2, 2, 2, 4),
c( 0, 2, 2, 4, 2))
rownames(Nij)<-c("A","B","C","D","E")
colnames(Nij)<-1:5
Nij
RBD(Nij)
RBD(Nij,SSeffects.max=2) #to give same level of detail as in book specify the maximum number of effects to display