scoreHistogramSingleSeq: Score histogram on a single sequence

Description Usage Arguments Value Examples

View source: R/score_wrapper.R

Description

This function computes the empirical score distribution by normalizing the observed score histogram for a given sequence.

Usage

1

Arguments

seq

A DNAString object

pfm

An R matrix that represents a position frequency matrix

bg

A Background object

Value

List containing

scores

Vector of scores

dist

Score distribution

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
# Load sequences
seqfile = system.file("extdata", "seq.fasta", package = "motifcounter")
seqs = Biostrings::readDNAStringSet(seqfile)

# Load background
bg = readBackground(seqs, 1)

# Load motif
motiffile = system.file("extdata", "x31.tab", package = "motifcounter")
motif = t(as.matrix(read.table(motiffile)))

# Compute the per-position and per-strand scores
motifcounter:::scoreHistogramSingleSeq(seqs[[1]], motif, bg)

motifcounter documentation built on Nov. 8, 2020, 5:44 p.m.