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

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

## no clustering by default
res.no <- simone(expr)
plot(res.no) ## "trop" de monde sur BIC / AIC
g.no <- getNetwork(res.no, 30)
plot(g.no)
plot(g.no, type = "cluster")

## try with clustering now
ctrl <- setOptions(clusters.crit=30)
res.cl <- simone(expr, clustering=TRUE, control=ctrl)
g.cl <- getNetwork(res.cl, 30)
plot(g.cl)
plot(g.cl, type = "circles")

## Let us compare the two networks
plot(g.no,g.cl)
plot(g.no,g.cl, type="overlap")