Sex-Specific Temporal Dynamics and Early Warning Signals of Adolescent Suicide Mortality in South Korea: A Nationwide Time-Series Study of Trends, Seasonality, Forecasting, and Lagged Associations With Suicidal Behaviors, 2014–2023

기술 요약: 한국 청소년 자살 사망의 성별 특이적 시계열 역동성 및 조기 경고 신호 (2014–2023)

문제 제기

자살 사고(ideation), 계획(planning), 시도(attempts)는 개념적으로 자살 연속체 상의 점진적인 단계로 이해되지만, 이러한 상위 행동들의 시계열적 변동이 인구 수준에서의 자살 사망 변화를 선행하는지에 대한 실증적 근거는 부족한 실정이다. 기존 연구들은 주로 개인 수준의 위험 요인이나 자살 역학의 고립된 측면(예: 예측 또는 계절적 변동 단독)에 집중해 왔으며, 이러한 구성 요소들을 통합하지 못했다. 또한, 전 세계적으로 여학생들의 심리적 고통과 비치명적 자살 행동이 급격히 증가하고 있음에도 불구하고, 자살 사망의 시계열적 역동성이 성별에 따라 어떻게 다른지, 그리고 특정 자살 전 단계가 성별 특이적인 조기 경고 지표로서 기능하는지는 여전히 불분명하다. 본 연구는 청소년 자살률이 지속적으로 높은 국가인 한국을 대상으로, 자살 행동과 청소년 자살 사망 간의 성별 특이적 시계열 역동성을 조사함으로써 이러한 공백을 메우고자 한다.

연구 방법론

본 연구는 2014년 1월부터 2023년 12월까지의 데이터를 분석하는 전국 단위 시계열 설계를 채택하였다.

• 데이터 출처:
  • 자살 사망: 통계청 마이크로데이터 통합서비스(MDIS)를 통해 13~19세 청소년의 월별 사망 데이터를 확보하였으며, ICD-10 코드 X60–X84를 통해 사망을 식별하였다.
  • 자살 행동: 한국 청소년 위험 행동 조사(KYRBS)로부터 자살 사고, 자살 계획, 자살 시도의 연간 유병률을 도출하였다. 도출된 연간 추정치는 조사 실시 시기의 중간값에 맞춰 선형 보간법을 통해 월별 시계열 데이터로 변환되었다.
• 통계적 프레임워크: 남성과 여성을 분리하여 다음의 통합된 시계열 기법들을 사용하여 분석을 수행하였다:
  • 일반화 가법 모델 (GAM): 비선형적 시간 궤적을 조사하기 위해 사용되었다.
  • STL 분해: 사망 시계열을 추세(trend), 계절성(seasonal), 잔차(residual) 성분으로 분리하였다.
  • 변곡점 분석 (Breakpoint Analysis): Bai–Perron 구조적 변화 프레임워크를 사용하여 유의미한 구조적 변화를 식별하였다.
  • 예측: 계절적 자기회귀 누적 이동평균(SARIMA) 모델을 사용하여 2030년까지의 사망 추세를 투영하였다.
  • 시차 연관성 (Lagged Associations): 교차 상관 함수(CCF)를 통해 최적의 시차 구조를 식별하였다. 외생 변수를 포함한 자기회귀 누적 이동평균(ARIMAX) 모델과 그랜저 인과관계 검정을 통해 자살 행동이 이후의 사망을 유의미하게 예측하는지 평가하였다.

주요 결과

1. 자살 사망의 시계열 역동성

• 추세: 양쪽 성별 모두 2014년부터 2023년까지 자살 사망의 지속적인 상승 궤적을 보였다. 그러나 여성은 2010년대 후반부터 더 큰 비선형성과 뚜렷한 가속화를 나타냈다.
• 구조적 변곡점: 성별 특이적인 변화가 식별되었다. 주요 구조적 변곡점은 여성의 경우 2019년 4월, 남성의 경우 2020년 4월에 발생하였다.
• 계절성: 양쪽 성별 모두 봄에서 초여여름(4월6월)에 정점을 찍고 겨울(12월1월)에 저점을 찍는 반복적인 12개월 주기 계절 패턴을 보였다.
• 예측: SARIMA 모델은 2030년까지 사망률이 계속 증가할 것으로 투영하였다. 현재 절대적인 비율은 남성이 더 높지만, 여성의 사망 궤적은 더 가파르게 상승하여 2030년까지 점진적으로 남성 수준에 근접할 것으로 예측된다.

2. 자살 행동과 사망 간의 연관성

• 여성: 자살 계획과 자살 사망 사이에 강하고 시간적으로 인접한 연관성이 발견되었다.
  • 교차 상관: 자살 계획은 1개월의 시차에서 가장 강력한 상관관계를 보였다 ($r = 0.40$).
  • ARIMAX: 1개월 시차의 자살 계획은 이후의 사망을 유의미하게 예측하였다 ($\beta = 0.41, p = .041$).
  • 그랜저 인과관계: 여성의 자살 계획은 자기회귀적 의존성을 넘어 이후의 사망 예측력을 유의미하게 향상시켰다 ($F = 9.37, p = .003$).
  • 2개월 시차의 자살 시도는 이후의 사망과 역의 상관관계를 보였다.
• 남성: 자살 사고, 계획 또는 시도와 이후의 자살 사망 사이에서 유의미한 단기적 시계열 연관성은 발견되지 않았다. 이러한 연관성은 여성에 비해 약하고 시간적으로 분산되어 있었다.

핵심 기여

• 인구 수준의 연결: 본 연구는 자살의 위험 요인을 개인 수준에서 식별하는 것을 넘어, 국가적 자살 사망 추세를 인구 수준의 자살 행동 지표와 직접 연결하였다.
• 성별 특이적 조기 경고 신호: 본 연구는 자살 계획을 여성의 청소년 자살 사망에 대한 구체적이고 시간적으로 인접한(1~2개월 창) 조기 경고 신호로 식별하였으며, 이는 남성에게서는 관찰되지 않은 관계이다.
• 통합적 시계열 프레임워크: 장기적 추세, 계절적 패턴, 구조적 변곡점 및 시차 연관성을 동시에 분석함으로써 자살 위험의 역동적 진행에 대한 포괄적인 관점을 제공한다.
• 예측 격차: 여성 사망률이 2030년까지 남성 수준에 수렴할 수 있다는 투영 결과는 타겟팅된 주의가 필요한 역학적 지형의 변화를 강조한다.

의의 및 주장

저자들은 연구 결과가 획일적인 자살 예방 전략에서 성별 특이적 시계열 접근법으로의 전환이 필요함을 시사한다고 주장한다:

• 여성 청소년 대상: 자살 계획과 사망 사이의 강한 12개월 시차는 임계 개입 창(intervention window)을 시사한다. 논문은 특히 자살 계획이 나타난 후 12개월 이내에 집중적인 스크리닝, 모니터링 및 신속한 개입을 수행할 것을 주장한다.
• 남성 청소년 대상: 보고된 행동으로부터의 단기 예측 신호가 부족하다는 점을 고려할 때, 예방 전략은 보다 광범적이고 지속적인 접근 방식에 집중해야 한다. 여기에는 단기적인 행동 지표에 의존하기보다 학교 및 또래 기반 지원 시스템 강화, 사회적 고립 감소, 그리고 급격한 위험 상승에 대한 지속적인 감시 체계 유지가 포함된다.
• 계절적 취약성: 양쪽 성별 모두에서 나타나는 봄-초여름의 정점은 불확실성(예: 새 학기)으로 특징지어지는 전환기가 성별과 관계없이 예방 정책에서 우선순위가 되어야 함을 나타낸다.

본 연구는 전통적으로 자살 시도가 주요 위험 예측 인자로 간주되어 왔으나, 자살 계획이 여성의 경우 임박한 사망에 대해 더 시간적으로 인접한 인구 수준의 지표 역할을 할 수 있으며, 이는 예방을 위한 임상적으로 의미 있는 창을 제공한다고 결론짓는다.