MDMix: EM fitting of the Multinomial Dirichlet MIMOSA model.
Description
This function fits the multinomial dirichelt MIMOSA model using EM. It can also be used to initialize the model parameters for the MCMC model.
Usage
MDMix(data = NULL, modelmatrix = NULL, alternative = "greater", initonly = FALSE)
Arguments
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.
Value
An object of class MDMixResult