# create the data: two binary relations among three elements
arr <- round( replace( array( runif(18), c(3 ,3, 2) ), array( runif(18),
c(3, 3, 2) ) > .5, 3 ) )
# remove isolates (if exist)
rm.isol(arr)
Run the code above in your browser using DataLab