CreateMotifMatrix

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 = StringToGRanges(rownames(atac_small), sep = c(":", "-")),
  pwm = pwm,
  genome = BSgenome.Hsapiens.UCSC.hg19,
  sep = c(":", "-")
)
}