# NOT RUN {
data("florentine_m")
P <- neighborhood_inclusion(florentine_m)
res <- exact_rank_prob(P)
mcmc <- mcmc_rank_prob(P,rp = vcount(g)^3)
# mean absolute error (expected ranks)
mean(abs(res$expected.rank-mcmc$expected.rank))
# }
Run the code above in your browser using DataLab