'Returns the Hessian matrix of the log partial likelihood
function evaluated at `params`.'

def hessian(self, params):

    if (self.ties == 'breslow'):
        return self.breslow_hessian(params)
    else:
        return self.efron_hessian(params)
