'Jacobian/Gradient of log-likelihood evaluated at params for each
observation.'

def score_obs(self, params, **kwds):

    kwds.setdefault('centered', True)
    return approx_fprime(params, self.loglikeobs, **kwds)
