Hybrid estimation of semivariogram parameters

Two widely used methods of semivariogram estimation are weighted least squares estimation and maximum likelihood estimation. The former have certain computational advantages, whereas the latter are more statistically efficient. We introduce and study a "hybrid" semivariogram estimation procedure that combines weighted least squares estimation of the range parameter with maximum likelihood estimation of the sill (and nugget) assuming known range, in such a way that the sill-to-range ratio in an exponential semivariogram is estimated consistently under an infill asymptotic regime. We show empirically that such a procedure is nearly as efficient computationally, and more efficient statistically for some parameters, than weighted least squares estimation of all of the semivariogram's parameters. Furthermore, we demonstrate that standard plug-in (or empirical) spatial predictors and prediction error variances, obtained by replacing the unknown semivariogram parameters with estimates in expressions for the ordinary kriging predictor and kriging variance, respectively, perform better when hybrid estimates are plugged in than when weighted least squares estimates are plugged in. In view of these results and the simplicity of computing the hybrid estimates from weighted least squares estimates, we suggest that software that currently estimates the semivariogram by weighted least squares methods be amended to include hybrid estimation as an option.