Question: LRT results DESeq2
gravatar for bekah
14 months ago by
bekah20 wrote:


​I am running analysis on RNA Seq data using DESeq2, I have two "treatments" for each sample each with two levels

status: sham or infected
temperature: high or low.

​I have run the LRT to look at whether I am make sure I have not made any assumptions based on experimental design on whether genes are significantly dif expressed. To do this I have used the LRT

dstemp_lrt<-DESeq(dstemp, test="LRT", reduced= ~status)

to give:
log2 fold change (MLE): temperature L vs H
LRT p-value: '~ status + temperature' vs '~ status'

​Am I correct in thinking that this gives the Pvalues of the genes significantly differentially expressed as a result of whether it was sham or infected in comparison to those genes sig. differentially expressed when temperature is also factored in (the interaction)?
Therefore I should also run with "reduced=~temperature" to test the same but with just temperature and then compare the two lists?

Best wishes,



deseq2 job • 247 views
ADD COMMENTlink modified 14 months ago • written 14 months ago by bekah20

Ah cheers! So this can be used to identify if the DEG I have found when looking at the four individual experimental groups (infected L, infected H, control L, control H) can be attributed to the difference in temperature or difference in status?


ADD REPLYlink written 14 months ago by bekah20

Basically, yes. You should see which have the differences.

ADD REPLYlink written 14 months ago by Michael Love24k
Answer: LRT results DESeq2
gravatar for Michael Love
14 months ago by
Michael Love24k
United States
Michael Love24k wrote:

The p-value you have there tests for differences in gene expression due to temperature, accounting for differences due to status.

ADD COMMENTlink written 14 months ago by Michael Love24k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 113 users visited in the last hour