Electronic Journal of Statistics

Maximum likelihood estimation for a bivariate Gaussian process under fixed domain asymptotics

Daira Velandia, François Bachoc, Moreno Bevilacqua, Xavier Gendre, and Jean-Michel Loubes

Full-text: Open access


We consider maximum likelihood estimation with data from a bivariate Gaussian process with a separable exponential covariance model under fixed domain asymptotics. We first characterize the equivalence of Gaussian measures under this model. Then consistency and asymptotic normality for the maximum likelihood estimator of the microergodic parameters are established. A simulation study is presented in order to compare the finite sample behavior of the maximum likelihood estimator with the given asymptotic distribution.

Article information

Electron. J. Statist. Volume 11, Number 2 (2017), 2978-3007.

Received: July 2016
First available in Project Euclid: 11 August 2017

Permanent link to this document

Digital Object Identifier

Zentralblatt MATH identifier

Bivariate exponential model equivalent Gaussian measures infill asymptotics microergodic parameters

Creative Commons Attribution 4.0 International License.


Velandia, Daira; Bachoc, François; Bevilacqua, Moreno; Gendre, Xavier; Loubes, Jean-Michel. Maximum likelihood estimation for a bivariate Gaussian process under fixed domain asymptotics. Electron. J. Statist. 11 (2017), no. 2, 2978--3007. doi:10.1214/17-EJS1298. https://projecteuclid.org/euclid.ejs/1502416821

Export citation


  • [1] Abt, M. and Welch, W. (1998). Fisher information and maximum-likelihood estimation of covariance parameters in Gaussian stochastic, processes.The Canadian Journal of Statistics, 26:127–137.
  • [2] Anderes, E. (2010). On the consistent separation of scale and variance for Gaussian random, fields.The Annals of Statistics, 38:870–893.
  • [3] Antognini, A. B. and Zagoraiou, M. (2010). Exact optimal designs for computer experiments via Kriging, metamodelling.Journal of Statistical Planning and Inference, 140:2607–2617.
  • [4] Bachoc, F. (2014). Asymptotic analysis of the role of spatial sampling for covariance parameter estimation of Gaussian, processes.Journal of Multivariate Analysis, 125:1–35.
  • [5] Bevilacqua, M., Vallejos, R., and Velandia, D. (2015). Assessing the significance of the correlation between the components of a bivariate Gaussian random, field.Environmetrics, 26:545–556.
  • [6] Chen, H., Simpson, D., and Ying, Z. (2000). Infill asymptotics for a stochastic process model with measurement, error.Statistica Sinica, 10:141–156.
  • [7] Cressie, N., (1993).Statistics for spatial data. J. Wiley.
  • [8] Cressie, N. and Lahiri, S. (1993). The asymptotic distribution of REML, estimators.Journal of Multivariate Analysis, 45:217–233.
  • [9] Cressie, N. and Lahiri, S. (1996). Asymptotics for REML estimation of spatial covariance, parameters.Journal of Statistical Planning and Inference, 50:327–341.
  • [10] Du, J., Zhang, H., and Mandrekar, V. S. (2009). Fixed-domain asymptotic properties of tapered maximum likelihood, estimators.The Annals of Statistics, 37:3330–3361.
  • [11] Furrer, R., Bachoc, F., and Du, J. (2016). Asymptotic properties of multivariate tapering for estimation and, prediction.Journal of Multivariate Analysis, 149:177–191.
  • [12] Furrer, R., Genton, M. G., and Nychka, D. (2006). Covariance tapering for interpolation of large spatial, datasets.Journal of Computational and Graphical Statistics, 15(3):502–523.
  • [13] Gneiting, T., Kleiber, W., and Schlather, M. (2010). Matérn cross-covariance functions for multivariate random, fields.Journal of the American Statistical Association, 105:1167–1177.
  • [14] Ibragimov, I. A. and Rozanov, Y. A., (1978).Gaussian Random Processes. Springer-Verlag New York.
  • [15] Kaufman, C. G., Schervish, M. J., and Nychka, D. W. (2008a). Covariance Tapering for Likelihood-Based Estimation in Large Spatial Data, Sets.Journal of the American Statistical Association, 103(484):1545–1555.
  • [16] Kaufman, C. G., Schervish, M. J., and Nychka, D. W. (2008b). Covariance tapering for likelihood-based estimation in large spatial data, sets.Journal of the American Statistical Association, 103:1545–1555.
  • [17] Kaufman, C. G. and Shaby, B. A. (2013). The role of the range parameter for estimation and prediction in, geostatistics.Biometrika, 100:473–484.
  • [18] Lehrke, S. G. and Ghorai, J. K. (2010). Large sample properties of ML estimator of the parameters of multivariate O–U random, fields.Communications in Statistics – Theory and Methods, 39(4):738–752.
  • [19] LeMone, M. A., Chen, F., Alfieri, J. G., Cuenca, R., Hagimoto, Y., Blanken, P. D., Niyogi, D., Kang, S., Davis, K., and Grossman, R. (2007). NCAR/CU surface vegetation observation network during the international H2O project 2002 field, campaign.Bull. Amer. Meteor. Soc., 88: 65–81.
  • [20] Loh, W. L. (2005). Fixed-domain asymptotics for a subclass of Matérn-type Gaussian random, fields.The Annals of Statistics, 33:2344–2394.
  • [21] Loh, W. L. and Lam, T. K. (2000). Estimating structured correlation matrices in smooth Gaussian random fields, models.The Annals of Statistics, 28:880–904.
  • [22] Mardia, K. V. and Marshall, R. J. (1984). Maximum likelihood estimation of models for residual covariance in spatial, regression.Biometrika, 71:135–146.
  • [23] Pascual, F. G. and Zhang, H. (2006). Estimation of linear correlation coefficient of two correlated spatial, processes.Sankhyā: The Indian Journal of Statistics, 68:307–325.
  • [24] Shaby, B. A. and Ruppert, D. (2012). Tapered covariance: Bayesian estimation and, asymptotics.Journal of Computational and Graphical Statistics, 21(2):433–452.
  • [25] Stassberg, D., LeMone, M. A., Warner, T., and Alfieri, J. G. (2008). Comparison of observed 10 m wind speeds to those based on Monin–Obukhov similarity theory using aircraft and surface data from the international H2O, project.Mon. Wea. Rev., 136:964–972.
  • [26] Stein, M. (1988). Asymptotically efficient prediction of a random field with a misspecified covariance, function.The Annals of Statistics, 16:55–63.
  • [27] Stein, M. (1990a). Bounds on the efficiency of linear predictions using an incorrect covariance, function.The Annals of Statistics, 18:1116–1138.
  • [28] Stein, M. (1990b). Uniform asymptotic optimality of linear predictions of a random field using an incorrect second-order, structure.The Annals of Statistics, 18:850–872.
  • [29] Stein, M. L., (1999).Interpolation of Spatial Data. Springer Series in Statistics. Springer-Verlag New York.
  • [30] van der Vaart, A. (1996). Maximum likelihood estimation under a spatial sampling, scheme.The Annals of Statistics, 5:2049–2057.
  • [31] Wang, D. and Loh, W. L. (2011). On fixed-domain asymptotics and covariance tapering in Gaussian random field, models.Electronic Journal of Statistics, 5:238–269.
  • [32] Weckwerth, T. M., Parsons, D. B., Koch, S. E., Moore, J. A., LeMone, M. A., Demoz, B. B., Flamant, C., Geerts, B., Wang, J., and Feltz, W. F. (2004). An overview of the international H2O project (IHOP 2002) and some preliminary, highlights.Bull. Amer. Meteor. Soc., 85:253–277.
  • [33] Ying, Z. (1991). Asymptotic properties of a maximum likelihood estimator with data from a Gaussian, process.Journal of Multivariate Analysis, 36:280–296.
  • [34] Ying, Z. (1993). Maximum likelihood estimation of parameters under a spatial sampling, scheme.The Annals of Statistics, 21:1567–1590.
  • [35] Zhang, H. (2004). Inconsistent estimation and asymptotically equivalent interpolations in model-based, geostatistics.Journal of the American Statistical Association, 99:250–261.
  • [36] Zhang, H. and Cai, W. (2015). When doesn’t cokriging outperform, Kriging?Statistical Science, 30:176–180.
  • [37] Zhang, H. and Zimmerman, D. L. (2005). Towards reconciling two asymptotic frameworks in spatial, statistics.Biometrika, 92:921–936.