[论文解读] On the Existence of the MLE for a Directed Random Graph Network Model with Reciprocation

A. Rinaldo, Sonja Petrović+1 more|arXiv (Cornell University)|Oct 5, 2010

Complex Network Analysis Techniques参考文献 25被引用 6

一句话总结

本文為具有互惠關係的有向網絡的Holland-Leinhardt $p_1$模型確立了最大似然估計量（MLE）存在的明確幾何條件。利用代數統計與多面體幾何，本文指出當觀察到的網絡數據在邊緣配置中表現出特定的零模式，違反了模型充分統計量的正性約束時，MLE會失效，特別是在這些模式不對應於零行或列邊際但仍位於模型指數族錐的相對內部時。

ABSTRACT

Holland and Leinhardt (1981) proposed a directed random graph model, the p1 model, to describe dyadic interactions in a social network. In previous work (Petrovic et al., 2010), we studied the algebraic properties of the p1 model and showed that it is a toric model specified by a multi-homogeneous ideal. We conducted an extensive study of the Markov bases for p1 that incorporate explicitly the constraint arising from multi-homogeneity. Here we consider the properties of the corresponding toric variety and relate them to the conditions for the existence of the maximum likelihood and extended maximum likelihood estimators or the model parameters. Our results are directly relevant to the estimation and conditional goodness-of-fit testing problems in p1 models.

研究动机与目标

• 解決有向網絡互惠$ p_1 $模型中最大似然估計量（MLE）存在性問題的長期懸而未決問題。
• 利用扭稜簇與 pointed 多面體錐，表徵$ p_1 $模型指數族的幾何結構。
• 識別觀察到的網絡數據中導致MLE不存在的特定零模式，即使行與列邊際非零。
• 基於設計矩陣及其相關錐結構，提供檢測MLE不存在性的算法標準。
• 將$ p_1 $模型中MLE存在性問題與擬獨立性與Rasch模型中的類似問題聯繫起來，揭示其共通的組合結構。

提出的方法

• 作者將$ p_1 $網絡建模為扭稜指數族，利用其充分統計量的多重齊次結構定義扭稜簇。
• 他們分析模型的設計矩陣$ A $，並提取對應於參數$\alpha_i$與$\beta_j$的子矩陣$ B $，該參數控制邊緣計數的邊際。
• 利用多面體幾何，他們研究 pointed 錐$\mathrm{cone}(B)$，並證明MLE的存在性取決於觀察數據是否位於該錐的相對內部。
• 他們識別出MLE不存在的三種不同情形：(1) 行或列邊際完全為0或$n-1$；(2) 無法由邊際捕獲的零模式；(3) 雖然複雜但仍違反正性的組合零模式。
• 他們使用計算代數軟件（polymake）分析小型網絡（$n \leq 10$），並推導出關於$\mathrm{cone}(B)$的面結構的猜想，該錐有$3n$個面：$2n$個來自零邊際，$n$個來自非平凡零模式。
• 模型的馬爾可夫基底與$B$的結構相關聯，從而實現了適合度檢驗與MLE存在性檢測的算法改進。

实验结果

研究问题

• RQ1 在具有互惠關係的有向網絡$ p_1 $模型中，MLE在何種條件下會不存在？
• RQ2 觀察到的網絡鄰接矩陣中的特定零模式如何影響MLE的存在性，而不受行與列邊際的影響？
• RQ3 $ p_1 $模型指數族的幾何結構是什麼？它與MLE存在性的關係為何？
• RQ4 是否存在零邊的組合模式，即使所有行與列和均為正且小於$n-1$，仍會導致MLE不存在？
• RQ5 $ p_1 $模型中的MLE存在性問題能否與已知的列聯表模型問題（如擬獨立性或Rasch模型）聯繫起來？

主要发现

• 當觀察到的網絡數據在阻止模型充分統計量空間中存在嚴格正向量的配置中包含零值時，即使行與列邊際非零，$ p_1 $模型的MLE也不存在。
• 對於$n=3$，MLE不存在的特定零模式如$\left[\begin{smallmatrix} \times & 0 & \cdot \\ 0 & \times & \cdot \\ \cdot & \cdot & \times \end{smallmatrix}\right]$，這些模式不對應於零或最大邊際。
• 對於$n=4$，存在四種不同的零模式導致MLE不存在，即使所有邊際均為正且小於3，表明存在非平凡的組合障礙。
• 作者猜想錐$\mathrm{cone}(B)$恰好有$3n$個面：$2n$個對應於零行或列邊際，$n$個對應於仍導致MLE失敗的非平凡零模式。
• 這些非平凡零模式與$2 \times 2$列聯表中擬獨立性模型以及項目與受試者數量相等的Rasch模型中導致MLE不存在的模式同構。
• 矩陣$B$（捕獲$\alpha_i$與$\beta_j$參數）被證明是決定模型馬爾可夫基底與MLE存在條件的核心。