'Return the gradient of the loglike at AB_mask.
Parameters
AB_mask : unknown values of A and B matrix concatenated
Notes
Return numerical gradient'

def score(self, AB_mask):

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