[BioC] Error in kooperberg procedure

michael watson (IAH-C) michael.watson at bbsrc.ac.uk
Fri Sep 17 15:49:04 CEST 2004


I ran this code and got an error:

>genepixFiles <- targets$Name
>read.series(genepixFiles, suffix=NULL,skip=26, sep="\t")
>layout <- RG$printer
>RGmodel <- kooperberg(genepixFiles, layout=layout) 
Error in integrate(.numeratorBayesianAdjustedFG, ifelse((fg - bg - 4 *
        a limit is missing

I guess this may be related to my earlier problem, where I have BG
estimates = zero for one of my arrays, but this is a little trickier to
solve as the data is in data.frames, and I cannot so easily add a small
value to Rb and Gb....


More information about the Bioconductor mailing list