Question: Heatmap of Differentially Expressed Genes
1
gravatar for ChIP
2.5 years ago by
ChIP480
Netherlands
ChIP480 wrote:

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

heatmap rna-seq R • 2.1k views
ADD COMMENTlink modified 2.5 years ago • written 2.5 years ago by ChIP480
1

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. 

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by mkulecka300
0
gravatar for EVR
2.5 years ago by
EVR460
Earth
EVR460 wrote:

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.

ADD COMMENTlink written 2.5 years ago by EVR460

Hi,

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

Could you please elaborate.

Thank you

ADD REPLYlink written 2.5 years ago by ChIP480
1

ANOVA analysis
 

ADD REPLYlink written 2.5 years ago by EVR460
1

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

ADD REPLYlink written 2.5 years ago by stolarek.ir550
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: 1984 users visited in the last hour