edgeR: the F-statistics changed after removing lowly expressed gene
1
0
Entering edit mode
Pei • 0
@e9de1a10
Last seen 20 months ago
United States

Hello everyone:

I used the glmQLFTest function by edgeR to detect differentially expressed genes. I have ONE dataset and two gene lists: (1) all genes (2) genes remained after removing lowly expressed genes by filterByExpr

In the result table, I found that logFC and logCPM for each gene remained unchanged between (1) and (2), which is as expected. However, the F statistic changed, tend to be decrease, for some genes.

As a result, genes identified as significant (FDR < 0.05) in (1) would become in-significant in (2).

Dose this make sense? Thanks in advance!

sessionInfo( )

```

glmQLFTest edgeR • 747 views
ADD COMMENT
2
Entering edit mode
@gordon-smyth
Last seen 43 minutes ago
WEHI, Melbourne, Australia

glmQLFTest() function is designed to be used after expression filtering. You should apply expression filtering, regardless of whether it results in more DE genes or not.

Filtering often increases the number of DE genes. If that is not what you are seeing, then the residual deviance may be getting underestimated by an overabundance of very low counts. The quasi-likelihood method is not designed to be used when the average count is much less than 1 for some genes, which can happen if filtering is not done.

ADD COMMENT
0
Entering edit mode

Thank you Prof. Smyth! In my case, I found that number of significant genes tend to decrease after filtering when using glmQLFTest but number of significant genes tend to increase after filtering when using LRT.

I am not sure what's going on but the result by LRT looked more reasonable.

ADD REPLY
1
Entering edit mode

There is nothing to be gained by doing analyses without filtering. We don't recommend it.

ADD REPLY

Login before adding your answer.

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