pcm-class {motifStack} | R Documentation

## Class `"pcm"`

### Description

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 from the Class

Objects can be created by calls of the form `new("pcm", mat, name, alphabet, color, background)`.

### Slots

`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.

`tags`

  a list to save metadata of motif.

`markers`

  a list of marker object.

### Methods

addBlank

  `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.

coerce

  `signature(from = "pcm", to = "matrix")`: convert object pcm to matrix

getIC

  `signature(x = "pcm",)` Calculate information content profile for position frequency matrix.

matrixReverseComplement

  `signature(x = "pcm")` get the reverse complement of position frequency matrix.

trimMotif

  `signature(x = "pcm", t= "numeric")` trim motif by information content.

plot

  `signature(x = "pcm")` Plots the sequence logo of the position count matrix.

### Examples

    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)

[Package motifStack version 1.30.0 Index]