HPV 6와 11 역학 모델링을 위한 적응형 마코프 체인 몬테카를로 전방 시뮬레이션

본 논문은 인간 유두종바이러스(HPV) 6형과 11형의 전파를 정량적으로 분석하기 위해, 연령·성별·성행위 위험군별로 구획된 미분방정식(ODE) 기반 전염병 모델을 설계하고, 이를 베이지안 통계 프레임워크 안에서 캘리브레이션하는 방법론을 제시한다. 1. **모델 설계** - 인구를 15‑59세 사이의 5년 구간 연령층으로 구분하고, 각 연령층을 네 단계의 성행위 위험군(활동 수준)으로 나눈다. - 각 구획은 감수성(S), 감염(I), 치료(T), 회복(P) 등 HPV 감염·치료·면역 과정을 나타내는 상태 변수로 구성된다. - 성적 접촉은 Garnett‑Anderson이 제안한 혼합 행렬을 기반으로 하며, 동질성(assortativity) 파라미터와 연령·위험군별 접촉 비율을 확률적 변수로 두어 데이터와 동시에 추정한다. 2. **베이지안 캘리브레이션** - 모델 파라미터(전파 확률, 잠복기, 감염 지속시간, 치료 기간, 면역 지속시간, 혈청전환 확률 등)과 혼합 행렬 파라미터에 대해 비정보적 사전분포를 설정한다. - 관측 데이터는 호주에서 수집된 성기 사마귀 발생률과 HPV 6·11 혈청역학 자료를 사용한다. - 사후분포는 복잡한 비선형 ODE와 높은 차원의 파라미터 공간 때문에 직접 계산이 불가능하므로, MCMC 샘플링이 필요하다. 3. **전방 투영 적응형 MCMC** - 기존 MCMC는 제안 분포의 공분산을 사전 튜닝해야 하지만, 고차원 ODE 상태 변수와 파라미터가 동시에 변동하면 튜닝 비용이 급증한다. - 저자들은 ‘전방 투영(Forward Projection)’ 단계에서 현재 파라미터 후보에 대해 ODE를 수치적으로 풀어 시간별 상태 궤적을 생성하고, 이를 관측값과 비교해 likelihood를 계산한다. - 적응형 메트로폴리스 알고리즘은 체인 진행 중에 제안 분포의 평균·공분산을 온라인으로 업데이트한다. 이 과정은 Diminishing Adaptation과 Bounded Convergence 조건을 만족하도록 설계되어, 체인의 에르고딕성을 이론적으로 보장한다. - 구현상, 적응 단계는 매 1000번의 샘플링마다 공분산을 재추정하고, 제안 스케일을 자동 조정한다. 4. **실험 및 검증** - **합성 데이터**: 알려진 파라미터값을 사용해 데이터를 생성하고, 제안된 방법으로 복원한다. 결과는 파라미터 추정 편향이 거의 없으며, 유효표본수(ESS)가 기존 고정 제안 MCMC 대비 2~3배 향상됨을 보여준다. - **실제 호주 데이터**: 성기 사마귀 연령별 발생률과 HPV 6·11 혈청역학을 이용해 모델을 캘리브레이션한다. 사후 예측 구간은 관측치를 잘 포함하고, 혼합 행렬의 동질성 파라미터가 높은 값을 가지며 연령·위험군 간 접촉이 제한적임을 시사한다. - 또한, 백신 도입 시나리오를 사후 예측 분포에 적용해 집단 면역(herd immunity) 효과를 정량화하였다. 5. **논의 및 결론** - 전방 투영을 포함한 적응형 MCMC는 고차원 비선형 역학 모델의 파라미터와 사회적 접촉 구조를 동시에 추정하는 데 효율적이며, 튜닝 비용을 크게 절감한다. - 제안된 베이지안 프레임워크는 파라미터 불확실성을 정량화하고, 정책 입안자가 백신 효과와 비용‑효과 분석에 활용할 수 있는 사후 예측을 제공한다. - 향후 연구에서는 남녀 비대칭 인구 구조, 이민 및 여행자 유입, 그리고 다형성(다중 HPV 유형) 모델을 확장하는 방향을 제시한다.