[논문 리뷰] Time reversal dualities for some random forests

Miraine Dávila Felipe, Amaury Lambert|arXiv (Cornell University)|2014. 09. 21.

Stochastic processes and statistical mechanics참고 문헌 24인용 수 6

한 줄 요약

이 논문은 i.i.d. 분할 수목으로 구성된 랜덤 숲에 대해 시간 역전 이중성(duality)를 수립하며, 초임계 분열 숲의 시간 흐름에 따른 인구 크기 과정이, 타의 생애율과 사망율을 바꾼 초부임계 숲의 시간 흐름에 따른 과정과 동일한 분포를 가짐을 보여준다. 이중성은 윤곽 과정과 레비 과정의 역전 가능성을 통해 유도되며, 전염병 전파의 역학 모델링에 적용 가능하다.

ABSTRACT

We consider a random forest $\mathcal{F}^*$, defined as a sequence of i.i.d. birth-death (BD) trees, each started at time 0 from a single ancestor, stopped at the first tree having survived up to a fixed time $T$. We denote by $\left(ξ^*_t,\ 0\leq t\leq T ight)$ the population size process associated to this forest, and we prove that if the BD trees are supercritical, then the time-reversed process $\left(ξ^*_{T-t},\ 0\leq t\leq T ight)$, has the same distribution as $\left(\widetildeξ^*_t,\ 0\leq t\leq T ight)$, the corresponding population size process of an equally defined forest $\widetilde{\mathcal{F}}^*$, but where the underlying BD trees are subcritical, obtained by swapping birth and death rates or equivalently, conditioning on ultimate extinction. We generalize this result to splitting trees (i.e. life durations of individuals are not necessarily exponential), provided that the i.i.d. lifetimes of the ancestors have a specific explicit distribution, different from that of their descendants. The results are based on an identity between the contour of these random forests truncated up to $T$ and the duality property of Lévy processes. This identity allows us to also derive other useful properties such as the distribution of the population size process conditional on the reconstructed tree of individuals alive at $T$, which has potential applications in epidemiology.

연구 동기 및 목표

• 이중성 결과가 알려진 타임-마이너스 과정으로부터 확장된, 분할 수목로 구성된 랜덤 숲에 대해 시간 역전 이중성을 수립하는 것.
• 초임계 숲의 인구 크기 과정이 시간을 거꾸로 뒤집었을 때, 생애율과 사망율을 바꾼 초부임계 숲의 정방향 과정과 일치하는 조건을 규명하는 것.
• 지속시간이 지수분포가 아닐 경우에도 이중성을 일반화하기 위해 조상 수명의 유일한 분포를 특정하여, 자손 수명과 다름을 명시하는 것.
• 이중성을 전염병 모델링에 적용하여, 생애를 전파로, 사망을 회복 또는 제거로 해석하는 것.
• 원래 수목 간선을 유지하면서 유전적 관계를 재해석하는 경로 기반의 이중 숲 구축 방법 제공

제안 방법

• 저자는 유전적 구조를 인코딩하기 위해 점프-크로노로지컬 윤곽 과정(JCCP)을 사용하며, 이 과정이 스펙트럴로 양의 레비 과정의 법칙을 따름을 보여준다. 이 과정은 흡수 시에 반사되고 종료된다.
• 원래 숲의 JCCP에 공간-시간 역전을 적용하며, 시간 역전 하에서 레비 과정의 알려진 이중성 성질을 활용한다.
• 이중성은 역전된 윤곽 과정이 하위 임계 분열 메커니즘을 가진 이중 숲의 JCCP와 동일한 법칙을 따른다는 것을 증명함으로써 수립된다.
• 핵심 기술적 단계는 경로 변환 하에서 국소 시간(레벨에 대한 방문 횟수)의 불변성을 보이는 것으로, 이는 인구 크기 과정에 대한 이중성이 유지됨을 보장한다.
• 강한 마코프 성질과 지수 기울임을 활용하여, 서로 다른 분열 메커니즘 하에서의 측도 간 관계를 규명한다.
• 비지수 수명을 위해, 조상 수명에 특정한 명시적 분포를 도입함으로써 구축을 확장한다. 이 분포는 자손 수명의 i.i.d. 분포와 다름을 명시한다.

실험 결과

연구 질문

• RQ1 지수 수명이 아닐 경우, 지수 수명 외의 분할 수목로 구성된 랜덤 숲에 대해 시간 역전 이중성을 수립할 수 있는가?
• RQ2 수명이 비지수일 경우, 이중성을 유지하기 위해 조상 수명에 어떤 조건이 필요한가?
• RQ3 원래 숲이 초임계이고 이중 숲이 초부임계일 경우, 이중성은 인구 크기 과정 측면에서 어떻게 나타나는가?
• RQ4 이중성은 인구 크기 과정을 넘어서 전체 숲의 구조에까지 확장될 수 있는가?
• RQ5 이 이중성은 전염병 전파 모델링에 어떤 함의를 갖는가?

주요 결과

• i.i.d. 분할 수목으로 구성된 초임계 랜덤 숲의 시간 역전 인구 크기 과정은, 생애율과 사망율을 바꾼 초부임계 숲의 정방향 인구 크기 과정과 동일한 분포를 가진다.
• 비지수 수명의 경우, 이중성은 조상 수명이 특정한 비지수 분포를 따라야만 성립하며, 이는 자손 수명의 i.i.d. 분포와 다름을 명시해야 한다.
• 이중성은 인구 크기 과정뿐만 아니라 전체 숲의 구조에 대해서도 성립하며, 원래 수목 간선을 유지하면서 유전적 관계를 재해석하는 경로 기반의 이중 숲 구축이 가능하다.
• 숲의 윤곽 과정이 시간을 거꾸로 뒤집었을 때, 이는 국소 시간의 불변성 덕분에 이중 숲의 윤곽 과정과 동일한 법칙을 따른다.
• 이중성은 초임계 과정에서 인구 크기 ≥T를 관측할 확률이, 종식 조건 하에서 초부임계 과정에서 최소 수준 ≤−T일 확률과 동일하다는 것을 시사한다.
• 이 결과는 전염병 전파 모델링을 위한 새로운 도구를 제공하며, 시간 반전 전파 역학이 다른, 초부임계 전파 과정에 대응함을 의미한다.