[논문 리뷰] Covid-19 spread: Reproduction of data and prediction using a SIR model on Euclidean network

Kathakali Biswas, Abdul Khaleque외 1명|arXiv (Cornell University)|2020. 03. 16.

COVID-19 epidemiological studies참고 문헌 21인용 수 69

한 줄 요약

해당 논문은 SIR 모델을 이용해 COVID-19 중국 데이터를 적합시켜 시간 및 거리 패턴을 재현하고, 지역적 유행 종료에 대한 결론을 예측한다.

ABSTRACT

We study the datafor the cumulative as well as daily number of cases in the Covid-19 outbreak in China. The cumulative data can be fit to an empirical form obtained from a Susceptible-Infected-Removed (SIR) model studied on an Euclidean network previously. Plotting the number of cases against the distance from the epicenter for both China and Italy, we find an approximate power law variation with an exponent $\sim 1.85$ showing strongly that the spatial dependence plays a key role, a factor included in the model. We report here that the SIR model on the Eucledean network can reproduce with a high accuracy the data for China for given parameter values, and can also predict when the epidemic, at least locally, can be expected to be over.

연구 동기 및 목표

• 유클리드 네트워크 위의 SIR 모델을 사용하여 중국 내 COVID-19의 공간적 및 시간적 확산을 이해한다.
• 유클리드 네트워크의 SIR에서 도출된 경험적 형태로 누적 확진 데이터를 맞춘다.
• 진원지로부터의 거리와 확진자 수의 관계를 탐색하여 공간 의존성을 평가한다.

제안 방법

• 거리 기반의 장거리 연결(P(ℓ) ∝ ℓ^{−δ})과 감염 확률 q를 갖는 유클리드 네트워크 위의 SIR 모델을 사용한다.
• 누적 감염 데이터 R(t)를 R(t) = a exp(t/T) / [1 + c exp(t/T)]에 맞춘다. 여기서 a, c, T는 적합 파라미터이다.
• R(t)를 미분하여 I(t)를 계산하고 피크 시간 t_p = T ln(1/c)을 찾는다.
• 가까운 이웃 링크와 임의의 긴 거리 링크 한 개를 포함한 2^12-노드 네트워크를 시뮬레이션하고, 100 구성을 평균하여 비교한다.
• 모델 출력을 실제 데이터의 포화에 맞게 재스케일하고 몬테카를로 스텝을 일수로 변환하여 비교한다.

실험 결과

연구 질문

• RQ1 유클리드 네트워크 위의 SIR 모델이 중국의 COVID-19 누적 확진 사례의 시간 전개를 재현할 수 있는가?
• RQ2 모델이 진원지로부터의 거리 증가에 따른 사례의 공간적 감소를 포착하는가, 그리고 이것이 접촉 범위 δ에 대해 무엇을 시사하는가?
• RQ3 데이터를 재현하는 매개변수 값(δ, q)은 무엇이며, 이것이 전파 역학과 유행 기간에 대해 무엇을 시사하는가?
• RQ4 모형에 따른 전염병 종료의 예상 기간은 어떻게 되는가?

주요 결과

• 중국의 누적 데이터는 R(t) = a exp(t/T) / [1 + c exp(t/T)]로 적합되며, a = 1.15698×10^−6, c = 0.0201529, T = 5.22 그리고 비교적 작은 매개변수 오차를 가진다.
• 거리 증가에 따른 감염율 감소는 멱 법칙 R(d) ∝ d^−γ를 따르며, χ γ = 1.85 ± 0.1은 중국(우한) 및 이탈리아(베르가모) 진원지에서의 값이다.
• δ ≈ 1.8 및 q ≈ 0.85인 유클리드 네트워크 위의 SIR 모델은 적절한 재스케일링(노출 분수 ρ ≈ 6.0783×10^−5) 후에 중국 데이터를 재현한다.
• 모형은 중국에서 2020년 1월 21일로부터 약 60–78일 사이에 유행이 종료될 수 있으며, 약 82 몬테카를로 스텝(약 78일)에서 감염자 제로를 예측한다.
• 더 큰 δ는 더 짧은 효과적 접촉 범위를 시사하고, 비교적 높은 q는 모델링 맥락에서 더 강한 전파를 반영하며, 결과는 이동성 없는 차단 해석과 일치한다.
• δ와 q의 작은 변화(±0.05)에 대해 적합도가 여전히 견고하며 데이터를 재현하기 위해 명시적 예방 요인을 필요로 하지 않는다.