# Generating multivariate normal data from a 'random' graph
data.sim <- bdgraph.sim( n = 80, p = 6, size = 6, vis = TRUE )
# Running algorithm based on GGMs
ssgraph.obj <- ssgraph( data = data.sim, iter = 1000 )
summary( ssgraph.obj )
# To compare the result with true graph
compare( pred = ssgraph.obj, actual = data.sim,
main = c( "Target", "ssgraph" ), vis = TRUE )
plotroc( pred = ssgraph.obj, actual = data.sim )
if (FALSE) {
# Running algorithm with starting points from previous run
ssgraph.obj2 <- ssgraph( data = data.sim, iter=5000, g.start = ssgraph.obj )
compare( pred = list( ssgraph.obj, ssgraph.obj2 ), actual = data.sim,
main = c( "Target", "Frist run", "Second run" ), vis = TRUE )
plotroc( pred = list ( ssgraph.obj, ssgraph.obj2 ), actual = data.sim,
label = c( "Frist run", "Second run" ) )
}
Run the code above in your browser using DataLab