DeSeq2 MA plot

0

Entering edit mode

7.5 years ago

niu.shengyong ▴ 70

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? Thanks! MA plot

ddsHTSeq <- DESeqDataSetFromHTSeqCount(sampleTable=sampleTable, directory = directory, 
design=~condition)
#ddsHTSeq
colData(ddsHTSeq)$condition<-factor(colData(ddsHTSeq)$condition, levels=c('control','test'))

dds<-DESeq(ddsHTSeq)
res<-results(dds)
res<-res[order(res$padj),]
head(res)
plotMA(dds, ylim=c(-2,2), main='M', alpha=0.1)

RNA-Seq DEseq2 MAplot • 6.0k views

ADD COMMENT • link 7.5 years ago by niu.shengyong ▴ 70

1

Entering edit mode

What version of DESeq 2 are you using? Note that the MA plots shown in the current DESeq2 vignette plot res or resLFC, and you are plotting dds.

Link to DESeq2 vignette:

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html

ADD REPLY • link 7.5 years ago by mastal511 ★ 2.1k

0

Entering edit mode

I use DESeq2 1.14.1. Is it too old?

ADD REPLY • link 7.5 years ago by niu.shengyong ▴ 70

0

Entering edit mode

Hello niu.shengyong!

It appears that your post has been cross-posted to another site: https://support.bioconductor.org/p/96645/

This is typically not recommended as it runs the risk of annoying people in both communities.

ADD REPLY • link 7.5 years ago by WouterDeCoster 47k

0

Entering edit mode

Thanks for your notice. Sorry I'm quite new in this area. I'll avoid this next time. Thanks again!

ADD REPLY • link 7.5 years ago by niu.shengyong ▴ 70

Login before adding your answer.