DESeq2: Different Dispersions for different sample groups
2
1
Entering edit mode
hhoeflin ▴ 10
@hhoeflin-7706
Last seen 8.4 years ago
Switzerland

Hi,

I am using DESeq2 for analyzing an RNASeq Experiment. The details of the setup are as follows:

I have two treatments for which I want to test for differential expression. However, when inspecting the data, the variability of the counts in one treatment group is considerably higher than the variability in the other group (for same mean).

If I understand the DESeq2 documentation correctly, the DESeq function fits one dispersion parameter per gene, not allowing for sample groups with different dispersion. 

a) Is this understanding correct?

b) Is there a way to allow for sample-group specific dispersions? In a linear model setting, this would usually be achieved by setting different weights for the samples in different groups.

Thanks!

Holger

deseq2 estimatedispersions • 2.5k views
ADD COMMENT
1
Entering edit mode
@steve-lianoglou-2771
Last seen 20 months ago
United States

Not sure if it has hit your radar, but the limma/voom folks have recently published a paper that describes the voomWithQualityWeights function, which allows you to achieve this behavior to analyze rna-seq data with sample level weights:

 

ADD COMMENT
0
Entering edit mode

Great thanks - that sounds very close to what I need. No, it hadn't hit my radar yet. Guess I will switch the workflow over to limma.

~Holger

 

ADD REPLY
0
Entering edit mode
@mikelove
Last seen 1 day ago
United States

hi Holger,

a) yes, you're right. we fit one dispersion estimate per gene, so locking the relationship Var = mean + dispersion * mean^2.

b) There is not a way to have group-specific dispersion values with DESeq2. Off the top of my head, I think baySeq offers different dispersion estimates per group: getPriors.NB(..., equalDispersions=FALSE)

ADD COMMENT

Login before adding your answer.

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