Abstract
A Markov chain is proposed that uses coupling from the past sampling algorithm for sampling m×n contingency tables. This method is an extension of the one proposed by Kijima and Matsui (Rand. Struct. Alg., 29:243–256, 2006). It is not polynomial, as it is based upon a recursion, and includes a rejection phase but can be used for practical purposes on small contingency tables as illustrated in a classical 4×4 example.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Chen, Y., Diaconis, P., Holmes, S., Liu, J.S.: Sequential Monte Carlo methods for statistical analysis of tables. J. Am. Stat. Assoc. 100, 109–120 (2005)
Cryan, M., Dyer, M., Goldberg, L.A., Jerrum, M., Martin, R.: Rapidly mixing Markov chains for sampling contingency tables with constant number of rows. In: Proceedings of the 43rd Annual Symposium on Foundations of Computer Science (FOCS 2002), pp. 711–720 (2002)
Diaconis, P., Efron, B.: Testing for independence in a two-way table: New interpretations of the chi-square statistic (with discussion). Ann. Stat. 13, 845–913 (1985)
Diaconis, P., Saloff-Coste, L.: Random walk on contingency tables with fixed row and column sums. Technical report, Department of Mathematics, Harvard University (1995)
Fisher, R.A.: The logic of inductive inference (with discussion). J. R. Stat. Soc. 98, 39–54 (1935)
Kijima, S., Matsui, T.: Polynomial time perfect sampling algorithm for two-rowed contingency tables. Rand. Struct. Algorithms 29, 243–256 (2006)
Propp, J., Wilson, D.: How to get a perfectly random sample from a generic Markov chain and generate a random spanning tree of a directed graph. Rand. Struct. Algorithms 9, 223–252 (1996)
Snee, R.D.: Graphical display of two-way contingency tables. Am. Stat. 38, 9–12 (1974)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Wicker, N. Perfect sampling algorithm for small m×n contingency tables. Stat Comput 20, 57–61 (2010). https://doi.org/10.1007/s11222-009-9115-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11222-009-9115-1