Markov chain Monte Carlo model determination for hierarchical and graphical log-linear models

Dellaportas, P. and Forster, J.J. (1999) Markov chain Monte Carlo model determination for hierarchical and graphical log-linear models. Biometrika, 86 (3), 615-633. (doi:10.1093/biomet/86.3.615).

Abstract

We use reversible jump Markov chain Monte Carlo methods (Green, 1995) to develop strategies for calculating posterior probabilities of hierarchical, graphical or decomposable log-linear models for high-dimensional contingency tables. Even for tables of moderate size, these sets of models may be very large. The choice of suitable prior distributions for model parameters is also discussed in detail, and two examples are presented. For the first example, a three-way table, the model probabilities calculated using our reversible jump approach are compared with model probabilities calculated exactly or by using an alternative approximation. The second example is a six-way contingency table for which exact methods are infeasible, because of the large number of possible models. We identify the most probable hierarchical, graphical and decomposable models, and compare the results with alternatives approaches.

This record has no associated files available for download.

Contributors

Author: J.J. Forster

Download statistics