[논문 리뷰] Markov branching processes with disasters: extinction, survival and duality to p-jump processes

Felix Hermann, Peter Pfaffelhuber|arXiv (Cornell University)|2018. 07. 31.

Stochastic processes and statistical mechanics참고 문헌 32인용 수 3

한 줄 요약

이 논문은 이항 재해를 동반한 마르코프 분열 과정과 p-점프 과정—[0,1]에서의 곱셈적 점프를 갖는 조각별 결정론적 마르코프 과정—사이의 이중성을 수립한다. 이 이중성과 p-점프 과정에 대한 극한 정리들을 활용하여, 재해가 있는 시간에 동일한 분열 과정과 시간에 비례하지 않는 분열 과정에서의 멸종 및 생존 확률에 대한 정확한 점근적 표현을 도출하며, 이는 감쇠 속도의 계단적 전이와 거의 확실한 멸종 또는 생존 조건을 드러낸다.

ABSTRACT

A $p$-jump process is a piecewise deterministic Markov process with jumps by a factor of $p$. We prove a limit theorem for such processes on the unit interval. Via duality with respect to probability generating functions, we deduce limiting results for the survival probabilities of time-homogeneous branching processes with arbitrary offspring distributions, underlying binomial disasters. Extending this method, we obtain corresponding results for time-inhomogeneous birth-death processes underlying time-dependent binomial disasters and continuous state branching processes with $p$-jumps.

연구 동기 및 목표

• 재해 발생 후 각 개체가 확률 p로 생존하는 조건에서 이항 재해를 동반한 마르코프 분열 과정의 장기적 행동을 분석하는 것.
• 시간에 비례하지 않는 비율과 연속 상태 분열 과정을 포함하여 기존의 멸종 및 생존 확률 결과를 확장하는 것.
• 확률 생성 함수를 활용한 이중성 프레임워크를 개발하여 p-점프 과정의 성질을 분열 과정의 결과로 변환하는 것.
• 일반적인 후손 및 재해 역학 하에서 생존 및 멸종 확률에 대한 정밀한 점근적 표현을 도출하는 것.
• 분열, 재해, 생존 파라미터 간의 상호작용에 기반한 생존 확률 감쇠 속도의 계단적 전이를 특성화하는 것.

제안 방법

• 재해 시점에 요인 p로 점프하는 [0,1]에서의 조각별 결정론적 마르코프 과정으로서 p-점프 과정을 도입한다.
• 확률 생성 함수를 통한 이중성으로 분열 과정과 p-점프 과정 사이의 관계를 수립하며, 분열 과정의 생성 함수를 확률적 미분 방정식의 해와 연결한다.
• 이중 과정에 대한 핵심 SDE 유도: $ \dot{X} = -\lambda(X - h(X)) $, 재해 빈도 $ \kappa $에서 요인 p로 점프가 포함된 형태로 수정한다.
• 데모보와 짐투니(1998)의 포아송 과정에 대한 대규모 탈락 이론을 활용하여 누적 점프 강도 및 생존 확률의 尾행동을 분석한다.
• 점프 강도 및 분열 강도에 대한 오목성 및 성장 조건 하에서 p-점프 과정에 대한 극한 정리들을 적용한다.
• 이중성 관계 $ \mathbb{E}[x^{Z_t} \mid Z_0 = z] = X_t^z $ 를 활용하여 $ X_t $ 의 수렴성을 $ Z_t $ 의 멸종/생존 확률로 번역한다.

실험 결과

연구 질문

• RQ1 이항 재해와 임의의 후손 분포를 가진 시간에 동일한 분열 과정에서 생존 확률의 점근적 행동은 무엇인가?
• RQ2 시간에 비례하지 않는 분열 및 재해 빈도는 인구의 장기적 생존 또는 멸종에 어떻게 영향을 미치는가?
• RQ3 어떤 조건에서 분열 과정이 재해를 동반하여 거의 확실하게 멸종하거나 멸종하지 않을 수 있으며, 생존 확률의 감쇠 속도는 어떻게 변화하는가?
• RQ4 p-점프 과정과 분열 과정 사이의 이중성은 연속 상태 분열 과정에 재해를 동반한 경우로 확장될 수 있는가?
• RQ5 재해 시간 기록 조건 하에서 생존 확률의 정확한 수렴 속도는 무엇인가?

주요 결과

• 시간에 동일한 분열 과정에서 재해가 있는 경우, 생존 확률은 지수적으로 감소하며, 감쇠 속도는 분열, 재해, 생존 파라미터 간의 균형에 따라 계단적 전이를 보인다.
• 적분 $ \int_0^t \log(1/p_s) \kappa_s \, ds $ 가 분열 강도보다 더 빠르게 증가하면 멸종은 거의 확실하며, 그렇지 않으면 생존 확률이 양수로 유지된다.
• 시간에 비례하지 않는 경우, 재해 시간 조건 하에서 생존 확률은 이중 p-점프 과정에 의해 결정되는 극한으로 수렴하며, 이는 $ \log \mathbb{P}(Z_t > 0 \mid \mathcal{D}_\infty) \sim \max\{ -L_t, \log \int_0^t e^{-L_s} b_s \, ds \} / h(t) $ 로 표현 가능한 수렴 속도를 가진다.
• 누적 드리프트 $ \ell(t) = \int_0^t (b_s - d_s - \kappa_s \log(1/p_s)) \, ds $ 가 상수로 수렴하면, 과정은 거의 확실하게 유한한 수의 출생을 가지며 비퇴화 극한 분포로 수렴한다.
• 연속 상태 분열 과정에 재해가 있는 경우, 멸종 확률은 동일한 이중성으로 특성화되며, 극한 행동은 이중 과정의 점근적 성장에 따라 달라진다.
• 적분 $ \int_0^t e^{-L_s} b_s \, ds $ 가 발산하면 생존 확률은 지수 감소보다 느리게 감소하며, 적분이 수렴하면 더 빠른 감쇠 또는 확률 1의 멸종이 발생한다.