[BioC] Illumina bead arrays

Ina Hoeschele inah at vt.edu
Mon Aug 13 16:46:24 CEST 2007


Hi,
   I am still struggling with my first try at working with Illumina bead 
array expression data, specifically getting the data read by the R 
package beadarray.

The data was created by an off-site collaborator using BeadStudio 
version  3.0.14. My collaborator sent two types of files to me:
(1) a number of .csv files (read by excel, 2 for each sample) with 
columns Illumicode, N, mean GRN, Dev GRN.
(2) a summary .txt file with columns TargetID  
MIN_Signal-1814647013_A    AVG_Signal-1814647013_A  
MAX_Signal-1814647013_A    NARRAYS-1814647013_A  
ARRAY_STDEV-1814647013_A    BEAD_STDEV-1814647013_A  
Avg_NBEADS-1814647013_A    Detection-1814647013_A, with the 
MIN_Signal-... to Detection-... columns repeated for the remaining samples.

I read the summary data(2) with the statement
readBeadSummaryData(dataFile, skip=8, columns = list(exprs = 
"AVG_Signal", BeadStDev = "BEAD_STDEV", NoBeads = "Avg_NBEADS"), sep="\t"),
but every time I do this R gets hung up (no error message).

I also tried reading the bead-level data (the .csv files) with the statement
BLData <- readIllumina(textType = ".csv")
but then I get the error message
Error in strtrim(x, width) : invalid 'width' argument

Any suggestions would be much appreciated.
Thank you.
Ina



More information about the Bioconductor mailing list