Estimating Space and Space-Time Covariance Functions for Large Data Sets: A Weighted Composite Likelihood Approach

In this article, we propose two methods for estimating space and space-time covariance functions from a Gaussian random field, based on the composite likelihood idea. The first method relies on the maximization of a weighted version of the composite likelihood function, while the second one is based on the solution of a weighted composite score equation. This last scheme is quite general and could be applied to any kind of composite likelihood. An information criterion for model selection based on the first estimation method is also introduced. The methods are useful for practitioners looking for a good balance between computational complexity and statistical efficiency. The effectiveness of the methods is illustrated through examples, simulation experiments, and by analyzing a dataset on ozone measurements.

[1]  M. Sherman,et al.  On the asymptotic joint distribution of sample space--time covariance estimators , 2008, 0803.2171.

[2]  A. V. Vecchia.  Estimation and model identification for continuous spatial processes , 1988 .

[3]  Xavier Guyon,et al.  Random fields on a network , 1995 .

[4]  Douglas W. Nychka,et al.  Statistical models for monitoring and regulating ground‐level ozone , 2005 .

[5]  C. Heyde,et al.  Multiple roots in general estimating equations , 1998 .

[6]  Patrick J. Heagerty,et al.  Window Subsampling of Estimating Functions with Application to Regression Models , 2000 .

[7]  Roger Woodard,et al.  Interpolation of Spatial Data: Some Theory for Kriging , 1999, Technometrics.

[8]  W Pan,et al.  Model Selection in Estimating Equations , 2001, Biometrics.

[9]  D. Zimmerman,et al.  Towards reconciling two asymptotic frameworks in spatial statistics , 2005 .

[10]  H. Rue,et al.  Fitting Gaussian Markov Random Fields to Gaussian Fields , 2002 .

[11]  Jennifer A Hoeting,et al.  Model selection for geostatistical models. , 2006, Ecological applications : a publication of the Ecological Society of America.

[12]  Montserrat Fuentes,et al.  Statistical assessment of geographic areas of compliance with air quality standards : Application of recent advances in space-time statistics to atmospheric data , 2003 .

[13]  Harry Joe,et al.  Composite Likelihood Methods , 2012 .

[14]  M. Palma,et al.  Prediction of particle pollution through spatio-temporal multivariate geostatistical analysis: spatial special issue , 2013 .

[15]  Douglas W. Nychka,et al.  Covariance Tapering for Likelihood-Based Estimation in Large Spatial Data Sets , 2008 .

[16]  P. Whittle,et al.  Topographic correlation, power-law covariance functions, and diffusion , 1962 .

[17]  C. Varin,et al.  A note on composite likelihood inference and model selection , 2005 .

[18]  N. Cressie,et al.  Fixed rank kriging for very large spatial data sets , 2008 .

[19]  Zhiyi Chi,et al.  Approximating likelihoods for large spatial data sets , 2004 .

[20]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[21]  M. Fuentes.  Approximate Likelihood for Large Irregularly Spaced Spatial Data , 2007, Journal of the American Statistical Association.

[22]  P. Doukhan.  Mixing: Properties and Examples , 1994 .

[23]  Michael L. Stein,et al.  A modeling approach for large spatial datasets , 2008 .

[24]  Hsin-Cheng Huang,et al.  Optimal Geostatistical Model Selection , 2007 .

[25]  A. Kuk.  A Hybrid Pairwise Likelihood Method , 2007 .

[26]  T. Gneiting.  Nonseparable, Stationary Covariance Functions for Space–Time Data , 2002 .

[27]  M. Sherman.  Variance Estimation for Statistics Computed from Spatial Lattice Data , 1996 .

[28]  Noel A Cressie,et al.  Statistics for Spatial Data, Revised Edition. , 1994 .

[29]  N. Cressie,et al.  Classes of nonseparable, spatio-temporal stationary covariance functions , 1999 .

[30]  S. Lahiri,et al.  Least squares variogram fitting by spatial subsampling , 2002 .

[31]  Subhash R. Lele,et al.  A composite likelihood approach to semivariogram estimation , 1999 .

[32]  Hao Zhang.  Inconsistent Estimation and Asymptotically Equal Interpolations in Model-Based Geostatistics , 2004 .

[33]  Jorge Mateu,et al.  Model comparison and selection for stationary space-time models , 2007, Comput. Stat. Data Anal..

[34]  P. Guttorp,et al.  Geostatistical Space-Time Models, Stationarity, Separability, and Full Symmetry , 2007 .

[35]  Xiongzhi Chen,et al.  Stat 598 Presentation on Nonseperable, Stationary Covariance Functions for Space-Time Data , 2010 .

[36]  Michael L. Stein,et al.  Statistical methods for regular monitoring data , 2005 .

[37]  E. Carlstein.  The Use of Subseries Values for Estimating the Variance of a General Statistic from a Stationary Sequence , 1986 .

[38]  P. Whittle.  ON STATIONARY PROCESSES IN THE PLANE , 1954 .

[39]  Noel A Cressie,et al.  On asymptotic distribution and asymptotic efficiency of least squares estimators of spatial variogram parameters , 2002 .

[40]  Richard L. Smith,et al.  Approximate Likelihoods for Spatial Processes , 2002 .

[41]  E. Andersen.  Composite likelihood and two-stage estimation in family studies. , 2004, Biostatistics.