Create a motif x feature matrix from a set of genomic ranges, the genome, and a set of position weight matrices.

CreateMotifMatrix(
  features,
  pwm,
  genome,
  score = FALSE,
  use.counts = FALSE,
  sep = c("-", "-"),
  ...
)

Arguments

features

A GRanges object containing a set of genomic features

pwm

A PFMatrixList or PWMatrixList object containing position weight/frequency matrices to use

genome

Any object compatible with the genome argument in matchMotifs

score

Record the motif match score, rather than presence/absence (default FALSE)

use.counts

Record motif counts per region. If FALSE (default), record presence/absence of motif. Only applicable if score=FALSE.

sep

A length-2 character vector containing the separators to be used when constructing matrix rownames from the GRanges

...

Additional arguments passed to matchMotifs

Value

Returns a sparse matrix

Details

Requires that motifmatchr is installed https://www.bioconductor.org/packages/motifmatchr/.

Examples

if (FALSE) { library(JASPAR2018) library(TFBSTools) library(BSgenome.Hsapiens.UCSC.hg19) pwm <- getMatrixSet( x = JASPAR2018, opts = list(species = 9606, all_versions = FALSE) ) motif.matrix <- CreateMotifMatrix( features = granges(atac_small), pwm = pwm, genome = BSgenome.Hsapiens.UCSC.hg19 ) }