Efficient sampling for Gaussian linear regression with arbitrary priors

This article develops a slice sampler for Bayesian linear regression models with arbitrary priors. The new sampler has two advantages over current approaches. One, it is faster than many custom implementations that rely on auxiliary latent variables, if the number of regressors is large. Two, it can...

ver descrição completa

Detalhes bibliográficos
Autores: Hahn, P. Richard, He, Jingyu, HEDIBERT FREITAS LOPES
Formato: artículo
Estado:Versión publicada
Fecha de publicación:2018
País:Brasil
Recursos:Instituição de Ensino Superior e de Pesquisa (INSPER)
Repositorio:Repositório Institucional da INSPER
Idioma:inglés
OAI Identifier:oai:repositorio.insper.edu.br:11224/4099
Acesso em linha:https://repositorio.insper.edu.br/handle/11224/4099
Access Level:acceso abierto
Palavra-chave:Bayesian computation
Linear regression
Shrinkage priors
Slice sampling
Descrição
Resumo:This article develops a slice sampler for Bayesian linear regression models with arbitrary priors. The new sampler has two advantages over current approaches. One, it is faster than many custom implementations that rely on auxiliary latent variables, if the number of regressors is large. Two, it can be used with any prior with a density function that can be evaluated up to a normalizing constant, making it ideal for investigating the properties of new shrinkage priors without having to develop custom sampling algorithms. The new sampler takes advantage of the special structure of the linear regression likelihood, allowing it to produce better effective sample size per second than common alternative approaches.