'Return the gradient of the loglikelihood at params.
Parameters
params : array-like
The parameter values at which to evaluate the score function.
Notes
Returns numerical gradient.'

def score(self, params):

    loglike = self.loglike
    return approx_fprime(params, loglike, epsilon=1e-08)
