Question: log fold change threshold in Deseq2 - using lfcthreshold
1
gravatar for lirongrossmann
23 months ago by
lirongrossmann40 wrote:

Hi Everyone, 

I am running a rna seq differential gene expression experiment with the following code:

dds <-DESeqDataSetFromMatrix(countData = ep,colData = cp,design =~Age+Gender+Response)
dds <- DESeq(dds)
res <- results(dds,contrast=c("Response","Yes","No"))

Looking at the results I get about thousands genes with absolute log fold change > 0.5.

However, when I run the following command:

resGA <- results(dds, lfcThreshold=0.5, altHypothesis="greaterAbs")

I only get 256 with absolute log fold change > 0.5.

Anyone can help me figure out what is the source of this discrepancy?

Thanks

deseq2 lfcthreshold • 1.3k views
ADD COMMENTlink modified 23 months ago by Michael Love26k • written 23 months ago by lirongrossmann40
1

Your first test is the following: pvalue is describing the probability that log fold change != 0, and the log fold change value is the most likely value given the data.

Your second test is: pvalue is describing the probability that the abs(log fold change value) < 0.5

You can therefore see how a gene with lfc of 0.51 may be significantly more than zero, but not significantly greater than 0.5

ADD REPLYlink written 23 months ago by kieran.mace20
Answer: log fold change threshold in Deseq2 - using lfcthreshold
4
gravatar for Michael Love
23 months ago by
Michael Love26k
United States
Michael Love26k wrote:

Using lfcThreshold is testing against an LFC threshold of 0.5, which is not the same as asking for which genes the estimated LFC is larger than 0.5. See the DESeq2 paper for discussion of this topic, under the heading, "Specifying minimum effect size":

https://genomebiology.biomedcentral.com/articles/10.1186/s13059-014-0550-8

ADD COMMENTlink written 23 months ago by Michael Love26k
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 16.09
Traffic: 150 users visited in the last hour