ENH add GLM.fit_constrained, almost all copy-paste not DRY
ENH GLM add llnull, fix starting mu to use predict (incl. offset, exposure)
TST/REF: return constraint result from fit_constrained function
         tests for GLM.fit_constrained and llnull