In addition to using acyclic graphs,
ontologyPlot lets you to visualise ontological annotation as grids. In order to do this for a given a
list of term sets (i.e. your annotation), one calls the function
plot_annotation_grid, passing parameters:
ontology, i.e. an
term_sets, the annotation list,
all_terms, a character vector of term IDs indicating which terms to include in the plot. This defaults to those terms which show only the most informative common ancestors for all subsets of annotated objects (otherwise there can be an ungainly number of columns).
We proceed by loading the package, and an example ontology - in our case, the Gene Ontology (GO). We then create a
list containing our annotation - here we use a list of 6 genes with the actual annotation downloaded from http://geneontology.org/ - and subsequently pass it to
library(ontologyIndex) library(ontologyPlot) data(go) genes <- list( A0A087WUJ7=c("GO:0004553", "GO:0005975"), CTAGE8=c("GO:0016021"), IFRD2=c("GO:0003674", "GO:0005515", "GO:0005634"), OTOR=c("GO:0001502", "GO:0005576", "GO:0007605"), TAMM41=c("GO:0004605", "GO:0016024", "GO:0031314", "GO:0032049"), ZZEF1=c("GO:0005509", "GO:0008270") ) plot_annotation_grid(go, term_sets=genes)
One can either use
annotation_grid to get a
logical matrix of term-inclusion, or plot directly using
?plot_annotation_grid for more details.