require(graph)
n.var <- 5 # number of variables
set.seed(123)
g <- randomEGraph(as.character(1:n.var), p=0.15)
Sigma <- qpG2Sigma(g, rho=0.5)
K <- solve(Sigma)
round(qpK2ParCor(K), digits=2)
as(g, "matrix")
Run the code above in your browser using DataLab