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]]
