Error in get_anno_data(ontology) : ontology not supported yet...
1
0
Entering edit mode
GODseeker • 0
@ae7374b1
Last seen 11 hours ago
China

kk <- enrichDO(gene =gene,ont = "DO", pvalueCutoff = 0.5, qvalueCutoff = 0.5,readable = T) Error in get_anno_data(ontology) : ontology not supported yet...

when i enrichDO but it says that...i have successfully installed this package but when i running the code ,it always appear .anyone could help me please...

enter image description here

EnrichDO • 30 views
ADD COMMENT
0
Entering edit mode
Guido Hooiveld ★ 4.1k
@guido-hooiveld-2020
Last seen 8 hours ago
Wageningen University, Wageningen, the …

Did you check the help page for the function enrichDO? Type ?enrichDO.

I am asking, because there it is stated that the argument ont should either be one of "HDO", "HPO" or "MPO". So "DO" (what you used) is indeed not a supported ontology!

BTW: HDO = Human Disease Ontology (link), HPO = Human Phenotype Ontology (link), MPO = Mouse Phenotype Ontology (link) )

Then it is working! Using the human example dataset (thus with the Human Disease Ontology [HDO]):

> library(DOSE)

DOSE v4.0.0 Learn more at https://yulab-smu.top/contribution-knowledge-mining/

Please cite:

Guangchuang Yu, Li-Gen Wang, Guang-Rong Yan, Qing-Yu He. DOSE: an
R/Bioconductor package for Disease Ontology Semantic and Enrichment
analysis. Bioinformatics. 2015, 31(4):608-609
> data(geneList)
> 
> gene = names(geneList)[geneList > 1]
> 
> kk <- enrichDO(gene = gene,
+                ont = "HDO",
+                pvalueCutoff = 0.5,
+                qvalueCutoff = 0.5,
+                readable = T) 

> 
> kk
#
# over-representation test
#
#...@organism    Homo sapiens 
#...@ontology    HDO 
#...@keytype     ENTREZID 
#...@gene        chr [1:525] "4312" "8318" "10874" "55143" "55388" "991" "6280" "2305" ...
#...pvalues adjusted by 'BH' with cutoff <0.5 
#...352 enriched terms found
'data.frame':   352 obs. of  12 variables:
 $ ID            : chr  "DOID:0060306" "DOID:934" "DOID:0080599" "DOID:11054" ...
 $ Description   : chr  "Meier-Gorlin syndrome" "viral infectious disease" "Coronavirus infectious disease" "urinary bladder cancer" ...
 $ GeneRatio     : chr  "8/302" "36/302" "20/302" "23/302" ...
 $ BgRatio       : chr  "10/7865" "336/7865" "138/7865" "180/7865" ...
 $ RichFactor    : num  0.8 0.107 0.145 0.128 0.087 ...
 $ FoldEnrichment: num  20.83 2.79 3.77 3.33 2.26 ...
 $ zScore        : num  12.54 6.7 6.57 6.31 5.44 ...
 $ pvalue        : num  1.82e-10 1.65e-08 2.58e-07 3.36e-07 1.52e-06 ...
 $ p.adjust      : num  1.43e-07 6.46e-06 6.59e-05 6.59e-05 1.63e-04 ...
 $ qvalue        : num  1.15e-07 5.20e-06 5.31e-05 5.31e-05 1.31e-04 ...
 $ geneID        : chr  "CDC45/CDT1/MCM5/CDC6/ORC6/ORC1/GINS3/GMNN" "MMP1/S100A8/CXCL10/CXCL9/GZMB/IDO1/CCL8/CXCL3/CCL20/CXCL8/CTLA4/CCL5/MKI67/LCN2/ACE2/MX1/TREM1/LCK/C1QB/PRF1/IL"| __truncated__ "MMP1/S100A8/CXCL10/CXCL9/GZMB/CCL8/CXCL3/CXCL8/CCL5/LCN2/ACE2/C1QB/IL2RA/CCR1/STAT1/CCL2/IL2RG/CCL4/SLPI/IL7" "MMP1/S100A9/S100A8/TOP2A/RRM2/CDK1/BIRC5/EZH2/AURKA/CDC25A/PLK1/MMP9/CCNE1/MKI67/FANCA/CA9/CDC25B/TYMS/MYC/E2F1"| __truncated__ ...
 $ Count         : int  8 36 20 23 38 38 9 34 14 9 ...
#...Citation
Guangchuang Yu, Li-Gen Wang, Guang-Rong Yan, Qing-Yu He. DOSE: an R/Bioconductor package for Disease Ontology Semantic and Enrichment analysis. Bioinformatics. 2015, 31(4):608-609 

> 
> as.data.frame(kk)[1:3,]
                       ID                    Description GeneRatio  BgRatio
DOID:0060306 DOID:0060306          Meier-Gorlin syndrome     8/302  10/7865
DOID:934         DOID:934       viral infectious disease    36/302 336/7865
DOID:0080599 DOID:0080599 Coronavirus infectious disease    20/302 138/7865
             RichFactor FoldEnrichment    zScore       pvalue     p.adjust
DOID:0060306  0.8000000      20.834437 12.540792 1.816539e-10 1.425983e-07
DOID:934      0.1071429       2.790326  6.702107 1.645690e-08 6.459332e-06
DOID:0080599  0.1449275       3.774355  6.570134 2.583957e-07 6.594132e-05
                   qvalue
DOID:0060306 1.147288e-07
DOID:934     5.196914e-06
DOID:0080599 5.305370e-05
                                                                                                                                                                                                           geneID
DOID:0060306                                                                                                                                                            CDC45/CDT1/MCM5/CDC6/ORC6/ORC1/GINS3/GMNN
DOID:934     MMP1/S100A8/CXCL10/CXCL9/GZMB/IDO1/CCL8/CXCL3/CCL20/CXCL8/CTLA4/CCL5/MKI67/LCN2/ACE2/MX1/TREM1/LCK/C1QB/PRF1/IL2RA/CA9/CCR1/STAT1/TLR8/CCL2/PLA2G2A/MYC/IL2RG/RSAD2/CCL4/SLPI/CD163/CDKN2A/CEBPB/IL7
DOID:0080599                                                                                         MMP1/S100A8/CXCL10/CXCL9/GZMB/CCL8/CXCL3/CXCL8/CCL5/LCN2/ACE2/C1QB/IL2RA/CCR1/STAT1/CCL2/IL2RG/CCL4/SLPI/IL7
             Count
DOID:0060306     8
DOID:934        36
DOID:0080599    20
> 
> 
> 
ADD COMMENT

Login before adding your answer.

Traffic: 690 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