[論文レビュー] Fractional Order Malaria Model With Temporary Immunity

Eric Okyere, Francis T. Oduro+3 more|arXiv (Cornell University)|Mar 21, 2016

Fractional Differential Equations Solutions参考文献 36被引用数 7

ひとこと要約

本稿では、一時的免疫を組み込んだ分数階微分方程式を用いたマラリア伝搬モデルを提案する。Caputo微分を用いて病気の動的変化における記憶効果を捉える。次世代行列法を用いて基本再生産数 $ R_0 $ を計算し、アダムス型予測・修正スキームによる数値シミュレーションにより、分数階 $ \alpha $ が小さくなるにつれて系に顕著な周期的挙動が現れることが判明した。これは、記憶効果が整数階モデルと比較して伝搬ダイナミクスに顕著な影響を及ぼしていることを示している。

ABSTRACT

In this paper, we propose and analyze a new fractional order malaria model with temporary immunity. We extend the existing malaria model to include fractional derivatives. Basic reproduction number is computed using the next generation matrix method. Model equilibria are determined and their local asymptotic stability analysis are considered. An efficient Adams-type predictor-corrector method is used to solve the initial value problem. Finally, phase plane portraits are generated to describe the dynamics between susceptible fractions of human hosts and infected fractions of mosquito vectors.

研究の動機と目的

• 既存の整数階マラリアモデルを拡張し、記憶効果を捉えるために分数階微分を導入する。
• 分数階フレームワークにおける無病均衡およびエンドミック均衡の局所的漸近的安定性を分析する。
• 分数階微積分の文脈において、次世代行列法を用いて基本再生産数 $ R_0 $ を計算する。
• 分数階微分方程式のためのアダムス型予測・修正スキームを用いて、モデルのダイナミクスを数値的にシミュレートする。
• 分数階 $ \alpha $ の変化が、系の周期的挙動および収束速度に与える影響を調査する。

提案手法

• モデルは、$ \alpha \in (0,1] $ のCaputo分数階微分を用いて、人間およびモスキートの個体群動態を記述する。
• モデルには、感受性、感染、回復した人間、および感受性・感染モスキートの状態を含む区分が含まれており、一時的免疫は回復と再感染のダイナミクスによってモデル化される。
• 基本再生産数 $ R_0 $ は次世代行列法を用いて導出され、病気閾値行動の解析的取り扱いを保証する。
• 初期値問題を解くために、アダムス型予測・修正スキームが実装され、係数 $ a_{j,k+1} $ および $ b_{j,k+1} $ は分数階二項係数によって定義される。
• 異なる $ \alpha $ 値における感受性人間と感染モスキートの相互作用を可視化するため、位相平面図が生成される。
• パラメータ値は Tumwiine et al. (2007a) から適応され、エンドミック状態を表すためにシミュレーションで $ R_0 = 1.5 $ が用いられる。

実験結果

リサーチクエスチョン

• RQ1 分数階微分の導入が、一時的免疫を有するマラリア伝搬モデルの安定性およびダイナミクスに与える影響は何か？
• RQ2 分数階マラリアモデルにおける基本再生産数 $ R_0 $ の値は何か？整数階の場合と比較するとどうなるか？
• RQ3 分数階 $ \alpha $ の異なる値が、系の周期的挙動および収束速度に与える影響は何か？
• RQ4 分数微分に埋め込まれた記憶効果が、古典的モデルと比較して予測される病気の経路にどの程度影響を及えるか？
• RQ5 感受性人間と感染モスキートの位相平面図が、$ \alpha $ を変化させた際にどのように変化するか？

主な発見

• 次世代行列法を用いて基本再生産数 $ R_0 = 1.5 $ が計算され、モデルにおける持続的伝搬可能性を示している。
• $ \alpha $ が 1 から 0.90 に減少するにつれて、解に顕著な振動が現れるようになり、記憶効果が病気の有病率における周期的挙動を強化していることが示唆された。
• $ \alpha = 1 $ のとき、分数階モデルは古典的整数階モデルに還元され、既存のダイナミクスと整合していることが確認された。
• アダムス型予測・修正スキームは、さまざまな $ \alpha $ における複雑な過渡ダイナミクス（振動や収束パターン）を的確に捉えた。
• 位相平面図では、感受性人間と感染モスキートの軌道パターンが $ \alpha $ に応じて変化しており、伝搬効率および持続性が変化していることが示された。
• 本モデルは、特にマラリア伝搬における長期記憶効果を捉える点で、整数階モデルよりも洗練されたダイナミカル構造を提供することを示した。