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]]


