Filtering genes with limma-voom
0
0
Entering edit mode
4.1 years ago
im ▴ 30

I am trying to follow this tutorial on TCGA-LUAD RNA-Seq data, but I am having trouble coming up with an appropriate cutoff to produce a "good" voom plot. I was doing some research on how they should be chosen, and since the filtering is generally done using cpm, the cutoff depends on library size. However, I am having some trouble finding how it should relate to library size - does anyone mind helping me out with what I should try?

limma voom limma-voom R RNA-Seq • 3.4k views
ADD COMMENT
0
Entering edit mode

edgeR has a function to perform automatic filtering, see ?filterByExpr. I guess filterByExpr is not used in that tutorial as it was (probably) introduced after the tutorial was published. There are several posts / tutorials incorporating filterByExpr() in their workflow, for example: RNA-seq analysis is easy as 1-2-3 with limma, Glimma and edgeR.

ADD REPLY

Login before adding your answer.

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