Dear all, Is there exist a function which provides an OPG (outer product of the gradient) estimate for variance - covariance matrix of the MLE? I tried to write my own function but it produces significantly different results comparing with empirical sd estimates computed from simulations and I cannot understand where is mistake. Thanks a lot, Vicky Landsman. [[alternative HTML version deleted]]