[BioC] Illumina beadlevel expression data

Ina Hoeschele inah at vt.edu
Wed Nov 14 23:24:16 CET 2007


Hi,
   I am new to the analysis of Illumina beadlevel expression data. I 
have data from a collaborator and have started to work with the data on 
the first Beadchip (with 12 arrays or strips).  These are single channel 
data. I am using the beadarray package.  I read the data in as follows:
BLData <- 
readIllumina(textType=".txt",useImages=TRUE,singleChannel=TRUE,imageManipulation="sharpen",
backgroundSize=17,storeXY=TRUE,metrics=FALSE,backgroundMethod="none",offset=0,normalizeMethod="none")
Note: above I'm doing NO background correction and note that I use 
backgroundSize=17.

Then I looked at the data:

an <- arrayNames(BLData)
BLData at beadData[[an[1]]]$G[1:10]
 [1]  0.0000000  2.5249867 -0.0359710  0.0000000  0.0000000 -6.8111267
 [7]  0.0000000 77.2345222 -0.8237067 -0.2324444
BLData at beadData[[an[1]]]$Gb[1:10]
 [1] 0 0 0 0 0 0 0 0 0 0

so I get quite a few G intensities that are zero and ALL of my 
background intensities are zero - why is that ?????

Then I created BLData with backgroundSize=4 and backgroundSize=40 (not 
knowing what value I should be using ?????). For backgroundSize=4, I 
find that
BLData at beadData[[an[1]]]$Gb[1:100]
  [1] -1  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  
0  0  0
 [26]  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  2  0  
0  0  0
 [51]  0  0 24  2  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  
0  0  0
 [76]  0  0  0  0  0  0  0  0  0  0 19  0  0  0  0  0  0  0  0  0  0  0  
0  0  0
so now I still have many zero Gb values, but also some nonzero ones and 
even negative (-1)???
But when I use backgroundSize=40, then R crashes (repeatedly).

And when I create BLData with backgroundMethod="minimum" and 
backgroundSize=17, then I find the following G intensities:
BLData at beadData[[an[1]]]$G[1:10]
 [1] 2.425286e-12 2.524987e+00 2.425286e-12 2.425286e-12 2.425286e-12 
2.425286e-12 2.425286e-12 7.723452e+01 2.425286e-12
[10] 2.425286e-12

These G and Gb values all look suspect to me - can someone please help?

Many thanks, Ina



More information about the Bioconductor mailing list