modOverlaps: Calculate overlaps of the modules

In tmod: Feature Set Enrichment Analysis for Metabolomics and Transcriptomics

Calculate overlaps of the modules

Description

Calculate overlaps of the modules

Usage

modOverlaps(modules = NULL, mset = NULL, stat = "jaccard")

Arguments

modules	either a character vector with module IDs from mset, or a list which contains the module members
mset	Which module set to use. Either a character vector ("LI", "DC" or "all", default: all) or an object of class tmod (see "Custom module definitions" below)
stat	Type of statistics to return. "jaccard": Jaccard index (default); "number": number of common genes; "soerensen": Soerensen-Dice coefficient; "overlap": Szymkiewicz-Simpson coefficient.

Details

For a set of modules (aka gene sets) determine the similarity between these. You can run modOverlaps either on a character vector of module IDs or on a list. In the first case, the module elements are taken from 'mset', and if that is NULL, from the default tmod module set.

Alternatively, you can provide a list in which each element is a character vector. In this case, the names of the list are the module IDs, and the character vectors contain the associated elements.

The different statistics available are: * "number": total number of common genes (size of the overlap) * "jaccard": Jaccard index, i.e. \frac{|A \cap B|}{|A \cup B|} (number of common elements divided by the total number of unique elements); * "soerensen": Soerensen-Dice coefficient, defined as \frac{2 \cdot |A \cap B|}{|A| + |B|} – number of common genes in relation to the total number of elements (divided by two, such that the maximum is 1) (number of common elements divided by the average size of both gene sets) * "overlap": Szymkiewicz-Simpson coefficient, defined as \frac{|A \cap B|}{\min(|A|, |B|)} – this is the number of common genes scaled by the size of the smaller of the two gene sets (number of common elements divided by the size of the smaller gene set)

tmod documentation built on March 31, 2023, 9 p.m.