Fused gaussian graphical model node pair path statistics
A simple wrapper for GGMpathStats.
GGMpathStats.fused(sparsePlist, ...)
sparsePlist |
A |
... |
Arguments passed to |
A list of path stats.
The function currently fails if no paths are present in one of the groups.
Anders E. Bilgrau, Carel F.W. Peeters <cf.peeters@vumc.nl>, Wessel N. van Wieringen
## Obtain some (high-dimensional) data set.seed(1) ns <- c(10, 11) Slist <- createS(ns, p = 7, topology = "banded") Tlist <- default.target.fused(Slist, ns) ## Obtain regularized precision and sparsify Plist <- ridgeP.fused(Slist, ns, Tlist, lambda = c(1, 1.6)) sparsePlist <- sparsify.fused(Plist, threshold = "absValue", absValueCut = 0.20) SPlist <- lapply(sparsePlist, "[[", "sparsePrecision") ## Obtain information on mediating and moderating paths between nodes 14 and 23 res <- GGMpathStats.fused(SPlist, node1 = 3, node2 = 4, graph = FALSE)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.