Added model selection: AIC, AICc, BIC, and BICMAP.
MAP estimation is done by regularizing the likelihood with a dirichlet prior.
Single step in EM procedure is altered to account for this.