Entering edit mode
Hi,
Thanks for reporting this bug. When using the .locs file each array
segment
is processes seperately to produce the imageplot. This bug manifests
itself
if the width of the segments in pixels differs modulo the squareSize
argument. I thought I had patched this, but clearly not.
My advice would be to increase the squareSize argument (the default
for
BeadChips is 50 pixels, so something larger than that) and try running
it
again.
Alternatively set the useLocs argument to FALSE and beadarray will
infer the
network of neighbouring beads itself. This issue won't be a problem
then,
but it has the downside of being both slower and not showing the
seperation
between individual array segments. Nonetheless it will still give a
good
indication of the intensities across the array.
I hope that is of some help
------------------------------------------------
Mike Smith
PhD Student
Computational Biology Group
Cambridge University
On Sat, Nov 6, 2010 at 1:52 AM, Ying W <daiyingw@gmail.com> wrote:
> Hi Dr Mike Smith,
> I have a problem with the imageplot() command in the latest
beadarray but
> it only occurs when I have a .locs file in the same directory.
> > library(beadarray)
> > BLData = readIllumina(dir = "4690176049", useImages = FALSE,
> illuminaAnnotation = "Humanv3", metricsFile = "Metrics.txt")
> ....
> > imageplot(BLData, array = 1, low = "lightgreen", high =
"darkgreen",
> squareSize = 25)
> Error in res[rowIdx, ] <- resList[[i]][j, ] :
> number of items to replace is not a multiple of replacement length
> > sessionInfo()
> R version 2.12.0 (2010-10-15)
> Platform: i486-pc-linux-gnu (32-bit)
>
> locale:
> [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
> [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
> [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8
> [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
> [9] LC_ADDRESS=C LC_TELEPHONE=C
> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> other attached packages:
> [1] beadarray_2.0.2 Biobase_2.10.0
>
> loaded via a namespace (and not attached):
> [1] BeadDataPackR_1.2.0 limma_3.6.6 tools_2.12.0
>
> I was trying to troubleshoot it for a bit and at the line with
"rowIdx <-
> 1"
> > dim(res)
> [1] 359 41
> some other information:
> length(resList) == 9
> nrow(resList[[i]]) == 39
> When I looked through resList[[i]][j, ] and brings up a list of
numbers of
> length 40 or 41 (with NA or NaN in 41st slot) so I'm guessing this
could be
> an issue w/stripping NA/NaNs out (not quite sure where tho).
>
> I hope this goes to the correct Mike Smith.
>
> Best,
> Ying Wu
>
[[alternative HTML version deleted]]