Computing p-values in conditional independence models for a contingency table

We present a Markov chain Monte Carlo (MCMC) method for generating Markov chains using Markov bases for conditional independence models for a four-way contingency table. We then describe a Markov basis characterized by Markov properties associated with a given conditional independence model and show how to use the Markov basis to generate random tables of a Markov chain. The estimates of exact p-values can be obtained from random tables generated by the MCMC method. Numerical experiments examine the performance of the proposed MCMC method in comparison with the χ2 approximation using large sparse contingency tables.

[1]  Ronald Christensen,et al.  Log-Linear Models and Logistic Regression , 1997 .

[2]  B. Everitt.  Log-linear models for contingency tables , 1992 .

[3]  Michel Mouchart,et al.  Discussion on "Conditional independence in statistitical theory" by A.P. Dawid , 1979 .

[4]  A. Dawid.  Conditional Independence in Statistical Theory , 1979 .

[5]  D. Edwards.  Introduction to graphical modelling , 1995 .

[6]  N. Jaffe,et al.  Weekly high-dose methotrexate and doxorubicin for osteosarcoma: the Dana-Farber Cancer Institute/the Children's Hospital--study III. , 1987, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[7]  Charles J. Geyer,et al.  Practical Markov Chain Monte Carlo , 1992 .

[8]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[9]  Jonathan J. Forster,et al.  Monte Carlo exact conditional tests for log-linear and logistic models , 1996 .

[10]  A Agresti,et al.  Exact inference for categorical data: recent advances and continuing controversies , 2001, Statistics in medicine.

[11]  S. Reiss,et al.  Data-swapping: A technique for disclosure control , 1982 .

[12]  A. Dobra.  Markov bases for decomposable graphical models , 2003 .

[13]  Michael I. Jordan.  Graphical Models , 2003 .

[14]  Xu Yu-bo,et al.  Reachability Checking of Finite Precision Timed Automata , 2006 .

[15]  P. Diaconis,et al.  Algebraic algorithms for sampling from conditional distributions , 1998 .

[16]  T. Speed,et al.  Markov Fields and Log-Linear Interaction Models for Contingency Tables , 1980 .

[17]  Brian S. Caffo,et al.  A Markov chain Monte Carlo Algorithm for Approximating Exact Conditional Probabilities , 2001 .

[18]  Brian Cao.  Exact Hypothesis Tests for Log-linear Models with exactLoglinTest , 2006 .

[19]  A. Agresti.  [A Survey of Exact Inference for Contingency Tables]: Rejoinder , 1992 .

[20]  Ronald W. Butler,et al.  An importance sampling algorithm for exact conditional tests in log-linear models , 1999 .

[21]  Brian S. Caffo,et al.  Exact Hypothesis Tests for Log-linear Models with exactLoglinTest , 2006 .