[BioC] GEOquery and parsing SOFT files

Wacek Kusnierczyk Waclaw.Marcin.Kusnierczyk at idi.ntnu.no
Mon May 25 23:00:08 CEST 2009

Wolfgang Huber wrote:
> Dear Wacek,
> thank you for the feedback and pointing this out. Two general remarks:

thank you for the response.

> 1. Please include a reproducible example (R script) for others to
> reproduce your experience, and subsequently the output of sessionInfo().

a reproducible (on my machine, at least) example is as follows:

    repeat {
       filename = getGEOfile('GSE13638')
       if (!is.null(filename)) break }
    system.time({gse13638 = getGEO(filename=filename)})
    # elapsed: 620

(so it was 'faster' this time, only 10 minutes...)

session info is as follows:

R version 2.10.0 Under development (unstable) (2009-05-25 r48607)

 [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           

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

other attached packages:
[1] GEOquery_2.9.0 RCurl_0.94-1   Biobase_2.5.2

> 2. Robert Gentleman's book "R Programming for Bioinformatics" (as well
> as many free sources on the web) describes how to profile R code in
> order to see in which functions the CPU time is spent. Based on this,
> you can investigate where to invest developer time for improving the
> code.

yes, i know;  but i thought someone could already have done this or have
a hint anyway, as this does not seem likely to be just my local
problem.  before i spend too much time on a potentially well-known
issue, it's surely fine to ask others, including the developer.

best regards,

More information about the Bioconductor mailing list