type = 'score' for multiple group as well
Avoid computing Hessian if not needed. Also fixed estimation bug when MHRM was used with #PLNRM models