[BioC] Problem reading Cel files - Oligo Package

Atul atulkakrana at outlook.com
Wed Aug 28 05:18:13 CEST 2013


Hi All,

I am trying to read four *.Cel files into oligo and getting this error:

 > celFiles <- list.celfiles()
 > celFiles
[1] "Iris.CEL" "Liv1.CEL" "Liv2.CEL" "Liv3.CEL"
 > AF_data = read.celfiles(celFiles)
All the CEL files must be of the same type.
Error: checkChipTypes(filenames, verbose, "affymetrix", TRUE) is not TRUE

Then I tried reading files separately (one by one) and found that one 
sample (Iris.CEL) shows annotation package as 'pd.huex.1.0.st.v1' while 
rest (Liv1,Liv2,Liv3) are 'pd.huex.1.0.st.v2'. I checked on GEO and 
found that though all the samples are from different studies but were 
generated using same chip - Human Exon 1.0 ST Arrays and the one which 
is giving error (Iris.cel )have 'HuEx-1_0-st-v2.r2.dt1.hg18.core.ps' 
mentioned under data processing description, that means it is also 
version2 of HuEx 1.0ST.

So I explicitly mentioned annotation package 'pd.huex.1.0.st.v2' instead 
of the one recognized by oligo ('pd.huex.1.0.st.v1') and file is read 
without any problem:

 > celFiles <- list.celfiles()
 > celFiles
[1] "Iris.CEL"
 > AF_data = read.celfiles(celFiles,pkgname='pd.huex.1.0.st.v2')
Platform design info loaded.
Reading in : Iris.CEL

But if I add other files and try same thing, than the error is back:
 > celFiles <- list.celfiles()
 > celFiles
[1] "Iris.CEL" "Liv1.CEL" "Liv2.CEL" "Liv3.CEL"
 > AF_data = read.celfiles(celFiles,pkgname='pd.huex.1.0.st.v2')
All the CEL files must be of the same type.
Error: checkChipTypes(filenames, verbose, "affymetrix", TRUE) is not TRUE


Can anybody please tell me why annotation package for Iris.cel which is 
from HuEx 1.0ST v2 (from NCBI GEO description) is recognized as 
'pd.huex.1.0.st.v1'? If explicitly mention package name 
pd.huex.1.0.st.v2 and try to read Iris.cel alone, it works. But if read 
with other cel files with same annotation (pd.huex.1.0.st.v2) it gives 
error??

NCBI GEO ID:
Iris.cel - GSM1008547
Liv1/2/3 - GSM486433/GSM486434/GSM486435

Awaiting help.

AK


Session Info:

 > sessionInfo()
R version 3.0.1 (2013-05-16)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C LC_TIME=en_US.UTF-8        
LC_COLLATE=en_US.UTF-8 LC_MONETARY=en_US.UTF-8
  [6] LC_MESSAGES=en_US.UTF-8    LC_PAPER=C LC_NAME=C                  
LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets methods   
base

other attached packages:
[1] pd.huex.1.0.st.v2_3.8.0 RSQLite_0.11.4 DBI_0.2-7               
oligo_1.24.2 Biobase_2.20.1          oligoClasses_1.22.0
[7] BiocGenerics_0.6.0

loaded via a namespace (and not attached):
  [1] affxparser_1.32.3     affyio_1.28.0 BiocInstaller_1.10.1  
Biostrings_2.28.0     bit_1.1-10 codetools_0.2-8
  [7] ff_2.2-11             foreach_1.4.0 GenomicRanges_1.12.4  
IRanges_1.18.1        iterators_1.0.6 preprocessCore_1.22.0
[13] splines_3.0.1         stats4_3.0.1 tools_3.0.1           
zlibbioc_1.6.0



More information about the Bioconductor mailing list