Skip to contents

Run hypergeometric enrichment a list of target genes over the gene ontology

rs_gse_geom_elim_list.Rd

[Experimental] This function implements a Rust version of the gene ontology enrichment with elimination: the starting point are the leaves of the ontology and hypergeometric tests will first be conducted there. Should the hypergeometric test p-value be below a certain threshold, the genes of that gene ontology term will be removed from all ancestors. This function is designed to leverage Rust-based threading for parallel processing of a list of target genes.

Usage

rs_gse_geom_elim_list(
  target_genes_list,
  levels,
  go_obj,
  gene_universe_length,
  min_genes,
  elim_threshold,
  min_overlap,
  fdr_threshold
)

Arguments

target_genes_list

A list of target genes against which to run the method.

levels

A character vector representing the levels to iterate through. The order will be the one the iterations are happening in.

go_obj

The GeneOntologyElim S7 class. See GeneOntologyElim().

gene_universe_length

The length of the gene universe.

min_genes

number of minimum genes for the gene ontology term to be tested.

elim_threshold

p-value below which the elimination procedure shall be applied to the ancestors.

min_overlap

Optional minimum overlap threshold.

fdr_threshold

Optional fdr threshold.

Value

A list containing (results of all target sets concatenated):

  • go_ids - The gene ontology identifier.

  • pvals - The calculated p-values.

  • fdr - The calculated FDRs.

  • odds_ratios - The calculated odds ratios.

  • hits - The size of the overlap.

  • gene_set_lengths - The length of the gene sets.

  • no_test - The number of tests for that target set that passed the thresholds.