How to cluster in heatmap.2

Question

[Björn]

How to cluster heatmap using different distance matrix (Manhattan or Euclidean) and split column or row-wise in following command

heatmap.2(lcpm, Colv=FALSE, dendrogram="row",labRow = y2$genes$genes, Rowv = TRUE,  col=rev(morecols(20)),trace="none",tracecol = "black",main="Top 20 variable genes across samples",cex.main=0.6, ColSideColors=group.col,scale="row",margins=c(10,9))

Thanks

Answer 1

[Kevin Blighe]

I give an answer here, that indirectly answers your question: A: Heatmap based with FPKM values

In a nutshell, just add the following as parameters to heatmap.2():

#Euclidean distance with Ward's linkage
heatmap.2(
  ...,
  distfun = function(x) dist(x, method="euclidean"),
  hclustfun = function(x) hclust(x, method="ward.D2"))

#1 minus Pearson correlation distance with average linkage
heatmap.2(
  ...,
  distfun = function(x) as.dist(1-cor(t(x))),
  hclustfun = function(x) hclust(x, method="average"))

et cetera

[Björn]

@Kevin Blighe, Thanks, It worked