Sampling contingency tables

From MaRDI portal





This article is devoted to the study of random sampling from a set of contingency tables. First, the authors show that the problem of determining the exact number of contingency tables is \#P-complete. They prove this result by establishing a relationship between their problem and the computation of the \((n-1)\)-dimensional volume of the polytope \(\sum_{j=1}^{n-1} a_jy_j\leq b\), where \(0\leq y_j\leq 1\) \((j=1,\dots, n-1)\) and \(a_1,\dots, a_{n-1},b\) are given positive integers. Then they reduce the problem of sampling from the discrete set of contingency tables to the problem of sampling with near-uniform density from a contingency polytope. They show that the problem of random sampling can be solved in polynomial time if the row and column sums are large enough. This provides, under the same hypothesis as above, a polynomial time algorithm to solve the problem of counting the number of contingency tables approximately. The algorithm that they present is based on a random walk approach and it is related to volume computation from log-concave distributions.




Cited in
(45)








This page was built for publication: Sampling contingency tables

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4345364)