[논문 리뷰] Moments of the frequency spectrum of a splitting tree with neutral Poissonian mutations

Nicolas Champagnat, Benoît Henry|arXiv (Cornell University)|2015. 09. 22.

Stochastic processes and statistical mechanics참고 문헌 17인용 수 9

한 줄 요약

이 논문은 중립적 희귀 돌연변이를 가진 초임계 분열 수목에서 유전자 빈도 스펙트럼의 계승 모멘트에 대한 명시적 재귀 공식을 유도한다. 새로운 공통점 과정의 구성과 일반화된 캄브렐 공식을 도입함으로써, 적절한 지수 시간 스케일링 후 빈도 스펙트럼의 거의 확실 수렴을 입증하며, 일반 분열 과정 이론에 의존하지 않는 새로운이고 간단한 증명을 제공한다.

ABSTRACT

We consider a branching population where individuals live and reproduce independently. Their lifetimes are i.i.d. and they give birth at a constant rate b. The genealogical tree spanned by this process is called a splitting tree, and the population counting process is a homogeneous, binary Crump-Mode-Jagers process. We suppose that mutations affect individuals independently at a constant rate $\\\ heta$ during their lifetimes, under the infinite-alleles assumption: each new mutation gives a new type, called allele, to his carrier. We study the allele frequency spectrum which is the numbers A(k, t) of types represented by k alive individuals in the population at time t. Thanks to a new construction of the coalescent point process describing the genealogy of individuals in the splitting tree, we are able to compute recursively all joint factorial moments of (A(k, t)) k$\\ge$1. These moments allow us to give an elementary proof of the almost sure convergence of the frequency spectrum in a supercritical splitting tree.

연구 동기 및 목표

• 분열 수목에서 중립적 희귀 돌연변이를 가진 경우의 빈도 스펙트럼의 계승 모멘트에 대한 명시적 재귀 공식을 유도한다.
• 지역 독립성 구조를 통해 모멘트 계산이 가능한 새로운 공통점 과정(CPP)의 구성법을 개발한다.
• 일반 분열 과정 이론에 의존하지 않고 초임계 영역에서 빈도 스펙트럼의 거의 확실 수렴을 증명하는 간단한 방법을 제공한다.
• 고차수 渐近 근사, 예를 들어 중심극한정리 유사 결과를 위한 정보를 제공할 수 있는 모멘트 공식을 수립한다.

제안 방법

• 공통점 과정을 서로 독립적인 동일분포를 가진 랜덤 변수의 수열로 표현함으로써, 재귀적 모멘트 계산이 가능하도록 한다.
• 지역 독립성 하에서 랜덤 측도에 대한 적분의 기대값을 위한 일반화된 캄브렐 공식을 개발한다.
• 공통점 과정의 공동 모멘트 구조를 활용하여 빈도 스펙트럼 $A(k,t)$의 계승 모멘트를 재귀적으로 계산한다.
• 모멘트 공식을 적용하여 $e^{-\alpha t} A(k,t)$가 $\mathcal{E}/\psi'(\alpha) \cdot c_k$ 비례하는 극한으로 거의 확실 수렴함을 증명한다. 여기서 $\mathcal{E}$는 지수 분포를 따르는 랜덤 변수이다.
• 독립 레비 과정을 이용한 윤곽 과정 구성법을 통해 잔류 수명을 특성화하고 공통점 과정에서의 독립성을 확립한다.
• 레비 과정의 척도 함수와 라플라스 지수를 활용하여 인구와 돌연변이 과정의 장기적 행동을 분석한다.

실험 결과

연구 질문

• RQ1 분열 수목에서 희귀 돌연변이를 가진 경우, 빈도 스펙트럼 $A(k,t)$의 공동 계승 모멘트는 어떻게 재귀적으로 계산할 수 있는가?
• RQ2 공통점 과정은 시간 $t$에서 인구의 계통수학적 구조를 어떻게 특성화하는가?
• RQ3 일반 분열 과정 이론을 활용하지 않고 모멘트 방법을 통해 빈도 스펙트럼의 거의 확실 수렴을 증명할 수 있는가?
• RQ4 초임계 경우에서 $e^{-\alpha t} A(k,t)$의 정확한 극한 분포는 무엇이며, 이는 인구의 마르타시안 파라미터와 어떻게 관련이 있는가?
• RQ5 시간 $t$에서 개인의 잔류 수명은 공통점 과정의 독립성 구조에 어떻게 기여하는가?

주요 결과

• 논문은 $A(k,t)$, 즉 시간 $t$에 $k$명의 개체가 생존해 있는 유전자 수를 나타내는 변수의 계승 모멘트에 대한 명시적 재귀 공식을 유도한다. 이는 더 낮은 순서의 모멘트에 기반한다.
• $e^{-\alpha t} A(k,t)$가 $\frac{\mathcal{E}}{\psi'(\alpha)} c_k$로 거의 확실 수렴함을 입증하였다. 여기서 $\mathcal{E}$는 지수 분포를 따르는 랜덤 변수이며, $c_k$는 모델 파라미터에 따라 결정되는 계산 가능한 상수이다.
• 수렴 속도가 지수적으로 빠르며, $\mathbb{E}_{t}[e^{-2\alpha t}(c_k N_t - A(k,t))^2] \sim \mathcal{C} e^{-\gamma t}$로 표현되며, 여기서 $\gamma = \theta$는 초임계/임계일 때이고, $\gamma = 2\alpha - \theta$는 초임계일 때이다.
• 시간 $t$에서 개인의 잔류 수명이 상호 독립임을 입증하였으며, 첫 번째 잔류 수명은 나머지와 독립적이다. 이는 독립 레비 과정을 이용한 윤곽 과정 구성법에 기반한다.
• 새로운 공통점 과정 구성법을 통해 일반 분열 과정 이론에 기반한 랜덤 특성 수세기 이론을 활용하지 않고도 직접적으로 모멘트 공식을 도출할 수 있다.
• 이 방법은 빈도 스펙트럼 수렴의 오차를 정량화함으로써, 기능 중심극한정리와 같은 고차수 渐近 결과를 도출하는 기초를 제공한다.