[논문 리뷰] On the variations of the principal eigenvalue with respect to a parameter in growth-fragmentation models

Fabien Campillo, Nicolas Champagnat외 1명|arXiv (Cornell University)|2016. 01. 11.

Evolution and Genetic Dynamics참고 문헌 24인용 수 10

한 줄 요약

이 논문은 성장-분열 모델에서 매개변수 $ S $ 에 대한 주 특이값의 단조성을 비확률적 개별 기반 모델을 사용하여 확립한다. 알려진 이원성에 의해 스토케스틱 과정의 생존 확률을 결정론적 특이값과 연결함으로써, 성장 및 분열 속도가 $ S $ 에 대해 단조일 경우 주 특이값 $ \Lambda_S $ 는 $ S $ 에 대해 비감소임을 증명한다.

ABSTRACT

We study the variations of the principal eigenvalue associated to a growth-fragmentation-death equation with respect to a parameter acting on growth and fragmentation. To this aim, we use the probabilistic individual-based interpretation of the model. We study the variations of the survival probability of the stochastic model, using a generation by generation approach. Then, making use of the link between the survival probability and the principal eigenvalue established in a previous work, we deduce the variations of the eigenvalue with respect to the parameter of the model.

연구 동기 및 목표

• 주 특이값 $ \Lambda_S $ 가 구조화된 인구의 점근적 성장률을 나타내며, 환경 매개변수 $ S $ 에 따라 어떻게 변화하는지 분석하는 것.
• 특히 비단조적이거나 복잡한 경우에 매개변수 변화에 따른 성장-분열 방정식에서 특이값 변화를 연구하기 위한 분석 도구의 부족을 해결하는 것.
• 생존 확률과 주 특이값 간의 연결을 활용하여 스토케스틱 개별 기반 모델과 결정론적 적분미분방정식 간 격차를 메우는 것.
• 스토케스틱 모델에서 생존 확률의 단조성을 바탕으로 특이값의 단조성을 유도하는 일반적인 방법을 제공하는 것.
• 연속적인 크기 구조와 시간 동질적 매개변수를 갖는 모델에 대해 적응 역학과 침입 적합도 개념의 적용 범위를 확장하는 것.

제안 방법

• 성장-분열-사망 방정식을 조각별 결정론적 마코프 과정(PDMP)으로 해석하는 확률적, 개별 기반 해석을 사용한다.
• 다른 세대 간의 겹침이 있음에도 불구하고, 생존 확률을 계산하기 위해 세대별 분석을 적용함으로써 단조성 분석이 가능해진다.
• 성장 및 분열 속도가 $ S $ 에 대해 단조일 것으로 가정하고, 서로 다른 매개변수 값 $ S_1 $ 과 $ S_2 $ 에서의 생존 확률을 비교하기 위해 쌍체 기법을 적용한다.
• 스토케스틱 모델의 생존 확률과 결정론적 모델의 주 특이값 $ \Lambda_S $ 를 연결하는 알려진 이원성 결과(정리 2.4)를 활용한다.
• 특이값 문제를 변형하고, 이완된 사망률 $ D' = D + \Lambda_{S^*} $ 를 통해 생존 확률 비교를 이용하여 특이값의 단조성을 도출한다.
• 특정 케이스, 예를 들어 성장 속도가 $ g(S,x) = \mu(S)\tilde{g}(x) $ 인 경우에 적용하여 $ \mu(S) $ 에 따라 $ \Lambda_S $ 가 증가함을 보인다.

실험 결과

연구 질문

• RQ1 성장 및 분열 속도를 조절하는 매개변수 $ S $ 에 따라 성장-분열 모델의 주 특이값 $ \Lambda_S $ 는 어떻게 변화하는가?
• RQ2 스토케스틱 개별 기반 모델에서 생존 확률의 단조성이 결정론적 주 특이값의 단조성을 유추하는 데 사용될 수 있는가?
• RQ3 성장 및 분열 속도에 대해 어떤 조건이면 특이값 $ \Lambda_S $ 가 $ S $ 에 따라 증가하는가?
• RQ4 스토케스틱 침입 적합도와 결정론적 특이값 간의 이원성을 어떻게 활용하여 특이값 변화에 대한 분석적 결과를 도출할 수 있는가?
• RQ5 이 방법은 적응 역학에 대해 어떤 함의를 가지며, 특히 환경 변화에 대한 돌연변이 침입 성공 예측에 어떤 영향을 미치는가?

주요 결과

• 성장 및 분열 속도가 $ S $ 에 대해 비감소일 경우 주 특이값 $ \Lambda_S $ 는 $ S $ 에 대해 비감소함을 확인하였으며, 이는 추론 3.13에 의해 입증된다.
• 만약 $ \Lambda_{S_1} > 0 $ 이면 모든 $ S_2 > S_1 $ 에 대해 $ \Lambda_{S_2} > 0 $ 이고, $ \Lambda_{S_1} \leq 0 $ 이면 모든 $ S_2 < S_1 $ 에 대해 $ \Lambda_{S_2} \leq 0 $ 이다. 이는 추론 3.12에 명시되어 있다.
• 성장 속도가 $ g(S,x) = \mu(S)\tilde{g}(x) $ 인 경우, 생존 확률과 특이값은 $ \mu(S) $ 에 대해 비감소하므로 $ \mu(S^1) \leq \mu(S^2) \Rightarrow \Lambda_{S^1} \leq \Lambda_{S^2} $ 가 성립한다. 이는 추론 3.14에 따라 성립한다.
• 생존 확률 비교 $ \mathbb{P}^{S_1,D} \geq \mathbb{P}^{S_2,D} $ 는 $ \Lambda^{D}_{S_1} \geq \Lambda^{D}_{S_2} $ 를 의미하며, 이는 생존과 특이값 단조성 간의 일반적인 이원성을 확립한다.
• 분열 핵함수 $ q $ 가 $ S $ 에 따라 달라지더라도, 생존 확률 비교가 성립한다면 이 방법은 여전히 유효하다.
• 겹침이 있는 세대 문제의 어려움을 피하기 위해 스토케스틱 모델의 세대 구조와 쌍체 기법을 활용함으로써, 특이값 문제에 대한 접근이 가능해진다.