SHARE
TWEET

Untitled

a guest Sep 8th, 2019 101 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. Y1 <- MASS::mvrnorm(500, rep(0, 16), Sigma = BGGM::ptsd_cor1, empirical = T)
  2. Y2 <- MASS::mvrnorm(500, rep(0, 16), Sigma = BGGM::ptsd_cor2, empirical = T)
  3.  
  4.  
  5. ggm_bf <- ggm_compare_bf(Y1, Y2,
  6.                          prior_sd = .35,
  7.                          iter = 5000,
  8.                          cores = 2)
  9.  
  10. # select with BF_cut = 3
  11. ggm_bf_sel <- select(ggm_bf, BF_cut = 3)
  12.  
  13. # summary
  14. summary(ggm_bf_sel)
  15.  
  16. # matrix for which there is evidence for a difference
  17. ggm_bf_sel$pcor_mat_10
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top