저감도 전염병 전파율 추정의 새로운 역방향 방법

초록

본 논문은 사례수가 적은 저감도 전염병 상황에서 전파율 β(t)을 안정적으로 추정하기 위해 로그 변환 후 B‑스플라인 보간을 수행하고, 이를 지수화하여 양의 값을 보장하는 ‘지수 B‑스플라인 역방향 방법’을 제안한다. 연속형 및 이산형 역방향 프레임워크에 모두 적용 가능하도록 설계했으며, 중국의 사람옹, 다중균주 인플루엔자, 연령구조 인플루엔자 데이터를 통해 기존 방법 대비 MAE가 크게 감소함을 실증하였다.

상세 분석

이 연구는 저감도 전염병 모델링에서 가장 빈번히 발생하는 ‘음수 전파율’ 문제를 근본적으로 해결한다는 점에서 의미가 크다. 기존 역방향 방법은 관측된 신규 감염 사례 yₖ가 매우 작거나 0에 가까울 때, 다항식이나 일반 B‑스플라인 보간을 적용하면 보간곡선이 음수 구간을 형성한다. 이는 SEIR 모델의 β(t) = … 식에 대입하면 비물리적인 음수 전파율을 초래하고, 수치적 불안정성까지 야기한다. 저자들은 이러한 문제를 ‘로그‑지수 변환’이라는 두 단계 절차로 회피한다. 먼저 yₖ에 로그를 취해 Yₖ = ln yₖ를 얻고, 이 로그값을 B‑스플라인으로 정확히 보간한다. B‑스플라인은 매끄러운 2차 이상 다항식 기반이므로, 로그값 자체는 부호가 보존된다. 이후 보간함수 f(t) = ∑Pᵢ Bᵢ,ₚ(t)를 지수화하여 ỹ(t) = exp f(t) 형태의 연속적인 신규 감염 추정치를 만든다. 이 과정에서 ỹ(t) ≥ 0가 보장되며, 특히 yₖ = 0인 경우에도 exp(−∞) → 0 형태로 자연스럽게 처리된다.

연속형 역방향에서는 ỹ(t)와 모델식 dE/dt = σ ỹ(t) 를 이용해 H(t) = dE/dt를 정의하고, 전체 인구 N(t), 감수성 S(t), 감염 I(t), 회복 R(t) 를 각각 미분방정식의 해로부터 구한다. β(t)는 H(t)와 기존 파라미터를 이용한 식 β(t) =