Non-Gaussian modeling of spatial data using scale mixing of a unified skew Gaussian process

JOURNAL OF MULTIVARIATE ANALYSIS

卷:114

Article

Zareifard, Hamid¹  Khaledi, Majid Jafari¹ 

[1] Tarbiat Modares Univ, Dept Stat, Tehran, Iran

关键词: Spatial modeling;    Scale mixing;    Unified skew Gaussian;    Random process;    Outlier;    EM algorithm;    Slice sampling;   

DOI  :  10.1016/j.jmva.2012.07.003

来源: Elsevier

【 摘 要 】

In this paper, we introduce a unified skew Gaussian-log Gaussian model and propose a general class of spatial sampling models that can account for both heavy tails and skewness. This class includes some models proposed previously in the literature. The likelihood function involves analytically intractable integrals and direct maximization of the marginal likelihood is numerically difficult. We obtain maximum likelihood estimates of the model parameters, using a stochastic approximation of the EM algorithm (SAEM). The predictive distribution at unsampled sites is approximated based on Markov chain Monte Carlo samples. The identifiability of the parameters and the performance of the proposed model is investigated by a simulation study. The usefulness of our methodology is demonstrated by analyzing a Pb data set in a region of north Iran. (c) 2012 Elsevier Inc. All rights reserved.

【 授权许可】

Free