Skip to content

Instantly share code, notes, and snippets.

@SachaEpskamp
Created June 13, 2019 14:50
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save SachaEpskamp/dbc0f5b330bd7452f0f23bbf29249915 to your computer and use it in GitHub Desktop.
Save SachaEpskamp/dbc0f5b330bd7452f0f23bbf29249915 to your computer and use it in GitHub Desktop.
library("qgraph")
library("bootnet")
set.seed(1)
# Number of nodes:
nNode <- 10
# random point estimates:
pointests <- runif(nNode,0,1)
# Add a CI Range:
error <- runif(nNode,0.05,0.2)
lowerbound <- pmax(0, pointests - error)
upperbound <- pmin(1, pointests + error)
# Some network:
net <- genGGM(nNode, p = 0.4, nei = 2)
# Plot:
qgraph(net, layout = "spring",
vsize = 10,
pieCImid = pointests,
pieCIlower = lowerbound,
pieCIupper = upperbound,
pieBorder = 0.3,
groups = sample(LETTERS[1:3],nNode,TRUE),
legend=FALSE,
theme = "colorblind"
)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment