Likelihood-based estimation for Gaussian MRFs