View source: R/rags2ridgesFused.R
sparsify.fused | R Documentation |
A simple wrapper for sparsify
which determines the support of
a list
of partial correlation/precision matrix by various methods and
returns the sparsified matrices.
sparsify.fused(Plist, ...)
Plist |
A |
... |
Arguments passed to |
A list
of the same length as Plist
with the output from
sparsify
.
Anders Ellern Bilgrau, Wessel N. van Wierigen, Carel F.W. Peeters <carel.peeters@wur.nl>
sparsify
ns <- c(10, 11)
Ylist <- createS(ns, p = 16, dataset = TRUE)
Slist <- lapply(Ylist, covML)
Tlist <- default.target.fused(Slist, ns)
# Obtain regularized precision under optimal penalty
opt <- optPenalty.fused.auto(Ylist, Tlist, cv.method = "aLOOCV",
maxit.ridgeP.fused = 1500)
# Use the optimal penalties
Plist <- ridgeP.fused(Slist, ns, lambda = opt$lambda, maxit = 1000)
# Determine support regularized (standardized) precision under optimal penalty
res <- sparsify.fused(Plist, threshold = "top", verbose = FALSE)
round(res[[1]]$sparsePrecision, 1)
round(res[[2]]$sparsePrecision, 1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.