Hi,

I was trying to draw a heatmap of the bicluster results. With the code given in the biclust package I can only get the heatmap for one cluster at a time (drawHeatmap function). Is there any way that I can get the heatmap for all the clusters at the same time?
The code that I am using (biclust documentation) is:
 #Random 100x50 matrix with a single, up-regulated 10x10 bicluster
       s2=matrix(rnorm(5000),100,50)
       s2[11:20,11:20]=rnorm(100,3,0.3)  # cluster 1
       s2[31:40,31:40]=rnorm(100,5,0.7) # cluster 2 
       set.seed(1)
       bics <- biclust(s2,BCPlaid(), back.fit = 2, shuffle = 3, fit.model = ~m + a + b,
       iter.startup = 5, iter.layer = 30,  verbose = TRUE)
       drawHeatmap(s2,bics,1) # heatmap for cluster 1
many thanks


      
	[[alternative HTML version deleted]]

