For each dyad, compute the posterior probability that the latent theta is positive (or exceeds a threshold)
Examples
# \donttest{
sim <- simulate_dynamic_dbn(n = 6, time = 5, seed = 1)
fit <- dbn(sim$Y, model = "dynamic", nscan = 200, burn = 100, verbose = FALSE)
ep <- edge_prob(fit)
# }