[BioC] Unable to load in Affymetrix data CEL files with readCelHeader
rgentlem at fhcrc.org
Tue May 13 21:31:03 CEST 2008
Matt K wrote:
> I am having problems reading in some publicly available chromosome X
> titration Nsp chip CEL files. The data are available from the Affymetrix
I did not see any file there that obviously contained the CEL files,
can you say what you downloaded (provided the questions below don't
solve your problem).
> I have not modified the data in anyway. Here is what happens when I try to
> read the data:
>> path <- "./rawData/3X/Mapping250K_Nsp/"
>> pathnames <- list.files(path=path, pattern="[.](cel|CEL)$",
>  "./rawData/3X/Mapping250K_Nsp//NA04626_NSP_R1.CEL"
>  "./rawData/3X/Mapping250K_Nsp//NA04626_NSP_R2.CEL"
>  "./rawData/3X/Mapping250K_Nsp//NA04626_NSP_R3.CEL"
>  "./rawData/3X/Mapping250K_Nsp//NA04626_NSP_R4.CEL"
>> hdr <- readCelHeader(pathnames)
> terminate called after throwing an instance of
That suggests that you may not have read permission on them. Did you
check and see if you could open those files with any other tool/editor?
You could just try and open them and read them using standard R tools
to see if that works (if they are binary CEL files then you will just
get junk from readLines, but that isn't the issue, you just want to know
if they can be opened from R).
Your R is out of date, and this sessionInfo is not correct, as you
should have had affxparser attached and it is not. Please don't do
that, mixing and matching error messages and sessionInfo output makes
life hard for anyone that wants to help. Step one is to update R and
> Process R aborted at Tue May 13 15:01:55 2008
> As you see R aborts. The same failure happens when I try to load in any of
> the other CEL files. My R session info is:
> R version 2.7.0 Under development (unstable) (2008-01-21 r44087)
> attached base packages:
>  stats graphics grDevices utils datasets methods base
> Thanks for any help.
> [[alternative HTML version deleted]]
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
Robert Gentleman, PhD
Program in Computational Biology
Division of Public Health Sciences
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N, M2-B876
PO Box 19024
Seattle, Washington 98109-1024
rgentlem at fhcrc.org
More information about the Bioconductor