Search
Question: In DESeq2, why LRT result's LFC doesn't match with any of Wald test result's LFC?
0
gravatar for gmhhope
16 months ago by
gmhhope0
gmhhope0 wrote:

Hi DESeq2,

I have a question with respect to DESeq2.

I have one factor with 4 levels (wild-type, mutant A, mutant B, mutantA&B; the wild-type is set as reference). Based on this, I want to do LRT test to get analysis just like ANOVA and then compare two by two among them Because there are no post-hoc test, so I also did multiple pairwise comparison using Wald test. I learned from the tutorial that the LRT test will give only one of the Log2FC among all others. So to confirm that, I try to match the Log2FC given by LRT test's result table with those other Log2FC given by Wald test. However, I just cannot find any of them fully matched. I found one which is very similar. I guessed it may be that; however, I want to know why will this happen? And I found difficult to understand how the output of coef() translate into LFC? Can anybody help explain to me? 

 

The R code is shown below:

dds <- DESeqDataSetFromMatrix(countData = countData, colData = colData, design = ~class)

dds <- dds[rowSums(counts(dds))>1,]

dds$class <- relevel(dds$class, ref="P.UF1")
### Perform DESeq function to get all MLE & bayes procedure (including MAP)
dds1 <- DESeq(dds,test="LRT",reduced=~1)
dds2 <- DESeq(dds)
### Perform results function to get the statistics and LFC
res_LRT <- results(dds1) 
summary(res_LRT)
res_DmmcAvsPUF1 <- results(dds2,contrast=c("class","DmmcA","P.UF1"))
res_DpduPvsPUF1 <- results(dds2,contrast=c("class","DpduP","P.UF1"))
res_DmmcA_pduPvsPUF1 <- results(dds2, contrast=c("class","DpduP_mmcA","P.UF1"))
res_DmmcAvsDpduP <- results(dds2, contrast=c("class","DmmcA","DpduP"))
res_DpduP_mmcAvsDmmcA <- results(dds2, contrast=c("class","DpduP_mmcA","DmmcA"))
res_DpduP_mmcAvsDpduP <- results(dds2, contrast=c("class","DpduP_mmcA","DpduP"))

 

The 2 columns representing similar log2FCs is shown below.

LRT_log2FoldChange MPvsPUF1_log2FoldChange
0.631377697 0.627936831
1.051816895 0.898240504
-0.253159341 -0.246927052
-0.010568974 -0.010054782
-0.788608064 -0.78341221
-0.092539348 -0.090469655
-0.692334811 -0.687365411
0.503294525 0.456056859
-0.441611727 -0.438624489

 

ADD COMMENTlink modified 16 months ago by Michael Love16k • written 16 months ago by gmhhope0
0
gravatar for Michael Love
16 months ago by
Michael Love16k
United States
Michael Love16k wrote:

You are showing the moderated LFC using the Wald test and the unmoderated LFC using the LRT. The LRT uses an unmoderated fold change, while it is your choice whether to moderate or not (betaPrior=TRUE or FALSE) for the Wald test.

This information (whether you produced MAP or MLE log fold change) is printed at the top of the results table if you print it:

> res
ADD COMMENTlink written 16 months ago by Michael Love16k
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.2.0
Traffic: 274 users visited in the last hour