감염병 인류 라는 책을 재미있게 읽는 중이다. 136페이지에는 기초감염재생산지수와 대유행의 가능성에 대한 간단한 수식이 나온다. \[\text{대유행의 가능성} = 1 - \frac{1}{R_0}\]
\(R_0\)는 기초감염재생산지수, 즉 한 사람의 감염자가 다른 모든 사람이 감수성자 (susceptible) 일때 감염시키는 평균 감염자수를 나타낸다. 위 수식은 제한적인 경우에만 적용된다. 감염병 유입 시 대유행의 가능성에 대해 좀 더 일반적으로 적용될 수 있는 방법에 대해 적어보려고 한다. 물론 아래 내용도 상당히 이상적인 상황에 대한 기술일 뿐이고 현실은 그 보다 훨씬 더 복잡할 것이다. 감염병 유입 시 대유행의 가능성에 대한 내용은 Niels G. Becker가 저술한 Modeling to Inform Infectious Disease Control의 제 2장에 자세하게 기술되어 있다.
동일한 사람들로 이루어진 인구 집단에서 감염병이 퍼져나가는 현상을 생각해 보자. 한 사람이 평균적으로 \(R_0\) 명을 감염시킨다고 할 때 \(R_0\)은 평균일뿐 후속 감염자수는 어떤 확률 분포를 가진다고 생각해볼 수 있다. 한 명의 감염자가 총 \(j\) 명의 후속 감염자를 만들어 내는 확률을 \(P(X=j) = p_j\)라고 할 때 \(R_0\)는 아래와 같이 표현할 수 있다.
\[R_0 = \sum_{k=0}^\infty j p_j.\]
대유행의 가능성은 역으로 생각하는 것이 유리하다. 즉 대유행이 아닌 소규모의 감염으로 막을 내리는 확률, \(\theta\), 을 구한 후 대유행의 확률은 \(1-\theta\)로 구하는 것이다. 소규모의 감염이 일어나기 위해서는 유입된 초기 환자 (index patient)가 아무도 감염시키지 않거나 혹은 몇 명을 감염시켰다고 할지라도 후속 감염자들이 추가적으로 일으키는 감염이 소규모일때만 가능할 것이다. 즉 소규모 감염의 확률, \(\theta\), 는 아래의 식을 만족한다.
따라서 \(\theta\)는 아래의 식을 계산하면 된다. 다만 해를 직접 구할 수는 없고 수치해석방법을 이용 해서 답을 구해야 한다.
\[\theta = e^{(1-\theta)R_0}.\]
위에서 언급한 \(\text{대유행의 가능성} = 1 - \frac{1}{R_0}\)는 후속 감염자수의 분포가 기하분포 (geometric distribution)를 따를때 성립한다. 기하분포는 아래와 같이 표현되 \[\text{Prob}(X=j) = (1-p)^j p\]
평균과 성공확률의 관게, \(R = \frac{1-p}{p}\), 이용하여 다시 표현하면 아래와 같다.
\[\text{Prob}(X=j) = (\frac{R_0}{1+R_0})^j \frac{1}{1+R_0}\] 위에서와 동일하게 계산하면 아래와 같다. \[\sum_{j=0}^\infty p_j \theta^{j} = \frac{1}{1 + (1-\theta)R_0}\] 따라서 아래의 식을 풀면 \(\theta\)를 구할 수 있다.
\[\theta = \frac{1}{1 + (1-\theta)R_0}.\]
\[\theta = \frac{1}{R_0}.\]
따라서 감염병 인류 책에서 언급된 것처럼 \[\text{대유행의 가능성} = 1 - \frac{1}{R_0}.\] 마지막으로 후속 감염자수가 이항분포 (negative binomial distribution)을 따른다고 가정해보자. 최근 연구들에서 빈번하게 언급되고 가장 현실에 가까운 가정인 듯 하다.
다양한 모수를 사용하여 표현할 수 있고 평균 \(R_0\)과 확산(dispersion; \(k\))를 이용하여 나타내면 아래와 같다. 이항 분포는 \(k=1\)일 때는 기하분포와 동일하고 \(k=\infty\) 일때는 푸아송분포와 동일하다. 코로나 19의 경우 \(k\) 값이 약 0.55이다.
\[\text{Prob}(X=j) = \frac{\Gamma(k+j)}{j!\Gamma(k)}\left(\frac{k}{k+R_0}\right)^k\left(\frac{R_0}{k+R_0}\right)^j \text{ for } k=0,1,2,...\] The probability of a minor outbreak, 위와 동일한 방법으로 계산하면 우변은 아래와 같다.