Performs HMRF for multivariate normal distribution. Accepts R data structures as inputs. Accepts a single beta.
Useful links