RNA seq differential expression analysis
1
0
Entering edit mode
@14ef1b09
Last seen 7 weeks ago
Egypt

I have a DGE object containing tumor and control RNA seq samples. Tumor n=126 and Control=390. I want to filter low expressed genes before starting limma analysis. I want to know how to set the threshold of filtering based CPM values and library size

The mean library size is 863989.

DEGseq RNASeqData limma • 730 views
ADD COMMENT
0
Entering edit mode

Just follow the limma-voom manual, it includes prefiltering.

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

Use filterByExpr(), for example as described here:

http://bioconductor.org/packages/release/workflows/vignettes/RNAseq123/inst/doc/limmaWorkflow.html

or in Chapter 15 of the limma User's Guide:

https://bioconductor.org/packages/release/bioc/vignettes/limma/inst/doc/usersguide.pdf

ADD COMMENT

Login before adding your answer.

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