# Make sure the observed network does not have isolated nodes
set.seed(27)
y <- network(7, directed = FALSE)
post.est <- bbeta(y ~ sociality(base = 0),
main.iters = 3000,
sigma.epsilon = diag(1.3, dim(y[,])[1]))
bbeta.output(post.est, plot = TRUE)
Run the code above in your browser using DataLab