A prognostic model and simplified risk score for predicting poor treatment response in children aged 2–59 months hospitalized with severe community-acquired pneumonia: Development and validation study in Hakim Gizaw Hospital, Ethiopia

기술 요약: 에티오피아 아동의 중증 폐렴에 대한 예후 모델 및 간소화된 위험 점수

문제 제기
폐렴은 특히 저소득 및 중저소득 국가에서 아동 사망의 주요 원인으로 남아 있다. 중증 폐렴은 주요한 임상적 관심사이지만, 임상적으로 매우 이질적이다. 동일한 진단 기준을 충족하는 아동이라 할지라도 치료 실패 위험은 매우 다를 수 있다. 기존의 소아 폐렴 예측 도구(예: RISC, RISC-Malawi)는 사례 구성, 결과 유병률, 임상 관리 방식의 차이로 인해 그 모델이 도출된 환경 이외의 지역에 적용될 때 가변적인 성능을 보였다. 따라서 입원 초기 단계에서 치료 반응이 좋지 않을 고위험군 아동을 식별하여 우선적인 모니터링과 보조적 관리를 용이하게 할 수 있는, 현지에 적합하고 내부 검증된 예후 모델이 절실히 필요하다.

방법론
본 연구는 에티오피아 하킴 기자우(Hakim Gizaw) 교육 병원에 입원한 중증 지역사회 획득 폐렴 환아 630명(연령 2~59개월)의 데이터를 사용하여 예후 모델 개발 및 내부 검증 분석을 수행하였다.

• 결과 정의: 일차 결과는 입원 후 72시간 이내의 "치료 반응 불량"으로 정의하였으며, 이는 입원 기간 연장(>72시간), 입원 중 사망, 및/또는 72시간 후 항생제 변경을 포함하는 복합 지표로 조작되었다.
• 예측 변수: 예측 변수는 예후적 유용성을 보장하기 위해 입원 시 또는 초기 임상 평가 시 사용 가능한 변수로 제한하였다. 여기에는 인구통계학적 특성, 예방 접종 상태, 활력 징후, 산소 포화도, 실험실 측정값 및 초기 항생제 처방이 포함되었다.
• 데이터 처리: 네 가지 실험실 변수(WBC, 호중구, 림프구, 헤모글로빈)에서 상당한 결측치(약 18%)가 발생하였다. 본 연구는 완전 사례 분석(complete-case analysis) 시 손실될 수 있는 정보를 보존하기 위해 예측 평균 매칭(predictive mean matching)을 이용한 연쇄 방정식에 의한 다중 대치법(MICE)을 사용하여 5개의 대치된 데이터 세트를 생성하였다.
• 모델 개발:
  • 스크리닝: 탐색적 단변량 분석을 통해 후보 변수를 선별하였다(p < 0.25 임계값).
  • 선택: 아이카이け 정보 기준(AIC)을 사용하여 다변량 로지스틱 회귀 모델을 비교하였다. 최종 모델은 16개의 예측 파라미터를 유지하였다.
  • 검증: 낙관주의 편향(overfitting)을 추정하고 교정하기 위해 1,000회의 부트스트랩 재표본 추출을 사용하여 내부 검증을 수행하였다.
• 성능 지표:
  • 판별력: 수신자 조작 특성 곡선 아래 면적(AUC) 및 Somers' Dxy를 통해 평가하였다.
  • 교정력: 교정 도표(calibration plots), 절편, 기울기, 최대 절대 오차(Emax), 브라이어 점수(Brier score)를 사용하여 평가하였다.
  • 간소화된 점수: 임상 현장에서의 사용을 용이하게 하기 위해 전체 모델으로부터 11개 구성 요소로 이루어진 정수형 위험 점수를 도출하였다.
  • 임상적 유용성: "전원 치료(treat-all)" 및 "미치료(treat-none)" 전략과 비교하여 다양한 임계 확률에 걸친 순 이익(net benefit)을 평가하기 위해 결정 곡선 분석(DCA)을 실시하였다.

주요 기여

1. 현지 모델 개발: 본 연구는 글로벌 도구를 현지 환경으로 전이할 때 발생하는 한계를 해결하기 위해 에티오피아 맥락에 특화된 다변량 로지스틱 회귀 모델을 제공한다.
2. 엄격한 통계적 접근: 본 연구는 결측 데이터 처리를 위한 다중 대치법, 내부 검증을 위한 부트스트랩 재표본 추출, 유용성 평가를 위한 결정 곡선 분석을 활용하여 현대적 보고 지침(예: TRIPOD)을 준수한다.
3. 이중 출력: 본 연구는 고정밀 위험 추정을 위한 전체 회귀 모델과 실질적인 임상 적용을 위한 11개 구성 요소의 간소화된 위험 점수를 모두 제공한다.
4. 투명한 보고: 저자들은 향후 외부 검증을 지원하기 위해 전체 회귀 계수, 오즈비(odds ratios) 및 상세한 성능 지표(낙관주의 교정값 포함)를 제공한다.

결과

• 결과 유병률: 치료 반응 불량은 267명의 아동(42.4%)에게서 나타났다.
• 전체 모델 성능:
  • 판별력: 겉보기 AUC는 0.875(95% CI 0.849–0.902)였다. 부트스트랩 교정 후, 낙관주의 교정 AUC는 약 0.857(Dxy 0.713)이었다.
  • 교정력: 모델은 교정된 기울기 0.927과 절편 0.021로 적절한 교정력을 보였다. 브라이어 점수는 0.153이었다.
  • 주요 예측 변수: 치료 반응 불량의 유의한 예측 변수로는 불완전한 예방 접종(AOR 3.79), 매우 낮은 산소 포화도(AOR 2.23), 12~23개월 연령(AOR 3.29), 암피실린-겐타마이신 치료(반코마이신 대비), 천명(stridor), 청색증, 섭식 불능, 의식 상태 변화, 구토, 설사 및 백혈구 수 증가가 포함되었다.
• 간소화된 점수 성능:
  • 점수 범위는 0에서 19까지였다.
  • 판별력: 점수의 AUC는 0.823으로, 전체 모델보다 통계적으로 유의하게 낮았으나(p < 0.001) 유용한 판별력을 유지하였다.
  • 위험 계층화: ≥8점 임계값에서, 이 점수는 저위험군(0–3점, 관찰된 치료 반응 불량률 7.7%)에 비해 현저히 높은 70.3%의 관찰된 치료 반응 불량률을 보이는 고위험군을 식별하였다.
  • 분류 지표: ≥8점 임계값에서 민감도는 76.4%, 특이도는 76.3%였다.
  • 유용성: 결정 곡선 분석 결과, 광범위한 임계 확률 범위에서 "미치료(treat-none)" 전략에 대해 양(+)의 순 이익을 나타냈다.

의의 및 주장
저자들은 개발된 모델이 이 특정 코호트에서 치료 반응 불량을 예측하는 데 있어 좋은 판별력과 수용 가능한 내부 교정력을 보여준다고 주장한다. 간소화된 점수는 전체 회귀 방정식보다 정확도는 낮지만, 아동들을 결과가 확연히 다른 그룹으로 분리하는 실용적인 위험 계층화 도구를 제공한다.

그러나 본 논문은 임상 적용에 대해 신중한 입장을 유지한다. 저자들은 다음과 같이 명시하였다:

• 결과는 내부 검증에 기반하며, 모델이 일상적인 임상 사용에 고려되기 위해서는 독립적인 코호트를 통한 외부 검증이 필요하다.
• 본 연구는 전향적 영향 평가를 수행하지 않았으므로, 결정 곡선 분석은 임상적 효과성이 입증된 것이 아니라 잠재적 유용성을 나타낸다.
• 결과는 사망 종점이 아닌, 치료 반응에 대한 복합적인 운영적 척도이다.
• 단일 센터 설계 및 데이터 기반 예측 변수 선택으로 인해, 모델의 타 병원 또는 타 인구 집단에 대한 전이 가능성은 아직 입증되지 않았다.

결론적으로, 본 연구는 향-연구를 위한 재현 가능하고 통계적으로 견고한 토대를 제공하며, 이러한 도구가 임상적 의사결정을 안내하기 위해서는 외부 검증과 전향적 영향 평가가 필수적인 전제 조건임을 강조한다.