Maximum Lilkelihood and Restricted Maximum Likelihood Estimation for a Class of Gaussian Markov Random Fields

Abstract

This work describes a Gaussian Markov random field model that includes several previously proposed models, and studies properties of their maximum likelihood (ML) and restricted maximum likelihood (REML) estimators in a special case. Specifically, for models where a particular relation holds between the regression and precision matrices of the model, we provide sufficient conditions for existence and uniqueness of ML and REML estimators of the covariance parameters, and provide a straightforward way to compute them. It is found that the ML estimator always exists while the REML estimator may not exist with positive probability. A numerical comparison suggests that for this model ML estimators of covariance parameters have, overall, better frequentist properties than REML estimators.Eigenvalues and eigenvectors; Profile likelihood; Restricted likelihood; Spatial data.