DESeq2 likehood ratio
Entering edit mode
SE • 0
Last seen 4 days ago
United Kingdom

I have RNAseq data from cells treated with 3 drugs , drugA, drugB, drugC and different combinations of those drugs so in total 8 arms, drugA, drugB, drugC, drugA+drugB, drugB+drugC, drugA+drugB+drugC, drugA+drugC and no treatment/control arms. I have 3-5 replicates per arm.

I am using DEseq2 for analysis but I am not sure what should be the design. Is a simple model with treatment as a single factor with 8 levels enough or a more complex model is more appropriate?.

Is LRT appropriate?. I tried to compare it to simple model using the code below and I found a lot more significantly expressed genes than using the simple model. Also on biological level, there mono therapy drugs are very different from the combinations with some resulting in additive effect and some combinations resulting in contradictory effect.

dds_LRT <- DESeq(ddsobj, test='LRT', reduced=~1)

res <- results(dds_LRT, lfcthreshold=0)

Can I obtain pair-wise comparisons from the dds_LRT to know the DEgenes are differentially expressed among which groups using the code below?.

res_LRT <- results(dds_LRT, contrast= list('drugA+drugB+drugC' , 'drugB+drugC'))

the aim of using the above code to study the effect of adding drugA to the combinations.

Is that correct?.

I am using DESeq2 version 1.32

Just to add

if dds_LRT is the ddsobject generated using LRT test and dds is the ddsobject generated by Wald test , there are a lot more Differentially expressed genes between the following 2 codes

res_LRT <- results(dds_LRT, contrast= c('Treatment', 'drugA+drugB+drugC' , 'drugB+drugC')) 

res <- results(dds, contrast= c('Treatment', 'drugA+drugB+drugC' , 'drugB+drugC'))


DESeq2 • 188 views
Entering edit mode

And further questions here:

DESeq2 likehood ratio

Entering edit mode
Last seen 8 hours ago
United States

For questions about how to design the statistical analysis and interpret the results, I recommend collaborating or consulting with a statistician or someone familiar with linear models in R. I have to reserve my time on the support site for software related questions.

To the other question, yes you can run pairwise (Wald) comparisons on the object after fitting an LRT.

Entering edit mode

Ok thanks for your reply.

Why there is a large difference between the differentially expressed genes between the pairwise comparisons obtained from results() using a dds_LRT object and using the same object dds_wald object?. which is correct?

dds_wald <- DESeq(ddsobj)

res <- results(dds_wald, lfcThreshold = 0, contrasts = c('B-control', 'A-control'))

dds_LRT <- DESeq(ddsobj, test="LRT", reduced = ~ 1)

res_LRT <- results(dds_LRT, lfcThreshold = 0, contrasts = c('B-control', 'A-control'))

The number of DE genes in res_LRT with adj p value <0.1 is 600 but the number of DE genes in res is only 40

Entering edit mode

You should explicitly specify in results() the test whether you want Wald (contrasts) or LRT p-values. Wald and LRT in this case are totally different tests with different meaning. You should certainly consult with a statistician on the analysis.


Login before adding your answer.

Traffic: 312 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6