ClusterProfiler KEGG Gene Set Enrichment Analysis
1
0
Entering edit mode
iidaL • 0
@b9e47314
Last seen 3.2 years ago
Finland

Dear all,

I try to use ClusterProfiler for Enrichment Analysis (KEGG). I face a problem right I try to create gseKEGG object. In the instructions there is said "organism KEGG Organism Code: The full list is here: https://www.genome.jp/kegg/catalog/org_list.html (need the 3 letter code)"

My problem is that my organism has four letter code "pfre" (Propionibacterium freudenreichii subsp. freudenreichii DSM 20271).

Is there anything what I could do to have the analysis for my organism?

kegg_organism = "pfre"
kk2 <- gseKEGG(geneList     = gene_list,
               organism     = pfre,
               nPerm        = 10000,
               minGSSize    = 3,
               maxGSSize    = 800,
               pvalueCutoff = 0.05,
               pAdjustMethod = "none",
               keyType       = "kegg")


Error in organismMapper(organism) : object 'pfre' not found
clusterProfiler KEGGGeneSetEnrichmentAnalysis • 2.0k views
ADD COMMENT
1
Entering edit mode
shepherl 3.8k
@lshep
Last seen 3 hours ago
United States

When you ran the function you used pfre without quotes or using your defined kegg_organism ... could you try running

kegg_organism = "pfre"
kk2 <- gseKEGG(geneList     = gene_list,
               organism     = "pfre",    # or  organism = kegg_organism 
               nPerm        = 10000,
               minGSSize    = 3,
               maxGSSize    = 800,
               pvalueCutoff = 0.05,
               pAdjustMethod = "none",
               keyType       = "kegg")
ADD COMMENT
0
Entering edit mode

Thank you. :)

ADD REPLY

Login before adding your answer.

Traffic: 772 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6