Search
Question: enrichGO from DOSE, clusterprofiler is returning empty result.
1
gravatar for kirannbishwa01
11 months ago by
kirannbishwa010 wrote:

I am trying to use enrichGO function from clusterProfiler package. All is good except that my results are empty.

# for significant genes

> sig.foldChangesAbv2 %>% length()
[1] 244

> sig.foldChanges %>% head(5)
[1] "NA1"       "AT1G64640" "AT1G64630" "AT1G64625" "NA3"

# for universal genes

> length(names(foldChanges))
[1] 1761

> foldChanges %>% head(5)
       NA1  AT1G64720  AT1G64710  AT1G64700  AT1G64690
-2.2955229  0.6501073  0.7451122  0.6378277 -0.9689359

 

ego <- enrichGO(gene = sig.foldChangesAbv2.5 ,
                universe = names(foldChanges),
                keyType = "TAIR",
                OrgDb = org.At.tair.db,
                ont = "BP",
                pAdjustMethod = "BH",
                pvalueCutoff  = 0.1,
                qvalueCutoff  = 0.5,
                readable = TRUE)

> head(ego)
[1] ID          Description GeneRatio   BgRatio     pvalue      p.adjust   
[7] qvalue      Count      
<0 rows> (or 0-length row.names)

Why is my result just empty? Is it because there are very few genes belonging to any GO category? or is the problem because I am using `Arabidopsis` database?

 

The `GO Classification` is working totally fine though.

> ggo <- groupGO(gene = names(foldChanges),
+                OrgDb = org.At.tair.db, keytype = "TAIR", ont = "BP",
+                level = 3, readable = TRUE)

> ggo
GO BP Profiles at level 3 of 1761 Arabidopsis thaliana genes


> head(ggo, 5)
                   ID                              Description Count GeneRatio
GO:0019953 GO:0019953                      sexual reproduction    20   20/1761
GO:0019954 GO:0019954                     asexual reproduction     0    0/1761
GO:0022414 GO:0022414                     reproductive process   147  147/1761

 

 

 

ADD COMMENTlink modified 11 months ago by Guangchuang Yu1.1k • written 11 months ago by kirannbishwa010
1
gravatar for Guangchuang Yu
11 months ago by
Guangchuang Yu1.1k
China/Guangzhou/Southern Medical University
Guangchuang Yu1.1k wrote:
ADD COMMENTlink written 11 months ago by Guangchuang Yu1.1k
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.2.0
Traffic: 242 users visited in the last hour