Question: How to add Log2FoldChange as an argument for plotMA in DESeq2
0
gravatar for gracechappell
3.6 years ago by
United States
gracechappell40 wrote:

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 • 1.7k views
ADD COMMENTlink modified 3.6 years ago • written 3.6 years ago by gracechappell40
0
gravatar for vivekbhr
3.6 years ago by
vivekbhr510
Germany
vivekbhr510 wrote:

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

ADD COMMENTlink written 3.6 years ago by vivekbhr510
0
gravatar for gracechappell
3.6 years ago by
United States
gracechappell40 wrote:

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. 

 

ADD COMMENTlink written 3.6 years ago by gracechappell40

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

ADD REPLYlink written 3.4 years ago by vivekbhr510
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 813 users visited in the last hour