DFA α1

DFA α1은 로그-로그 직선의 기울기로, 심장 박동 시계열의 흔들림이 관찰 창의 길이에 따라 어떻게 커지는지를 4~16박동 창에서 측정한 값입니다. 단위는 없습니다.

이 지표가 만들어진 이유

Peng, Havlin, Stanley, Goldberger는 1995년에 표준편차로는 답할 수 없는 질문을 던지기 위해 추세 제거 추세 제거 추세 제거 스펙트럼을 계산하기 전에 박동 간격 시계열에서 느린 드리프트를 없애는 작업입니다. 드리프트가 가장 낮은 주파수에 가짜 파워로 쌓이는 것을 막습니다. 추세 제거 더 알아보기 변동 분석(DFA)을 소개했습니다. 박동이 서로를 기억하는가 하는 질문입니다. 순수한 잡음에서는 각 박동이 독립적이어서 시계열의 누적합이 술 취한 사람의 걸음처럼 헤매고, 장거리 상관이 있는 시계열에서는 규모가 커질수록 이 헤맴이 더 빠르게 커집니다. 지수 $\alpha$가 이 차이를 담아냅니다. 백색 잡음은 0.5, 1/f 잡음은 1.0, 랜덤 워크는 1.5입니다. 건강한 심장 박동 시계열에는 장거리 상관이 있었고 이것이 중증 심부전에서는 무너졌습니다. 같은 논문은 교차점도 보고했는데, 짧은 창의 기울기(α1)가 긴 창의 기울기(α2)와 다르며 그 경계는 약 16박동입니다. Sassi 연구진의 2015년 입장문은 DFA를 신호에 대한 이해를 높인 프랙탈 방법의 하나로 꼽으면서, 지금까지 임상에서의 활용은 다소 제한적이라고 봅니다.

공식

$$y(k) = \sum_{i=1}^{k}\left(\mathrm{RR}_i - \mathrm{RR}_{\mathrm{mean}}\right),\qquad F(n) = \sqrt{\frac{1}{N}\sum_{k=1}^{N}\left(y(k) - y_n(k)\right)^2},\qquad F(n) \propto n^{\alpha}$$

시계열을 적분합니다. $y(k)$는 간격이 평균에서 벗어난 값의 누적합입니다. $y$를 $n$박동 크기의 박스로 자르고 각 박스에 직선 $y_n$을 맞춘 뒤, 남은 값의 제곱평균제곱근 $F(n)$을 구합니다. 모든 박스 크기에 대해 이를 반복하고 $\log F(n)$을 $\log n$에 대해 그립니다. 416박동에서의 기울기가 α1이고 1664박동에서의 기울기가 α2입니다(Kubios는 기본으로 412와 1364를 쓰고, Sassi 연구진은 많은 임상 연구가 α1을 약 4~11박동에서 구한다고 지적합니다).

계산 예시

리듬은 과거를 기억할까?

교육용 기록 5분

박스 크기(박동 수)에 대한 변동 F(n), 로그–로그 α1 α2

4 16 64

α1(4~16박동)

0.49

α2(16~64박동)

0.65

충분한 휴식 기록은 합성한 것입니다. 네 박동마다 깨끗한 호흡 파동이 있고 박동 사이에 무작위 잡음이 섞여 있어서, 4~16박동에서는 거의 잡음처럼 보이고(α1 0.49), 더 느린 파동 때문에 16~64박동에서는 어느 정도 기억이 생깁니다(α2 0.65). 안정 시의 실제 기록은 보통 α1이 약 1 이상입니다.

충분한 휴식 백색 잡음 1/f 잡음 이소성 박동 있음

섞은 백색 잡음 없음

JavaScript를 켜면 조작 도구가 작동합니다.

박스 크기 4~64박동에 대한 변동을 나타낸 로그–로그 그래프로, 박스 크기마다 점이 있고 맞춘 직선이 두 개 있습니다. 하나는 4~16박동 구간의 α1용이고, 다른 하나는 16~64박동 구간의 α2용입니다. 수치 표시에는 α1과 α2가 나오고, 그 뒤에 이를 설명하는 문장이 이어집니다. 버튼 네 개로 충분한 휴식 기록, 백색 잡음, 1/f 잡음, 이소성 박동 있는 기록 중에서 고르고, 슬라이더로 선택한 시계열에 백색 잡음을 섞습니다.

교육용 기록의 오 분 구간에서는 α1 = 0.49, α2 = 0.65입니다. 이 α1이 낮은 것은 의도한 결과입니다. 이 기록은 합성한 것으로, 네 박동에 한 번꼴의 깨끗한 호흡 파동에 박동 사이의 무작위 잡음이 섞여 있어 짧은 창에서는 거의 무작위처럼 보이고, 더 느린 혈압 파동과 드리프트가 긴 창에 어느 정도 기억을 줍니다. 안정 시의 실제 기록에서는 보통 α1이 약 1이거나 그 이상입니다. 백색 잡음으로 바꾸면 두 기울기가 모두 0.5 근처에 놓이고(α1 0.48, α2 0.50), 1/f 잡음으로 바꾸면 둘 다 1 근처에 놓입니다(α1 1.04, α2 1.03). 어떤 시계열이든 잡음을 섞으면 α1이 0.5 쪽으로 미끄러집니다.

DFA α1을 움직이는 요인

교육용 기록에서 보듯, 안정 시에 크고 규칙적인 호흡 흔들림이 있으면 α1이 내려갈 수 있습니다. 운동 강도도 값을 움직여, 가벼운 운동에서 약 1이던 값이 힘든 운동에서는 0.5 쪽으로 내려갑니다. Rogers 연구진은 2021년에 러너 15명에서 α1이 0.75를 가로지르는 지점이 제1 환기 역치(VT1)와 일치함을 확인했고, 이것이 일부 앱이 제공하는 DFA α1 역치 방법의 바탕입니다. 나이와 심장 질환도 값을 바꿉니다. Shaffer와 Ginsberg는 DFA가 몇 시간에 걸친 기록을 대상으로 설계되었다고 지적하며, α1의 초단기 최소 기록 시간으로 제안된 값으로 120초를 듭니다.

오해하기 쉬운 경우

조기 박동과 아티팩트는 기억이 없는 도약을 더해 α1을 잡음 쪽으로 끌어당기므로 시계열을 먼저 정리해야 합니다. 박스 범위는 프로그램마다 달라서 Kubios의 α1과 Peng의 α1이 같지는 않습니다. 짧은 기록에서는 박스가 너무 적어서 α2를 신뢰하기 어렵습니다. 그리고 1에 가까운 값이 ‘더 좋은’ 것은 아닙니다. 하나의 패턴일 뿐이며, 상황에 따라 더 높은 값도 더 낮은 값도 더 나쁜 결과와 함께 나타났습니다.

어디에서 볼 수 있나

Kubios HRV는 비선형 지표에서 DFA α1과 α2를 보여 줍니다. 내 가슴 스트랩 기록으로 α1을 계산하려면 RR 변환기로 Garmin, Polar, Movesense 내보내기 파일에서 박동 간 간격을 꺼낼 수 있습니다. 일부 훈련 앱은 운동 중에 α1을 실시간으로 보여 주어 유산소 역치를 찾는 데 씁니다. Task Force는 이 지표를 다루지 않으며, Sassi 연구진의 입장문이 근거를 검토합니다. 비선형 HRV 지표란 무엇이고 필요할까?에서 웨어러블 기기를 쓰는 독자에게 이런 내용이 중요한지 따져 봅니다.