[Bioperl-l] Re: [BioC] Error in kooperberg procedure

[Sean Davis]

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.

Sean

On Sep 17, 2004, at 9:49 AM, michael watson (IAH-C) wrote:

> Hello
>
> 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....
>
> Mick
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor