Fit GLM with LEP-based penalized maximum likelihood