시간 시계열 지표로 가족계획 조사 데이터와 모델 추정치 탐색

본 논문은 전 세계 가족계획 모니터링에서 핵심적인 현대 피임법 사용률을 다루며, 설문 조사 간격이 길고 데이터가 부족한 현실을 극복하기 위해 UN 인구부가 개발한 Family Planning Estimation Model(FPEM)을 활용한다. FPEM은 베이지안 계층시계열 모델로, 국가를 하위 지역, 지역, 전 세계 수준의 계층 구조에 배치해 정보 공유를 통해 데이터가 희박한 국가에도 신뢰성 있는 연간 추정치를 제공한다. 그러나 이러한 모델은 본질적으로 평활화(smoothing)된 결과를 산출하므로, 실제 설문 관측치와의 일치성을 검증하는 것이 필수적이다. 연구는 두 주요 데이터 소스를 준비한다. 첫째, UNPD가 관리하는 ‘World Contraceptive Use’ 데이터베이스에서 197개 국가·1970~2021년 기간의 현대 피임법 사용률 설문 데이터를 수집한다. 설문 유형은 DHS, MICS, PMA, 국가 조사 등으로 구분되며, 동일 연도·국가에 다중 관측이 존재할 경우 신뢰도가 높은 순서(DHS → MICS → PMA → 국가 조사)로 하나의 값으로 정제한다. 둘째, 동일 기간에 대한 FPEM 연간 추정치를 fpemlocal 패키지를 통해 확보한다. 두 데이터셋은 국가‑연도 형태의 패널 구조로 정렬되어 wdiexplorer 패키지에 입력된다. wdiexplorer는 세 단계 워크플로우를 제공한다. (1) 데이터 완전성·결측 검토 단계에서는 국가·연도별 관측 빈도를 시각화해 하위 지역별 데이터 격차를 파악한다. (2) 진단 지표 계산 단계에서는 총 10개의 시계열 특성 지표를 제공하지만, 본 연구에서는 실루엣 폭(silhouette width)과 추세 강도(trend strength)를 선택한다. 실루엣 폭은 사전 정의된 하위 지역 그룹 내에서 각 국가의 시계열 패턴이 얼마나 일관된지를 측정한다. 값이 +1에 가까울수록 해당 국가가 그룹 특성을 잘 따르고, −1에 가까울수록 다른 그룹과 더 유사함을 의미한다. 추세 강도는 시계열 분해에서 트렌드와 잔차의 분산 비율로 정의되어, 0에 가까우면 잡음이 많고 1에 가까우면 부드러운 장기 추세를 나타낸다. (3) 시각화·해석 단계에서는 설문 관측치와 모델 추정치를 동일 그래프에 겹쳐 보여, 두 시계열이 시간적으로 어떻게 일치하거나 차이가 나는지를 직관적으로 확인한다. 분석 결과, 데이터가 풍부한 국가(예: 케냐, 인도네시아)에서는 실루엣 폭이 0.7~0.9 수준으로 높고, 추세 강도도 0.8 이상으로 나타나 모델이 실제 변화를 잘 포착한다는 것을 확인했다. 반면, 설문 데이터가 드문 사하라 이남 아프리카나 일부 중남미 국가에서는 실루엣 폭이 0.2~0.4로 낮고, 추세 강도가 0.5 이하로 감소한다. 이는 모델이 하위 지역 평균에 과도히 끌려가 실제 국가별 변동성을 충분히 반영하지 못함을 의미한다. 특히, 특정 국가에서 설문이 급격히 상승하거나 하락하는 경우(예: 탄자니아 2018년 급증) 모델은 완만한 로지스틱 성장 곡선으로 평활화해 차이를 보인다. 모델 적합성을 추가로 검증하기 위해 잔차(관측치 − 추정치)에 대해 선형성(linearity)과 곡률(curvature) 지표를 적용한다. 잔차 시계열을 시간에 대한 1차·2차 정규 다항식으로 회귀시켜 β1(선형 계수)와 β2(곡률 계수)를 추정한다. 이상적인 경우 β1≈0, β2≈0이어야 하지만, 사하라 이남 아프리카 일부 국가에서는 β2가 0.15~0.25 수준으로 양의 곡률을 보이며, 모델이 비선형 상승을 충분히 포착하지 못함을 시사한다. 이러한 국가에서는 추가적인 지역 특화 사전 분포나 외부 변수(예: 보건 투자, 교육 수준)를 모델에 포함시키는 것이 개선 방안으로 제시된다. 논문은 또한 wdiexplorer의 적용 가능성을 논의한다. 원래는 세계 개발 지표(WDI) 패널 분석을 위해 설계된 패키지이지만, 가족계획과 같이 관측 간격이 불규칙하고 데이터 양이 제한적인 분야에도 유연하게 적용될 수 있음을 보여준다. 다만, 실루엣 폭과 추세 강도는 연속적인 연도 데이터가 전제이므로, 설문 간격이 3~5년인 경우 보간법이나 시계열 보정이 필요할 수 있다. 결론적으로, 본 연구는 (1) FPEM이 전반적으로 설문 데이터와 높은 일치성을 보이지만, 데이터가 희박한 국가에서는 하위 지역 평균에 과도히 의존하는 경향, (2) 실루엣 폭과 추세 강도라는 시계열 진단 지표가 모델-관측 일치를 정량적으로 평가하는 데 유용함, (3) 잔차의 선형성·곡률 분석을 통해 모델이 놓치고 있는 비선형 패턴을 식별할 수 있음을 입증한다. 향후 연구에서는 더 정교한 베이지안 사전 설정, 국가별 동적 공변량 도입, 그리고 시계열 보간·다중 관측 통합 방법을 개발해 모델 정확도를 향상시키는 방향을 제시한다.