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!
Prior answer from @Kevin: A: How can I generate Heat Map with dendograms, and PCA analysis in "R Programming"
Also see this collection of tutorials: Heatmaps in R
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.
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.