[BioC] PANTHER in R.... Gene lists in sigPathway?
Ovokeraye Achinike-Oduaran
ovokeraye at gmail.com
Wed May 23 09:09:35 CEST 2012
Thanks a bunch, Felipe.
Regards,
Avoks
On Tue, May 22, 2012 at 6:45 PM, Felipe Riveroll Aguirre
<friveroll at gmail.com> wrote:
> May be you can get the limma logFC associated values, for associated
> GWAS genes, in the example from SPIA, I started only with entrezID's
> from KEGG link, and grab the logFC values from limma in sig_genes
> vector
>
> sig_genes <- subset(top, adj.P.Val<0.01)$logFC
> names(sig_genes) <- subset(top, adj.P.Val<0.01)$ENTREZ
>
> with this code
>
> genes.both <- sig_genes[which(names(sig_genes) %in% genes.via.sel)]
>
> Where genes.via.sel, was a gene list set of entrezID's to get limma
> associated values.
>
> So you can do an approach like this to get associated logFC values
> from associated GWAS genes
>
> Also there is another server I known for gene lists
>
> http://www.bioprofiling.de/R_spider.html
>
> Here you get an analysis combining signaling and metabolic pathways
> from Reactome and KEGG databases.
>
> So, R only helps to get the required inputs, for external web appz, or
> you can download pantherdb database and write an script for automate
> the whole process, from R output
>
> Best,
>
> Felipe
>
> 2012/5/22 Ovokeraye Achinike-Oduaran <ovokeraye at gmail.com>:
>>
>> Hi Felipe,
>>
>> Thanks. The thing is that I have genes from both GWAS and DE analysis from Limma in one list. So they don't all have logFC values, hence my question about gene lists and not necessarily gene sets. I ran my list through the PANTHER web interface and got some results but I just wanted to know if it had an R implementation.
>>
>> DAVIDQuery() might work if it gives the same options as is web interface of choosing the PANTHER pathway option but I didn't find the instruction manual very straightforward.
>>
>> Thanks again.
>>
>> -Avoks
>>
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list