Perform gene set overrepresentation analysis for each bicluster and annotate cells based on the best match.

per_cluster_goa loads the required gene set, formats it and performs gene overrepresentation analysis for each bicluster in the cabic object.

Usage

per_cluster_goa(
  cabic,
  universe,
  org,
  set = "CellMarker",
  min_size = 10,
  max_size = 500,
  verbose = TRUE
)

Arguments

cabic

A biclustering object of class "caclust" as obtained from `caclust`.

universe

All genes in data set.

org

Short name of the organism. "mm" for mouse, "hs" for human.

set

Name of the gene set. Currently only supports "CellMarker"

min_size

Min. number of genes in the gene set.

max_size

Max number of genes in the gene set. Set to Inf if you want to keep all genes.

verbose

Toggles verbosity of warnings.

Value

A list contain the goa results for each cluster.