GGMpathStats.fused {rags2ridges}

Fused gaussian graphical model node pair path statistics

Description

A simple wrapper for GGMpathStats.

Usage

GGMpathStats.fused(sparsePlist, ...)

Arguments

sparsePlist

A list of sparsified precision matrices.

...

Arguments passed to GGMpathStats.

Value

A list of path stats.

Note

The function currently fails if no paths are present in one of the groups.

Author(s)

Anders E. Bilgrau, Carel F.W. Peeters <[e-mail]>, Wessel N. van Wieringen

See Also

GGMpathStats

Examples

## 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)