### Abstract

The Bayes factor has become an important tool for model selection. The marginal likelihoods are also important because they can be used to rank the models. In fact, the Bayes factor is the ratio of the marginal likelihoods for two models with proper prior densities. We discuss the marginal likelihood for a class of generalized linear models used in small area estimation for mortality data analysis. Computation in these models is intensive and requires the implementation of Markov chain Monte Carlo (MCMC) methods. A sophisticated method for computing the marginal likelihoods for generalized linear models using reduced Metropolis-Hastings (M-H) samplers has recently been introduced. Also, a much simpler method that uses the Laplace approximation has been proposed. Our method lies between these two in simplicity, and it uses importance sampling via a simple output analysis from a MCMC sampler. We also show that the new method can be approximated without using the MCMC sampler, We illustrate our methods for two Poisson regression models which have been used for mapping mortality data, and we compare them with the methods based on reduced M-H samplers and the Laplace approximation. The method based on the Laplace approximation gives smaller marginal log-likelihoods than the other three methods, but they select the same model, We use a simulation study in which the marginal likelihood is used to investigate the improvement in the goodness of fit of one hierarchical Poisson regression model over the other.

Original language | English |
---|---|

Pages (from-to) | 319-340 |

Number of pages | 22 |

Journal | Journal of Statistical Computation and Simulation |

Volume | 72 |

Issue number | 4 |

DOIs | |

Publication status | Published - 2002 Dec 1 |

### Fingerprint

### All Science Journal Classification (ASJC) codes

- Statistics and Probability
- Modelling and Simulation
- Statistics, Probability and Uncertainty
- Applied Mathematics

### Cite this

}

*Journal of Statistical Computation and Simulation*, vol. 72, no. 4, pp. 319-340. https://doi.org/10.1080/00949650212842

**Marginal likelihood for a class of Bayesian generalized linear models.** / Nandram, Balgobin; Kim, Hyunjoong.

Research output: Contribution to journal › Article

TY - JOUR

T1 - Marginal likelihood for a class of Bayesian generalized linear models

AU - Nandram, Balgobin

AU - Kim, Hyunjoong

PY - 2002/12/1

Y1 - 2002/12/1

N2 - The Bayes factor has become an important tool for model selection. The marginal likelihoods are also important because they can be used to rank the models. In fact, the Bayes factor is the ratio of the marginal likelihoods for two models with proper prior densities. We discuss the marginal likelihood for a class of generalized linear models used in small area estimation for mortality data analysis. Computation in these models is intensive and requires the implementation of Markov chain Monte Carlo (MCMC) methods. A sophisticated method for computing the marginal likelihoods for generalized linear models using reduced Metropolis-Hastings (M-H) samplers has recently been introduced. Also, a much simpler method that uses the Laplace approximation has been proposed. Our method lies between these two in simplicity, and it uses importance sampling via a simple output analysis from a MCMC sampler. We also show that the new method can be approximated without using the MCMC sampler, We illustrate our methods for two Poisson regression models which have been used for mapping mortality data, and we compare them with the methods based on reduced M-H samplers and the Laplace approximation. The method based on the Laplace approximation gives smaller marginal log-likelihoods than the other three methods, but they select the same model, We use a simulation study in which the marginal likelihood is used to investigate the improvement in the goodness of fit of one hierarchical Poisson regression model over the other.

AB - The Bayes factor has become an important tool for model selection. The marginal likelihoods are also important because they can be used to rank the models. In fact, the Bayes factor is the ratio of the marginal likelihoods for two models with proper prior densities. We discuss the marginal likelihood for a class of generalized linear models used in small area estimation for mortality data analysis. Computation in these models is intensive and requires the implementation of Markov chain Monte Carlo (MCMC) methods. A sophisticated method for computing the marginal likelihoods for generalized linear models using reduced Metropolis-Hastings (M-H) samplers has recently been introduced. Also, a much simpler method that uses the Laplace approximation has been proposed. Our method lies between these two in simplicity, and it uses importance sampling via a simple output analysis from a MCMC sampler. We also show that the new method can be approximated without using the MCMC sampler, We illustrate our methods for two Poisson regression models which have been used for mapping mortality data, and we compare them with the methods based on reduced M-H samplers and the Laplace approximation. The method based on the Laplace approximation gives smaller marginal log-likelihoods than the other three methods, but they select the same model, We use a simulation study in which the marginal likelihood is used to investigate the improvement in the goodness of fit of one hierarchical Poisson regression model over the other.

UR - http://www.scopus.com/inward/record.url?scp=0038826034&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0038826034&partnerID=8YFLogxK

U2 - 10.1080/00949650212842

DO - 10.1080/00949650212842

M3 - Article

VL - 72

SP - 319

EP - 340

JO - Journal of Statistical Computation and Simulation

JF - Journal of Statistical Computation and Simulation

SN - 0094-9655

IS - 4

ER -