Group comparison using DESeq
1
0
Entering edit mode
amir_b_e • 0
@amir_b_e-15705
Last seen 4.3 years ago

Hello all,

I have two group of samples (control vs drug treated) and within each group, I have three cell lines in biological triplicates. I am interested in the comparison of control vs drug treated. There are two ways to do this differential expression analysis.

1. Assume we have 9 samples in each group (three cell line in biological triplicates) and consider them as the replicate of each other. Then we can simply compare 9 vs 9 samples. However in this way we don't consider biological heterogeneity between cell lines.

2.we have three cell line and three replicates for each and we want to do group comparison like the following

((control - cell line1) and (control - cell line 2) and (control - cell line 3) Vs. (treated - cell line 1) and (treated - cell line 2) and (treated - cell line 3))

So then instead of 9 vs 9 comparison, we do (3+3+3 vs 3+3+3) which consider each cell line as a group and then we compare three groups vs three groups. 

I want to run the second type of analysis but have no idea how to this. Could you please introduce a tutorial for that? Any help would be highly appreciated.

Thanks! 

deseq2 multiple groups • 528 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 9 hours ago
United States

There is a huge (HUGE!) vignette that covers pretty much everything you need to know about DESeq2, including how to fit a model with blocking factors. If you plan to use the package, you should read the vignette first.

ADD COMMENT
0
Entering edit mode

What I could find was the first way of analysis. Could you please help me with the second type? In which part of vignette I can find some info about the second way?

ADD REPLY
0
Entering edit mode

What you are describing is a simple model where you block on cell line. This is analogous to fitting a paired model (where you block on subject), or really any other blocking.

 

ADD REPLY

Login before adding your answer.

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