# NOT RUN {
data("Townships")
## original data
pimage(Townships)
criterion(Townships)
## seriated data
order <- seriate(Townships, method = "BEA", control = list(rep = 5))
pimage(Townships, order)
criterion(Townships, order)
# }
Run the code above in your browser using DataCamp Workspace