selectBatchHVG: *[Experimental]* Batch-aware highly variable gene selection

In MacoskoLab/liger: Linked Inference of Genomic Experimental Relationships

Batch-aware highly variable gene selection

Description

Method to select HVGs based on mean dispersions of genes that are highly variable genes in all batches. Using a the top target_genes per batch by average normalize dispersion. If target genes still hasn't been reached, then HVGs in all but one batches are used to fill up. This is continued until HVGs in a single batch are considered.

This is an rliger implementation of the method originally published in SCIB. We found the potential that it can improve integration under some circumstances, and is currently testing it.

This function currently only works for shared features across all datasets. For selection from only part of the datasets and selection for dataset-specific unshared features, please use selectGenes().

Usage

selectBatchHVG(object, ...)

## S3 method for class 'liger'
selectBatchHVG(
  object,
  nGenes = 2000,
  verbose = getOption("ligerVerbose", TRUE),
  ...
)

## S3 method for class 'ligerDataset'
selectBatchHVG(
  object,
  nGenes = 2000,
  features = NULL,
  verbose = getOption("ligerVerbose", TRUE),
  ...
)

## S3 method for class 'dgCMatrix'
selectBatchHVG(
  object,
  nGenes = 2000,
  returnStats = FALSE,
  verbose = getOption("ligerVerbose", TRUE),
  ...
)

Arguments

object	A liger object, ligerDataset object or a sparse/dense matrix. The liger objects must have raw counts available. A direct matrix input is preferably log-1p transformed from CPM normalized counts in cell per column orientation.
...	Arguments passed to S3 methods.
nGenes	Integer number of target genes to select. Default 2000.
verbose	Logical. Whether to show a progress bar. Default getOption("ligerVerbose") or TRUE if users have not set.
features	For ligerDataset method, the feature subset to limit the selection to, due to liger downstream requires non-zero features to be considered. Default NULL selects from all features.
returnStats	Logical, for dgCMatrix-method, whether to return a data frame of statistics for all features, or by default FALSE just return a character vector of selected features.

Value

• liger-method: Returns the input liger object with the selected genes updated in varFeatures slot, which can be accessed with varFeatures(object). Additionally, the statistics are updated in the featureMeta slot of each ligerDataset object within the datasets slot of the object.

• ligerDataset-method: Returns the input ligerDataset object with the statistics updated in the featureMeta slot.

• dgCMatrix-method: By default returns a character vector of selected variable features. If returnStats = TRUE, returns a data.frame of the statistics.

References

Luecken, M.D., Büttner, M., Chaichoompu, K. et al. (2022), Benchmarking atlas-level data integration in single-cell genomics. Nat Methods, 19, 41–50. https://doi.org/10.1038/s41592-021-01336-8.

See Also

selectGenes()

Examples

pbmc <- selectBatchHVG(pbmc, nGenes = 10)
varFeatures(pbmc)

MacoskoLab/liger documentation built on Jan. 19, 2025, 8:22 a.m.