[BioC] GO question

[Gustavo Fernández Bayón]

Hi Cristobal,  

thank you very much for the answer. I'll write it down in case my current workflow decides to work no more. For now, it does work just by loading the GOstats library inside the scope of the inner foreach. Maybe I have the multithreaded version of SQlite, I don't know.  

I was wondering why Gostats seems too slow when compared with the DAVID web tool. Is it just a matter of hardware (I do not know what is running at DAVID's backyard), or are there more efficient implementations? Is topGO a more efficient alternative? I currently have more than 100 groups of genes on which I want to do a GO analysis, that is why I am experimenting with parallel computing for it.

Thank you again for your answer.

Regards,
Gus


---------------------------
Enviado con Sparrow (http://www.sparrowmailapp.com/?sig)


El jueves 15 de noviembre de 2012 a las 16:55, Cristobal Fresno Rodríguez escribió:

> HI Gus,
>  
> The same problem here but using parallel. The problem lies in sqlite threadsafe mode (Single-thread, Multi-thread, Serialized). As far as I know, in windows the default binary comes with Serialized (thread safe) and in unix no, thus you have to compile it from the source. But, if you are using fork to parallelize, as in parallel or multicore library it stills breaks the database conection. I don't know if for foreach works. Maybe you should give it a try.  
>  
> At present, the workaround I am using is to manualy split HyperGTest into two functions: one to access annotation packages and other for the actual hypergeometric test. So, in the code first secuentially access all the annotation package/s to build the GO graphs and then run the tests in parallel. This workaround is pretty much what you have been doing so far.  
>  
> Regards,
>  
> Cristobal
>  
>  
> 2012/11/15 Gustavo Fernández Bayón <gbayon at gmail.com (mailto:gbayon at gmail.com)>
> > Hi everybody.
> >  
> > A simple question: is there any way I can perform a GO enrichment analysis without using the annotation packages?
> >  
> > Problem is, I am trying to perform a series of GO analyses in parallel (with foreach), and I am experiencing problems with every call trying to access the same SqLite database. For now, I have solved it by putting "library(GOstats)" inside the inner foreach, but I was wondering if there is a better way.
> >  
> > Regards,
> > Gus
> >  
> > ---------------------------
> > Enviado con Sparrow (http://www.sparrowmailapp.com/?sig)
> >  
> > _______________________________________________
> > Bioconductor mailing list
> > Bioconductor at r-project.org (mailto:Bioconductor at r-project.org)
> > https://stat.ethz.ch/mailman/listinfo/bioconductor
> > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
>