[BioC] An error in oligo: All the CEL files must be of the same type??
Setsuko Sahara
sahara at salk.edu
Thu Jun 9 18:47:07 CEST 2011
Hello all,
I'm a new to bioconductor and encounter an error when I use oligo for
the analysis of affymetrix mouse Gene 1.0 ST array.
A month ago when I set up oligo with GEO data, it worked perfectly,
but started having a problem when the function of read.celfiles.
No matter the files are ( mine or GEO one that had worked fine), now I
always see an error as
"All the CEL files must be of the same type.
Error: checkChipTypes(filenames, verbose, "affymetrix", TRUE) is not
TRUE"....
So I'm wondering if anyone has any idea what would be wrong?
Here is a script.
Thank you very much for your comments in advance!
Best,
> library(pdInfoBuilder)
> baseDir <- "/Users/setsukosahara/Documents/data/Affy/oligo/
pdInfoVignette/AffyGene"
> (pgf <- list.files(baseDir, pattern = ".pgf",
+ full.names = TRUE))
[1] "/Users/setsukosahara/Documents/data/Affy/oligo/pdInfoVignette/
AffyGene/MoGene-1_0-st-v1.r4.pgf"
> (clf <- list.files(baseDir, pattern = ".clf",
+ full.names = TRUE))
[1] "/Users/setsukosahara/Documents/data/Affy/oligo/pdInfoVignette/
AffyGene/MoGene-1_0-st-v1.r4.clf"
> (prob <- list.files(baseDir, pattern = ".probeset.csv",
+ full.names = TRUE))
[1] "/Users/setsukosahara/Documents/data/Affy/oligo/pdInfoVignette/
AffyGene/MoGene-1_0-st-v1.na31.mm9.probeset.csv"
> (transcript <- list.files(baseDir, pattern = ".transcript.csv",
+ full.names = TRUE))
[1] "/Users/setsukosahara/Documents/data/Affy/oligo/pdInfoVignette/
AffyGene/MoGene-1_0-st-v1.na31.mm9.transcript.csv"
> (coreMps <- list.files(baseDir, pattern = ".mps",
+ full.names = TRUE))
[1] "/Users/setsukosahara/Documents/data/Affy/oligo/pdInfoVignette/
AffyGene/MoGene-1_0-st-v1.r4.mps"
> seed <- new("AffyGenePDInfoPkgSeed",
+ pgfFile = pgf, clfFile = clf,
+ probeFile = prob, transFile = transcript, coreMps = coreMps, author
= "setsu sahara",
+ email = "sahara at salk.edu",
+ biocViews = "AnnotationData",
+ organism = "Mouse", species = "mouse",
+ url = "http://www.salk.edu")
> makePdInfoPackage(seed, destDir = ".")
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
========================================================================
Building annotation package for Affymetrix Gene ST Array
PGF.........: MoGene-1_0-st-v1.r4.pgf
CLF.........: MoGene-1_0-st-v1.r4.clf
Probeset....: MoGene-1_0-st-v1.na31.mm9.probeset.csv
Transcript..: MoGene-1_0-st-v1.na31.mm9.transcript.csv
Core MPS....: MoGene-1_0-st-v1.r4.mps
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
========================================================================
Parsing file: MoGene-1_0-st-v1.r4.pgf... OK
Parsing file: MoGene-1_0-st-v1.r4.clf... OK
Creating initial table for probes... OK
Creating dictionaries... OK
Parsing file: MoGene-1_0-st-v1.na31.mm9.probeset.csv... OK
Parsing file: MoGene-1_0-st-v1.r4.mps... OK
Creating package in ./pd.mogene.1.0.st.v1
Inserting 44 rows into table chrom_dict... OK
Inserting 5 rows into table level_dict... OK
Inserting 8 rows into table type_dict... OK
Inserting 241576 rows into table core_mps... OK
Inserting 241576 rows into table featureSet... OK
Inserting 899636 rows into table pmfeature... OK
Counting rows in chrom_dict
Counting rows in core_mps
Counting rows in featureSet
Counting rows in level_dict
Counting rows in pmfeature
Counting rows in type_dict
Creating index idx_pmfsetid on pmfeature... OK
Creating index idx_pmfid on pmfeature... OK
Creating index idx_fsfsetid on featureSet... OK
Creating index idx_core_meta_fsetid on core_mps... OK
Creating index idx_core_fsetid on core_mps... OK
Saving DataFrame object for PM.
Saving NetAffx Annotation... OK
Done.
Warning messages:
1: In is.na(x) : is.na() applied to non-(list or vector) of type 'NULL'
2: In is.na(x) : is.na() applied to non-(list or vector) of type 'NULL'
> library(oligo)
> library(pdInfoBuilder)
> geneCELs <- list.celfiles("Users/setsukosahara/data/Affy/oligo/CEL
files/Fgf10 cell files", full.names = TRUE)
> affyGeneFS <- read.celfiles(geneCELs)
All the CEL files must be of the same type.
Error: checkChipTypes(filenames, verbose, "affymetrix", TRUE) is not
TRUE
>
Setsuko Sahara
Senior Research Associate
Molecular Neurobiology Laboratory (MNL-O)
The Salk Institute for Biological Studies
10010 North Torrey Pines Road
La Jolla, CA 92037 USA
Tel:(858) 453-4100 ext.1449
Fax:(858) 558-6207
More information about the Bioconductor
mailing list