[论文解读] Exponential convergence to a steady-state for a population genetics model with sexual reproduction and selection

Gaël Raoul|arXiv (Cornell University)|Apr 13, 2021

Mathematical Biology Tumor Growth参考文献 31被引用 8

一句话总结

该论文在具有弱选择的有性生殖种群遗传模型中，利用Wasserstein距离不等式和尾部估计，建立了稳态的指数收敛性。当选择强度α较小时，证明了解以指数速度收敛至局部唯一的稳态，利用了繁殖算子的压缩性质及选择扰动下的稳定性分析。

ABSTRACT

We are interested in the dynamics of a population structured by a phenotypic trait. Individuals reproduce sexually, which is represented by a non-linear integral operator. This operator is combined to a multiplicative operator representing selection. When the strength of selection is small, we show that the dynamics of the population is governed by a simple macroscopic differential equation, and that solutions converge exponentially to steady-states that are locally unique. The analysis is based on Wasserstein distance inequalities using a uniform lower bound on distributions. These inequalities are coupled to tail estimates to show the stability of the steady-states.

研究动机与目标

• 分析在有性生殖和选择作用下，按表型特征分层的种群的长期动力学行为。
• 在选择强度α较小时，建立解对稳态的指数收敛性。
• 使用Wasserstein距离而非基于正则性的方法，对稳态进行稳定性分析。
• 平衡选择的破坏性效应与有性生殖算子的压缩效应。
• 通过一致下界和尾部估计，证明在小α下稳态的局部唯一性及对稳态的指数收敛性。

提出的方法

• 使用非线性积分算子建模有性生殖，其中后代特征从以亲代特征均值为中心的正态分布中抽取。
• 将繁殖算子与乘法选择项(1 + αa(x))结合，以建模适应度差异。
• 应用Wasserstein距离W₂来度量随时间推移的特征分布之间概率测度的距离。
• 利用Raoul（2017）推导出的繁殖算子的Wasserstein压缩不等式，控制解的演化。
• 通过分布的统一下界和尾部估计，控制选择项的破坏性影响。
• 通过质心和选择加权矩等宏观量分析动力学，利用函数不等式估计其差异。

实验结果

研究问题

• RQ1 在具有弱选择的有性生殖种群中，种群分布如何随时间演化？
• RQ2 在何种条件下系统会指数收敛至稳态？
• RQ3 选择的破坏性效应能否被有性生殖算子的压缩效应所抵消？
• RQ4 Wasserstein距离在此类非局部、非线性模型中在量化对稳态收敛性方面起什么作用？
• RQ5 稳态的局部唯一性与选择强度α之间有何关系？

主要发现

• 当α > 0较小且F′(Z̄) < 0时，动能方程的解以指数速度收敛至稳态，收敛速率由F′(Z̄)α决定。
• 稳态是局部唯一的，由引理2.2确立，并通过压缩性论证得到验证。
• 通过Lyapunov型论证证明指数收敛性，其中量√α X(t) + |Zₙ∘φₙ(t) − Zₘ∘φₘ(t)|在时间t₀ = −C ln α / α后指数衰减。
• 两个解之间的Wasserstein距离满足W₂(𝒩(t,·), 𝒪(t,·)) ≤ C exp((F′(Z̄) + 𝒪(√α))αt)，表明对稳态的指数衰减。
• 尾部估计和分布的统一下界对于控制选择项的破坏性影响至关重要。
• 该方法避免依赖正则性估计（如U = −σ log n的导数），转而使用基于输运的度量，为这类模型提供了新的分析框架。