Contrast Specific Dispersion Estimate
0
0
Entering edit mode
tk • 0
@tk-12362
Last seen 7.2 years ago

We have RNASeq data from 15 conditions, for which many contrasts but not all pairwise contrasts are of interest.  

When we run the dispersion estimate across all 15 conditions, perform the fit and extract one of contrasts of interest, we obtain much less significant DE genes than if we are first subset the data by the 2 condition use for the contrast and perform an dispersion estimate based solely on the data used. 

Question is the following:

1-Wouldn't it make sense to estimate the dispersion on a per contrasts (or even per condition) basis? Granted we lose samples in the dispersion estimate, but given that the conditions are so different, we don't see why we would estimate it altogether.  

Dispersion estimate using all 15 conditions:

Inline image 2

 

Inline image 3

 

Dispersion estimate using just 2 conditions for the contrast of interest:

Inline image 1

Inline image 4

edgeR design and contrast matrix • 742 views
ADD COMMENT
0
Entering edit mode

Your pictures aren't showing up.

ADD REPLY

Login before adding your answer.

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