scoreHistogram: Score histogram

Description Usage Arguments Details Value See Also Examples

View source: R/score_wrapper.R

Description

This function computes the empirical score distribution for a given set of DNA sequences.

Usage

1
scoreHistogram(seqs, pfm, bg)

Arguments

seqs

A DNAStringSet or DNAString object

pfm

An R matrix that represents a position frequency matrix

bg

A Background object

Details

It can be used to compare the empirical score distribution against the theoretical one (see scoreDist).

Value

List containing

scores

Vector of scores

dist

Score distribution

See Also

scoreDist

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 empirical score histogram
scoreHistogram(seqs, motif, bg)

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