Question: Different results with heatmap.2 and pheatmap ?
0
gravatar for Björn
4 months ago by
Björn30
Björn30 wrote:

I got two completely different results (gene numbers) using heatmap.2 and pheatmap commands while doing DE analysis in edgeR

heatmap.2 command:

heatmap.2(highly_variable_lcpm,labRow = rownames(y$genes),Rowv = TRUE,  col=rev(morecols(20)),trace="column",tracecol = "black",main="Top 20 most variable genes across samples",ColSideColors=group.col,scale="row",margins=c(10,5))

Graph:![heatmap_graph][1] [1]: https://ibb.co/h8BqoH

pheatmap command:

pheatmap(highly_variable_lcpm,labRow = rownames(y$genes),Rowv = TRUE,  col=rev(morecols(20)),trace="column",tracecol = "black",main="Top 20 most variable genes across samples",ColSideColors=group.col,scale="row",margins=c(10,5), cutree_rows = 4)

Graph![enter image description here][1] [1]: https://ibb.co/gTgjFx

Additionally, How to convert the row-names to real gene names. The "labRow" don't seem to work ! Thanks Björn

heatmap.2 edger pheatmap • 332 views
ADD COMMENTlink written 4 months ago by Björn30

Can you show us the results of head(y) and head(highly_variable_lcpm) ?

ADD REPLYlink written 4 months ago by Carlo Yague4.1k

The problem probably stems from the fact that you are using two different objects for the measurements and gene names. You should transfer gene names (as row names) to the highly_variable_lcpm object before selecting / ordering the subset of genes you want to plot, and then plot as

heatmap.2( highly_variable_lcpm, labRow = rownames(highly_variable_lcpm), ...

If you show us your code, we will be able to help more specifically.

ADD REPLYlink written 4 months ago by h.mon16k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1256 users visited in the last hour