파트너 내·외부 HIV 전파율 추정: 성별 차이를 반영한 확률 모델

초록

본 연구는 이성 파트너십 내에서 발생하는 HIV 감염을 ‘내부 전파(τ)’와 외부로부터의 ‘외부 전파(λ)’ 두 가지 요인으로 구분하고, 이를 확률적 SI 페어 모델에 수식화한다. 남·여 각각의 외부 감염률(λₘ, λ𝒻)과 성별 특이 내부 전파율(τₘ→𝒻, τ𝒻→ₘ)을 포함한 확장 모델을 제시하고, 탄탄한 폐쇄형 해와 다항식 likelihood를 이용해 2년간 1,802쌍의 코호트 데이터를 기반으로 최대우도 추정(MLE)과 불확실성(신뢰구간)을 계산한다. 결과는 외부 전파가 상대적으로 낮지만 내부 전파가 파트너 내 감염을 주도한다는 점을 보여주며, 성별에 따라 전파 강도가 다름을 확인한다.

상세 분석

이 논문은 HIV 전파를 “파트너 내 전파(τ)”와 “파트너 외 전파(λ)”라는 두 개의 강도로 단순화한 SI(감수성‑감염) 페어 모델을 구축한다. 기본 모델은 세 가지 페어 상태 SS, SI, II를 정의하고, ODE 시스템
dP_SS/dt = –2λP_SS,
dP_SI/dt = 2λP_SS – (λ+τ)P_SI,
dP_II/dt = (λ+τ)P_SI
을 제시한다. 이 시스템은 폐쇄형 해를 갖으며, P_SS(t)=P_SS(0)e^{-2λt} 등으로 명시된다. 여기서 τ와 λ은 각각 내부와 외부 감염률이며, τ는 λ와 결합된 형태(τ=λ+θ)로 재파라미터화해 식의 안정성을 확보한다.

성별 차이를 반영하기 위해 모델을 확장한다. 남성·여성 각각에 대한 외부 감염률 λₘ, λ𝒻와 성별 전파율 τₘ→𝒻, τ𝒻→ₘ을 도입하고, 네 개의 상태(SS, IₘS𝒻, SₘI𝒻, IₘI𝒻)를 추적한다. 이 확장 모델 역시 폐쇄형 해를 제공하며, λₘ+λ𝒻가 전체 외부 감염 압력을, τₘ→𝒻와 τ𝒻→ₘ이 파트너 내 전파 비대칭성을 설명한다.

데이터는 탄자니아 무안자에서 2년간 추적한 1,802쌍의 이성 파트너 코호트이며, 초기와 2년 후의 SS, SI, II 상태를 관측했다. 관측값을 다항분포와 연결해 로그-우도 함수를 구성하고, λ와 τ에 대한 편미분을 통해 MLE를 도출한다. 특히 λ는 λ̂ = (1/2T)·ln(N₀_SS/N_T_SS) 로 간단히 계산되며, τ는 재파라미터화된 ϕ를 이용해 근사식으로 추정한다. 성별 모델에서는 λ을 λ̂·q, λ̂·(1–q) 형태로 분해하고, θₘ, θ𝒻를 관측된 성별 SI 변화를 통해 추정한다.

모델 식별성은 로그-우도 곡면의 볼록성 및 Hessian 행렬의 역을 통해 검증한다. 파라미터 추정 결과는 외부 전파율 λ≈0.003 yr⁻¹(≈0.3%/년)이며, 내부 전파율 τ는 약 0.02 yr⁻¹(≈2%/년)으로, 내부 전파가 외부 전파보다 약 6~7배 강함을 시사한다. 성별 분석에서는 남성→여성 전파가 여성→남성보다 약 1.5배 높게 추정돼, 남성 감염자가 파트너에게 전파할 위험이 더 크다는 공중보건적 함의를 제공한다.

통계적 불확실성은 로그-우도 곡면의 2차 미분(Hessian)으로부터 얻은 공분산 행렬을 이용해 95% 신뢰구간을 계산한다. 결과는 λ와 τ 모두 통계적으로 유의미하지만, 특히 τ의 신뢰구간이 넓어 데이터의 제한성(관측 시점이 두 개뿐)과 파트너 내 재감염을 무시한 단순화가 원인일 수 있다.

이 논문의 주요 기여는 (1) 폐쇄형 해를 갖는 간단하면서도 실용적인 파트너 기반 SI 모델을 제시, (2) 다항식 likelihood와 MLE를 결합해 파라미터를 직접 추정, (3) 성별 차이를 정량화해 정책 설계에 활용 가능한 정보를 제공한다는 점이다. 한계로는 파트너 지속 기간, 항레트로바이러스 치료(ART) 효과, 프리엑스포저 프로필(PrEP) 등 시간에 따라 변하는 요인을 고려하지 않았으며, 관측 시점이 두 개뿐이라 모델 검증이 제한적이다. 향후 연구에서는 연속적인 관찰 데이터와 동적 파트너 형성·해체 과정을 포함한 확장 모델이 필요하다.