if (FALSE) {
# Load the florentine marriage network
data(florentine)
# Posterior parameter estimation:
p.flo <- bergm(flomarriage ~ edges + kstar(2),
burn.in = 50,
aux.iters = 500,
main.iters = 1000,
gamma = 1.2)
# MCMC diagnostics plots:
plot(p.flo)
}
Run the code above in your browser using DataLab