Easily extract a network from one of the nodes in a networktree object
getnetwork(tree, id = 1L, transform = "detect", verbose = FALSE, ...)a networktree object
the node in the tree to extract. Use summary(tree) to see id numbers for each split
should stored correlation matrices be transformed to partial correlations or graphical lasso? Can be set to "cor", "pcor", or "glasso". Defaults to automatic detection
should warnings and messages from transformation functions (qgraph) be printed?
arguments passed to qgraph (e.g., "tuning", "threshold")
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)
#> y1 y2 y3
#> y1 0.0000000 0.1421727 0.4671353
#> y2 0.1421727 0.0000000 0.4563241
#> y3 0.4671353 0.4563241 0.0000000