Computationally efficient restricted maximum likelihood estimation of generalized covariance functions

Computational aspects of the estimation of generalized covariance functions by the method of restricted maximum likelihood (REML) are considered in detail. In general, REML estimation is computationally intensive, but significant computational savings are available in important special cases. The approach taken here restricts attention to data whose spatial configuration is a regular lattice, but makes no restrictions on the number of parameters involved in the generalized covariance nor (with the exception of one result) on the nature of the generalized covariance function's dependence on those parameters. Thus, this approach complements the recent work of L. G. Barendregt (1987), who considered computational aspects of REML estimation in the context of arbitrary spatial data configurations, but restricted attention to generalized covariances which are linear functions of only two parameters.