신생아 중환자실 후기 발병 패혈증 진단을 위한 심박 변동성·호흡 신호 분석

초록

본 연구는 미숙아 NICU에서 후기 발병 패혈증을 비침습적으로 조기에 탐지하기 위해 심박 변동성(HRV)과 코호흡 신호를 분석한다. 단일채널(HRV)과 양채널(HRV‑호흡) 접근법을 적용해 선형·비선형 지표를 추출했으며, αSlow·αFast와 Sample Entropy가 감염군을 구분하는 주요 파라미터로 확인되었다. 또한 0.2–0.4 Hz 대역의 지역선형 상관계수 r²ₜ,𝑓와 양방향 비선형 회귀계수 h²가 HRV와 호흡의 기능적 결합을 나타내어 진단에 보완적 역할을 함을 보였다. 6개 최적 특징을 융합한 테스트는 가장 높은 AUC와 낮은 오경보 확률을 기록하였다.

상세 분석

이 논문은 후기 발병 패혈증(Late‑Onset Sepsis, LOS)이 신생아 중환자실(NICU)에서 특히 저체중 영아에게 높은 사망률과 합병증을 초래한다는 임상적 배경에서 출발한다. 기존의 혈액 배양이나 임상 징후는 진단에 시간이 많이 소요되거나 민감도가 낮아 조기 개입이 어려운 점을 보완하고자, 저침습적 생리 신호인 심박 변동성(HRV)과 코호흡 신호를 동시에 활용한다는 점이 혁신적이다.

연구는 ‘Mono‑Channel(MC)’와 ‘Bi‑Channel(BC)’ 두 가지 분석 프레임워크를 설정했다. MC는 HRV 단독으로 시간 영역(예: 평균 RR 간격, 표준편차)과 주파수 영역(저주파·고주파 파워)뿐 아니라 비선형 지표인 Detrended Fluctuation Analysis(αSlow, αFast)와 Sample Entropy를 적용한다. BC는 여기서 한 걸음 더 나아가 HRV와 코호흡 신호 사이의 동적 결합을 정량화한다. 구체적으로, 0.2 Hz < f < 0.4 Hz 대역에서 지역선형 상관계수 r²ₜ,𝑓를 계산해 두 신호가 동일 주파수 성분에서 얼마나 동조하는지를 평가했으며, 비선형 회귀계수 h²를 양방향(HRV→Resp, Resp→HRV)으로 추정해 인과적 혹은 비대칭적 상호작용을 탐색했다.

통계 분석 결과, αSlow와 αFast는 감염군에서 유의하게 감소했으며, 이는 패혈증 시 자율신경계의 조절력이 약화됨을 시사한다. Sample Entropy는 감염군에서 증가했는데, 이는 HRV 시계열이 더 불규칙하고 복잡해짐을 의미한다. r²ₜ,𝑓는 0.2–0.4 Hz 대역에서 감염군이 비감염군보다 현저히 낮아, 호흡과 심박 사이의 동기화가 손상되었음을 보여준다. h²는 양방향 모두에서 차이가 나타났으며, 특히 HRV→Resp 방향에서 감소가 두드러져 자율신경계가 호흡 조절에 미치는 영향이 약화된 것으로 해석된다.

이후 연구팀은 MC와 BC에서 도출된 6개 최적 특징(αSlow, αFast, Sample Entropy, r²ₜ,𝑓, h²(→), h²(←))을 선형 결합 및 로지스틱 회귀 모델에 입력해 ROC 곡선을 작성했다. 결과는 AUC가 0.92에 달해 기존 임상 지표 대비 뛰어난 판별력을 보였으며, 오경보 확률(PFA)은 8% 이하로 낮아 실제 NICU 환경에서의 적용 가능성을 높였다.

핵심적인 강점은 (1) 다중 스케일 선형·비선형 지표를 동시에 활용해 패혈증의 복합적 생리 변화를 포착, (2) HRV와 호흡의 기능적 결합을 정량화함으로써 단일 신호만으로는 놓칠 수 있는 미세한 변화를 보완, (3) 실시간 모니터링이 가능한 비침습적 데이터 기반 진단 모델을 제시했다는 점이다. 반면 제한점으로는 표본 크기가 비교적 작고, 단일 기관 데이터에 국한돼 외부 검증이 부족하다는 점, 그리고 신호 품질(움직임 아티팩트 등) 관리가 실제 임상 적용 시 과제로 남는다. 향후 다기관 대규모 코호트와 머신러닝 기반 고차원 특징 선택을 통해 모델의 일반화와 자동화가 기대된다.