Question: Reading CEL files
0
gravatar for libya.tahani
3.2 years ago by
Libyan Arab Jamahiriya
libya.tahani0 wrote:

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

 

affy cel • 621 views
ADD COMMENTlink modified 3.2 years ago by James W. MacDonald52k • written 3.2 years ago by libya.tahani0

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 REPLYlink written 3.2 years ago by Wolfgang Huber13k

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 REPLYlink written 3.2 years ago by libya.tahani0
Answer: Reading CEL files
1
gravatar for James W. MacDonald
3.2 years ago by
United States
James W. MacDonald52k wrote:

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 COMMENTlink written 3.2 years ago by James W. MacDonald52k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 157 users visited in the last hour