dmmEM-methods: dmmEM method
In ge11232002/TFBSTools: Software Package for Transcription Factor Binding Site (TFBS) Analysis

Description Usage Arguments Details Value Methods Author(s) See Also Examples

This function trains the Dirichlet multinomial mixture models parameters for a set of profile matrices.

1	dmmEM(x, K=6, alg=c("C", "R"))

x

x can be a matrix, PFMatrixList or JASPAR2014 to be trained.

K

The maximal number of components to test in the mixture model when alg is "C". Then an optimal number of components between 1 and K will be chosen based on the fitness of the model.

The fixed number of components to use when alg is "R". The default is 6.

alg

The algorithm to use. "C" uses the implementation from DirichletMultinomial package which has more advanced feature and performance. "R" uses our own implemention in R.

When using the implementation from DirichletMultinomial package, the final number of components can be 1:K. An internal selection will be made based on the maximum likelihood.

When using the implementation of R, the number of component is fixed to K.

A list of trainned alpha0, pmix and likelihood during the training.

signature(x = "ANY")
signature(x = "matrix")
signature(x = "PFMatrixList")

Ge Tan

rPWMDmm

  
    data(MA0003.2)
    data(MA0004.1)
    pfmList <- PFMatrixList(pfm1=MA0003.2, pfm2=MA0004.1, use.names=TRUE)
    dmmParameters <- dmmEM(pfmList, K=6, alg="C")