[논문 리뷰] Epidemiological dynamics of the 2009 Influenza A(H1N1)v outbreak in India

T. Jesan, Gautam I. Menon외 1명|arXiv (Cornell University)|2010. 06. 03.

COVID-19 epidemiological studies참고 문헌 6인용 수 7

한 줄 요약

이 연구는 2009년 6월 1일부터 9월 30일까지 인도의 랩터리 확진자 시간열 자료를 사용하여 2009년 H1N1 패ニック의 기본 재생수 $ R_0 \approx 1.45 $ 를 추정한다. 부드럽게 처리된 발생률 자료에 대해 기하급수적 성장 모델링과 부트스트랩 재표집을 적용한 결과, 평균 이하의 전파성으로 나타나 전국적 도시 확산에도 불구하고 지역 간 전파 역학의 이질성이 존재함을 시사한다.

ABSTRACT

We analyze the time-series data for the onset of A(H1N1)v influenza pandemic in India during the period June 1- September 30, 2009. Using a variety of statistical fitting procedures, we obtain a robust estimate of the exponential growth rate $\langle λ angle \simeq 0.15$. This corresponds to a basic reproductive number $R_0 \simeq 1.45$ for influenza A(H1N1)v in India, a value which lies towards the lower end of the range of values reported for different countries affected by the pandemic.

연구 동기 및 목표

• 2009년 H1N1 패ニック의 초기 유행기 동안 인도에서 기본 재생수 $ R_0 $ 를 추정하기 위해.
• 매일의 확진자 발생률 자료에 대한 통계적 피팅을 통해 A(H1N1)v의 전파성을 평가하기 위해.
• 부트스트랩 재표집과 이동 평균을 사용하여 추정된 기하급수적 성장률 $ \lambda $ 의 탄력성 평가하기 위해.
• 특히 푸네와 챔라이 사이의 전파에 영향을 미칠 수 있는 기후적 또는 인구통계적 요인이 지역 간 전파 차이를 설명하는지 조사하기 위해.
• 다른 국가에서 보고된 값들과 비교하여 인도의 $ R_0 $ 가 패닉 확산의 지역적 다양성에 어떻게 기여하는지 평가하기 위해.

제안 방법

• 2009년 6월 1일부터 9월 30일까지 인도 보건복지부에서 제공한 매일의 랩터리 확진자 자료를 사용하였다.
• 잡음이 많은 발생률 자료의 노이즈를 줄이고 일일 변동성을 감소시키기 위해 5일 이동 평균을 적용하였다.
• 부드러운 자료에 기하급수적 성장 모델 $ I(t) \propto e^{\lambda t} $ 를 피팅하여 기하급수적 성장률 $ \lambda $ 를 추정하였다.
• 평균 세대 간격 $ \tau = 3 $ 일로 설정하여 $ R_0 = 1 + \lambda \tau $ 를 계산하였으며, 결과적으로 $ R_0 \approx 1.45 $ 를 도출하였다.
• 1000개의 표본을 사용한 부트스트랩 재표집을 통해 $ \lambda $ 의 불확실성 추정을 수행하였으며, 재표집 기울기 분포로부터 신뢰구간을 유도하였다.
• 2009년 7월 말경 주요 대도시들 사이에서 감염 수가 동시에 증가하는 패턴을 관찰하여 지역 전파의 시작 시점을 특정하였다.

실험 결과

연구 질문

• RQ1 2009년 H1N1 패닉의 기본 재생수 $ R_0 $ 는 얼마이며, 다른 국가에서 보고된 값들과 비교해 볼 때 어떻게 다를까?
• RQ2 자료 노이즈와 보고 변동성으로 인해, 인도에서 추정된 기하급수적 성장률 $ \lambda $ 의 탄력성은 어느 정도일까?
• RQ3 H1N1 바이러스는 어느 시점에서 수입 사례에서 지속적인 지역 전파로 전환되었는가?
• RQ4 기후적 또는 지리적 요인(예: 고도 또는 도시 간 연결성)이 푸네와 챔라이와 같은 도시 간 발생률 차이를 설명하는 데 어느 정도 기여하는가?
• RQ5 9월 30일 이후의 통제 조치 및 보고 방식 변화는 추정된 $ R_0 $ 의 신뢰성에 어떤 영향을 미치는가?

주요 결과

• 2009년 6월 1일부터 9월 30일까지 부드러운 발생률 자료를 기반으로, 인도의 유행에 대한 기하급수적 성장률 $ \lambda $ 는 $ \langle\lambda\rangle \simeq 0.15 $ 매일로 추정되었다.
• 세대 간격 $ \tau = 3 $ 일로 설정하여 기본 재생수를 $ R_0 \simeq 1.45 $ 로 계산하였으며, 이는 글로벌 범위에서 관찰된 값들 중 낮은 편에 속한다.
• 2009년 7월 31일부터 8월 15일까지의 기간에 대해 1000개의 표본을 사용한 부트스트랩 재표집 결과, 평균 기울기 $ \langle\lambda'\rangle = 0.166 $ 이며 표준편차는 0.024였으며, 이는 추정치의 탄력성을 확인한다.
• 유행은 2009년 7월 말경에 지역 전파로 전환되었으며, 이는 6개의 주요 대도시에서 감염 수가 동시에 증가함으로써 확인되었다.
• 낮은 인구 밀도와 낮은 연결성에도 불구하고 푸네는 1238건의 사례를 기록한 반면, 쿠타카타는 113건에 그쳐, 기후적 요인이 전파에 영향을 미쳤을 가능성이 있다.
• Chennai는 주로 수입 사례를 보였고, 푸네는 지역 감염 비율이 높아, 더 시원한 기후가 푸네에서의 지속적 전파를 촉진했을 가능성이 있다.