생체표지자 변화와 사건 위험 연계 분석 소아 HIV/AIDS 장기 추적 연구

초록

본 논문은 반복 측정되는 바이오마커의 변화와 사건 위험을 동시에 추정하는 새로운 공동 종단‑생존 모델을 제안한다. 종단 부분에 cubic B‑spline를 적용해 기울기와 적분을 손쉽게 계산하고, MCMC로 파라미터를 추정한다. 우간다 소아 HIV 감염 영아 데이터를 이용해 모델을 적용하고, DIC, CPO, ROC 곡선 등을 통해 모델 선택과 예측력을 평가한다.

상세 분석

본 연구는 바이오마커의 시간적 변화를 사건 위험과 연결짓는 공동 모델링에 대한 이론적·실용적 기여가 크다. 종단 모델에 cubic B‑spline 기반의 베이시안 프레임워크를 도입함으로써, 관측 시점이 불규칙하고 누락이 존재하는 실제 임상 데이터에서도 매끄러운 연속 곡선을 추정한다. B‑spline의 계수를 통해 손쉽게 1차 미분(즉, 순간 변화율)과 구간 적분(누적 노출량)을 구할 수 있어, 위험 함수에 두 가지 형태의 시간‑의존적 요인을 동시에 포함시킬 수 있다. 위험 모델은 일반적인 Cox 비례위험 구조를 확장해, 선형 predictor에 바이오마커의 현재값, 기울기, 적분값을 각각 가중치와 함께 삽입한다. 이는 “현재 상태”와 “변화 속도”, “누적 노출”이 사건 발생에 미치는 독립적·복합적 영향을 정량화한다는 점에서 의미가 크다. 파라미터 추정은 Gibbs sampler와 Metropolis‑Hastings를 결합한 MCMC 알고리즘으로 수행되며, 사전분포는 비정보적 정규·역감마를 사용해 데이터에 의한 후향적 업데이트를 강조한다. 모델 적합도와 예측력을 평가하기 위해 Deviance Information Criterion(DIC)과 Conditional Predictive Ordinate(CPO)를 계산하고, 시간별 ROC 곡선을 도출해 민감도·특이도 균형을 시각화한다. 특히, CPO는 사후 예측 분포를 기반으로 각 관측치의 교차 검증 가능성을 제공해, 과적합 위험을 정량적으로 판단하게 한다. 실증 분석에서는 우간다의 장기 추적 코호트(생후 0~24개월)에서 CD4%와 바이러스 부하를 주요 바이오마커로 설정하고, 사망·중증 감염을 사건으로 정의한다. 결과는 바이오마커의 급격한 감소(음의 기울기)와 높은 누적 부하가 사망 위험을 유의하게 증가시킴을 보여, 임상적 의사결정에 직접 활용 가능한 위험 예측 지표를 제공한다. 전반적으로 이 모델은 복합적인 시간‑의존 효과를 포괄적으로 다루면서도 베이시안 추정의 유연성을 유지해, 다양한 의료·역학 연구에 적용 가능성을 넓힌다.