MDMix: EM fitting of the Multinomial Dirichlet MIMOSA model.
In RGLab/MIMOSA: Mixture Models for Single-Cell Assays

This function fits the multinomial dirichelt MIMOSA model using EM. It can also be used to initialize the model parameters for the MCMC model.

MDMix(
  data = NULL,
  modelmatrix = NULL,
  alternative = "greater",
  initonly = FALSE
)

`data`	The observed data
`modelmatrix`	a model matrix specifying which components should be computed
`alternative`	either 'greater' or 'not equal' to fit the one-sided or two-sided model.
`initonly`	`TRUE` or `FALSE` to return just the initialization parameters.