[논문 리뷰] A simple computational approach to the Susceptible-Infected-Recovered (SIR) epidemic model via the Laplace-Adomian Decomposition Method

Tiberiu Harko, M. K. Mak|arXiv (Cornell University)|2020. 06. 12.

Fractional Differential Equations Solutions참고 문헌 25인용 수 4

한 줄 요약

이 논문은 SIR 전염병 모델에 대해 단순하고 높은 정확도를 가진 급수 해를 도출하기 위해 라플라스-아도미안 분해 방법(Laplace-Adomian Decomposition Method, LADM)을 제안한다. 이는 효율적인 매개변수 피팅과 역학적 분석을 가능하게 하며, 오직 세 항으로만 구성된 지수 급수 형태의 반해석적 해를 제공한다. 이는 정확한 매개변수 해에 비해 계산 비용이 훨씬 낮은 대안을 제공한다.

ABSTRACT

The Susceptible-Infected-Recovered (SIR) epidemic model is extensively used for the study of the spread of infectious diseases. Even that the exact solution of the model can be obtained in an exact parametric form, in order to perform the comparison with the epidemiological data a simple but highly accurate representation of the time evolution of the SIR compartments would be very useful. In the present paper we obtain a series representation of the solution of the SIR model by using the Laplace-Adomian Decomposition Method to solve the basic evolution equation of the model. The solutions are expressed in the form of infinite series. The series representations of the time evolution of the SIR compartments are compared with the exact numerical solutions of the model. We find that there is a good agreement between the Laplace-Adomian semianalytical solutions containing only three terms, and the numerical results.

연구 동기 및 목표

• 실제 역학적 데이터와의 비교를 더 쉽게 하기 위해 SIR 모델의 시간 진화를 단순하면서도 정밀한 해석적 표현으로 제공하는 것.
• 정확한 매개변수 해를 사용할 경우 계산 부담이 크므로, 급수 기반의 반해석적 대안을 제공함으로써 이를 극복하는 것.
• 간편한 함수 형태를 통해 전파($\beta$) 및 회복($\gamma$) 속도의 효율적 추정을 가능하게 하기 위한 데이터 피팅을 가능하게 하는 것.
• 단 세 항으로 잘라낸 급수 해가 다양한 매개변수 영역에서 SIR 모델의 역학을 정확하게 포착할 수 있는지 보여주는 것.

제안 방법

• SIR 모델의 회복자 집단 $z(t)$를 기술하는 기본 비선형 상미분방정식에 라플라스-아도미안 분해 방법(LADM)을 적용한다.
• 해는 무한급수 $z(t) = \sum_{i=0}^{\infty} z_i(t)$ 형태로 표현되며, 아도미안 다항식을 사용해 순차적으로 계산된다.
• 비선형 미분방정식을 선형화하기 위해 라플라스 변환을 사용하여 급수의 각 항을 순차적으로 구한다.
• 최종적으로 도출된 급수 해는 세 항으로 잘라내어 $x(t)$, $y(t)$, $z(t)$에 대한 실용적인 반해석적 근사치를 생성한다.
• 이 방법은 지수 함수 형태의 해를 제공하여 실질적 데이터에 직접 피팅이 가능하다.
• 여러 매개변수 조합에 대해 잘라낸 세 항 급수 해와 정확한 수치 해를 비교하여 방법의 타당성을 검증한다.

실험 결과

연구 질문

• RQ1 LADM를 통해 유도된 잘라낸 급수 해가 최소한의 계산 비용으로 SIR 모델의 각 집단의 시간 진화를 정확하게 표현할 수 있는가?
• RQ2 세 항으로 구성된 LADM 급수 근사치가 $\beta$ 및 $\gamma$ 값이 다양한 경우 정확한 수치 해와 얼마나 잘 일치하는가?
• RQ3 반해석적 급수 해가 실제 역학적 데이터에 대한 $\beta$ 및 $\gamma$ 피팅을 얼마나 단순화할 수 있는가?
• RQ4 LADM 기반의 급수 해가 피크 감염 및 회복 경향과 같은 SIR 모델의 정성적 역학을 유지하는가?

주요 결과

• 세 항으로 잘라낸 라플라스-아도미안 급수 해는 테스트된 모든 매개변수 조합에서 SIR 모델 역학을 매우 정확하게 근사한다. 정확한 수치 해와 거의 일치한다.
• 오직 세 항으로도 정확한 수치 결과와 양호한 일치를 보이며, 실용적 응용에 있어 강력한 수렴성과 신뢰성을 입증한다.
• 반해석적 해는 지수 함수의 합으로 표현되며, 이는 실질적 역학적 데이터에 직접 피팅하기에 매우 용이하다.
• 높은 전파 속도나 빠른 회복 속도를 가진 경우를 포함하여 $\beta$ 및 $\gamma$의 다양한 값에서도 정확성이 유지된다.
• 두 항 근사치 역시 합리적인 일치를 보이며, 이는 추가 계산 비용을 최소한으로 증가시켜도 정확도가 향상됨을 시사한다.
• 정확한 매개변수 해는 각 단계에서 시간 적분을 수치적으로 평가해야 하므로 계산 비용이 크지만, 본 방법은 이를 훨씬 효율적인 대안으로 제공한다.