[논문 리뷰] A note on Exact solution of SIR and SIS epidemic models

Ghulam Shabbir, Haroon Ahmed Khan외 1명|arXiv (Cornell University)|2010. 12. 22.

Mathematical and Theoretical Epidemiology and Ecology Models참고 문헌 12인용 수 26

한 줄 요약

이 논문은 일정한 인구 수를 가진 특정 경우의 SIR 및 SIS 전염병 모델에 대해 직접 통합 기법을 사용하여 정확한 해석적 해를 제시한다. 저자들은 제약 조건이 없는 파라미터 조건 하에서 두 모델에 대해 닫힌 형태의 해를 유도하며, 일반적인 초기 조건 하에서 전염병 역학의 정확한 분석을 가능하게 함으로써 수학적 전염병학 분야에서 중요한 발전을 이룬다.

ABSTRACT

In this article we have successfully obtained an exact solution of a particular case of SIR and SIS epidemic models given by Kermack and Mckendrick [1] for constant population, which are described by coupled nonlinear differential equations. Our result has no limiting conditions for any parameter involved in the given models. In epidemiology many researchers believe that it is very hard to get an exact solution for such models. We hope this solution will be an opening window and good addition in the area of epidemiology.

연구 동기 및 목표

• 일정한 인구를 가진 SIR 및 SIS 전염병 모델에 대해 정확한 해석적 해를 도출함으로써, 오랫동안 남아있던 수학적 전염병학의 과제를 해결한다.
• 질병 전파를 기술하는 비선형 연립 미분방정식의 일반적인 비통합성 문제를 극복한다.
• 모델 파라미터에 대한 제한 조건이 없는 해 프레임워크를 제공함으로써 모델의 적용 가능성과 해석 가능성을 향상시킨다.
• 정확한 함수 형태를 통해 질병 역학의 더 깊은 정성적 분석을 위한 기반을 제공한다.

제안 방법

• SIS 모델의 경우, 보존법칙 s(t) + i(t) = k를 사용하여 시스템을 i(t)에 대한 단일 비선형 상미분방정식으로 감소시키며, 이후 y = i⁻¹의 치환을 통해 해결한다.
• SIS 해는 감소된 미분방정식에 y = i⁻¹을 대입한 후 직접 통합을 통해 유도되며, 통합 상수 C를 포함한 지수 형태로 표현된다.
• SIR 모델의 경우, 방정식을 합산함으로써 보존법칙 s(t) + i(t) = 1 + C e⁻ᵘᵗ를 도출하며, 이를 통해 시스템을 i(t)에 대한 단일 미분방정식으로 감소시킨다.
• SIR 해는 z = i⁻¹의 치환을 적용하고, e⁻ᵘᵗ에 대한 급수 전개를 사용하여 통합을 단순화한다. 고차항은 무시한다.
• 두 모델 모두 초기 조건을 활용하여 통합 상수를 구하고, s(t) 및 i(t)에 대한 닫힌 형태의 표현식을 도출한다.
• 이 방법은 수치적 또는 섭동 접근 방식이 아닌 직접 통합과 치환 기법에 의존하므로 정확성을 보장한다.

실험 결과

연구 질문

• RQ1 제약 조건이 없는 파라미터 조건 하에서 일정한 인구를 가진 SIR 및 SIS 모델에 대해 정확한 해석적 해를 도출할 수 있는가?
• RQ2 이러한 고전적 전염병 모델의 비선형 연립 미분방정식은 어떻게 정확한 해석적 통합을 통해 정확하게 해결할 수 있는가?
• RQ3 보존법칙 s(t) + i(t) = k (또는 s(t) + i(t) = 1 + C e⁻ᵘᵗ)는 정확한 해를 가능하게 하는 데 어떤 역할을 하는가?
• RQ4 정확한 해가 전염병 역학의 이해 및 모델 검증에 있어 어느 정도 향상시킬 수 있는가?

주요 결과

• SIS 모델의 정확한 해는 i(t) = β / (r + βC e⁻ᵝᵗ)로 표현되며, 여기서 β = rk - α이며, C는 초기 조건 i(0) = i₀로 결정된다.
• SIS 해는 s(t) + i(t) = k를 정확히 만족함으로써 총 인구의 보존을 확인한다.
• SIR 모델의 경우, 해는 i(t) = λ / (β + λD e⁻ᵘᵗ e^(βC/μ))로 표현되며, 여기서 λ = β - μ + βC 이고, C = s₀ + i₀ - 1이다.
• SIR 해는 초기 조건 s(0) = s₀ 및 i(0) = i₀를 활용하여 통합 상수 D 및 C를 완전히 결정한다.
• SIS의 경우 근사치를 사용하지 않지만, SIR의 경우 단순화를 위해 급수 전개(e⁻ᵘᵗ ≈ 1 - μt)를 사용하며 고차항은 무시한다.
• 유도된 해는 r, α, β, μ 및 초기 조건 s₀, i₀ > 0의 모든 양수 값에 대해 유효하며, 제한 조건 없이 적용 가능하다.