if (FALSE) {
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