# get golub data
library(vsn)
library(golubEsets)
data(Golub_Merge)
# use 10% most variable genes
e <- exprs(Golub_Merge)
vars <- apply(e, 1, var)
e <- e[vars > quantile(vars,0.9),]
# use diana2means to get splits and scores
diana2means(e)
diana2means(e, return.cut=TRUE)
Run the code above in your browser using DataLab