Search
Question: enrichGO from DOSE, clusterprofiler is returning empty result.
1
gravatar for kirannbishwa01
27 days 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 27 days ago by Guangchuang Yu800 • written 27 days ago by kirannbishwa010
1
gravatar for Guangchuang Yu
27 days ago by
Hong Kong
Guangchuang Yu800 wrote:
ADD COMMENTlink written 27 days ago by Guangchuang Yu800
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: 419 users visited in the last hour