selex.kmerPSFM: Construct a PSFM from a K-mer table

Description Usage Arguments Details Value See Also Examples

View source: R/SELEX.R

Description

A function used to calculate and return the PFSM (Position Specific Frequency Matrix) for a K-mer table of length k from sample. If an offset value is provided, K-mer counting takes place at a fixed position within the variable region.

Usage

1

Arguments

sample

A sample handle to the dataset on which K-mer counting should be perfomed.

k

K-mer length(s) to be counted.

offset

Location of window for which K-mers should be counted for. If not provided, K-mers are counted across all windows. offset starts from 0.

Details

A K-mer table will be constructed for the specified sample, length k, and offset if it does not already exist.

The output can be used by the seqLogo package to create a sequence logo.

Value

selex.kmerPSFM returns a matrix containing the frequences for each base at every position.

See Also

selex.counts, selex.fastqPSFM

Examples

1
2
3
4
5
6
# Build the PSFM
psfm1 = selex.kmerPSFM(sample=r0, k=8, 0)

# Can make sequence logos using the seqLogo package:
#library(seqLogo)
#seqLogo(makePWM(t(psfm1)))

Bioconductor-mirror/SELEX documentation built on June 1, 2017, 1:14 a.m.