DeSeq2 MA plot
1
0
Entering edit mode
@niushengyong-12910
Last seen 6.8 years ago

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 :https://drive.google.com/file/d/0ByJ9Aiypgye9aHNWSHlDZEk0dEk/view?usp=sharing

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)

 

 

 

deseq2 MA • 2.4k views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 10 hours ago
United States

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.

ADD COMMENT

Login before adding your answer.

Traffic: 960 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6