Exact Goodness‐of‐Fit Testing for the Ising Model

The Ising model is one of the simplest and most famous models of interacting systems. It was originally proposed to model ferromagnetic interactions in statistical physics and is now widely used to model spatial processes in many areas such as ecology, sociology, and genetics, usually without testing its goodness of fit. Here, we propose various test statistics and an exact goodness‐of‐fit test for the finite‐lattice Ising model. The theory of Markov bases has been developed in algebraic statistics for exact goodness‐of‐fit testing using a Monte Carlo approach. However, finding a Markov basis is often computationally intractable. Thus, we develop a Monte Carlo method for exact goodness‐of‐fit testing for the Ising model that avoids computing a Markov basis and also leads to a better connectivity of the Markov chain and hence to a faster convergence. We show how this method can be applied to analyze the spatial organization of receptors on the cell membrane.

[1]  E. Ising.  Beitrag zur Theorie des Ferromagnetismus , 1925 .

[2]  S. Brush.  History of the Lenz-Ising Model , 1967 .

[3]  D. Vere-Jones.  Markov Chains , 1972, Nature.

[4]  J. Besag.  Nearest‐Neighbour Systems and the Auto‐Logistic Model for Binary Data , 1972 .

[5]  J. Besag.  Spatial Interaction and the Statistical Analysis of Lattice Systems , 1974 .

[6]  S. Galam,et al.  Sociophysics: A new approach of sociological collective behaviour. I. mean‐behaviour description of a strike , 1982, 2211.07041.

[7]  J. Besag.  On the Statistical Analysis of Dirty Pictures , 1986 .

[8]  Donald Geman,et al.  Bayesian Image Analysis , 1986 .

[9]  J. Besag,et al.  Generalized Monte Carlo significance tests , 1989 .

[10]  S. Galam,et al.  Towards a theory of collective phenomena: Consensus and attitude changes in groups , 1991 .

[11]  이기수,et al.  II. , 1992 .

[12]  Julian Besag,et al.  Towards Bayesian image analysis , 1993 .

[13]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .

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

[15]  J. Ott,et al.  The Ising model in physics and statistical genetics. , 2001, American journal of human genetics.

[16]  Fabio Rapallo.  Algebraic Markov Bases and MCMC for Two‐Way Contingency Tables , 2003 .

[17]  Michael J. Berry,et al.  Ising models for networks of real neurons , 2006, q-bio/0611072.

[18]  Michael J. Berry,et al.  Weak pairwise correlations imply strongly correlated network states in a neural population , 2005, Nature.

[19]  Seth Sullivant,et al.  A finiteness theorem for Markov bases of hierarchical models , 2007, J. Comb. Theory, Ser. A.

[20]  Seth Sullivant,et al.  Lectures on Algebraic Statistics , 2008 .

[21]  Stephen E. Fienberg,et al.  Algebraic Statistics for a Directed Random Graph Model with Reciprocation , 2009, 0909.0073.

[22]  Ruriko Yoshida,et al.  Algebraic and Geometric Methods in Statistics: Markov chains, quotient ideals and connectivity with positive margins , 2009 .

[23]  Ruriko Yoshida,et al.  Markov bases and subbases for bounded contingency tables , 2009, 0905.4841.

[24]  Giovanni Pistone,et al.  The algebra of reversible Markov chains , 2010, 1007.4282.

[25]  Giovanni Pistone,et al.  Toric Statistical Models: Ising and Markov , 2011, 1110.2544.

[26]  Adrian Dobra,et al.  Dynamic Markov Bases , 2011, 1103.4891.

[27]  Hisayuki Hara,et al.  Graver basis for an undirected graph and its application to testing the beta model of random graphs , 2011, 1102.2583.

[28]  青木 敏,et al.  Markov bases in algebraic statistics , 2012 .

[29]  Michael J. Berry,et al.  Searching for Collective Behavior in a Large Network of Sensory Neurons , 2013, PLoS Comput. Biol..

[30]  A. Takemura,et al.  Markov degree of configurations defined by fibers of a configuration , 2014, 1405.2676.

[31]  Elizabeth Gross,et al.  Goodness of fit for log-linear network models: dynamic Markov bases using hypergraphs , 2014, 1401.4896.

[32]  Jianhua Guo,et al.  A simplification of computing Markov bases for graphical models whose underlying graphs are suspensions of graphs , 2014 .

[33]  A. Slavkovic,et al.  Fibers of multi-way contingency tables given conditionals: relation to marginals, cell bounds and Markov bases , 2014, Annals of the Institute of Statistical Mathematics.

[34]  Seth Sullivant,et al.  Lifting Markov bases and higher codimension toric fiber products , 2014, J. Symb. Comput..