여성 1명이 낳는 자녀의 수가 평균 6명을 넘던 1960년대 우리나라는 인구 증가 억제 정책을 도입하였다. 그 이후 20년에 걸쳐 아이 적게 낳기, 둘 낳기, 하나 낳기와 같은 산아제한 정책으로 1980년대에 이르러 여성 1명이 평생 낳을 것으로 예상되는 평균 출생아 수를 의미하는 지표인 합계출산율이 2명으로 감소되었다. 1990년대 말 외환위기와 맞물려 산아 제한 정책을 공식 폐지하였고 2000년대 중반이 되어서야 저출산 문제를 인식하고 출산율 회복을 위한 방향으로 인구정책을 변경하였다[1,2].

현재 우리나라는 합계출산율이 1명도 되지 않은 심각한 저출산 문제와 함께 고령화 사회에 직면하게 되었다. 출산율의 저하는 생산가능 인구의 감소로 미래의 노동력 부족 현상뿐만 아니라 보건 의료분야에도 영향을 미칠 것으로 예상한다. 실제로 산부인과와 소아청소년과 분야의 위축으로 인하여 전문의 서비스보다는 일차진료 성격의 서비스 제공이 증가하고 있는 것이 현실이다. 아동의 구강건강을 다루는 소아치과도 예외는 아닐 것으로 생각되며 저출산 정책이 시행된 2005년부터 지속적으로 출산율 현황, 소아치과의 환자 분포 및 진료 현황, 소아치과 전문의원 분포, 보호자의 인식 변화 등과 같은 연구들이 진행되어 왔다[3-5].

이 연구는 최근 10년간의 광주광역시 출생아 수와 전남대학교 치과병원 소아치과 내원 환자 수라는 시계열 자료에 대해 이론적 우수성과 실용성이 증명된 예측기법인 시계열 분석 방법을 적용하였다. 시계열 분석 방법은 시간의 흐름에 따라 관측된 자료를 바탕으로 확률적 또는 통계적 성질들을 파악하고, 과거 시계열의 형태가 미래에도 같은 형태로 반복될 것이라는 가정 하에 모형을 구축하여 미래에 대해 예측을 하는 방법이다. 시계열 분석을 실시하고자 할 때, 기본적인 3가지 절차가 있다. 먼저 시계열 그래프를 작성하여, 시계열 자료가 안정적 시계열인지, 비안정적 시계열 인지를 파악하는 과정이며, 두 번째는 시계열 자료가 비안정적일 경우에 차분과 분산의 안정화 변환을 통하여 안정화시킨다. 세 번째 절차는 여러 가지 모형 식별 통계량을 이용하여 적합한 모형을 선택하는 것이다[6].

이러한 시계열 분석을 통하여 출생아 수와 내원 환자 수에 대한 모형을 구축하고 향후 1년을 예측하여 실제 관측값과 비교하여 모형 안정성을 확인하였다. 소아치과 분야와 시계열 분석을 접목하여 감소하는 출산율에 대한 소아치과 영역의 수요 경향을 파악하고 미래를 대비하고 예측하는 데 이용할 수 있는 모형 구축 방법을 제시하고자 하였다.

이 연구는 전남대학교의 생명윤리심의위원회(Institutional Review Board, IRB)의 심의를 통과한 후 시행되었다(IRB No: CNUDH-EXP-2021-008).

이 연구는 통계청의 월별 인구동향조사로부터 얻어진 광주광역시의 출생아 수와 전남대학교 치과병원 소아치과의 월별 내원 환자 수에 대하여 120개월에 걸친 시계열 자료를 분석하고자 하였다. 이 연구 결과에서 출생아수는 지속적으로 감소하고 소아치과의 내원 환자 수는 유지됨을 확인하였고, 저출산이 가져올 소아치과 분야의 변화에 대하여 대비할 수 있는 기초자료를 제공하고자 하였다.

2020년 통계청 발표에 따르면 출생통계를 작성하기 시작한 1970년대 합계출산율이 4.53명인 이래로 2020년 합계출산율은 역대 최저치인 0.84명이었다[7,8]. 임신 가능한 연령대의 인구 감소, 만혼과 초산 연령의 상승, 혼인율의 하락 등을 저출산의 주요 원인으로 보고되고 있다[9]. 출산율 저하는 사회 경제적으로 노동력 부족 현상뿐 아니라 의료 및 기타 산업에도 영향을 미치고 있다. 산부인과는 저출산과 직결된 분야로 저출산이 특히 산과의 의료서비스 수요에 직접적인 영향을 미쳐 전문인력의 급격한 감소를 초래했지만, 초산 연령의 증가로 인한 부인과 서비스의 수요 증가에 따라 산부인과 의료서비스 이용이 다시 안정세에 오르고 있다고 하였다. 출생아 수가 감소함에 따라 전체적인 아동 인구가 감소하게 되었고, 아동의 건강 문제를 다루는 소아과 진료의 수요도 감소하였다. 소아과 서비스의 수요 감소는 새로 배출되는 소아과 전문의의 감소에도 영향을 미치지만 이미 배출된 소아과 전문의의 진료영역을 아동의 특수 분야가 아닌 다른 연령층의 건강문제를 다루는 일차 진료의사로서의 기능으로 확대되고 있다고 하였다[10].

이에 소아의 구강건강을 다루는 소아치과 영역도 예외는 아닐 것으로 생각되며 출산율의 변화와 함께 실제 소아치과 의료서비스 수요의 변화를 확인하는 것이 큰 의미가 있을 것으로 생각되었다. 2008년 치과전문의 제도가 시행됨에 따라 소아치과 전문의 배출로 소아전문 치과의 증가와 의료 정보에 대한 접근성 향상으로 보호자들의 인식이 많이 변화하고 있다[5]. Son 등[11]과 Kim 등[12]의 연구에서 구강위생에 대한 인식이 높아지면서 소아치과에 초진 내원하는 환자의 평균연령이 감소한다고 하였다. Kim 등의 연구에 따르면 전남대학교 치과병원 소아치과의 진료현황이 2011년에 수복치료, 예방치료, 외과적 치료, 치수치료, 교정치료 순에서 2020년에는 예방치료, 수복치료, 외과적 치료, 교정치료, 치수치료 순으로 변화하여 치수치료가 감소하고 예방치료가 높은 비율을 차지하게 되었다고 하였다. 출산율 감소와 함께 소아치과에 내원한 환자 수에 관한 연구도 2000년대 이후 지속적으로 이루어지고 있다. Ko 등[13]의 연구에서도 1998년부터 2005년까지 전북대학교 소아치과에 내원한 환자의 수는 지속적으로 증가하다가 2005년에 다소 감소하였는데 이는 소아전문 치과의 개원과 함께 전라북도 지역의 지속적인 출생률 감소의 영향이라고 하였다. Kim 등[14]의 연구에 따르면 2005년부터 2010년까지 조선대학교 소아치과의 내원 환자 수와 광주에 있는 소아전문 치과의 내원 환자 수의 연구에 따르면 내원 환자 수는 지속적으로 증가하였으며 월별 분포의 경우 개인 치과의 경우에는 방학기간과 내원 사이에 특정 상관관계가 없었으나 대학병원에서는 방학기간 내원 환자의 수가 증가한다고 하였으며 이는 이 연구에서 1월과 8월에 내원 환자 수가 가장 많았던 결과와 일치하였다. Kang 등[15]의 연구에서 2008년부터 2012년까지 연세대학교 소아치과에 내원한 환자의 수는 뚜렷한 변화가 없이 비교적 일정하게 유지되었다고 하였다. 월별 분포의 경우 방학기간에 내원 환자 수가 증가하였으며 이는 이 연구와 일치한 결과이다. 최근 출생률 감소로 인해 소아치과 수요 감소를 우려하였으나 우려와는 달리 지역별로 차이가 있을 수 있지만 일정 수준 꾸준히 유지되고 있음을 알 수 있었다. 이는 출생률의 감소와 더불어 보호자의 인식 변화로 자녀에 대한 구강관리에 관심이 증가하고, 전문성이 강조되는 사회변화와 IT 발달로 정보 교환이 용이해지는 등 소아치과에 내원할 수 있는 접근성 향상의 결과로 볼 수 있다.

Woo 등[16]의 연구에 따르면 시계열 분석이 우리나라에서 학술적인 목적으로 이용된 시기는 1960년대부터이며, 금융, 환경, 공학, 교육, 정책 등의 여러 분야에서 다양한 방법으로 이용되어 왔다. 국내 보건복지 분야를 대상으로 한 논문은 1987년 이래로 100건 정도이며 전체 논문 수는 많지 않으나 점차 증가하는 추세라고 하였다. 국민 건강을 주제로 질병 발생 및 사망, 자살의 원인과 그 영향요인을 다루거나 환자 수와 의료 이용량, 진료비 변화와 같은 의료서비스 이용에 대한 내용이 주를 이루었고 치의학과 관련된 연구는 전무한 실정이다. 주로 보건 복지 분야에서 사용되는 시계열 분석 모형은 AR 모형, MA 모형, 누적 모형, 이들의 조합으로 제안된 ARMA 모형, ARIMA 모형이 사용된다. 시계열 분석의 결과를 신뢰하기 위해서 일반적으로 50개 이상의 자료가 있어야 한다고 알려져 있으며 모형 적합을 통하여 미래를 예측하는 경우 향후 1개월 - 3년까지 예측하는 것이 어느 정도 신뢰성을 확보할 수 있고, 6개월에서 1년의 기간을 예측하는 것이 일반적이라고 하였다[17,18]. 이 연구도 120개월의 데이터를 토대로 기본적인 분석절차에 따라 수행하였으며 적합한 모형을 통하여 향후 1년을 예측하였다. 광주광역시의 출생아 수의 경우 10년간의 월별 시계열도가 계절성을 가지며 급격하게 감소하는 추세를 보이고 1월, 3월의 출생아 수가 많고 11월, 12월이 출생아 수가 적었으며 자기회귀오차모형이 최종모형으로 선택되었다. 전남대학교 치과병원 소아치과의 내원 환자 수의 경우에도 계절성을 가지며 약간 상승하는 추세를 보였으며 8월, 1월의 내원 환자 수가 많고, 6월, 3월의 내원 환자 수가 적었으며 자기회귀오차 모형이 최종모형으로 선택되었다. 출생아 수와 내원 환자 수에 대해 향후 1년을 예측하고 실제 관측값과 비교하였는데 코로나 19라는 특수한 상황으로 인하여 예측치보다 다소 낮게 확인되었다. 출생아 수에 대한 내원 환자 수의 비율은 지속적으로 증가하였으나 향후 1년의 실제 비율값이 예측값보다 더 큰 이유는 출생아 수가 예측했던 수치보다 적었기 때문으로 생각된다. 특히 연구에 포함된 내원 환자 수에 출생률의 영향을 받지 않는 환자도 포함되어 있을 수 있어 비율 값으로 결과를 해석하기에는 한계가 존재한다. 이는 장기간에 걸쳐 관측이 이루어지는 시계열 자료의 특성상 정책의 변화, 환경적 재해 등과 같은 제어할 수 없는 사건들에 의하여 기존의 관측된 시계열 자료의 패턴과 달라질 수 있으므로 추가 연구에서 충분히 고려되어야 할 사항이다[17].

이 연구의 한계점은 연구 대상이 광주지역의 출생아 수와 전남대학교 치과병원 소아치과에 내원한 환자만을 대상으로 하였다는 것이다. 특히 내원 환자 수를 연령, 초진, 재진과 같은 기준이 없이 모든 연구대상으로 포함시켰고, 3차 의료기관인 대학병원의 특수성을 고려하지 않았다는 점에서 이 연구의 결과를 일반화하기는 어렵다고 할 수 있다. 추후 전국적인 출생아 수, 소아 전문치과의원 및 타 대학병원 소아치과 내원 환자의 자료 수집을 통한 추가적인 연구가 필요할 것으로 생각된다.

이 연구에서는 시계열 분석이라는 방법을 통하여 최근 10년간의 광주지역의 출생아 수와 전남대학교 소아치과 내원 환자 수에 대하여 분석하고 향후 1년을 예측하여 저출산이 실제 소아치과의 수요에 영향을 미치는지 알아보고자 하였다. 최근 10년간 광주지역 출생아 수의 급격한 감소에도 불구하고 소아치과에 내원한 환자의 수는 일정 수준 유지되며 다소 증가할 것으로 예측되었다. 하지만 COVID-19이라는 특수한 상황으로 인하여 실제 출생아 수와 내원 환자 수가 예측치보다 다소 낮게 확인되었으나 예측 범위 내에 포함됨을 확인하였다. 시계열 분석 모형은 저출산 시대에 소아치과 영역에 대한 수요를 유지하고 미래를 대비하기 위한 기초 도구로써 활용될 수 있을 것이다.