pcm-class {motifStack} | R Documentation |
"pcm"
An object of class "pcm"
represents the position count
matrix of a DNA/RNA/amino-acid sequence motif. The entry stores a matrix,
which in row i, column j gives the
counts of observing nucleotide/or amino acid i in
position j of the motif.
Objects can be created by calls of the form new("pcm", mat, name, alphabet, color, background)
.
mat
Object of class "matrix"
The position count matrix
name
code"character" The motif name
alphabet
"character"
The sequence
alphabet. "DNA", "RNA", "AA" or "others".
color
a "character"
vector. The color setting for each symbol
background
a "numeric"
vector. The background frequency.
signature(x="pcm", n="numeric", b="logical")
add space
into the position count matrix for alignment. b is a bool value, if TRUE, add
space to the 3' end, else add space to the 5' end. n indicates how many
spaces should be added.
signature(from = "pcm", to = "matrix")
: convert object pcm to matrix
signature(x = "pcm",)
Calculate information
content profile for position frequency matrix.
signature(x = "pcm")
get the
reverse complement of position frequency matrix.
signature(x = "pcm", t= "numeric")
trim motif by information content.
signature(x = "pcm")
Plots the
sequence logo of the position count matrix.
pcm <- read.table(file.path(find.package("motifStack"), "extdata", "bin_SOLEXA.pcm")) pcm <- pcm[,3:ncol(pcm)] rownames(pcm) <- c("A","C","G","T") motif <- new("pcm", mat=as.matrix(pcm), name="bin_SOLEXA") plot(motif)