Question: Using Plant GO Slim in topGO
0
gravatar for Sreeraj Thamban
6 months ago by
Indian Institute of Science Education and Research
Sreeraj Thamban80 wrote:

Hi I am using topGO for GO enrichment analysis, I extracted the GO terms using biomart package and it worked perfectly. Now I wanted to do the GO analysis using Plant GO slim to get a broad overview of the ontology content without the detail of the specific fine grained terms.Is there any way to extract Plant GO Slim data and use that to create the geneID2GO object. I am attaching the code below. Thanks

library("biomaRt")
mart <- biomaRt::useMart(biomart = "plants_mart",`dataset = "athaliana_eg_gene",
                     host = 'plants.ensembl.org')
#Get ensembl gene ids and GO terms
GTOGO <- biomaRt::getBM(attributes = c( "ensembl_gene_id",
                                 "go_id"), mart = mart)
head (GTOGO)
#Remove blank entries
GTOGO <- GTOGO[GTOGO$go_id != '',]
# convert from table format to list format
geneID2GO <- by(GTOGO$go_id,
            GTOGO$ensembl_gene_id,
            function(x) as.character(x))
#examine result


head (geneID2GO)
ADD COMMENTlink written 6 months ago by Sreeraj Thamban80

I have found some old script from Thomas Girke link. It looks like it is what you are looking for. For example, in the variable slimv the slim GO terms are listed.

HTH

ADD REPLYlink written 6 months ago by e.rempel540
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1183 users visited in the last hour