HI All,
thanks very much for the help, doing local fit fixes the issue
I a not sure I understand what is the impact of doing this in terms of the
accuracy of the estimates
can someone explain?
will it affect greatly the number of differential calls at a particular
cutoff?

thanks
Lucia



On Wed, Apr 17, 2013 at 1:41 PM, Simon Anders <anders@embl.de> wrote:

> Hi
>
>
> On 17/04/13 17:32, Lucia Peixoto wrote:
>
>> Error in parametricDispersionFit(means, disps) :
>>    Parametric dispersion fit failed. Try a local fit and/or a pooled
>> estimation. (See '?estimateDispersions')
>>
> [...]
>
>  PS: This is not "traditional" RNA-seq data, each row is a window in the
>> genome and the values on the matrix are mapped reads in each sample for
>> each window
>>
>
> The parametric fit fails sometimes, especially on non-standard data sets.
> This is why we offer the local fit as a fall-back.
>
>   Simon
>
>
> ______________________________**_________________
> Bioconductor mailing list
> Bioconductor@r-project.org
> https://stat.ethz.ch/mailman/**listinfo/bioconductor<https://stat.ethz.ch/mailman/listinfo/bioconductor>
> Search the archives: http://news.gmane.org/gmane.**
> science.biology.informatics.**conductor<http://news.gmane.org/gmane.science.biology.informatics.conductor>
>



-- 
Lucia Peixoto PhD
Postdoctoral Research Fellow
Laboratory of Dr. Ted Abel
Department of Biology
School of Arts and Sciences
University of Pennsylvania

"Think boldly, don't be afraid of making mistakes, don't miss small
details, keep your eyes open, and be modest in everything except your
aims."
Albert Szent-Gyorgyi

	[[alternative HTML version deleted]]

