# NOT RUN {
p <- 20
n <- 500
true_net <- gen_net(p = p, edge_prob = 0.25)
y <- MASS::mvrnorm(n = n,
mu = rep(0, p),
Sigma = true_net$cors)
# default
fit_atan <- ggmncv(R = cor(y),
n = nrow(y),
penalty = "atan",
progress = FALSE)
# lasso
fit_l1 <- ggmncv(R = cor(y),
n = nrow(y),
penalty = "lasso",
progress = FALSE)
# atan scores
score_binary(estimate = true_net$adj,
true = fit_atan$adj,
model_name = "atan")
score_binary(estimate = fit_l1$adj,
true = true_net$adj,
model_name = "lasso")
# }
Run the code above in your browser using DataLab