[논문 리뷰] The D model for deaths by COVID-19

J. E. Amaro|arXiv (Cornell University)|2020. 03. 30.

COVID-19 epidemiological studies참고 문헌 2인용 수 4

한 줄 요약

이 논문은 SI 전염병학 모델에서 유도된 간단한 분석 프레임워크인 D 모델을 제안한다. 이 모델은 회복이 없다는 가정 하에 감염자 수의 지연되고 스케일링된 함수로 사망자를 기술하고 예측한다. 중국, 스페인, 이탈리아의 데이터에 대해 세 가지 조정 가능한 매개변수를 사용하여 모델을 적합시켰으며, 최소한의 가정에도 불구하고 단기 예측에서 패안의 정점 시점과 정점 수준을 정확하게 예측할 수 있다.

ABSTRACT

We present a simple analytical model to describe the fast increase of deaths produced by the corona virus (COVID-19) infections. The 'D' (deaths) model comes from a simplified version of the SIR (susceptible-infected-recovered) model known as SI model. It assumes that there is no recovery. In that case the dynamical equations can be solved analytically and the result is extended to describe the D-function that depends on three parameters that we can fit to the data. Results for the data from Spain, Italy and China are presented. The model is validated by comparing with the data of deaths in China, which are well described. This allows to make predictions for the development of the disease in Spain and Italy.

연구 동기 및 목표

• 복잡한 수치 시뮬레이션에 의존하지 않고도 코로나19 사망 경로를 신속하게 예측할 수 있는 분석 모델을 개발하는 것.
• 회복이 없다는 가정 하에 SI 프레임워크를 기반으로 한 단순화된 모델을 실제 사망 데이터에 적합시켜 검증하는 것.
• 현재 유행 중인 국가에서 초기 팬데믹 예측을 위한 저매개변수, 해석 가능한 도구를 제공하는 것.
• 중국에서는 팬데믹이 이미 정점에 도달한 상태였고, 스페인과 이탈리아에서는 유행이 활발히 진행 중이었기 때문에, 중국의 데이터를 기반으로 한 모델의 예측 능력을 평가하는 것.

제안 방법

• 모델은 SI 방정식에서 유도되며, 회복이 없다는 가정(R=0) 하에 감염자 수 I(t)에 대한 단일 미분방정식으로 단순화된다.
• SI 방정식의 해는 I(t)에 대한 로지스틱 유형의 함수로 변환되며, 이는 사망자 수 D(t)를 I(t)의 지연되고 스케일링된 버전으로 모델링하는 데 사용된다.
• D(t) 함수는 D(t) = a * exp((t-t₀)/b) / (1 + c * exp((t-t₀)/b))로 표현되며, 세 가지 피팅 매개변수인 a(진폭), b(시간 상수), c(초기 조건 스케일링)를 포함한다.
• a, b, c 매개변수는 전파율이나 회복률과 같은 기초 생물학적 매개변수를 알 필요 없이 중국, 스페인, 이탈리아의 실제 사망 데이터에 비추어 비선형 회귀를 통해 피팅된다.
• 감염과 사망 사이의 고정된 지연 τ와 일정한 사망률 m은 모두 피팅된 매개변수에 통합된다.
• 모델은 팬데믹 곡선이 정점에 도달한 중국의 데이터를 바탕으로 검증된 후, 스페인과 이탈리아에서 정점 시점과 최종 사망자 수를 예측하는 데 사용된다.

실험 결과

연구 질문

• RQ1 SI 프레임워크를 기반으로 한 단순화된 분석 모델이 여러 국가에서 코로나19 사망 경로를 정확하게 기술할 수 있는가?
• RQ2 세 매개변수만을 사용하는 D 모델이 중국, 스페인, 이탈리아의 실제 사망 데이터에 얼마나 잘 맞는가?
• RQ3 유행의 상승 단계에 있는 국가에서 모델이 팬데믹 정점과 최종 사망자 수를 신뢰성 있게 예측할 수 있는가?
• RQ4 회복이 없다는 가정이 모델의 예측 정확도와 적용 가능성에 미치는 영향은 무엇인가?

주요 결과

• D 모델은 중국의 사망 데이터를 놀랍게 잘 맞추었으며, 이는 다른 국가에 대한 예측 도구로의 활용 가능성을 검증한다.
• 스페인의 경우 모델은 팬데믹이 유행 시작 후 약 21일째에 정점에 도달할 것으로 예측했으며, 2020년 4월 중순까지 정점에 도달할 것으로 예상했다.
• 이탈리아의 경우 모델은 데이터 수집 시작일로부터 약 20일 만에 팬데믹 정점에 도달할 것으로 예측했으며, 스페인과 유사한 경향을 보였다.
• D 모델의 세 매개변수 함수는 모든 세 국가에서 사망 곡선의 기하급수적 증가와 이후 정점에 도달하는 곡선을 매우 정밀하게 반영하였다.
• a, b, c 매개변수만으로도 전파율이나 회복 시간과 같은 기초 생물학적 매개변수를 추정할 필요 없이 전체 사망 곡선을 기술하는 데 충분하였다.
• 중국의 데이터만을 기반으로도 스페인과 이탈리아와 같이 팬데믹 초기 단계에 있는 국가에 대해 모델이 강력한 예측 능력을 보였다.