# NOT RUN {
# Load the florentine marriage network
data(florentine)
# Create missing data
set.seed(22101992)
missNode <- sample(1:16,1)
flomarriage[missNode,] <- NA
flomarriage[,missNode] <- NA
# Posterior parameter estimation:
m.flo <- missBergm(flomarriage ~ edges + kstar(2),
burn.in = 50,
aux.iters = 500,
main.iters = 500,
gamma = 1,
nImp = 5)
# Posterior summaries:
bergm.output(m.flo)
# Bayesian goodness-of-fit test:
bgof(m.flo,
aux.iters = 500,
sample.size = 50,
n.deg = 10,
n.dist = 9,
n.esp = 6)
# }
Run the code above in your browser using DataLab