filterbyexp() in edgeR
2
0
Entering edit mode
@edc5aa0a
Last seen 4 months ago
France

Dear R users, I want to perform a GSEA on RNA-seq data, I use the filterbyexp() function in the DGELlist in edgeR. I am not sure about hte arguments I use : I would like to keep all genes with an expression of at least 10 cpm in 10% of samples. Is it filterbyexp(DGElistobject, min.count = 10, min.prop = 0.1) ?

It keeps almost the same number of gene sthan the default filter which is I think min.count = 10, min.prop = 0.7

Best

Simon

edgeR • 236 views
0
Entering edit mode
Yunshun Chen ▴ 770
@yunshun-chen-5451
Last seen 1 day ago
Australia

min.count is a threshold for actual counts rather than cpm. If you want to filter genes based on their cpm values, you may need to take into account the library sizes and get an equivalent cut-off for min.count. E.g., if the average library size of all the samples is about 20 million, then 10 cpm would be equivalent to min.count = 200.

0
Entering edit mode
@gordon-smyth
Last seen 7 hours ago
WEHI, Melbourne, Australia

The arguments are already set to optimal values for differential expression analyses so you don't need to change them at all.

The only thing you need to do is to make sure that you specify the design of your experiment, either by setting the group variable of your DGEList or specifying the design or group arguments to filterByExpr().