## Create some "high-dimensional" data
set.seed(1)
p <- 10
ns <- c(5, 6)
Slist <- createS(ns, p)
## Obtain sparsified partial correlation matrix
Plist <- ridgeP.fused(Slist, ns, lambda = c(5.2, 1.3), verbose = FALSE)
PCsparse <- sparsify.fused(Plist , threshold = "absValue", absValueCut = 0.2)
SPlist <- lapply(PCsparse, "[[", "sparsePrecision") # Get sparse precisions
## Calculate GGM network statistics in each class
if (FALSE) GGMnetworkStats.fused(SPlist)
Run the code above in your browser using DataLab