Creating a heatmap from DGE data
0
0
Entering edit mode
2.1 years ago

Hi there, Apologies if this question is a repeat, but I couldn't find the answer I'm looking for on here or youtube. Perhaps someone could point me in the right direction?

I have the DGE RNA-seq data for a gene knock down experiment versus a control (with 3 repeats and 3 control repeats). I have 718 genes and have the log2 fold change (LFC) for all of the genes.

In my .txt file that I read into R: I have gene names in the row, and LFC in a single column. Can I plot a heatmap from this data? I want to be able to visualise which genes are up-regulated or down-regulated and how they cluster. All the other heatmap tutorials I have seen have plotted data for 2 or more columns.

I am unsure how to proceed now I have read this file in.

Best wishes!

heatmap • 1.1k views
ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Also see this collection of tutorials: Heatmaps in R

ADD REPLY
0
Entering edit mode
  1. You can plot heatmap with 1 column. But you may need some knowledge of writing the script to customize you heatmap.
  2. It doesn't make sense to me to cluster with just 1 column. Isn't it just a sorting?
ADD REPLY
0
Entering edit mode

Right! I am struggling to find a guide that will help me with the code for this. I keep getting an error in R saying I need more than one column.

ADD REPLY
0
Entering edit mode

Instead of plotting the average LFC you could plot the read counts in each replicate and each treatment status (total 6) and then plot that heat map. You can do this only for De genes and see how they cluster.

ADD REPLY

Login before adding your answer.

Traffic: 2352 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