[論文レビュー] A study of gravity-linked metapopulation models for the spatial spread of dengue fever

Marta Sarzynska, Oyita Udiani+1 more|arXiv (Cornell University)|Aug 20, 2013

COVID-19 epidemiological studies参考文献 19被引用数 3

ひとこと要約

本研究では、人口規模と距離を用いて局所間の伝播重みを定義することで、ペルーにおける空間的デングア熱拡散を模擬する重力連携メタポピュレーションモデルを開発した。重力ベースのリンクは、均一なリンク重みと比較して、地域間でより現実的で非同期な流行ピークを生み出すことが示された。一方で、不規則な発生数と限られた気象データのため、実データへの適合は困難であったが、年間の流行パターンをうまく捉えることができ、さらに精度を高めるには多系統型および確率的モデリングの統合が不可欠であることが明らかになった。

ABSTRACT

Metapopulation (multipatch) models are widely used to study the patterns of spatial spread of epidemics. In this paper we study the impact of inter-patch connection weights on the predictions of these models. We contrast arbitrary, uniform link weights with link weights predicted using a gravity model based on patch populations and distance. In a synthetic system with one large driver city and many small follower cities, we show that under uniform link weights, epidemics in the follower regions are perfectly synchronized. In contrast, gravity-based links allow a more realistic, less synchronized distribution of epidemic peaks in the follower regions. We then fit a three-patch metapopulation model to regional dengue fever data from Peru -- a country experiencing yearly, spatially defined epidemics. We use data for 2002-2008 (studying the seasonal disease patterns in the country and the yearly reinfection patterns from jungle to the coast) and 2000-2001 (one large epidemic of a new disease strain across the country). We present numerical results.

研究の動機と目的

• メタポピュレーションモデルにおける局所間接続重みが空間的疾患拡散予測に与える影響を調査すること。
• 一中枢都市と99の周辺都市からなる合成ネットワークにおいて、均一なリンク重みと重力ベースのリンク重みを比較し、デングア熱流行を模擬すること。
• ペルーの実際のデングア熱発生データ（2000–2008年）に、決定論的で重力ベースの3ピッチメタポピュレーションモデルを適合させ、季節的および系統特異的伝播ダイナミクスに焦点を当てること。
• データの不規則さや適合不良によるモデルの限界を評価し、気象データの統合、確率的要因、複数のデングアウイルス系統の導入といった改善方向を同定すること。
• 人間移動に基づく伝播予測のための代替モデル（例：放射モデル）の可能性を検討すること。

提案手法

• 本研究では、異なる地理的領域を表す複数のピッチにわたり、決定論的SIR型の区分モデルを用いる。
• 局所間の伝播重みは重力モデルで定式化される：$ P_{ij} = \theta \frac{n_i^\alpha n_j^\beta}{d_{ij}^\gamma} $、ここで$ n_i $および$ n_j $は人口、$ d_{ij} $は距離、$ \alpha, \beta, \gamma, \theta $はフィッティングされたパラメータである。
• モデル適合には最小二乗法を用い、モデル予測と観察されたデングア熱発生データとの差の二乗和を最小化する。
• まず、1つの主要都市と99の小規模都市からなる合成系を用いて、均一なリンク重みと重力ベースの接続性を比較する。
• 次に、ペルーの実際のデングア熱データに、当初7ピッチで適用したが、データ品質と不規則性のため3ピッチに縮小した。
• さらに、2000–2001年の流行中にデングア熱発生が確認された全49の州を表す49ピッチモデルでもテストを行い、予測の形状と大きさを実データと比較した。

実験結果

リサーチクエスチョン

• RQ1 均一な局所間伝播重みと重力ベースの重みの違いが、メタポピュレーションモデルにおける流行ピークの同期性にどのように影響するか？
• RQ2 重力ベースのメタポピュレーションモデルは、ペルーにおける観察されたデングア熱の空間的・時間的パターンを正確に再現できるか？
• RQ3 モデルは一般に発生ピークの形状を捉えているにもかかわらず、なぜ流行の規模を過大評価してしまうのか？
• RQ4 発生数の不規則さや気象データの欠落といったデータ制限が、モデルの適合と予測精度にどのように影響するか？
• RQ5 多系統型ダイナミクス、確率的要因、または代替移動モデルの導入といった、モデル性能を向上させるための改善策は何か？

主な発見

• 合成モデルでは、均一なリンク重みはすべてのフォロワーカンプで完全な同期を生じたが、重力ベースのリンクはより現実的で非同期な拡散を生み出した。
• 3ピッチの重力ベースモデルは、ペルーにおける年間デングア熱発生ピークの一般形状を捉えたが、観察値より高い流行規模を予測した。
• 49ピッチモデルでも良好な適合は得られず、最小二乗法の目的関数の性質により、形状はやや良いが振幅が低い曲線が最良の適合とされた。
• 49ピッチモデルの初期適合では、$ \beta $および$ \gamma $の値が高めに推定され、3ピッチスケールでは重力モデルの効果を十分に捉えられていない可能性が示唆された。
• モデルが複数のデングアウイルス系統を考慮できないこと（各系統が別々の年間流行を引き起こす）が、正確性を制限しており、1つの系統に対する免疫は他の系統に対しては保護効果を持たない。
• 本研究は、気象データを伝播パラメータに統合し、確率的要因を組み込み、放射モデルのような代替移動モデルの検討を進めることが、より良いモデルの構築に不可欠であると結論づけた。