[論文レビュー] Non-reversible Markov chain Monte Carlo for sampling of districting maps

Gregory Herschlag, Jonathan C. Mattingly+2 more|arXiv (Cornell University)|Aug 18, 2020

Markov Chains and Monte Carlo Methods参考文献 25被引用数 5

ひとこと要約

本稿では、重心移動フローと地域間遷移に基づく新しい非可逆ダイナミクスを用いて、非偏見の赤字計画を効率的にサンプリングする非可逆マルコフ連鎖モンテカルロ（MCMC）手法を導入する。複数の運動量変数を用いた一般化された歪詳細釣合条件を用いることで、可逆な手法と比較して著しく高速な混合（mixing）を達成し、50,000ステップで相関が低下するのに対し、標準的手法では数万ステップ以上を要する。

ABSTRACT

Evaluating the degree of partisan districting (Gerrymandering) in a statistical framework typically requires an ensemble of districting plans which are drawn from a prescribed probability distribution that adheres to a realistic and non-partisan criteria. In this article we introduce novel non-reversible Markov chain Monte-Carlo (MCMC) methods for the sampling of such districting plans which have improved mixing properties in comparison to previously used (reversible) MCMC algorithms. In doing so we extend the current framework for construction of non-reversible Markov chains on discrete sampling spaces by considering a generalization of skew detailed balance. We provide a detailed description of the proposed algorithms and evaluate their performance in numerical experiments.

研究の動機と目的

• ギャングリーマンディングを検出するための、効率的で非党派的な赤字計画のサンプリングのニーズに対処すること。
• 非可逆ダイナミクスを活用することで、赤字計画のMCMCサンプラーにおける混合性能を向上させること。
• 一般化された歪詳細釣合条件を用いた非可逆MCMCの数学的枠組みを、離散空間上で構築すること。
• 相関減衰とハミング距離の指標を用いて、実世界の赤字計画問題における性能を評価すること。
• ギャングリーマンディング分析におけるスケーラブルで原理的根拠のある非可逆MCMCの基盤を提供すること。

提案手法

• 複数の運動量変数を用いた歪詳細釣合条件の一般化により、離散状態空間上での非可逆MCMCを構築する。
• 一般化された歪詳細釣合条件の下で、エルゴディック性と正しい定常分布を保証するため、混合歪メトロポリス＝ハスティングス法を提案する。
• 地域再割り当てを効果的に進めるために、ベクトル場に基づく重心移動フローを採用し、ランダムウォーク的挙動を低減することで混合性能を向上させる。
• 個々のノードの反転に制限されない、全地域間の遷移を可能にする地域間フロー機構を設計する。
• 単位運動量を用いた状態空間の拡張により、非可逆ダイナミクスを埋め込みつつ、目的測度を保存する。
• 1000万ステップのチェーンを用いたブートストラップサンプリングにより、相関減衰（G(t)）を推定し、混合速度を評価する。

実験結果

リサーチクエスチョン

• RQ1 非可逆MCMC手法は、可逆メトロポリス＝ハスティングス法と比較して、赤字計画のサンプリングにおける混合性能を著しく向上させることができるか？
• RQ2 赤字計画グラフのような離散的・構造的空間上に、非可逆ダイナミクスを体系的に構築する方法は何か？
• RQ3 ベクトル場の選定と運動量の割り当てが、混合効率に果たす役割は何か？
• RQ4 複数の運動量を用いた一般化された歪詳細釣合条件を用いて、組合せ的サンプリングに適した非可逆MCMCアルゴリズムを設計できるか？
• RQ5 重心移動に基づく非可逆フローは、従来のノード単位の反転法やマージ・スプリット法と比較して、相関減衰の観点でどのように異なるか？

主な発見

• 重心移動フロー法では、50,000ステップ未満で相関が低下し、他の手法と比較して著しく高速である。
• 地域間遷移法と温度調整付きノード単位反転法は、ほぼ同じ速度で相関が低下し、いずれも標準的なノード単位反転法よりも高速である。
• 標準的なノード単位反転アルゴリズムは、独立性に達するまでに著しく多くのステップを要し、最も遅い相関減衰を示す。
• 4つの手法（重心移動フローを含む）すべてが50,000ステップ以内に相関が低下しており、効果的な混合が達成されていることが示された。
• 提案された非可逆MCMCフレームワークにより、赤字計画アンサンブルの収束が高速化され、サンプリング効率が向上した。
• 混合歪メトロポリス＝ハスティングス法は、一般化された歪詳細釣合条件の下で、エルゴディック性と正しい定常分布を保証する。