# NOT RUN {
# group 1
Y1 <- MASS::mvrnorm(500, mu = rep(0,16),
Sigma = BGGM::ptsd_cor2,
empirical = FALSE)
# group 2
Y2 <- MASS::mvrnorm(500, mu = rep(0,16),
Sigma = BGGM::ptsd_cor2,
empirical = FALSE)
# fit model
fit <- ggm_compare_bf(Y1, Y2,
prior_sd = 0.20,
iter = 50,
cores = 2)
# select E
E <- select(fit, BF_cut = 3)
# plot E
plot(E)
# }
Run the code above in your browser using DataLab