The classical linear model of coregionalization and its simpler counterpart known as the proportional covariance model, or intrinsic correlation model, have become standard tools in many areas of application for the analysis of multivariate spatial data. Despite the merits of this model, it guarantees optimal predictions only in the case of Gaussian data and can lead to erroneous conclusions in all other circumstances, in particular in the presence of skew data. To deal with multivariate geostatistical data showing some degree of skewness, this article proposes a latent spatial factor model in which all finite-dimensional marginal distributions are multivariate unified skew-normal. For this model, we can write the log-likelihood function of the data and implement a maximum likelihood estimation procedure which enables the simultaneous estimation of all parameters of the model. Moreover, we also show how the computational burden involved in the nonlinear mapping of the latent factors can be substantially reduced by exploiting a linearity property of the predictions. The sampling performances of the inferential procedures are investigated in some thorough simulation studies, and an application to radioactive contamination data is presented to show the flexibility of the model. Detailed derivations of our results are available as Supplementary Material.

A unified skew-normal geostatistical factor model

Minozzo, M;Bagnato, L
2021-01-01

Abstract

The classical linear model of coregionalization and its simpler counterpart known as the proportional covariance model, or intrinsic correlation model, have become standard tools in many areas of application for the analysis of multivariate spatial data. Despite the merits of this model, it guarantees optimal predictions only in the case of Gaussian data and can lead to erroneous conclusions in all other circumstances, in particular in the presence of skew data. To deal with multivariate geostatistical data showing some degree of skewness, this article proposes a latent spatial factor model in which all finite-dimensional marginal distributions are multivariate unified skew-normal. For this model, we can write the log-likelihood function of the data and implement a maximum likelihood estimation procedure which enables the simultaneous estimation of all parameters of the model. Moreover, we also show how the computational burden involved in the nonlinear mapping of the latent factors can be substantially reduced by exploiting a linearity property of the predictions. The sampling performances of the inferential procedures are investigated in some thorough simulation studies, and an application to radioactive contamination data is presented to show the flexibility of the model. Detailed derivations of our results are available as Supplementary Material.
factor model
multivariate geostatistics
spatial prediction
unified skew normal distribution
File in questo prodotto:
File Dimensione Formato  
env.2672.pdf

non disponibili

Descrizione: Articolo principale
Tipologia: Versione dell'editore
Licenza: Accesso ristretto
Dimensione 1.55 MB
Formato Adobe PDF
1.55 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
env2672-sup-0001-supinfo.pdf

non disponibili

Descrizione: Supplementary material
Tipologia: Versione dell'editore
Licenza: Accesso ristretto
Dimensione 96.82 kB
Formato Adobe PDF
96.82 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11562/1038779
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 1
social impact