[논문 리뷰] On the Existence of the MLE for a Directed Random Graph Network Model with Reciprocation

A. Rinaldo, Sonja Petrović외 1명|arXiv (Cornell University)|2010. 10. 05.

Complex Network Analysis Techniques참고 문헌 25인용 수 6

한 줄 요약

이 논문은 방향성 네트워크에서 상호작용을 고려한 Holland-Leinhardt $p_1$ 모델에서 최대우도추정량(MLE)의 존재를 위한 명시적인 기하 조건을 수립한다. 대수적 통계와 다면체 기하학을 활용하여, 관측된 네트워크 데이터가 모델의 충분통계량의 양성 제약을 위반하는 특정한 0-패턴을 보일 경우, MLE가 존재하지 않음을 규명한다. 특히 이러한 패턴이 행 또는 열의 마진이 0이 아니더라도 모델의 지수족 코너의 상대 내부에 위치할 수 있음을 밝혀낸다.

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)이 존재하는지를 오랫동안 미해결된 문제로 해결한다.
• 토닉 다양체와 지향 다면체 코너를 활용하여 $p_1$ 모델의 지수족의 기하적 구조를 특성화한다.
• 행과 열의 마진이 0이 아니더라도 MLE가 존재하지 않는다는 것을 보여주는 관측된 네트워크 데이터의 특정 0-패턴을 규명한다.
• 설계 행렬과 그에 관련된 코너 구조를 기반으로 MLE 존재 불가를 탐지하기 위한 알고리즘 기준을 제시한다.
• 이론적 연결: $p_1$ 모델에서의 MLE 존재 문제를 준독립성 및 Rasch 모델과 유사한 문제와 연결하여 공통의 조합적 구조를 드러낸다.

제안 방법

• 저자들은 $p_1$ 네트워크를 토닉 지수족으로 모델링하고, 그 충분통계량의 다중동차 구조를 활용하여 토닉 다양체를 정의한다.
• 모델의 설계 행렬 $A$를 분석하고, 모델의 행과 열 마진을 제어하는 매개변수 $\alpha_i$와 $\beta_j$에 대응하는 부분행렬 $B$를 분리한다.
• 다면체 기하학을 활용하여 지향 코너 $\mathrm{cone}(B)$를 연구하고, MLE 존재 여부가 관측 데이터가 이 코너의 상대 내부에 있는지 여부에 달려 있음을 보여준다.
• MLE 존재 불가의 경우를 세 가지로 분류한다: (1) 행 또는 열 마진이 0 또는 $n-1$인 경우, (2) 마진으로 설명되지 않는 0 패턴, (3) 여전히 양성 제약을 위반하는 조합적으로 복잡한 0 패턴.
• 소규모 네트워크($n \leq 10$)를 대상으로 계산기반 대수 소프트웨어(polymake)를 활용하여 분석하고, $\mathrm{cone}(B)$의 면 구조에 대한 추측을 도출한다. 이 코너는 정확히 $3n$개의 면을 가진다: $2n$개는 0 마진에 해당하고, $n$개는 비자명한 0 패턴에 해당한다.
• 모델의 마르코프 기저는 행렬 $B$의 구조와 연결되며, 이는 적합도 검정 및 MLE 존재 여부 확인의 알고리즘적 개선을 가능하게 한다.

실험 결과

연구 질문

• RQ1 방향성 네트워크에서 상호작용을 고려한 $p_1$ 모델에서 최대우도추정량이 존재하지 않는 조건은 무엇인가?
• RQ2 관측된 네트워크 인접행렬의 특정 0-패턴이 행과 열의 마진과는 무관하게 MLE 존재에 어떤 영향을 미치는가?
• RQ3 $p_1$ 모델의 지수족의 기하적 구조는 무엇이며, MLE 존재와 어떻게 관련되어 있는가?
• RQ4 모든 행과 열의 합이 양수이면서 $n-1$ 이하일 때도 MLE가 존재하지 않는 조합적 패턴의 0 간선이 존재하는가?
• RQ5 $p_1$ 모델에서의 MLE 존재 문제는 준독립성 또는 Rasch 모델과 같은 연관 표본 모델의 알려진 문제와 연결될 수 있는가?

주요 결과

• MLE가 존재하지 않는 조건은 관측된 네트워크 데이터에 모델의 충분통계량 공간에서 엄격히 양수인 벡터가 존재하지 않는 구성요소가 포함되어 있을 경우 발생한다. 이는 행과 열의 마진이 0이 아니더라도 해당된다.
• $n=3$일 경우, $\left[\begin{smallmatrix} \times & 0 & \cdot \\ 0 & \times & \cdot \\ \cdot & \cdot & \times \end{smallmatrix}\right]$와 같은 특정 0 패턴이 존재할 경우 MLE가 존재하지 않으며, 이는 0 또는 최대 마진과 대응되지 않는다.
• $n=4$일 경우, 모든 마진이 양수이면서 3 미만임에도 불구하고 MLE가 존재하지 않는 네 가지의 서로 다른 0 패턴이 존재하여, 비자명한 조합적 장애가 있음을 시사한다.
• 저자들은 $\mathrm{cone}(B)$가 정확히 $3n$개의 면을 가진다고 추측한다: $2n$개는 0 행 또는 열 마진에 해당하고, $n$개는 여전히 MLE 실패를 유도하는 비자명한 0 패턴에 해당한다.
• 이 비자명한 0 패턴들은 $2 \times 2$ 연도표 모델에서의 준독립성 모델과 동일한 MLE 존재 불가 원인과 동형이며, 동일한 수의 항목과 피험자가 있는 Rasch 모델에서도 동일한 패턴을 보인다.
• 모델의 마르코프 기저와 MLE 존재 조건을 결정하는 데 핵심적인 역할을 하는 행렬 $B$는 $\alpha_i$와 $\beta_j$ 매개변수를 캡처하고 있음을 입증하였다.