Error in kooperberg procedure
1
0
Entering edit mode
@michael-watson-iah-c-378
Last seen 9.7 years ago
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
• 610 views
ADD COMMENT
0
Entering edit mode
Matthew Ritchie ▴ 1000
@matthew-ritchie-650
Last seen 17 days ago
Australia
Dear Mick, Yes, a 0 value for the background intensity will cause this error. I have modified the code in the latest version of limma (1.7.7) to handle 0's, so hopefully this is no longer a problem. Best wishes, Matt Ritchie 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 > [[alternative HTML version deleted]]
ADD COMMENT

Login before adding your answer.

Traffic: 386 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6