[BioC] DESeq error with glm.fit during estimateDispersions

Simon Anders anders at embl.de
Sat May 26 11:24:06 CEST 2012


On 2012-05-24 22:15, chris_utah wrote:
> cdsFull<- estimateDispersions( cdsFull, method="pooled-CR", modelFormula = count ~ sample+allele+allele:diet)
> Error in glm.fit(mm, y, family = MASS::negative.binomial(initialGuess),  :
>    NA/NaN/Inf in foreign function call (arg 1)

As just discussed off-list with Chris: This was due to a missing "try" 
in the code to 'estimateDispersion'.

I have just fixed the issue in DESeq version 1.8.3 (release) and 1.9.7 


More information about the Bioconductor mailing list