Statistics and Its Interface
Volume 13 (2020)
Sampling high dimensional tables with applications to assessing linkage disequilibrium
Pages: 157 – 166
We propose a sequential importance sampling strategy to sample high dimensional tables with fixed one way margins. The proposal distribution for the method is constructed by adapting an approximation to the total number of tables available in the literature. We apply the method to estimating the total number of tables and assessing linkage disequilibrium in multimarker genetic data with the table representing haplotype count data. We demonstrate efficient and accurate performance in these practical, real-world examples. The method may be applied in any situation in which uniformly sampling high dimensional tables with fixed one way margins is of interest. Detailed derivations are provided in the appendix.
counting problem, exact test, high dimensional table, linkage disequilibrium, Monte Carlo method, sequential importance sampling
Yuguo Chen was partially supported by the NSF grant DMS-1406455.
Received 12 April 2019
Received revised 27 July 2019
Accepted 19 September 2019