'Second derivative of the link function g\'\'(p)
implemented through numerical differentiation'

def deriv2(self, p):

    from statsmodels.tools.numdiff import approx_fprime_cs
    return np.diag(approx_fprime_cs(p, self.deriv))
