Bayesian inference for partial orders from random linear extensions: power relations from 12th Century Royal Acta

技術的サマリー：ランダムな線形拡張からの部分順序に対するベイズ推論

問題定義
本論文は、基礎となる構造が必ずしも完全（全）順序ではないランク順序データから、時間とともに変化する社会的ヒエラルキーを推論するという統計的課題に取り組む。著者らは、12 世紀のイングランドおよびノルマンの王領 Acta（法文書）から得られた 371 の目撃者リストを分析し、67 人の司教が記載された順序を記録した。歴史家は、この順序が地位、年功、政治的影響力に基づく社会的ヒエラルキーを反映していると仮説を立てている。しかし、データはノイズを含み、不完全（各リストには司教のサブセットのみが現れる）であり、基礎となるヒエラルキーは単一の線形ランキングではなく、比較不可能な行動者対を含む部分順序（partial order）である可能性がある。マロー（Mallows）やプラケット・ルー（Plackett-Luce）などの既存のランクデータモデルは、通常、完全順序または文脈に依存しない確率を仮定しており、推移性が期待されるが順序が部分的であるというこの歴史的文脈には不適切である可能性がある。

手法
著者らは、潜在状態が進化する部分順序（poset）であり、観測データがその poset のランダムな線形拡張であるベイズ隠れマルコフモデル（HMM）を提案する。

1. 潜在変数のパラメータ化:

   • 時刻 $t$ における poset $h(t)$ は、潜在の $K$ 次元ステータス行列 $Z(t)$ から生成される。
   • 各司教 $j$ について、潜在の権威ベクトル $U^{(t)}_j$ は、1 次のベクトル自己回帰（VAR）として時間とともに進化する。
   • 年功に基づく共変量効果 $\beta_{s_{t,j}}$ が権威ベクトルに追加され、ステータスベクトルを形成する：$Z^{(t)}_j = U^{(t)}_j + \mathbf{1}_K \beta_{s_{t,j}}$。
   • poset $h(t)$ は、$Z(t)$ の列から導かれる $K$ つの全順序の交差によって定義される。司教 $i$ が poset において $j$ に先行するのは、すべての $k=1,\dots,K$ に対して $Z^{(t)}_{i,k} > Z^{(t)}_{j,k}$ である場合に限られる。
   • 次元 $K$ は $\lfloor D/2 \rfloor$（$D$ は活動中の司教の最大数）として選択され、モデルが任意の可能な poset を表現できることを保証する（普遍表現性）。

2. 観測モデル（キュー・ジャンプ）:

   • 観測されたリストは、真の poset の「ノイズの多い」線形拡張としてモデル化される。
   • 著者らは「キュー・ジャンプ」メカニズムを導入する：確率 $p$ で、リストの次の行動者は一様にランダムに選択され（poset の制約を無視する）；そうでない場合、残りの部分順序のランダムな線形拡張に従って行動者が選択される。
   • これにより、観測された順序が基礎となるヒエラルキーと矛盾する「ノイズ」を許容しつつ、有効な線形拡張の方が確率的に高いという制約を維持する。
   • 本モデルは、存在する行動者によって誘発される部分順序に条件付けることで、不完全なリスト（司教のサブセット）を処理する。

3. 事前分布と整合性:

   • poset 過程に対する事前分布は、周辺整合的に構成される。つまり、部分順序（行動者のサブセットによって誘発される）の分布は、完全順序の分布と整合的である。これは、潜在ベクトル $U$ と相関パラメータ $\rho$ に対する特定の事前分布を使用することで達成される。
   • 超パラメータ $\rho$ は poset の「深さ」（最長の鎖の長さ）を制御する。事前分布は深さに関して無情報となるように設計されており、すべての poset に対する一様事前分布で見られる深さ 3 の poset への集中を回避する。
   • 年功効果 $\beta$ は、年功が地位と相関するという歴史的期待を反映するため、順序付けられた制約（$\beta_1 > \beta_2 > \dots$）を課す。ただし、制約のないモデルもテストされている。

4. 推論:

   • 事後推論は、メトロポリス・ヘイスティングス更新を用いたマルコフ連鎖モンテカルロ（MCMC）によって行われる。
   • 尤度は線形拡張の数を数えることを含み、これは #P 完全問題である。著者らは再帰的数え上げアルゴリズム（特に lecount パッケージ）を利用しており、このデータセットにおける行動者の数が少ない（$m \le 22$）場合には実行可能であるが、より大規模な問題へのスケーラビリティを制限する。

主要な貢献

• 新規モデルクラス: 本論文は、基礎となる構造が部分順序である時系列ランクデータをモデル化するための一般的な枠組みを導入し、マローやプラケット・ルーなどの既存モデルを一般化する。
• 周辺整合性: 著者らは、poset に対する潜在変数事前分布が周辺整合的であることを導き、証明した。これは、部分順序に対するアドホックな事前分布でしばしば侵害される性質である。
• 普遍表現性: 潜在次元 $K \ge \lfloor D/2 \rfloor$ により、モデルが任意の可能な poset 時系列を表現できることを示し、複雑なヒエラルキーに対する構造的バイアスを排除する。
• 歴史的応用: これは、12 世紀の王領 Acta 目撃者リストを統計的に分析し、社会的ヒエラルキーを明示的に部分順序としてモデル化して年功効果と個人的権威を分離する最初の研究である。

結果

• モデル適合度: poset ベースのモデルは、適合度（Leave-One-Out クロスバリデーションとベイズ因子で測定）の点で、より単純な代替手段（頂点系列並列順序、バケット順序）およびプラケット・ルー時系列混合モデルを上回る。
• 年功対権威: 分析は、地位が年功に強く依存することを確認し、年功の効果はランク番号が増加するにつれて減少する。しかし、モデルは「権威」（年功からの逸脱）を成功裏に分離する。
• 歴史的洞察:
  • ノルマン人司教: 12 世紀初頭から、アヴランシュ、ベイユー、エヴルー出身のノルマン人司教の地位がイングランド人司教に対して低下していることをモデルが検出する。これは以前に定量化されていなかった傾向である。
  • 政治的影響力: ヘンリー・ド・ブロワ（スティーブン王の弟）やニジェル・オブ・イリーなどの特定の司教の地位は、年功とは無関係に、政治的恩寵または失脚と相関する形で変動する。
  • ヒエラルキー構造: 再構成されたコンセンサス poset は比較的浅く、通常、少数の行動者が上下に位置する中位ランクの司教のいくつかのグループから構成され、深く硬直的なヒエラルキーではない。

意義と主張
本論文は、その主な貢献として、ノイズの多い不完全なランクデータから変化する部分順序を推論するための統計的に厳密な手法の開発であると主張する。この手法を歴史的 Acta に適用することで、著者らは 12 世紀の社会的ヒエラルキーがどのように機能したかについての定量的証拠を提供し、構造的な年功と流動的な政治的権威を区別する。

著者らは、その手法の計算的なスケーラビリティについては謙虚であり、線形拡張の数を数える #P 完全な性質が、モデルを行動者の数が少ない（約 20 人程度）データセットに限定することを認めている。より大規模なデータセットについては、VSP やバケット順序などの制限されたクラスの poset が必要な妥協案である可能性を示唆している。また、モデルはデータによく適合するが、「ノイズ」モデルは誤りが一方方向（キュー・ジャンプ）で発生すると仮定しており、将来の研究では双方向のノイズモデルを検討できる可能性にも言及している。

最終的に、本論文は、標準的なランキングモデルのような完全順序を仮定することは、この種の社会的データに対して正当化されず、部分順序モデルの方が歴史的社会的構造のより正確で解釈可能な表現を提供すると論じている。