donaldRwilliams

# estimate covariance matrix with regression
library(GGMnonreg)

p = 20

pcors = GGMnonreg::gen_net(p = p)

y = MASS::mvrnorm(n = 500, 
                  mu = rep(0, p), 
                  Sigma = pcors$cors)