Add sufficient statistics for Hessian of observed-data logl.
Factor out variational code into separate file.
Fix covariance estimation bug, and include updated documentation.