SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 12-23-2013, 03:41 AM   #1
giu.facchetti
Junior Member
 
Location: Germany

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

hello,
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)?

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

Join Date: Aug 2009
Posts: 109
Default

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
__________________
Wolfgang Huber
EMBL
Wolfgang Huber is offline   Reply With Quote
Old 01-17-2014, 01:07 AM   #3
giu.facchetti
Junior Member
 
Location: Germany

Join Date: Dec 2013
Posts: 2
Default

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
Reply

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 03:15 AM.


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