How to filter glmTreat results by a false discovery rate (FDR) of 0.05
1
0
Entering edit mode
mictadlo ▴ 10
@mictadlo-10885
Last seen 4.2 years ago

Hello, Currenty, I call glmTreat method as following:

tr <- glmTreat(fit, contrast=B.LvsP, coef = 2, lfc=log2(1.5))
topTags(tr)

I discovered this paper and In the "Differential gene and microRNA expression" section they described that the differentially expressed genes (DEGs) were filtered using the glmTreat function in edgeR, which tests for significant differences relative to a set log2-fold change (logFC) cutoff of 0.5. Results were then filtered by a false discovery rate (FDR) of 0.05.

How did they performed the above task?

Thank you in advance.

edger • 1.3k views
ADD COMMENT
2
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 58 minutes ago
The city by the bay

Seems pretty simple, just do lfc=0.5.

ADD COMMENT

Login before adding your answer.

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