# NOT RUN {
# Quick toy example
library(MASS)
p <- 5 # number of variables
data <- mvrnorm(n=100, mu=rep(0, p), Sigma = diag(p))
out <- OLS(data = data)
# point estimates + CIs; show 3 largest effects only
plot(out, labels = colnames(ptsd_data),
order=TRUE, subset = 1:3)
# }
# NOT RUN {
# Fit GGM to empirical PTSD data
set.seed(1)
out <- lasso_dsp(data = ptsd_data)
# Plot first 20 edges
plot(out, labels = colnames(ptsd_data),
order=TRUE, subset = 1:20)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab