Entering edit mode
7.0 years ago
12021560-002 ▴ 30
I am reading microarray data (CEL files) through R. While reading GSE35021 based on platform Affymetrix Human Exon 1.0 ST Array [CDF: HuEx-1_0-st-v2,coreR3,A20071112,EP]. I am getting this error when I give command to read these files. let me clear I am using oligo pkg to read celfiles for this platform. Could please help how to remove this error? Is there is some other pkg to read files for this package? if yes then kindly share the name of pkg and command as well.
rawData <- read.celfiles(celfiles) Error in read.celfiles.header (x): Is GSE35021/CEL/GPL9793_gene_assignment
GSE35021/CEL/GPL9793_gene_assignmentis not a CEL file. It should end by .CEL or .CEL.gz
this is the script I am using. Am i giving write commands please pick up the error
gse30521dat = getGEO('GSE30521')
tarfile<- grep(".tar$", rownames(gse30521raw), value = TRUE)
untar("GSE30521/GSE30521_RAW.tar", exdir = "GSE30521/CEL")
celfiles<- list.files("GSE30521/CEL", full = TRUE)
Sometimes you need to see/read the error properly. Apart from *.CEL.gz files, your folder has 3 more files:
Either move them to other folder or do
what does 4:26 mean?
It is vector index in R. It means pick the elements from 4 to 26. I would highly recommend you to go thru some R programming basics before you start with microarray analysis.