|
|
Bachoc, F., Porcu, E., Bevilacqua, M., Furrer, R., & Faouzi, T. (2022). Asymptotically equivalent prediction in multivariate geostatistics. Bernoulli, 28(4), 2518–2545.
Abstract: Cokriging is the common method of spatial interpolation (best linear unbiased prediction) in multivariate geo-statistics. While best linear prediction has been well understood in univariate spatial statistics, the literature for the multivariate case has been elusive so far. The new challenges provided by modern spatial datasets, being typ-ically multivariate, call for a deeper study of cokriging. In particular, we deal with the problem of misspecified cokriging prediction within the framework of fixed domain asymptotics. Specifically, we provide conditions for equivalence of measures associated with multivariate Gaussian random fields, with index set in a compact set of a d-dimensional Euclidean space. Such conditions have been elusive for over about 50 years of spatial statistics. We then focus on the multivariate Matern and Generalized Wendland classes of matrix valued covariance functions, that have been very popular for having parameters that are crucial to spatial interpolation, and that control the mean square differentiability of the associated Gaussian process. We provide sufficient conditions, for equivalence of Gaussian measures, relying on the covariance parameters of these two classes. This enables to identify the parameters that are crucial to asymptotically equivalent interpolation in multivariate geostatistics. Our findings are then illustrated through simulation studies.
|
|
|
|
Bevilacqua, M., Camano-Carrillo, C., & Porcu, E. (2022). Unifying compactly supported and Matern covariance functions in spatial statistics. J. Multivar. Anal., 189, 104949.
Abstract: The Matern family of covariance functions has played a central role in spatial statistics for decades, being a flexible parametric class with one parameter determining the smoothness of the paths of the underlying spatial field. This paper proposes a family of spatial covariance functions, which stems from a reparameterization of the generalized Wendland family. As for the Matern case, the proposed family allows for a continuous parameterization of the smoothness of the underlying Gaussian random field, being additionally compactly supported.
More importantly, we show that the proposed covariance family generalizes the Matern model which is attained as a special limit case. This implies that the (reparametrized) Generalized Wendland model is more flexible than the Matern model with an extra-parameter that allows for switching from compactly to globally supported covariance functions.
Our numerical experiments elucidate the speed of convergence of the proposed model to the Matern model. We also inspect the asymptotic distribution of the maximum likelihood method when estimating the parameters of the proposed covariance models under both increasing and fixed domain asymptotics. The effectiveness of our proposal is illustrated by analyzing a georeferenced dataset of mean temperatures over a region of French, and performing a re-analysis of a large spatial point referenced dataset of yearly total precipitation anomalies.
|
|
|
|
Faouzi, T., Porcu, E., & Bevilacqua, M. (2022). SPACE-TIME ESTIMATION AND PREDICTION UNDER FIXED-DOMAIN ASYMPTOTICS WITH COMPACTLY SUPPORTED COVARIANCE FUNCTIONS. Stat. Sin., 32(3), 1187–1203.
Abstract: We study the estimation and prediction of Gaussian processes with spacetime covariance models belonging to the dynamical generalized Wendland (DGW) family, under fixed-domain asymptotics. Such a class is nonseparable, has dynamical compact supports, and parameterizes differentiability at the origin similarly to the space-time Matern class.
Our results are presented in two parts. First, we establish the strong consistency and asymptotic normality for the maximum likelihood estimator of the microergodic parameter associated with the DGW covariance model, under fixed-domain asymptotics. The second part focuses on optimal kriging prediction under the DGW model and an asymptotically correct estimation of the mean squared error using a misspecified model. Our theoretical results are, in turn, based on the equivalence of Gaussian measures under some given families of space-time covariance functions, where both space or time are compact. The technical results are provided in the online Supplementary material.
|
|
|
|
Faouzi, T., Porcu, E., Kondrashuk, I., & Bevilacqua, M. (2023). Convergence arguments to bridge cauchy and matern covariance functions. Stat. Pap., Early Access.
Abstract: The Matern and the Generalized Cauchy families of covariance functions have a prominent role in spatial statistics as well as in a wealth of statistical applications. The Matern family is crucial to index mean-square differentiability of the associated Gaussian random field; the Cauchy family is a decoupler of the fractal dimension and Hurst effect for Gaussian random fields that are not self-similar. Our effort is devoted to prove that a scale-dependent family of covariance functions, obtained as a reparameterization of the Generalized Cauchy family, converges to a particular case of the Matern family, providing a somewhat surprising bridge between covariance models with light tails and covariance models that allow for long memory effect.
|
|
|
|
Khosravi, M., Leiva, V., Jamalizadeh, A., & Porcu, E. (2016). On a nonlinear Birnbaum-Saunders model based on a bivariate construction and its characteristics. Commun. Stat.-Theory Methods, 45(3), 772–793.
Abstract: The Birnbaum-Saunders (BS) distribution is an asymmetric probability model that is receiving considerable attention. In this article, we propose a methodology based on a new class of BS models generated from the Student-t distribution. We obtain a recurrence relationship for a BS distribution based on a nonlinear skew-t distribution. Model parameters estimators are obtained by means of the maximum likelihood method, which are evaluated by Monte Carlo simulations. We illustrate the obtained results by analyzing two real data sets. These data analyses allow the adequacy of the proposed model to be shown and discussed by applying model selection tools.
|
|
