getnetwork {networktree}

getnetwork

Description

Easily extract a network from one of the nodes in a networktree object

Usage

getnetwork(tree, id = 1L, transform = "detect", verbose = FALSE, ...)

Arguments

tree

a networktree object

id

the node in the tree to extract. Use summary(tree) to see id numbers for each split

transform

should stored correlation matrices be transformed to partial correlations or graphical lasso? Can be set to "cor", "pcor", or "glasso". Defaults to automatic detection

verbose

should warnings and messages from transformation functions (qgraph) be printed?

...

arguments passed to qgraph (e.g., "tuning", "threshold")

Examples

set.seed(1)
d <- data.frame(trend = 1:200, foo = runif(200, -1, 1))
d <- cbind(d, rbind(
  mvtnorm::rmvnorm(100, mean = c(0, 0, 0),
          sigma = matrix(c(1, 0.5, 0.5, 0.5, 1, 0.5, 0.5, 0.5, 1), ncol = 3)),
  mvtnorm::rmvnorm(100, mean = c(0, 0, 0),
          sigma = matrix(c(1, 0, 0.5, 0, 1, 0.5, 0.5, 0.5, 1), ncol = 3))
))
colnames(d)[3:5] <- paste0("y", 1:3)

## Now use the function
tree1 <- networktree(nodevars=d[,3:5], splitvars=d[,1:2])

getnetwork(tree1, id=1)