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()
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.

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()
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

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.