Signac (version 1.1.0)

GeneActivity: Create gene activity matrix

Description

Compute counts per cell in gene body and promoter region.

Usage

GeneActivity(
  object,
  assay = NULL,
  features = NULL,
  extend.upstream = 2000,
  extend.downstream = 0,
  verbose = TRUE,
  ...
)

Arguments

object

A Seurat object

assay

Name of assay to use. If NULL, use the default assay

features

Genes to include. If NULL, use all protein-coding genes in the annotations stored in the object

extend.upstream

Number of bases to extend upstream of the TSS

extend.downstream

Number of bases to extend downstream of the TTS

verbose

Display messages

...

Additional options passed to FeatureMatrix

Examples

# NOT RUN {
fpath <- system.file("extdata", "fragments.tsv.gz", package="Signac")
fragments <- CreateFragmentObject(
  path = fpath,
  cells = colnames(atac_small),
  validate.fragments = FALSE
)
Fragments(atac_small) <- fragments
GeneActivity(atac_small)
# }