MARKOV CHAIN MONTE CARLO AND EXACT CONDITIONAL TESTS WITH THREE-WAY CONTINGENCY TABLES
Atkinson, Michael P.
MetadataShow full item record
We propose an algorithm modifying a popular exact conditional test involving the goodness-of-fit of contingency tables. This study focuses on improving the efficiency of Markov chain Monte Carlo (MCMC) when sampling three-way contingency tables--defined as log-linear models with three discrete random categorical variables consisting of finite levels--under the no-three-way interaction model. Standard to MCMC, we approximate the null distribution by sampling tables from the conditional distribution. However, our proposal involves expanding the conditional state space to include tables with cell count values of -1. We apply the proposed methodology, described in full detail, to randomly generated sparse and non-sparse data sets. Our results show that traditional asymptotic methods on sparse contingency tables yield inaccurate results. We also prove mathematically that a Markov chain with our proposed method is connected (i.e., ergodic) on the conditional state space for 3x3xK, with K >= 3. The output from applying the proposed methodology provides conclusive evidence that the distribution of the test statistics for sparse data sets does not resemble the asymptotic distribution.
Approved for public release. distribution is unlimited
Showing items related by title, author, creator and subject.
HYBRID SIS AND MARKOV CHAIN MONTE CARLO SAMPLING METHODOLOGY FOR GOODNESS-OF-FIT TESTS ON CONTINGENCY TABLES Saluke, Patrick M. (Monterey, CA; Naval Postgraduate School, 2018-09);Logistic regression is one of the most popular means of modeling contingency table data due to its ease of use. Simple asymptotic inference (like a χ2 approximation) for evaluating goodness-of-fit tests, however, may not ...
Stone, Aaron (Monterey, CA; Naval Postgraduate School, 2018-06);We develop an algorithm blending Sequential Importance Sampling (SIS) and Markov Chain Monte Carlo (MCMC) to conduct goodness of fit testing on three-way contingency tables under the no-three-way interaction model. Unlike ...
Kahle, David; Yoshida, Ruriko; Garcia-Puente, Luis (Springer, 2017-09);Exact conditional goodness-of-fit tests for discret eexponential family models can be conducted via Monte Carlo estimation of p values by sampling from the conditional distribution of multiway contingency tables. The two ...