library(simone)
data(cancer)
str(cancer, max.level=1)

attach(cancer)
boxplot(expr, las=3, cex.axis=0.6)
table(status)

out <- simone(expr, tasks=status)

plot(out)

glist <- getNetwork(out, "AIC")
plot(glist[[1]],glist[[2]])
glist <- getNetwork(out, "BIC")
plot(glist[[1]],glist[[2]])
glist <- getNetwork(out, 65)
plot(glist[[1]],glist[[2]])
plot(glist[[1]],glist[[2]], type="overlap")

detach(cancer)