Abstract

Generalized linear models are routinely used in many environment statistics problems such as earthquake magnitudes prediction. Hu et al. proposed Pareto regression with spatial random effects for earthquake magnitudes. In this paper, we propose Bayesian spatial variable selection for Pareto regression based on Bradley et al. and Hu et al. to tackle variable selection issue in generalized linear regression models with spatial random effects. A Bayesian hierarchical latent multivariate log gamma model framework is applied to account for spatial random effects to capture spatial dependence. We use two Bayesian model assessment criteria for variable selection including Conditional Predictive Ordinate (CPO) and Deviance Information Criterion (DIC). Furthermore, we show that these two Bayesian criteria have analytic connections with conditional AIC under the linear mixed model setting. We examine empirical performance of the proposed method via a simulation study and further demonstrate the applicability of the proposed method in an analysis of the earthquake data obtained from the United States Geological Survey (USGS).

Highlights

  • The earthquake magnitude data has become increasingly popular over the last decade

  • Conditional Predictive Ordinate (CPO) and Deviance Information Criterion (DIC). We show that these two Bayesian criteria have analytic connections with conditional Akaike’s information criterion (AIC) under the linear mixed model setting

  • In order to have more explicit understanding of dependent covariates of earthquake magnitudes, variable selection approaches should be considered in a Pareto regression model with spatial random effects

Read more

Summary

Introduction

The earthquake magnitude data has become increasingly popular over the last decade. Statistical models for earthquake have been proposed since 1800s. Hu and Bradley [5] proposed using the Pareto regression with spatial random effects for earthquake magnitudes, but they did not consider the model selection problems. In order to have more explicit understanding of dependent covariates of earthquake magnitudes, variable selection approaches should be considered in a Pareto regression model with spatial random effects. We use CPO and DIC criteria to carry out Bayesian variable selection for Pareto regression models due to the performance of the conjugate priors (see [16], for a discussion). Both CPO and DIC are criteria-based methods and they have some advantage over other criteria.

Pareto Regression with Spatial Random Effects
Bayesian Model Assessment Criteria
MCMC Scheme
Simulation Study
Simulation for Estimation Performance
Simulation for Model Selection
Simulation for Model Comparison
Data Description
Analysis
Discussion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call