The support.bioconductor.org editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: DeSeq2 MA plot
0
20 months ago by
niu.shengyong0 wrote:

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!

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)

plotMA(dds, ylim=c(-2,2), main='M', alpha=0.1)

deseq2 ma • 1.2k views
modified 20 months ago by Michael Love22k • written 20 months ago by niu.shengyong0
Answer: DeSeq2 MA plot
1
20 months ago by
Michael Love22k
United States
Michael Love22k wrote:

The visual artifacts are typically when one group (here control) is sequenced at much lower depth than the treated group in a two group comparison. The patterns are unavoidable. It is not great to have large differences in sequence depth across the conditions, but I don't have a decision boundary or anything at which it's too much of a problem. I would zoom out on the plot (e.g. to [-6,6]) so you can see which genes are being called DE if any.