Filtering the basemean value in volcano plots
I wanted to know if there was a way to filter the basemean value in volcano plots. I tried using baseMeanThreshold() function but it produced an error saying unused argument. I am trying to create volcano plots for Deseq2 results using the EnhancedVolcano() function. Could you please help with that as I am trying to filter out lowly expressed genes?
Thank you!
• 1,559 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Look, all these plotting packages are essentially just wrappers around ggplot, meaning they accept whatever you give them given it's in the expected format. If you want to only plot genes above a certain baseMean then filter your DESeq2 results table for that baseMean and then pass the output to the plotter. By the way, the baseMean is not necessarily a filter for expression level as it is the average for a gene across all samples. A baseMean can be low because 3 groups do not express a gene but one group expresses or moderately. DESeq2 is aware of that while calculating stats so best would be to either use all genes or exclude genes with low power as e.g. identified by running
edgeR::filterByExpr()before the DE analysis. It is not standard to filter for baseMean for what's worth.Okay got it. Thank you so much for your help!