library(qgraph)
library(bootnet)

set.seed(1337)

x1<-rnorm(200)
x2<-rnorm(200)+0.5*x1
x3<-rnorm(200)+0.3*x2

set.seed(1337)

y1<-rnorm(200)
y2<-rnorm(200)+0.9*y1
y3<-rnorm(200)+0.1*y2

data1<-as.data.frame(cbind(x1,x2,x3))
data2<-as.data.frame(cbind(y1,y2,y3))

g1 <- estimateNetwork(data1, default="EBICglasso")
g2 <- estimateNetwork(data2, default="EBICglasso")

plot(g1, layout="spring", cut=0)
plot(g2, layout="spring", cut=0)

g1mat <- getWmat(g1)
g2mat <- getWmat(g2)

L <- averageLayout(g1mat, g2mat)

plot(g1, layout=L, cut=0)
plot(g2, layout=L, cut=0)