fit_mdm: Fit a mixture of dirichlet-multinomial distributions using...

In dwinter/mdm: Mixtures of Dirichelet-Multinomial Distributions

Description

Fit a mixture of dirichlet-multinomial distributions using maximum likelihood

Usage

fit_mdm(x, f = 1L, phi = NULL, p = NULL, phTol = 100, cycles = 1000,
  cyclesInner = NULL, phAcc = 40, traceLevel = 0L, fixStart = FALSE)

Arguments

x	A matrix of integers from which to estimate the DM-model

Examples

genotypes <- matrix( 0.01/3, nrow=4, ncol=4)
diag(genotypes) <- 0.99
nfreq <- c(0.4, 0.1, 0.1, 0.4)
mixed <- rmdm(1000, 40, phi=0.01, p=genotypes, f=nfreq)
fit_mdm(mixed)

dwinter/mdm documentation built on May 15, 2019, 4:54 p.m.