[BioC] GEOquery error when parsing Agilent annotation (?)

axel.klenk at actelion.com axel.klenk at actelion.com
Tue May 25 17:51:51 CEST 2010


Dear BioC list and esp. Sean,

I'm trying to download GSE15197 using GEOquery and get the following 
error:

> gse15197 <- getGEO(GEO = "GSE15197", destdir = ".", AnnotGPL = TRUE)
Found 1 file(s)
GSE15197_series_matrix.txt.gz
trying URL 
'ftp://ftp.ncbi.nih.gov/pub/geo/DATA/SeriesMatrix/GSE15197/GSE15197_series_matrix.txt.gz'
ftp data connection made, file length 6566671 bytes
opened URL
==================================================
downloaded 6.3 Mb

File stored at: 
/tmp/RtmpDBAlMH/GPL6480.soft
Error in `row.names<-.data.frame`(`*tmp*`, value = c("A_23_P100001", 
"A_23_P100011",  : 
  missing values in 'row.names' are not allowed
> traceback()
7: stop("missing values in 'row.names' are not allowed")
6: `row.names<-.data.frame`(`*tmp*`, value = c("A_23_P100001", 
"A_23_P100011", 
   "A_23_P100022", "A_23_P100056", "A_23_P100074", "A_23_P100092", 

[... lots of Agilent oligo IDs ...]

   "ETG10_13482", "ETG10_195139", "ETG10_234183", "ETG10_236652", 
   "GE_BrightCorner", NA, "RC1", "RC10", "RC11", "RC12", "RC2", 
   "RC3", "RC4", "RC5", "RC6", "RC7", "RC8", "RC9"))
5: `row.names<-`(`*tmp*`, value = c("A_23_P100001", "A_23_P100011", 
   "A_23_P100022", "A_23_P100056", "A_23_P100074", "A_23_P100092", 

[... lots of Agilent oligo IDs ...]

   "ETG10_13482", "ETG10_195139", "ETG10_234183", "ETG10_236652", 
   "GE_BrightCorner", NA, "RC1", "RC10", "RC11", "RC12", "RC2", 
   "RC3", "RC4", "RC5", "RC6", "RC7", "RC8", "RC9"))
4: `rownames<-`(`*tmp*`, value = c("A_23_P100001", "A_23_P100011", 
   "A_23_P100022", "A_23_P100056", "A_23_P100074", "A_23_P100092", 

[... lots of Agilent oligo IDs ...]

   "ETG10_13482", "ETG10_195139", "ETG10_234183", "ETG10_236652", 
   "GE_BrightCorner", NA, "RC1", "RC10", "RC11", "RC12", "RC2", 
   "RC3", "RC4", "RC5", "RC6", "RC7", "RC8", "RC9"))
3: parseGSEMatrix(file.path(tmp, b[i]))
2: getAndParseGSEMatrices(GEO)
1: getGEO(GEO = "GSE15197", destdir = ".", AnnotGPL = TRUE)


which apparently occurs after successful download of the data files when
trying to parse GPL6480 annotation (?) and according to the error message
the NA row.name near the end of the value list may be the culprit. Setting
AnnotGPL to FALSE gives the same result. 

Since the identical procedure usually works, is there a problem with this
particular data set or am I missing something?

Thanks a lot in advance,

 - axel


> sessionInfo()
R version 2.11.0 (2010-04-22) 
x86_64-pc-linux-gnu 

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C 
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8 
 [5] LC_MONETARY=C              LC_MESSAGES=en_US.UTF-8 
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C 
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C 

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base 

other attached packages:
 [1] GEOquery_2.12.0      RCurl_1.4-2          bitops_1.0-4.1 
 [4] limma_3.4.0          sfsmisc_1.0-11       geneplotter_1.26.0 
 [7] lattice_0.18-8       annotate_1.26.0      AnnotationDbi_1.10.0
[10] Biobase_2.8.0        MASS_7.3-5           gmodels_2.15.0 
[13] fortunes_1.3-7 

loaded via a namespace (and not attached):
[1] DBI_0.2-5          gdata_2.8.0        grid_2.11.0        gtools_2.6.2  
 
[5] RColorBrewer_1.0-2 RSQLite_0.9-0      xtable_1.5-6 


Axel Klenk
Research Informatician
Actelion Pharmaceuticals Ltd / Gewerbestrasse 16 / CH-4123 Allschwil / 
Switzerland



The information of this email and in any file transmitted with it is strictly confidential and may be legally privileged.
It is intended solely for the addressee. If you are not the intended recipient, any copying, distribution or any other use of this email is prohibited and may be unlawful. In such case, you should please notify the sender immediately and destroy this email.
The content of this email is not legally binding unless confirmed by letter.
Any views expressed in this message are those of the individual sender, except where the message states otherwise and the sender is authorised to state them to be the views of the sender's company. For further information about Actelion please see our website at http://www.actelion.com



More information about the Bioconductor mailing list