Last data update: 2014.03.03
By Ravi Varadhan. Constrained optimisation for likelihood function, returning hessian matrix used for parameter standard errors.