This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to add Log2FoldChange as an argument for plotMA in DESeq2

Hello,

I would like to add an argument to the plotMA function in DESeq2. Specifically, I would like to highlight points on the MA plot that are not only alpha<0.1, but also log2FoldChange < -1 or >1. Is there a way to achieve this? I tried the following "ifelse" statement but it did not work.

plotMA(res,alpha=0.1,col=ifesle(log2FoldChange>0.75 | log2FoldChange<-0.75, "red", "grey"), cex=2, ylim=c(-2,2))

Thank you!

deseq2
• 4,720 views • link
updated by Ram Ram • written by gracechappell • gracechappell

1 answer

You can subset the DESeq result object by abs(logFC) > 1, then plot by plotMA. That's the fastest..

• 0 views • link
written by vivekbhr vivekbhr

Thanks for the reply, Vivek. I'd like to retain the non-significant points in my plot as grey. By subsetting I only plot those that meet my criteria.

• 0 views • link
updated by Ram Ram • written by gracechappell • gracechappell

You can manually make the MA plot with all the genes, then color the genes of interest.

• 0 views • link
updated by Ram Ram • written by vivekbhr vivekbhr

Hello, how are you?

Excuse me, I'm new to Deseq. I don't understand how to do it? I also have that question, how to change the ranges of Log2FoldChange in plotMA?

I thank you if you explain to me

:D

• 0 views • link
updated by Ram Ram • written by Johan Largo • Johan Largo

Log in to answer this question.

More posts like this
  • Change log2FoldChange range - plotMA
    written by Johan Largo • Johan Largo
    Hi all I am trying to make a graph with the plotMA tool and it throws me the graph below, now I would like from …
  • Generating a heatmap in DEG analysis
    written by Pranathi • Pranathi
    After performing DESEQ2 on my data, I could able to plot MA, PCA, and EnhancedVolcano plots including a (.xlsx) file consisting of log fold change …
  • DESeq2: How do one label specific points (gene) on the DESEq2 plot?
    written by jaqx008 jaqx008
    Hello everyone. I am doing a DESeq2 analysis and wanting to identify all the genes (points) colored in red. I found this somewhere (see below) …
  • add cut off lines to volcano plot
    written by Sam Sam
    Dear all I got the volcano plot of deseq2 output with this code : with(res, plot(log2FoldChange, -log10(pvalue), pch=20, main="Volcano plot", xlim=c(-8,8))) with(subset(res, padj<0.1 ), points(log2FoldChange, …
  • Volcano_plot using R
    written by Shina • Shina
    using the following R Script i want to generate Volcano plot.but my volcano plot is not correct because their is no connecting point between -logpvalue …
  • Volcano Plot from DEseq2
    written by 1769mkc 1769mkc
    Im using this code to make based on log2foldchange and padj value ,im getting the plot but i want those value for my reference how …
  • DESeq2 Heatmap issue
    written by dazhudou1122 dazhudou1122
    Dear Community, I stole the codes from Stephen Turner to analyze RNAseq data using DESeq2 (out put of featureCount). Everything worked well, except that I …
  • How to add single gene name in volcano plot ?
    written by Mike Mike
    Hi all, I am using R script from Getting Genetics Done ( http://www.gettinggeneticsdone.com/2014/05/r-volcano-plots-to-visualize-rnaseq-microarray.html) to generating Volcano plots but my problem is how to show particular …
  • DeSeq2 MA plot
    written by niu.shengyong • niu.shengyong
    I use the following codes and generate this MA plot. It looks quite weird for me. Is there any problem I should take care of? …
  • Ma plot labelling
    written by op263 • op263
    Hi there, I'm trying to generate a graph with the results from my first RNA seq experiment. I just used R for the first time …