Implemented constrained poisson model.
Not sure if correct.