[Bioperl-l] Re: [BioC] Error in kooperberg procedure
sdavis2 at mail.nih.gov
Fri Sep 17 09:57:25 EDT 2004
Actually, you can add constants to your matrices.
RGmodel$Rb <- RGmodel$Rb + 0.001
I'm not sure this will solve your problem directly, though. Try it and
see, but someone with a better understanding of the function may be
able to comment in the context of your last post.
On Sep 17, 2004, at 9:49 AM, michael watson (IAH-C) wrote:
> 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....
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
More information about the Bioperl-l