[BioC] estimating size factors in DESeq produces warnings

Assa Yeroslaviz frymor at gmail.com
Tue Jan 14 18:55:29 CET 2014


I am using Deseq to analyze my data. I have three replica for each of the
two conditions.

When I am running the script in the normal way:

cds = newCountDataSet( Counts_set, condition )
cds = estimateSizeFactors( cds ) # uses the median as location function
cds = estimateDispersions( cds )
I get the following error massage:
There were 16 warnings (use warnings() to see them)
> warnings()
Warning messages:
1: In log(ifelse(y == 0, 1, y/mu)) : NaNs produced
2: step size truncated due to divergence
The dispersion plot looks like that: DispersionPlot.png

but with fitType="local" ot looks like that: DispersionPlotLocal.png

When I am doing the same analysis, but with fitType="local" the warnings
I read that this happens ,When One is using the analysis without replica,
but I have three for each condition.
Is it preferable to run the analysis with local fit type?

What is the meaning of these warnings?

-------------- next part --------------
A non-text attachment was scrubbed...
Name: DispersionPlot.png
Type: image/png
Size: 51440 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/bioconductor/attachments/20140114/2f903926/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: DispersionPlotLocal.png
Type: image/png
Size: 50242 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/bioconductor/attachments/20140114/2f903926/attachment-0001.png>

More information about the Bioconductor mailing list