model.NB function

Function to perform an EM algorithm (Gaussian naive Bayes (parametric, independent features)) for a 2-component mixture where each observation ii has a p-value pip_i, and an auxiliary feature vector viv_i.