R/plot.owlents.R

In vjcitn/ontoProc: processing of ontologies of anatomy, cell lines, and so on

#' visualize ontology selection via onto_plot2, based on owlents
#' @rawNamespace S3method(plot,owlents)
#' @param x owlents instance
#' @param y character() vector of entries in x$clnames
#' @param \dots passed to onto_plot2
#' @param dropThing logical(1) defaults to TRUE; if "Thing" is
#' present in terms to display, it is removed
#' @examples
#' cl3k = c("CL:0000492", "CL:0001054", "CL:0000236", 
#'   "CL:0000625", "CL:0000576", 
#'   "CL:0000623", "CL:0000451", "CL:0000556")
#' cl3k = gsub(":", "_", cl3k)
#' clont_path = owl2cache(url="http://purl.obolibrary.org/obo/cl.owl")
#' o2 = try(reticulate::import("owlready2"), silent=TRUE)
#' if (!inherits(o2, "try-error")) {
#'  clont = setup_entities(clont_path)
#'  plot(clont,cl3k)
#' }
#' @export
plot.owlents = function(x, y, ..., dropThing=TRUE) {
  if (missing(y)) stop("must provide classes for plotting as 'y'")
  oky = intersect(y, x$clnames)
  if (length(oky)<2) stop("not enough selected names present in x$clnames")
  litc = x[oky]
  alitc = ancestors(litc)
  anames = ancestors_names(alitc)
  alluse = unique(unlist(anames))
  if (dropThing) alluse = setdiff(alluse, "Thing")
  red = x[alluse]
  suppressWarnings({  # will complain if Thing or other un-parented term is present
    myo = ontologyIndex::ontology_index(parents=parents(red),
         name=labels(red))
    })
  if (length(unlist(myo$parents))==0) stop("terms chosen lack parents in specified ontology")
  onto_plot2(myo, oky, ...)
}