clustering in theta + lrr working with fit.em() function. beta.em() is broken.