Question: DeSeq2 MA plot
0
gravatar for niu.shengyong
23 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!

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 • 1.3k views
ADD COMMENTlink modified 23 months ago by Michael Love23k • written 23 months ago by niu.shengyong0
Answer: DeSeq2 MA plot
1
gravatar for Michael Love
23 months ago by
Michael Love23k
United States
Michael Love23k 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.

ADD COMMENTlink written 23 months ago by Michael Love23k
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 16.09
Traffic: 307 users visited in the last hour