Reading CEL files
1
0
Entering edit mode
@libyatahani-9206
Last seen 4.9 years ago
Libyan Arab Jamahiriya

Hi all,

I have a problem during I try to read CEL files ReadAffy() , my code was:

> data<-ReadAffy()
Error in read.celfile.header(as.character(filenames[[1]])) : 
  Is C:/Users/SONY/Documents/myproject/CEL/~$M187272.CEL really a CEL file? tried reading as text, gzipped text, binary, gzipped binary, command console and gzipped command console formats

And 

> sessionInfo() 

R version 3.2.2 (2015-08-14)

Platform: x86_64-w64-mingw32/x64 (64-bit)

Running under: Windows 7 x64 (build 7601) Service Pack 1



locale:

[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252    LC_MONETARY=English_United States.1252 LC_NUMERIC=C                          

[5] LC_TIME=English_United States.1252    



attached base packages:

[1] parallel  stats     graphics  grDevices utils     datasets  methods   base     



other attached packages:

[1] affy_1.48.0          Biobase_2.30.0       BiocGenerics_0.16.1  BiocInstaller_1.20.3



loaded via a namespace (and not attached):

[1] zlibbioc_1.16.0       tools_3.2.2           affyio_1.40.0         preprocessCore_1.32.0
 
Any advice please?
Regards.
 
Tahani.H

 

cel affy • 762 views
ADD COMMENT
0
Entering edit mode

Have a look at the error message, it points to a file ~$M187272.CEL whose name looks like a temporary file generated by a Windows application. Try removing it, and generally please clean up the directory so that it only contains valid CEL files.

You can also give the names of the CEL files to ReadAffy() directly, this a more robust way to call that function.

ADD REPLY
0
Entering edit mode

Hi  Wolfgang Huber , thank you for your answer . I really did what you  did explain above , but it still giving an error !!

> data<-ReadAffy()
Error in read.celfile.header(as.character(filenames[[1]])) : 
  Is C:/Users/SONY/Documents/myproject/CEL1/GSM187268.CEL really a CEL file? tried reading as text, gzipped text, binary, gzipped binary, command console and gzipped command console formats
 

Even I put the CEL files in a separated folder and remove all other files !!!

> getwd()
[1] "C:/Users/SONY/Documents/myproject/CEL1"
> list.celfiles()
 [1] "GSM187268.CEL" "GSM187269.CEL" "GSM187271.CEL" "GSM187272.CEL" "GSM187273.CEL" "GSM187274.CEL" "GSM187275.CEL" "GSM187276.CEL" "GSM187277.CEL" "GSM187278.CEL"

 

How can I give the names of the CEL files to ReadAffy() directly according to my CEL files above , please?

Regards.

Tahani.H

 

 

ADD REPLY
1
Entering edit mode
@james-w-macdonald-5106
Last seen 9 hours ago
United States

Most likely you have one or more corrupted files, which you could figure out by trying to read each one individually. But do note that the authors who submitted those data used RMA in R, so unless you want to do some other normalization/summarization (and there is very little profit in doing so), you can just use getGEO in the GEOquery package to get the summarized data directly, without having to bother with the raw data.

ADD COMMENT

Login before adding your answer.

Traffic: 426 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6