p=30
n=30
# simulate graph
eta=0.11
Gr <- simulateGraph(p,eta)
# simulate data
X <- rmvnorm(n, mean=rep(0,p), sigma=Gr$C)
# estimate graph
if (FALSE) GRest <- selectFast(X)
# plot result
if (FALSE) library(network)
if (FALSE) par(mfrow=c(1,2))
if (FALSE) gV <- network(Gr$G)
if (FALSE) plot(gV,jitter=TRUE, usearrows = FALSE, label=1:p,displaylabels=TRUE)
if (FALSE) g <- network(GRest$EW$G)
if (FALSE) plot(g, jitter=TRUE, usearrows = FALSE, label=1:p,displaylabels=TRUE)
Run the code above in your browser using DataLab