Computing engine for binomial mixture models

Description

A function used by binomixEstimate for fitting binomial mixture models.

Usage

1
binomixMachine(y,K,core.detect.prob=1.0)

Arguments

y

Vector with data, the number of gene clusters observed in 1,2,...,N genomes.

K

The number of binomial components to use.

core.detect.prob

The fixed detection probability of the core genes.

Details

This is the function that performs the truncated maximum likelihood estimation for binomial mixture models, and not normally something you would call directly. Use binomixEstimate to fit binomial mixture models. This function is only of interest for those who want to explore the computational details of the method.

For details on the computations, see Snipen et al (2009).

Value

The output is similar to that of binomixEstimate.

Author(s)

Lars Snipen and Kristian Hovde Liland.

References

Snipen, L., Almoy, T., Ussery, D.W. (2009). Microbial comparative pan-genomics using binomial mixture models. BMC Genomics, 10:385.

See Also

binomixEstimate, negTruncLogLike.

Examples

1
# You should never run this function explicitly, use binomixEstimate().

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.