Deriving the second partial derivatives of the log likelihood function of the log-binomial model in bsw() (Hessian matrix)
bsw()
Useful links