A Markov basis for conditional test of common diagonal effect in quasi-independence model for square contingency tables

In two-way contingency tables we sometimes find that frequencies along the diagonal cells are relatively larger (or smaller) compared to off-diagonal cells, particularly in square tables with the common categories for the rows and the columns. In this case the quasi-independence model with an additional parameter for each of the diagonal cells is usually fitted to the data. A simpler model than the quasi-independence model is to assume a common additional parameter for all the diagonal cells. We consider testing the goodness of fit of the common diagonal effect by the Markov chain Monte Carlo (MCMC) method. We derive an explicit form of a Markov basis for performing the conditional test of the common diagonal effect. Once a Markov basis is given, MCMC procedure can be easily implemented by techniques of algebraic statistics. We illustrate the procedure with some real data sets.

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

[2]  Fabio Rapallo,et al.  Markov bases and structural zeros , 2006, J. Symb. Comput..

[3]  Akimichi Takemura,et al.  MATHEMATICAL ENGINEERING TECHNICAL REPORTS Markov Bases for Two-way Subtable Sum Problems , 2007, 0708.2312.

[4]  A. Dobra,et al.  Monte Carlo Algorithms for Hardy–Weinberg Proportions , 2006, Biometrics.

[5]  Alan Agresti,et al.  Categorical Data Analysis , 1991, International Encyclopedia of Statistical Science.

[6]  Seth Sullivant,et al.  A Divide-and-Conquer Algorithm for Generating Markov Bases of Multi-way Tables , 2004, Comput. Stat..

[7]  Seth Sullivant,et al.  Gröbner Bases and Polyhedral Geometry of Reducible and Cyclic Models , 2002, J. Comb. Theory, Ser. A.

[8]  Satoshi Kuriki,et al.  Testing marginal homogeneity against stochastically ordered marginals for r × r contingency tables , 2006 .

[9]  E. Lehmann Testing Statistical Hypotheses , 1960 .

[10]  Satoshi Aoki,et al.  Minimal and minimal invariant Markov bases of decomposable models for contingency tables , 2010 .

[11]  Satoshi Aoki,et al.  Markov chain Monte Carlo exact tests for incomplete two-way contingency tables , 2005 .

[12]  Shelby J. Haberman,et al.  A Warning on the Use of Chi-Squared Statistics with Frequency Tables with Small Expected Cell Counts , 1988 .

[13]  D. Geiger,et al.  On the toric algebra of graphical models , 2006, math/0608054.

[14]  A. Takemura,et al.  Minimal Basis for a Connected Markov Chain over 3 × 3 ×K Contingency Tables with Fixed Two‐Dimensional Marginals , 2003 .

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

[16]  B. Sturmfels Gröbner bases and convex polytopes , 1995 .

[17]  Satoshi Aoki,et al.  Fibers of sample size two of hierarchical models and Markov bases of decomposable models for contingency tables , 2007 .

[18]  Stephen E. Fienberg,et al.  Discrete Multivariate Analysis: Theory and Practice , 1976 .

[19]  Satoshi Aoki,et al.  Distance-reducing Markov bases for sampling from a discrete sample space , 2005 .