# NOT RUN {
data(marty)
##Clustering on 50 most variant genes amongst 500 first
mv.genes<-genes.selection(marty[1:500,], thres.num=50)
c.sample<-clustering(marty[mv.genes,], metric="pearson", metho="ward")
clustering.plot(c.sample, lab=marty.type.cl, title="H.Clustering\nPearson-Ward")
c.gene<-clustering(data=t(marty[mv.genes,]), metric="pearson",method="ward")
##Two-ways clustering
clustering.plot(tree=c.sample, tree.sup=c.gene, data=marty[mv.genes,], trim.heatmap=0.99)
# }
Run the code above in your browser using DataLab