Generate the posterior predictive distribution of counts for every pairwise interaction.
Examples
if (FALSE) { # interactive()
data(web)
dt <- prepare_data(mat = web, sampl.eff = rep(20, nrow(web)))
fit <- fit_model(dt, refresh = 0)
predict_counts(fit, dt)
}