DESeq MA-plot diagnosis

Dear Biostarers,

I recently did RNA-seq analysis with DESeq with the counts from HTSeq, and the MA plot looks different from what I have before:

thousands of differentially expressed genes were identified. one thing confused me is that the red spots ( adjust pvalue <0.1) spread around. for lowly expressed genes, the log2 fold needs to be bigger to get it more significant right?

This is a simple experimental design: condition<- c("control", "control", "treatment1")

I do not have duplicate for the treatment1. I thought it was due to this, but I have duplicates for treatment2.

condition2<- c("control","control","treatment2","treatment2")

the MA plot is similar:

previously, I have another MA plot below which I think is more "right".

Any insights? Thank you so much!

Tommy

rna-seq

• 7,278 views • link

updated 4.7 years ago by Ram

• written 12.1 years ago by Ming Tommy Tang

1 answer

This is likely just a bug in how the colors are assigned. See here: Deseq's plotMA color-coding at random?

• 1 views • link

written 12.1 years ago by Devon Ryan

Thank you Devon Ryan, I saw that post and when I reordered the res dataframe and plot with plotMA, it looked fine, but why sometimes plotMA woks fine while sometimes it works not as expected?...

• 0 views • link

written 12.1 years ago by Ming Tommy Tang

More posts like this

PCA of expression data what PC to use?
written 10.6 years ago by Floris Brenk

Hi all, I'm looking for covariates in my expression dataset (countdata CAGEseq). I have 120 samples and about 20.000 expression values. I found this tutorial …
Which program, tool, or strategy do you use to visualize genomic rearrangements?
written 10.9 years ago by stianlagstad

**Which program, tool, or strategy do you use to visualize genomic rearrangements**? In relation to my master thesis I'm working on tools to visualize fusion …
MA plots DESeq: strange MA plots
written 11.1 years ago by VHahaut

Hello! I am currently analyzing with DESeq2 a group 32 ovine tumors against 3 healthy controls. 80% of the tumor samples are technical replicates. I …
RNAseq heatmap.2 log2FC clustering
written 11.7 years ago by Ming Tommy Tang

Hi there, I have a RNA-seq experiment with four samples. control, treatment1, treatment2, treatment3. each with duplicates, I used DEseq to determined the differentially expressed …
estimation of size factors in DESeq2 analysis
written 12.1 years ago by Assa Yeroslaviz

Hi all, in relation to a mail from January this yearin the r-help community, I followed Simon's advice to do my analyses in DESeq2 instead …
HiSeq2000 RNA-seq QC
written 12.3 years ago by Ming Tommy Tang

Hi biostarers, I am relatively new to RNA-seq, and have some question about the HiSeq 2000 RNA-seq 100bp pair-end data. I got the data from …
Edger Gof Plots
written 12.5 years ago by seqsearch •

Dear all, I'm doing an EdgeR analysis and made some GOF plots after estimating geneewise dispersion and fiting my data. I noticed that my values …
Plot Positions On Chromosomes
written 12.8 years ago by Nicolas Rosewick

Hi, How can I plot positions on chromosomes using R. I put a little example to better explain. So I've a dataframe for each sample …
Extract Dendrogram Information From Heatmap Generated By Heatmap.2
written 13.0 years ago by Ming Tommy Tang

Hi all, I am having a hard time to extract the corresponding subgroup members in the heatmap generated by heatmap.2 in gplots. I have looked …
Wgcna -Dilemma In Choosing Soft Threshold
written 13.1 years ago by pixie@bioinfo

I have two expression matrices (normalized, using dChip), with the same data except that in one of them , the expression values are log-transformed. I …