Bayesian Methods for Generalized Linear Models

  • Paul E. Green (Associate Professor, Department of Statistical Informtion, Catholic University) ;
  • Kim, Dae-Hak (Associate Professor, Department of Statistical Informtion, Catholic University)
  • Published : 1999.08.01

Abstract

Generalized linear models have various applications for data arising from many kinds of statistical studies. Although the response variable is generally assumed to be generated from a wide class of probability distributions we focus on count data that are most often analyzed using binomial models for proportions or poisson models for rates. The methods and results presented here also apply to many other categorical data models in general due to the relationship between multinomial and poisson sampling. The novelty of the approach suggested here is that all conditional distribution s can be specified directly so that staraightforward Gibbs sampling is possible. The prior distribution consists of two stages. We rely on a normal nonconjugate prior at the first stage and a vague prior for hyperparameters at the second stage. The methods are demonstrated with an illustrative example using data collected by Rosenkranz and raftery(1994) concerning the number of hospital admissions due to back pain in Washington state.

Keywords

References

  1. Categorical Data Analysis Agresti, A.
  2. J. Amer. Statist. Assoc. v.85 Sampling-based approaches to calculating marginal densities Gelfand, A.E.;Smith A.F.M.
  3. Bayesian Data Analysis Gelman, A.;Carlin, J.B.;Stern, H.S.;Rubin, D.B.
  4. Markov Chain Monte Carlo in Practice Gilks, W.R.;Richardson, S.;Spiegelhalter, D.J.
  5. Generalized Linear Models McCullagh, P.;Nelder, J.A.
  6. Appl. Statist. v.31 Applications of a method for the efficient computation of posterior distributions Naylor, J.C.;Smith, A.F>M.
  7. Technical report no.268 Covariate selection in hierarchical models of hospital admission counts: a Bayes factor approach Rosenkranz, S.L.;Raftery, A.E.
  8. Technical report General hit-and-run Monte Carlo sampling for evaluating multidimensional integrals Schmeiser, B.;Chen, M.H.
  9. BUGS: Bayesian inference Using Gibbs Sampling, Version 0.5 Spiegelhalter, D.J.;Thomas, A.;Best, N.G.lGilks, W.R.
  10. Ann. Statist v.22 Markov chains for exploring posterior distributions (with Discussion) Tierney, L.
  11. J. Amer. Statist. Assoc. v.81 Accurate approximations for posterior moments and margins Tierney, L.;Kadane, J.
  12. Biometrika v.76 Approximate marginal densities for nonlinear functions Tierney, L.;Kass, R.E.;Kadane, J.B.