Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
count v dispersion plot stephenhart Bioinformatics 1 07-18-2013 05:04 AM
In DESeq, should one estimate sizeFactors and Dispersion in a pairwise fashion? carmeyeii Bioinformatics 4 06-19-2013 11:56 AM
data inclusion for dispersion estimates in DEseq? amcloon Bioinformatics 0 06-19-2013 01:44 AM
DESeq high withing group dispersion wdt RNA Sequencing 0 07-16-2012 04:54 PM
DESeq Estimate Dispersion mattia Bioinformatics 1 02-09-2012 02:42 AM

Thread Tools
Old 12-23-2013, 04:41 AM   #1
Junior Member
Location: Germany

Join Date: Dec 2013
Posts: 2
Default dispersion from DESeq

I have started working with DESeq for RNASeq data from plants.
Among the first plots I have got, there is the dispersion (D) vs counts mean (M). This plot is fitted with the curve D=a+b/M (which in a log-log scale is a decreasing curve that reaches a plateau for high value of M, see attachment).
However, my data show a linear correlation (still in the log-log plot).
Is my analysis wrong? can I change the fitting expression? why the fitting procedure does not set a=0 which would give a linear relationship between Log(D) and Log(M)?

Attached Images
File Type: jpg dispersion.jpg (25.8 KB, 25 views)
giu.facchetti is offline   Reply With Quote
Old 12-26-2013, 01:06 AM   #2
Wolfgang Huber
Senior Member
Location: Heidelberg, Germany

Join Date: Aug 2009
Posts: 109

Dear Giu.facchetti,

thank you for the feedback. Some suggestions:
1. The attached jpeg image came out very small (100x115x) making it hard to impossible to read (at least for me).
2. Such a strange-looking plot may (or may not) be indicative of a data quality problem, e.g. large batch effect or sample-swap. Have you explored that? (The DESeq2 vignette gives some suggestions on how to do that.)
3. Have you tried the argument fitType = "local" to estimateDispersions?
4. Since you are just starting, the authors of the package recommend you use DESeq2 now, since it offers many improvements.

Kind regards
Wolfgang Huber
Wolfgang Huber is offline   Reply With Quote
Old 01-17-2014, 02:07 AM   #3
Junior Member
Location: Germany

Join Date: Dec 2013
Posts: 2

Thank you for the suggestions!
Unfortunately the option "fit=local" did not improve the fitting (the fitting line was even much curved). Same results with DESeq2.
At the moment I am doing again the alignment and reads counts, trying to improve the sensitivity and accuracy.
giu.facchetti is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:48 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO