#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 RBD(Nij) #gives default level of detail RBD(Nij,SSeffects.max=2) #to give same level of detail as in book specify the maximum number of effects to display