[论文解读] Phase transition in random contingency tables with non-uniform margins

Sam Dittmer, Hanbaek Lyu+1 more|arXiv (Cornell University)|Mar 20, 2019

Markov Chains and Monte Carlo Methods参考文献 26被引用 4

一句话总结

该论文在具有非均匀边缘的随机列联表中建立了极限分布的尖锐相变，表明当 B < B_c 时，左上角条目 𝔼[X₁₁] 的期望值保持一致有界，而当 B > B_c 时，其增长量级为 Θ(n¹⁻⁵)，其中临界阈值为 B_c = 1 + √(1 + 1/C)。相变源于高边缘块与全局约束之间的平衡，在亚临界区域表现出类似几何的特性，而在超临界区域则呈现类似正态的波动。

ABSTRACT

For parameters $n,δ,B,$ and $C$, let $X=(X_{k\ell})$ be the random uniform contingency table whose first $\lfloor n^δ floor $ rows and columns have margin $\lfloor BCn floor$ and the last $n$ rows and columns have margin $\lfloor Cn floor$. For every $0B_{c}$. Our main result shows that $\mathbb{E}[X_{11}]$ is uniformly bounded for $BB_{c}$. We also establish a strong law of large numbers for the row sums in top right and top left blocks.

研究动机与目标

• 分析具有非均匀行和列边缘的随机列联表中条目的渐近行为。
• 识别个体条目分布（特别是 𝔼[X₁₁]）随边缘比参数 B 变化的尖锐相变。
• 确定临界阈值 B_c = 1 + √(1 + 1/C) 如何控制左上角条目期望值从有界到增长的转变。
• 将结果扩展至非均匀边缘情况，特别是 δ > 1/2 的情形，并对 δ ≤ 1/2 的情形提出猜想。
• 为列联表 MCMC 模拟中观察到的相变现象提供严格的理论依据。

提出的方法

• 作者将列联表建模为 𝒫(𝐫,𝐜) 中的均匀随机矩阵，其中前 ⌊n^δ⌋ 行和列的边缘值为 ⌊BCn⌋，其余 n 行和列的边缘值为 ⌊Cn⌋。
• 他们应用概率与组合技术，包括集中不等式和 Borel–Cantelli 论证，分析顶部左上块与右上块的行和波动。
• 分析依赖于表格的分块结构分解，以及条目分布与几何或正态分布之间总变差距离的渐近分析。
• 关键估计涉及在亚临界区域（B < B_c）中，对 X₁₁ 的分布与几何分布之间总变差距离的上界估计。
• 在超临界区域（B > B_c），他们通过尺度变换论证推导出 𝔼[X₁₁] 的渐近界，并证明其增长量级为 C(B − B_c)n¹⁻⁵。
• 证明使用了类似 Borel–Cantelli 的论证，表明归一化行和几乎必然收敛于零，从而支持块和的强大数律。

实验结果

研究问题

• RQ1 在具有非均匀边缘的随机列联表中，当 n → ∞ 时，左上角条目 𝔼[X₁₁] 的期望值如何渐近表现？
• RQ2 X₁₁ 的分布在极限行为上经历尖锐相变的临界值 B_c 是什么？
• RQ3 在亚临界区域（B < B_c）中，X₁₁ 的分布是否在总变差距离下收敛于几何分布？在超临界区域（B > B_c）中，其分布是否呈现类似正态的波动，特别是当 δ > 1/2 时？
• RQ4 该相变如何依赖于控制高边缘块大小的参数 δ？
• RQ5 结果能否推广至 δ = 0（即单个特殊行和列）的情形，如经验模拟和 Barvinok 的猜想所提示？

主要发现

• 当 1/2 < δ < 1 且 B < B_c 时，𝔼[X₁₁] 在 n 上一致有界，表明行为稳定且类似几何分布。
• 当 B > B_c 时，𝔼[X₁₁] 渐近增长为 C(B − B_c)n¹⁻⁵，表明期望值从有界到线性增长的尖锐转变。
• 临界阈值精确为 B_c = 1 + √(1 + 1/C)，其依赖于边缘比参数 C。
• 在亚临界区域（B < B_c），X₁₁ 的分布在总变差距离下收敛于几何分布。
• 在超临界区域（B > B_c），X₁₁ 的分布转变为类似正态的波动，本文猜想当 δ < 1/2 时存在正态极限。
• 对于顶部左上与右上块的归一化行和，强大数律成立，其归一化和几乎必然收敛于零。