networktree (version 0.2.1)

getnetwork: getnetwork

Description

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

Usage

getnetwork(tree, id = 1L, type = "detect", ...)

Arguments

tree

a networktree object

id

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

type

"cor", "pcor", or "glasso". Defaults to automatic detection

...

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

Examples

# NOT RUN {
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)

# }