show the names of the genes on the Y axis on the heatmap Deeptools
1
0
Entering edit mode
4.6 years ago
Bioiris • 0

Hello, I made the heatmap of épégénétiques marks around TSS using Deeptools (ComputeMatrix then plotHeatmap) but the problem that the latest tools (plotHeatmap) gives me a heatmap with the region (+ -1kb TSS and 1KB) in X axis.

is there any way to show the names of the genes on the Y axis on the heatmap?

thank you

RNA-Seq ChIP-Seq deepTools • 1.7k views
ADD COMMENT
1
Entering edit mode
ADD REPLY
0
Entering edit mode

Is there a way to just get the top 100 genes from the heatmap as a list?

ADD REPLY
1
Entering edit mode
4.6 years ago

There's no option in deepTools to display the gene/transcript/region names on the heatmap, normally there are so many of them that you wouldn't be able to see them. If you really want the names displayed then you'll need to tweak the code to do that (if you're interested in doing so, just mention that in a comment and I'll look through the code for where you'll likely need to modify things on Monday).

ADD COMMENT
0
Entering edit mode

yes I'm very interested if you can do it I will be very grateful

ADD REPLY
0
Entering edit mode

This will take a bit of playing around to get things correct, but the following is roughly how it'll work:

  1. You'll need to delete this line to be something more like ax.axes.set_yticks(sub_matrix['matrix'].shape[0])
  2. You'll then need to add new ticks and labels, which is something like ax.axes.yticks(np.arange(sub_matrix['regions'].shape[0]), labels).
  3. The labels will have to be created and is in hm.matrix.regions I think.

Hopefully you only have one group, since it'll be tougher to do correctly with more.

ADD REPLY
0
Entering edit mode

Hello,

I tried this but I still do not get the gene names on the y-axis. Is there a way I can just obtain the list of the top 100 genes? I am interested in obtaining the top genes where the peaks are near the TSS. Thank you.

ADD REPLY

Login before adding your answer.

Traffic: 2410 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6