FeatureMatrix: Feature Matrix

View source: R/quantification.R

FeatureMatrixR Documentation

Feature Matrix

Description

Construct a feature x cell matrix from a genomic fragments file

Usage

FeatureMatrix(
  fragments,
  features,
  cells = NULL,
  process_n = 2000,
  sep = c("-", "-"),
  verbose = TRUE
)

Arguments

fragments

A list of Fragment objects. Note that if setting the cells parameter, the requested cells should be present in the supplied Fragment objects. However, if the cells information in the fragment object is not set (Cells(fragments) is NULL), then the fragment object will still be searched.

features

A GRanges object containing a set of genomic intervals. These will form the rows of the matrix, with each entry recording the number of unique reads falling in the genomic region for each cell. If a genomic region provided is on a chromosome that is not present in the fragment file, it will not be included in the returned matrix.

cells

Vector of cells to include. If NULL, include all cells found in the fragments file

process_n

Number of regions to load into memory at a time, per thread. Processing more regions at once can be faster but uses more memory.

sep

Vector of separators to use for genomic string. First element is used to separate chromosome and coordinates, second separator is used to separate start and end coordinates.

verbose

Display messages

Value

Returns a sparse matrix

Examples

fpath <- system.file("extdata", "fragments.tsv.gz", package="Signac")
fragments <- CreateFragmentObject(fpath)
FeatureMatrix(
  fragments = fragments,
  features = granges(atac_small)
)

Signac documentation built on Sept. 11, 2024, 9:30 p.m.