Question: Pairwise comparisons in Differential Expression
0
gravatar for Ron
17 months ago by
Ron920
United States
Ron920 wrote:

Hi all,

I have 4 RNAseq samples, both at Day 0 and Day 3. So total 8 samples.

I want to do differential Expression analysis between Day 0 and Day 3,but pairwise.

I wanted to confirm if this is The best way to do this is using multi-factor design ? from DEseq2 package. The condition of interest should go at the end of the design formula, e.g. ~ subject + condition. Does my metadata look fine for multi-factor design ?

Sample  Subject Condition
Sample1_day0    Day0    Rep1
Sample2_day0    Day0    Rep2
Sample3_day0    Day0    Rep3
Sample4_day0    Day0    Rep4
Sample1_day3    Day3    Rep1
Sample2_day3    Day3    Rep2
Sample3_day3    Day3    Rep3
Sample4_day3    Day3    Rep4

Also, When we do differential expression of 1 group vs the other normally (without multicondition), we get results in the form of 1 log foldchange column ,1 p value column etc. How would we get the results in this ?

Thanks,

Ron

rna-seq deseq • 963 views
ADD COMMENTlink written 17 months ago by Ron920

Your metadata does not look like a multi-factor design. It looks like an experiment with a single factor (Day) with two levels (Day 0 and Day 3). Does "Rep" stand for replicate, or is it a designation of paired samples between the days? The model you mention (~ subject + condition) would I believe not be fittable given you have one observation per subject:condition combination. I guess what you want is a test between days, paired by your "Condition" column. If using a t.test, that could be accomplished in R by (for example) t.test(..., paired=TRUE).

ADD REPLYlink written 17 months ago by Ahill1.5k

Rep1 is an annotation I added to make sure that Sample1_day0 and Sample1_day3 are replicates since they need to be paired, while doing pairwise comparison.

ADD REPLYlink written 17 months ago by Ron920

OK, in that case the model (~subject + condition) will fit a two factor model for each gene. Interpretation would depend on if your condition column is technical factor you are looking to normalize away, or if it is a biological factor of interest that you are trying to quantify. What is your main interest, the differences between Day0 and Day3, or the differences among Rep1, Rep2, Rep3, and Rep4? You would extract the corresponding effects (logFC) and p-values from the two factor model. A good idea is to consult a statistician if you are not clear on the best way to model the data.

ADD REPLYlink modified 17 months ago • written 17 months ago by Ahill1.5k

I was looking at this post's answer and think this is what I want too:

Deseq2 pairwise comparision

ADD REPLYlink written 17 months ago by Ron920
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: 1857 users visited in the last hour