Y1 <- MASS::mvrnorm(500, rep(0, 16), Sigma = BGGM::ptsd_cor1, empirical = T)
Y2 <- MASS::mvrnorm(500, rep(0, 16), Sigma = BGGM::ptsd_cor2, empirical = T)


ggm_bf <- ggm_compare_bf(Y1, Y2,
                         prior_sd = .35,
                         iter = 5000,
                         cores = 2)

# select with BF_cut = 3
ggm_bf_sel <- select(ggm_bf, BF_cut = 3)

# summary
summary(ggm_bf_sel)

# matrix for which there is evidence for a difference
ggm_bf_sel$pcor_mat_10