# NOT RUN {
# Load the florentine marriage network
data(florentine)
# Calibrated pseudo-posterior:
cpp.flo <- calibrate.bergm(flomarriage ~ edges + kstar(2),
aux.iters = 3000,
mcmc = 10000,
burnin = 500,
tunePL = 2.5)
# MCMC diagnostics and posterior summaries:
bergm.output(cpp.flo)
# Bayesian goodness-of-fit test:
bgof(cpp.flo,
aux.iters = 500,
sample.size = 50,
n.deg = 10,
n.dist = 9,
n.esp = 6)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab