This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Heatmap of Differentially Expressed Genes

Hi,

I am bit lost with the heatmap in R.

I have 4 groups of sample, group 1 (untreated) has 5 sample, group2 (treatment1) has 15 samples, group3 (treatment2) has 11 samples and group4 (treatment3) has 21 samples.

The question that I am trying to answer is to generate a heatmap, that can point the DEG across these groups. The values are in log2.

I know, how to do this in R with two samples (Control Vs Treatment) but not with so many samples with different types of treatments.

A figure like this will be helpful: https://www.google.be/url?sa=i&rct=j&q=&esrc=s&source=images&cd=&cad=rja&uact=8&ved=0ahUKEwiXq6iO19PJAhVFuBoKHdJQATsQjRwIBw&url=http%3A%2F%2Fwww.mdpi.com%2F1422-0067%2F16%2F10%2F24555%2Fhtm&psig=AFQjCNFUuwoZ-A9VIBWPow7E81Izai33tw&ust=1449918885050234

Or even better representation is welcomed.

Also, How could one find DEG across all groups with some p-value associated with it

Can anyone help?

Thank you

rna-seq r heatmap
• 6,120 views • link
updated by Ram Ram • written by ChIP ChIP

I think you should consider mapping your patients groups to different colours and plotting this alongside your heatmap with ColSideColors option. You can have a look at this post: http://www2.warwick.ac.uk/fac/sci/moac/people/students/peter_cock/r/heatmap/. If you have more than two groups, you just have to modify this function:

color.map <- function(mol.biol) { if (mol.biol=="ALL1/AF4") "#FF0000" else "#0000FF" }

to something like this:

color.map <- function(group) { if (group==my_first_group) "#FF0000" else if (group==my_second_group) "#0000FF" else if (group==my_third_group) "#006400" }

Also, I think most of the tools for DGE handle more than two groups: DESeq2 and edgeR for sure do.

• 1 views • link
updated by Ram Ram • written by mkulecka mkulecka

1 answer

Hi,

I would find the DEGs across all the groups and take their fold change or normalised counts and generate a heatmap juts like how you make for Control vs treatment.

• 1 views • link
written by EVR EVR

Hi,

How would you find DEG across all groups with some p-value associated with it?

Could you please elaborate.

Thank you

• 1 views • link
written by ChIP ChIP

ANOVA analysis

• 1 views • link
updated by Ram Ram • written by EVR EVR

ANOVA tells you that one or more are different, but doesn't indicate which ones

• 1 views • link
updated by Ram Ram • written by stolarek.ir stolarek.ir

Log in to answer this question.

More posts like this