[논문 리뷰] Reactive Social distancing in a SIR model of epidemics such as COVID-19

Rose Baker|arXiv (Cornell University)|2020. 03. 18.

COVID-19 epidemiological studies참고 문헌 1인용 수 12

한 줄 요약

이 논문은 SIR 프레임워크 내에서 감염 비율이 감염 비율이 증가함에 따라 감소하는 반응형 사회적 거리두기 모델을 제안한다. 이는 힐 함수를 통해 모델링되며, 최종 유행 규모와 정점 감염 비율에 대한 해석적 해를 유도한다. 극단적인 반응형 거리두기가 감염 수를 최소화하면서 유행 기간을 연장함으로써, 코로나19와 같은 유행 동안 건강과 경제적 영향을 균형 잡는 데 도움이 되는 정책 수립자들을 위한 도구를 제공한다.

ABSTRACT

A model of reactive social distancing in epidemics is proposed, in which the infection rate changes with the number infected. The final-size equation for the total number that the epidemic will infect can be derived analytically, as can the peak infection proportion. This model could assist planners during for example the COVID-19 pandemic.

연구 동기 및 목표

• 상승하는 감염 수준에 대응해 자율적 또는 정부가 명령하는 거리두기가 전파율을 어떻게 낮출 수 있는지 모델링하기.
• 이러한 반응형 행동 하에서 최종 유행 규모와 정점 감염 비율에 대한 해석적 표현을 도출하기.
• 유행 발생 시 정책 결정을 지원하기 위해 유행 결과를 추정하고 지침을 내릴 수 있는 실용적인 모델을 제공하기.
• 고정된 접촉률 감소 모델과 달리, 거리두기 조치 해제 후 유행 재발을 방지하는 반응형 거리두기의 효과를 보여주기.

제안 방법

• 상수 전파율 β를 β(t) = β / (1 + γi(t))로 대체하여 접촉률을 감염 비율에 반비례하도록 설정.
• 피드백 메커니즘을 모델링하기 위해 힐 유형 함수를 사용하여 해석적 해법이 가능하면서도 현실적인 행동 반응을 반영.
• 시간을 제거하고 감염 동역학을 저항성 s에 따라 기술함으로써 i(s)에 대한 미분방정식 유도.
• 표준 방법(보조 함수 및 특수 적분)을 사용해 유도된 ODE를 해석하고, 경계 조건을 s=1, i=0로 설정.
• 최종 크기 방정식을 도출하기 위해 i(∞) = 0로 설정하여 s∞에 대한 비선형 방정식을 유도하고, 뉴턴-라프슨 방법을 통해 수치적으로 해를 구함.
• di/ds = 0으로 설정하여 정점 감염 비율을 찾고, 이를 통해 비선형 시스템을 해결하여 s_max 및 i_max에 대한 닫힌 형식의 표현을 도출.

실험 결과

연구 질문

• RQ1 감염 비율이 상승함에 따라 접촉률이 감소하는 반응형 거리두기가 SIR 모델에서 최종 유행 규모에 어떤 영향을 미치는가?
• RQ2 이러한 피드백 메커니즘 하에서 정점 감염 비율과 최종 규모에 대한 해석적 표현을 도출할 수 있는가?
• RQ3 사회적 거리두기가 극도로 반응적일 경우(γ → ∞), 유행의 점차적 행동은 어떻게 되는가?
• RQ4 고정된 감소 모델과 비교해 볼 때, 거리두기 조치가 해제된 후 유행 재발을 방지하는 데 모델이 얼마나 효과적인가?
• RQ5 감염 수를 최소화하면서도 경제적 및 의료 체계의 부담을 고려할 때 최적의 반응형 거리두기 수준은 무엇인가?

주요 결과

• 최종 크기 방정식은 s∞^(γ/R₀) = [1/γ + s∞/(1 - γ/R₀)] / [1/γ + 1/(1 - γ/R₀)] 로 유도되며, s∞에 대해 뉴턴-라프슨 방법을 통해 수치적으로 해를 구할 수 있다.
• γ → ∞ 일 때, 감염되지 않은 비율의 최종 비율은 1/R₀에 수렴하며, 이는 집단면역의 이론적 최소값이다.
• 정점 감염 비율 i_max 는 i_max = (R₀/γ) × [ (1 - γ/R₀ + γ)/R₀ ]^(1/(1 - γ/R₀)) - 1/γ 으로 주어진다.
• γ = R₀ 일 때, 해는 i(s) = (s - 1)/R₀ - s ln s 로 단순화되며, 일반 해의 극한과 일치한다.
• γ → 0 일 때, s_max → 1/R₀ 이고 i_max → 0 이며, 표준 SIR 모델의 행동을 회복한다.
• 정점 감염 시점에서 효과적 복제수 R₀ / (1 + γ i_max) → 1 이 되며, γ → ∞ 일 때 유행이 정점에서 겨우 유지됨을 나타낸다.