[DESeq2] setting LFC threshold
1
0
Entering edit mode
Shail • 0
@dba6855a
Last seen 20 months ago
Hong Kong

Hi,

I want to set the p-value and LFC threshold to DESeq2 output results

The way i understood is

independentFiltering = TRUE,alpha = 0.001 #will filter results for p-value 0.001.

However, how to set an LFC threshold of 1.5. Is it lfcThreshold = 1.5?

Also, how to identify upregulated and downregulated genes.

The summary list by the p-value of DESeq2 run is as follows

res <- res[order(res$padj),]
head(res)

results

In this, I assume the genes with negative LFC are downregulated. But in which group LD1 or control?

DESeq2 • 1.3k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 12 hours ago
United States

However, how to set an LFC threshold of 1.5. Is it lfcThreshold = 1.5?

This depends entirely on what is a relevant effect size for your experiment. 2^1.5 is a pretty large change. I sometimes recommend log2(1.5) ~= .585 to start, as this covers when expression goes up or down by 1.5.

which group LD1 or control?

See the vignette or workflow for details like this. Note that the answer is printed above the table: LD1 vs control.

ADD COMMENT
0
Entering edit mode

Thank you for the help

ADD REPLY

Login before adding your answer.

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