A seemingly working version

As EM progress, the marginal of observed increases, which is a good
sign that EM is work.

Needs to fix the 'glmnet' error returned by Fortran routine, which
indicate the problem is numerically unstable.