[논문 리뷰] Reciprocal Maximum Likelihood Degrees of Brownian Motion Tree Models

Tobias Boege, Jane Ivy Coons외 3명|arXiv (Cornell University)|2020. 09. 24.

Polynomial and algebraic computation참고 문헌 12인용 수 8

한 줄 요약

이 논문은 토닉 모델과의 연결을 통해 트로닉 피어 프로덕트 분해를 이용하여 브라운 운동 트리 모델의 역최대우도도수(rmld)에 대한 명시적 공식을 유도한다. 주요 결과는 모든 내부 정점 v에 대해 (2^{outdeg(v)} - outdeg(v) - 1)의 곱과 같다는 것으로, 이는 계통발생 모델에서 이러한 기하학적 불변량에 대한 폐쇄형 해를 제공한다.

ABSTRACT

We give an explicit formula for the reciprocal maximum likelihood degree of Brownian motion tree models. To achieve this, we connect them to certain toric (or log-linear) models, and express the Brownian motion tree model of an arbitrary tree as a toric fiber product of star tree models.

연구 동기 및 목표

• 통계적 추정에서의 핵심 불변량인 브라운 운동 트리 모델의 역최대우도도수(rmld)를 결정하는 것.
• 브라운 운동 트리 모델과 토닉(log-linear) 모델 간의 연결을 수립하여 대수기하학적 도구를 활용하는 것.
• 일반적인 트리 모델을 별개의 별형 트리 모델의 토닉 피어 프로덕트로 분해하여 rmld의 재귀적 계산을 가능하게 하는 것.
• 트리의 내부 정점의 출도수에만 의존하는 폐쇄형 공식을 제공하는 것.

제안 방법

• 저자는 브라운 운동 트리 모델을 별형 트리 모델의 토닉 피어 프로덕트로 표현하여 rmld의 재귀적 계산을 가능하게 한다.
• 브라운 운동 모델에서의 역최대우도추정 문제(reciprocal MLE problem)가 해당하는 토닉 모델의 표준 최대우도추정 문제와 동치임을 보인다.
• homogenization과 사영기하학을 이용하여, 토닉 아이디얼과 선형 공간의 교차로 정의된 사영대안의 차수를 통해 rmld를 계산한다.
• 내부 정점의 수에 대한 귀납법을 사용하며, 트리를 별형 부분트리와 축소된 트리로 분해하는 방식으로 증명이 진행된다.
• 주요 대수기하학적 도구로는 농도 행렬, 선형부분공간, 그리고 정의부호가 양인 콘에 있는 대칭행렬의 기하학적 성질을 사용한다.
• 저자는 특히 토닉 피어 프로덕트 이론과 아핀 및 사영대안의 차수 이론을 포함한 대수통계 이론의 결과를 활용한다.

실험 결과

연구 질문

• RQ1 임의의 루트가 있는 트리에 대해 브라운 운동 트리 모델의 역최대우도도수는 무엇인가?
• RQ2 브라운 운동 모델에서의 역최대우도추정 문제는 어떻게 표준 최대우도추정 문제와 관련이 있는가?
• RQ3 트리를 간단한 구성요소로 분해하여 역최대우도도수를 재귀적으로 계산할 수 있는가?
• RQ4 이 모델에서의 역우도임계점의 대수적 구조는 무엇인가?
• RQ5 rmld는 트리의 정점의 차수에 따라 폐쇄형 표현식을 갖는가?

주요 결과

• 루트가 있는 트리 T에 대한 브라운 운동 트리 모델의 역최대우도도수는 모든 내부 정점 v에 대해 (2^{outdeg(v)} - outdeg(v) - 1)의 곱으로 주어진다.
• n개의 잎을 가진 별형 트리의 경우(출도수 n), rmld는 2^n - n - 1이며, 이는 관련 토닉 대안의 차수와 일치한다.
• 그림 1의 트리 모델에서 rmld는 16이다. 두 내부 정점의 출도수가 모두 3이므로 (2^3 - 3 - 1)^2 = 16이 된다.
• 역최대우도추정 문제는 사영대안의 교차가 무한원점에서 점을 가지지 않아 정의가 잘 되어 있으며, 이는 rmld와 일치한다.
• 역최대우도추정 문제의 해는 해당하는 토닉 모델의 표준 최대우도추정 문제와 동치이며, 이는 대수기하학적 도구의 활용을 가능하게 한다.
• 모델의 구조는 토닉 피어 프로덕트를 통한 분해를 허용하며, rmld는 이러한 분해에 대해 곱셈적 성질을 가진다.
• 토닉 아이디얼의 사영적 폐포와 데이터로 정의된 선형공간의 교차는 무한원점에서 점을 가지지 않으며, 이는 차수가 잘 정의져져 있고 rmld와 일치함을 보장한다.