'Compute the score function at params.
Parameters
params : array_like
Array of parameters at which to evaluate the score
function.
transformed : boolean, optional
Whether or not `params` is already transformed. Default is True.'

def score(self, params, transformed=True):

    params = np.array(params, ndmin=1)
    return approx_fprime_cs(params, self.loglike, args=(transformed,))
