서울시 주거지의 근린환경특성이 개인의 신체활동과 건강 수준에 미치는 영향 분석 : 국민건강영양조사 (2007~2012) 자료를 활용한 다수준 분석의 적용

대부분의 선행 연구는 물리적 환경과 신체활동 그리고 건강 수준 사이의 통합적인 관점 보다는 단순화한 모형의 분석에 치중하고 있다. 먼저 커뮤니티의 물리적 환경과 지역주민의 보행시간에 관한 연구로 이경환·안건혁(2007)은 2005년 국민건강영양조사 자료를 분석하여 토지이용혼합도, 가로의 연결도, 근린공원 접근성 등이 보행시간과 양(+)의 상관관계에 있는 것으로 보고하였다. 도시계획시설인 공원을 주제로 문혜식·김수봉(2009)은 대구시 20세 이상 성인남녀를 대상으로 신체활동 공간으로써 공원이용실태를 분석하였다. 394개의 설문조사자료를 이용하여 빈도 및 교차분석을 진행하였으며 많은 사람이 신체활동을 위해 주거지 주변의 공원을 이용하고 있음을 확인하였다. 이는 국민생활체육참여실태조사(2012)에서 발표한 결과¹⁾와도 일치하며 상기 연구는 도시계획시설인 공원의 효용을 검증하였다는 점에서 의의가 있다. 같은 유형의 연구로 이슬기 외 4인(2013)은 창원시의 근린생활권을 대상으로 물리적 환경이 신체활동 목적의 공원 이용에 미치는 영향을 분석하였다. 설문조사자료를 이용하여 개인의 주소지와 운동 목적의 공원이용 여부 등을 변수로 만들어 로지스틱회귀분석 하였다. 분석결과, 근린단위에서의 토지이용혼합도, 도로비율, 교차로밀도, 공원까지의 최단거리와 공원의 수, 공원 출입구 수 등이 영향을 미치는 요인임을 확인하였다.

해외연구로 Kaczynski et al.(2008)은 공원 크기, 공원의 구성요소들과 신체활동참여자의 주거지로부터 공원의 거리가 신체활동을 위한 공원을 이용하는 데 있어 주요한 지표가 되는지를 연구하였다. 33개 공원의 28가지의 특정 구성요소를 자료로 활용하였고 공원 주변의 성인 거주자 380명을 대상으로 7일간의 관찰기록을 사용하였다. 로지스틱 회귀분석 결과, 공원의 구성요소 또는 관련시설이 증가할 때 신체활동 목적의 공원이용 가능성이 상대적으로 높고 공원의 구성요소 또는 관련 시설 등이 쾌적성보다 더욱 높은 지표임을 실증하였다. 특히 산책로가 많을 경우, 신체활동을 위해 공원을 더욱 이용한다는 강한 상관성을 확인하였다.

최근 발표된 성현곤 외 3인(2014)의 실증연구는 물리적 환경과 보행 간의 관계를 다룬 선행연구의 결과와 비교 및 재검증하는 형식으로 진행하였다. 2007~2011년 국민건강영양조사자료 서울시 149개 행정동, 4,361개의 설문표본을 활용하여 다수준 회귀분석을 실행하였다. 종속변수는 하루 보행시간, 일주일 걷기일 수, 일주일 보행시간 등을 사용하였다. 연구 결과, 지하철 접근성은 일주일간 보행시간과의 관계에서 유의하나 나머지 물리적 환경 특성은 유의하지 않은 것으로 나타났고 행정동 단위 근린환경의 물리적 특성과 보행활동 간의 연관성이 떨어지는 것으로 보고하였다. 이 연구는 근린의 물리적 환경과 보행활동 간의 관계에만 집중하였다는 데 있어 본 연구와 차이점이 있다.

해외 연구로 Giles-Corti and Donovan(2002)은 호주의 노동자, 주부 등을 대상으로 물리적 환경이 신체활동에 미치는 영향을 분석하였다. 로지스틱 회귀분석결과 여가시설의 인접은 시설 이용 빈도의 증가와 밀접한 연관성이 있는 것으로 나타났다. 사람들은 가로, 공개공지, 해변 순으로 이용 빈도가 높으며 충분히 물리적 환경, 개인 환경, 사회적 환경 등이 신체활동에 영향을 미침을 확인하였다. Frank et al.(2005)은 애틀랜타 대도시권 내 13개 카운티 거주자 357명을 대상으로 거주지 주변부 토지이용과 개발 밀도, 교차로 밀도 등이 신체활동 실측치와 어떠한 연관성이 있는지 분석하였다. 연구 결과 토지이용 특성과 주변 밀도, 교차로 밀도 등의 특성은 모두 신체활동량과 양(+)의 관계에 있는 것으로 보고하였다. 이어진 Frank et al.(2007) 연구에서는 개인의 선호선택을 고려하여 물리적 환경과 개인의 통행행태 간의 관계를 분석하였다. 미국 2,056명의 개인 통행 자료를 이용하여 로지스틱 회귀분석하였다. 분석결과 보행환경의 질이 좋을수록 신체활동량이 증가하는데 영향을 미치나 신체활동 자체의 선호수준이 낮은 경우 물리적 환경의 개선 효과가 나타나지 않아 실질적으로 개인의 의식수준에 따라 신체활동이 결정됨을 보고하였다.

신체활동량을 직접 측정하여 물리적 환경과의 관계를 살펴본 Troped et al.(2010)의 연구는 주거지와 직장을 중심으로 각각 50m, 1km 버퍼를 이용하여 개인의 신체활동량과 버퍼 영역별 물리적 환경 간의 관계를 분석하였다. 분석 결과, 교차로의 밀도, 토지이용 혼합도, 인구 밀도, 주거 밀도 등은 신체활동량 증가에 영향을 미치는 요소임을 확인하였다. 비슷한 방법으로 진행된 Ding et al.(2012)의 연구는 미국의 2개 지역에서 2,199명의 성인의 신체활동 수치를 측정한 자료를 이용하여 건조 환경과 인간의 심리 간의 상호작용이 신체활동참여에 미치는 영향을 분석하였다. 연구 결과, 자기효능감이 낮은 성인일수록 건조환경과 여가보행의 상관성이 강하게 나타나 상호작용 효과가 있음을 검증하였고 심리적으로 활동적이지 못한 경향이 있는 성인일수록 건조 환경에 따라 그 영향이 크게 좌우될 수 있음을 강조하였다.

Kelly et al.(2014)은 미국 인디아나 폴리스와 세인트루이스를 대상으로 건조 환경과 신체활동의 관계를 살펴보기 위해 지역적 특성과 신체활동 간의 연관성을 분석하였다. 인구분포와 사회경제적인 특성을 기준으로 근린단위의 가로구간의 건조 환경을 측정하였고 신체활동 행태를 관찰하였다. 분석결과, 복합적인 토지이용, 비 주거지역, 보행 관련 시설과 대중교통이 있는 가로구간에서 신체활동을 하는 개인의 총량이 매우 높게 나타남을 확인하였다. 이 연구는 미시적 단위의 특성과 신체활동 간의 관련성을 밝혀냈다는 점에서 의미가 있다.

물리적 환경과 건강 수준 간의 관계는 국내의 경우 최근까지 다양한 연구가 진행되었다. 초기 연구로 이경환·안건혁(2008)은 근린환경과 개인의 건강 수준이 통계적으로 유의한 관계에 있음을 실증적으로 분석하였다. 이와 유사한 주제로 비만도(BMI지수)와 물리적인 환경과의 관계를 살펴본 김은정·강민규(2011a)의 연구는 수도권의 53개 시군구별 비만 인구 비율과 건강인구 비율에 미치는 영향을 분석하였다. 국민건강영양조사자료(2005)를 이용하여 공간계량(자기상관)분석을 이용하였다. 연구의 결과, 비만인구 비율은 공간적인 자기 상관관계가 있고 도시의 물리적 환경은 개인의 건강에 영향을 미치는 요소이며 특히 비만인구에 더 많은 영향을 미치고 있음을 보고하였다. 국민건강영양조사 2008년 자료를 활용하여 전반적인 건강 수준인 EQ-5D, 삶의 질을 나타내는 EQ-VAS, 스트레스, 우울증 등을 종속변수로 사용한 성현곤(2011)은 다수준 회귀분석을 통해 주거지의 근린환경인 개발밀도, 토지이용복합도, 공원과 지하철 접근도 등 근린의 물리적 환경이 전반적인 건강 수준에 긍정적인 영향을 미치고 있음을 주장하였다.

해외 연구로 Ewing(2003)은 미국의 83개 대도시권, 448개 카운티를 대상으로 스프롤지수(용도의 혼합, 주거 밀도, 커뮤니티시설 밀집도 등을 활용한 지수)와 신체활동, 비만, 사망률 간의 관계를 분석하였다. 다층회귀분석 결과, 도시 스프롤지수, 보행, 비만은 서로가 영향관계에 놓여있음을 확인하였다. 비슷한 연구로 Kelly-Schwartz(2004)는 미국 29개 대도시권 지역을 대상으로 스프롤 지수와 주거밀도, 토지이용혼합도, 도로망 연결도와 체질량 지수, 만성질환과의 관계를 분석하였다. 다층회귀분석을 통해 스프롤과 신체활동, 건강과의 상호관계는 일정하지 않으며 도로망 연결도가 높은 주거지와 건강 간의 관계가 양(+)의 관계를 갖는 것으로 보고하였다. 많은 표본을 이용하여 도시의 스프롤과 비만 간의 관계를 분석한 Eid et al.(2008)은 도시의 스프롤과 토지이용 혼합도 등은 비만과 연관성이 낮으나 개인의 선호선택 비중이 높아 이 부분을 고려해야 함을 주장하였다. 비만과 관련된 지표 이외의 건강지표를 사용한 연구로 이경환(2012)의 연구는 지방 중소도시를 대상으로 물리적 환경과 개인의 건강상태 간의 관계를 분석하였다. 연구 결과, 인구밀도와 토지이용혼합도 등이 높으면 비만도는 낮으며 건강수준은 높은 것으로 보고 하였다.

그러나 최근 발표된 최태규·김흥순(2013)의 연구에서는 다른 결과를 제시하고 있다. 이 연구는 울산을 제외한 5개 광역시를 대상으로 토지이용 특성, 교통 환경 특성과 건강지표 간의 관계를 분석하였다. 단독주택 비율과 버스정류장의 수가 건강 수준에 긍정적인 영향을 미치는 것으로 나타났으나 전체적으로는 명확한 관계가 드러나지 않아 물리적 환경과 개인의 건강 수준 간의 관계가 거의 없음을 보고하였다. 그러나 서울을 제외한 지방 도시를 대상으로 연구를 진행한 상위 두 연구는 충분히 변수를 고려하지 못하였고 분석에 사용한 자료의 수가 상대적으로 적다는 한계를 가지고 있다.

선행연구의 결과를 다른 시각에서 바라본 김욱진·김태연(2013)의 연구는 개인의 건강상태를 물리적, 사회적 환경이라는 측정된 변수를 이용하여 예측하고 개인이 인지한 물리적, 사회적 환경이 개인의 건강과 실제적 환경 간의 관계를 매개하는 과정을 분석하였다. 서울서베이(2011)자료의 45,600개 표본을 이용하여 다수준 분석하였으며 개인의 건강은 객관적인 측정을 통한 실체적 지표이자 주관적인 특성이 있는 구성물이기 때문에 물리적인 환경은 개인의 건강에 영향을 미치지만, 당사자가 어떻게 느끼는지에 따라 그 효과가 다를 수 있음을 밝혔다. 결국, 지역적인 효과는 개인의 인식 속에서 상쇄되거나 소멸할 수 있는 것으로 보았다.

지금까지 정리한 연구들을 종합하면, 해외에서 물리적 환경과 개인의 신체활동 참여, 건강 수준 등과의 관계를 지속해서 연구해왔으며 우리나라에서도 일부 연구자를 중심으로 연구가 진행되는 추세에 있다. 그러나 아직은 해외에서 진행되어온 다양한 연구에 비해 부족한 실정이다. 국내 연구의 경우 근린의 물리적 환경과 신체활동이나 건강 수준과의 관계는 서로 일관성 있는 결과를 도출하고 있지 않다. 이러한 현상은 분석 자료의 유형과 분석단위, 분석변수, 분석 방법론의 차이에서 나타날 수 있으므로 다양한 측면에서 후속연구가 필요한 실정이다. 그리고 선행된 연구에서 주로 사용하는 토지이용 밀도와 혼합도, 가로망 연결도 등 단순한 물리적 변수에서 나아가 개인의 신체활동과 건강에 영향을 미칠 수 있는 다양한 변수를 종합적으로 고려할 필요가 있다. 예를 들어 개인의 인구 및 사회·경제적 지표뿐만 아니라 개인의 건강상태, 음주와 흡연 같은 생활습관, 아파트나 단독주택 같은 주거형태, 운동시설과의 접근성 등도 신체활동과 건강에 큰 영향을 미칠 수 있는 지표이다.

본 연구에서는 전술한 바와 같이 국민건강영양조사가 정기적으로 실시된 제4기 1차 년도(2007년)부터 제5기 3차연도(2012년) 까지 총 6년간의 자료를 누적하여 사용하였다. 상기 자료에서는 성별, 연령, 학력, 직업, 소득 수준, 거주하고 있는 주택의 유형 등 개인 및 가구 특성과 흡연, 음주, 신체활동, 자가건강도, 건강 수준, 삶의 질, 비만도 등의 항목을 독립변수와 종속변수로 재가공하여 사용하였다. 분석에서는 만 19세 이상 성인의 자료만 사용하였다.

근린단위의 물리적 환경은 국민건강영양조사(2007~2012), 서울시 새주소사업DB(2010), 서울 열린 데이터 광장(2011), BIZ-GIS XSDB(학교)(2011), 전국생활체육시설데이터(2012) 등을 이용하여 분석 자료를 구축하였다. 물리적 환경 변수는 ArcGIS 10.0을 이용하여 구축하였으며 필요한 경우 외부자료를 연결하여 사용하였다.

물리적 환경 변수 중 토지이용 및 공원과 관련된 변수는 새주소사업DB를 이용하여 구축하였다. GIS 프로그램을 이용하여 통계청을 통해 받은 행정동 자료와 건축물 자료를 공간조인(spatial join) 한 후, 건물의 주거, 상업, 업무 용도별 연면적을 계산하였다. 이후 계산된 건물의 연면적은 동별로 다시 합산하여 전체 동별, 용도별 면적자료를 구축하였다. 이 경우 건축물의 연면적을 계산하여 토지이용 특성을 분석한 이지은·이명훈(2011)의 기준²⁾을 준용하여 건물의 용도군을 설정하였다. 본 연구는 주거지 근린을 대상으로 하므로 공업적 토지이용을 제외한 6개의 토지이용 용도로 구분하여 토지이용 특성을 도출하였다. 공원변수의 경우「도시공원 및 녹지 등에 관한 법률」에 따라 근린공원, 어린이공원, 소공원의 세 가지 유형으로 구분하였고 GIS를 이용하여 각각의 면적을 측정하였다.

접근성 지표는 정확한 주소가 없는 자료의 한계로 인하여 행정동의 모든 주거용 건물(주택)을 점(point) 자료로 변환하여 각 주택에서 체육시설, 버스정류장, 지하철역 등 유형별 시설까지의 평균 거리를 계산하였다. 추가로 체육 활동을 할 수 있는 운동장의 대리변수로 학교 수를 사용하였다. 이와 더불어 행정동별 자전거도로 연장 등을 구축하여 사용하였다.

분석모형에 사용한 변수는 선행연구의 검토와 실증분석 이전의 기초적 분석을 통하여 선정하였다. 개인의 신체활동 참여와 건강 수준은 여러 가지 환경의 영향을 받는다는 가정을 통계적으로 확인하기 위하여 신체활동 지표로 일일 보행시간, 일일 중강도 신체활동 참여시간, 일일 고강도 신체활동 참여시간, 개인의 건강 수준지표로 EQ-5D³⁾, 삶의 질을 나타내는 EQ-VAS⁴⁾, 비만도(BMI) 총 6가지 종속변수를 사용하여 각각의 모형을 정립하였다. 이처럼 여러 가지 종속변수를 선정한 이유는 신체활동과 건강 수준이 보행이나 전반적인 건강 수준을 나타내는 EQ-5D와 같은 변수로만 설명하기에는 다차원적인 측면을 가지고 있기 때문이다. 특히 신체활동의 경우 대부분의 선행연구에서 저강도 신체활동인 보행만을 분석하는 경우가 많지만 저강도 신체활동과 중강도 또는 고강도 신체활동과 비교하여 분석하는 것이 필요하다.

신체활동과 건강 수준 종속변수와 더불어 종속변수에 영향을 미치는 구체적으로 설명변수를 선정하였다. 독립변수항목은 개인특성과 물리적 환경특성 두 가지로 나눌 수 있다. 많은 선행연구에서 다룬 개인 및 가구특성 변수는 대부분 명목척도로 이루어져 있어 별도로 <표 2>와 같이 정리하였다. 기본적인 특성으로 개인의 성별, 나이, 교육수준, 직업, 결혼 여부, 월간 음주 수준, 하루 흡연수준, 경제적 지표로 가구소득 등이 포함되었다. 특히 선행연구에서 잘 다루어지지 않았던 ‘자가건강인식수준’ 변수를 추가하였다. 이는 자신의 주관적인 건강인식 수준이 신체활동이나 건강 수준에 영향을 미칠 가능성이 크기 때문이다.

Table 2.

Individual and household characteristics

근린의 물리적 환경변수는 다시 개인이 거주하는 주택의 유형, 주변 근린의 특성을 나타내는 밀도와 토지이용, 공원 및 체육시설, 교통시설까지의 접근성, 일반적인 시설특성 등으로 나눌 수 있다. 주택유형의 경우 ‘국민건강영양조사’ 자료에 포함된 개인의 거주주택유형을 재분류한 것으로 단독주택, 다가구·다세대·연립주택, 아파트, 기타 주거형태의 네 가지 유형으로 구분하였다. 김승연 외 2인(2013)의 연구는 주거유형을 포함하는 주거환경이 건강 수준에 영향을 미치는 변수로 제시하고 있다.

토지이용변수는 면적을 사용하는 지표의 특성상 비교가 가능하도록 동별 면적으로 표준화하여 비율 변수로 사용하였다. 주거지역의 경우 단독주택 연면적, 다가구·다세대·연립주택 연면적, 아파트 연면적과 상업시설, 업무시설 각각의 건축물 연면적의 합계를 변수로 사용하였다. 토지이용혼합도는 선행연구에서 많이 사용한 토지이용 혼합도인 엔트로피지수(Frank and Pivo, 1995)를 활용하였다. 주거·상업·업무의 세 가지 용도가 잘 혼합된 경우는 1, 그렇지 못한 경우는 0에 가까워지는 지표로써 선행연구에서 보편적으로 사용하고 있는 변수이다.

공원지표의 경우 근린공원, 어린이공원, 소공원의 세 가지 공원 유형별 면적을 산출하였으며 토지이용지표와 마찬가지로 동별 면적으로 표준화를 하여 사용하였다. 본 연구와 유사한 대부분의 선행연구는 공원의 위계나 유형의 분류 없이 단순히 공원의 면적이나 접근성 변수를 사용하고 있다. 그러나 일정 규모 이상의 근린공원과 어린이공원은 규모와 사용대상에서 큰 차이가 있으므로 유형별로 나누어 살펴보는 것이 더 타당하다고 할 수 있다.

체육시설의 경우 민간 체육시설은 배제하고 운동장이나 축구장, 야구장 등의 공공 야외체육시설과 무료 또는 저렴한 가격으로 시설을 이용할 수 있는 국민체육센터 등을 포함하였다. 지하철역과 버스정류장까지의 거리지표는 사람들이 이용하고 있는 대중교통시설의 접근성 효과와 함께 지역의 특성을 함께 살펴볼 수 있는 지표로 분석에 사용하였다. 시설특성에서 학교 수는 학교운동장을 나타내는 대리변수로써 국민생활체육참여실태조사(2012)에서 사람들이 걷기나 조깅 등을 위해 가장 많이 찾는 곳 중 하나로 나타나 신체활동과 높은 관련성이 있을 것으로 판단되어 분석에 사용하였다. 자전거도로 연장의 경우 자전거를 활용한 신체활동 역시 활발히 일어나고 있으므로 자전거를 이용할 수 있는 물리적 환경인 자전거도로의 구축현황과 신체활동, 건강 간의 관계를 살펴볼 필요가 있다고 판단하여 변수로 사용하였다. 앞서 언급한 종속변수와 독립변수 중 근린의 물리적 환경 변수는 <표 3>과 같이 요약할 수 있다.

Table 3.

Descriptive statistics for dependent variables and neighborhood environment

본 연구에서는 통계적 분석을 위해 구축한 자료의 구조와 내용의 특성상 가장 적합하다 할 수 있는 다수준 분석 모형(Multilevel analysis model)을 사용하였다. 신체활동과 건강관련 선행연구에서도 자료의 구조에 따라 다수준 분석을 이용하여 연구의 결과를 도출하고 있다(이경환·안건혁, 2008; 최태규·김흥순, 2013; 성현곤 외 3인, 2014). 응답자 개인은 가구원으로서 가구의 특성을 공유하며 자신이 거주하고 있는 근린지역인 행정동의 물리적 환경 속성을 공유한다고 할 수 있다. 이는 같은 동네에 사는 사람의 경우 그 동네의 토지이용, 공원, 체육시설, 대중교통 시설 등의 영향을 받을 수 있기 때문이다. 따라서 분석 자료의 위계가 있는 경우 전통적인 다중회귀분석보다는 위계 수준의 차이를 고려할 수 있는 다수준 분석 모형이 더 적합하다(이희연·노승철, 2013). 또한, 행정동보다 넓은 기초자치단체(시, 군, 구)의 경우 행정동 보다 공간적인 범위 자체가 월등히 넓으므로 같은 시 또는 구에 거주하는 개인이라 하더라도 같은 공간적인 특성을 공유하는 관계가 굉장히 미약한 것으로 이해할 수 있다. 따라서 본 연구에서는 자치구 수준보다 더 작은 공간 단위로써 근린을 대표할 수 있는 행정동과 가구 수준 그리고 자료의 시간적 특성을 고려한 다수준 분석모형을 정립하여 통계 분석결과를 도출하였다.

신체활동 강도별로 신체활동 참여시간에 미치는 영향을 분석하기 위해 다수준 분석에서 독립변수를 제외한 상태에서 위계 수준만 고려하는 무제약 모형(unconditional model)을 수행하였다. 무제약 모형은 가구, 행정동, 시간(연도) 수준을 2수준, 3수준, 4수준 모두 고려하여 분석하였으며, 분석결과는 <표 4>에 제시하였다. 신체활동 강도별 모든 무제약 모형에 사용한 위계 수준은 통계적 유의성을 가지고 있는 것으로 나타났으며, 이는 다수준 분석에 적합한 것을 의미한다. 그리고 모형의 적합도를 나타내는 AIC 또는 BIC 값을 통해 4수준 모형이 2수준이나 3수준보다 모형의 적합도가 더 높게 나타나는 것을 확인하였다. 그리고 전체적으로 가구수준이 행정동이나 시간수준 보다는 ICC값이 높아 그룹 분산의 설명력이 높은 것으로 나타났다.

주거지 근린의 물리적 환경 특성이 저강도 신체활동(보행)과 중강도 및 고강도 신체활동 참여시간에 미치는 영향을 모든 독립변수를 고려한 임의절편모형(random intercept model)을 통해 분석하였다. 분석의 결과는 <표 5>와 같이 나타낼 수 있으며 결과의 해석은 다음과 같다. 우선, 전체적인 모형은 충분한 통계적 유의성을 확보한 것을 볼 수 있다. 일반적인 선형회귀분석과의 모형적합도 비교를 할 수 있는 LR test 결과도 다수준 모형을 사용하는 것이 더 효과적임을 나타내고 있다. 모형의 적합도(goodness of fit)를 나타내는 AIC, BIC 값 또한 무제약 모형보다 적합도가 더 증가한 것을 보여주고 있다. 분석모형 1의 분석 자료의 수준별 설명력을 나타내는 ICC(intra-class correlation) 값을 살펴보면, 시간(1.10%)과 시간-행정동(2.29%)의 설명력은 3% 이하로 나타났지만, 시간-행정동-가구 수준으로 구분할 경우 전체 6.62%의 설명력을 가지는 것으로 나타났다. 이는 시간-행정동-가구 수준만 가지고 응답자의 저강도 신체활동인 보행시간의 분산을 6.62% 설명할 수 있는 것을 의미한다.

Table 4.

Multi-level unconditional model test results for physical activity levels

Table 5.

Multi-level analysis for physical activity levels

개인 및 가구특성을 살펴보면 남성은 여성보다 저강도(보행), 중강도, 고강도 신체활동에 더 오래 참여하는 것으로 나타났다. 나이가 많을수록 보행시간이 더 긴 것으로 나타났고, 중학교 졸업자와 비교했을 때 상대적으로 높은 교육수준을 가진 사람일수록 저강도, 중강도 신체활동 참여시간이 모두 높은 것을 알 수 있다. 그리고 결혼한 사람은 결혼하지 않은 사람보다 보행시간과 고강도 신체활동 참여시간이 짧은 것으로 나타나 결혼한 사람은 대체로 신체활동참여에 소극적인 것을 알 수 있다.

또 다른 개인특성으로 자신이 인식한 건강수준을 ‘매우 나쁨’ 으로 응답한 사람과 비교할 때 ‘매우 좋음’으로 응답한 사람이 저강도(보행) 신체활동 참여시간이 높은 것으로 나타났다. 중강도와 고강도에서는 통계적으로 유의한 영향은 없는 것으로 나타났다. 경제적 지표로 가구소득은 중강도 및 고강도 신체활동 참여와 유의미한 관계에 있는 것으로 나타났는데 중강도 및 고강도 신체활동은 보행보다 소득수준과 높은 연관성이 있는 변수로 볼 수 있다. 이는 중, 고강도 신체활동이 보행활동보다 요구되는 시설, 장비, 이용료 등이 부가적으로 발생하기 때문이다. 또한, 보행참여 시간이 긴 사람은 중·고강도 신체활동에도 적극적으로 참여하고 있는 것으로 나타나 보행참여시간과 중강도, 고강도 신체활동 간의 상관성이 있음을 확인할 수 있다.

물리적 환경지표를 살펴보면, 고강도 신체활동의 경우 행정동의 사업체 밀도가 높은 지역에 사는 사람이 고강도 신체활동 참여시간이 긴 것으로 나타났다. 이는 고강도 신체활동이 직업적 활동을 포함하고 있으므로 높은 강도의 직업 활동을 많이 하는 사람이기 때문이라고 볼 수 있다. 또 다른 측면에서 업무시설의 밀도가 높은 도심, 부도심 등의 업무중심지에 거주하는 사람일 수 있다. 도심과 부도심의 외부환경은 달리기 등 고강도 신체활동에 적합하지 않지만, 실내 체육시설들이 자리 잡고 있어 나타날 수 있는 결과로 볼 수 있다.

저강도 신체활동인 보행활동에 영향을 미치는 물리적 환경 변수는 토지이용 변수 중에서 상업시설의 면적 비율이 높을수록, 업무시설 면적 비율이 낮을수록, 근린공원의 면적 비율이 높을수록, 그리고 학교시설에 대한 접근성이 좋을수록 보행활동 참여시간이 높은 것으로 나타났다. 근린의 상업시설은 보행을 유발하는 용도이기 때문에 적절한 결과로 판단되며, 업무중심지는 업무를 위한 주요시설과 함께 많은 편의 시설이 있어 다른 주거지역보다 접근성이 좋아 상대적으로 보행시간이 짧을 수 있다. 근린공원과 운동장이 있는 학교시설은 저강도 신체활동을 할 수 있는 장소이기 때문에 적절한 결과라고 볼 수 있다.

그러나 소규모 공원을 제외한 대부분의 토지이용 변수는 중강도 또는 고강도 신체활동과는 거의 관련이 없는 것으로 나타났다. 고강도 신체활동과 소공원 비율의 경우 음(-)의 관계를 보이는 것으로 나타났는데, 이는 고강도 신체활동이 많이 일어나는 도심 주거지, 아파트 밀집지역의 특성과 소공원을 포함한 자투리 녹지공간이 많은 중저층(다가구, 다세대, 연립주택 등) 주거지의 특성이 서로 다른 데서 나타난 결과로 보여진다.

주거·상업·업무 토지이용혼합도 변수 역시 저강도, 중강도, 고강도 신체활동과는 통계적으로 유의한 관계를 보이지 않았다. 이러한 결과는 최근 성현곤 외 3인(2014)의 연구결과와 유사하지만, 토지이용혼합도와 보행활동과의 유의미한 관계를 보인 이경환·안건혁(2007)의 연구결과와 차이가 있는 내용이다. 그리고 단독주택, 다가구·다세대·연립주택, 아파트와 같이 거주하고 있는 주택의 유형과도 신체활동 참여시간은 관계가 없는 것으로 나타났다.

중강도 신체활동의 경우 대부분의 물리적 환경 변수는 유의하지 않게 나타났지만, 자전거도로의 연장과는 양(+)의 관계에 놓인 것으로 나타났다. 이는 자전거도로 연장이 긴 지역에 사는 사람은 중강도 신체활동 참여시간이 긴 것을 의미한다. 일차적으로 자전거도로를 이용하는 자전거 이용자들이 주로 자전거도로를 이용하면서 신체활동을 참여했기 때문으로 나타난 결과로 받아들일 수 있다. 그러나 이차적으로 자전거도로는 자전거의 통행과 함께 주변부의 보행로, 작은 체육시설, 공원 등이 함께 조성된 경우가 많아 다양한 활동들이 수반될 수 있으므로 자전거도로와 중강도 신체활동과의 관계는 여러 의미로 해석할 수 있다. 이러한 결과는 서울시의 자전거도로 인프라 구축이 중강도 신체활동 증진에 효과가 있음을 의미한다.

이처럼 신체활동 유형별 분석결과는 물리적 환경과 신체활동과의 관계를 연구할 때, 저강도 신체활동인 보행활동과 운동이나 여가활동으로 접근하는 중강도, 고강도 신체활동은 별도로 분리하여 분석할 필요가 있음을 의미한다.

신체활동모형에 이어 물리적 환경과 개인의 건강 수준 간의 관계에 대해 분석하였다. 우선 건강상태변수에 대하여 독립변수를 제외한 위계 수준별 다수준 무제약 모형을 실행하였다(표 6 참조). 앞에서 설명한 신체활동 모형과 마찬가지로 대부분 모형에서 다수준 분석이 타당함을 보였다. ICC값을 비교해보면 가구 수준의 설명력이 가장 높고 그다음으로 행정동과 시간 수준으로 나타났다. 시간 수준의 경우 모형6 비만도 모형에서는 설명력이 매우 낮은 것으로 나타나 시간적인 차이는 크지 않은 것으로 나타났다. 그리고 2수준이나 3수준 모형보다는 시간-행정동-가구 수준을 고려한 4수준 모형의 적합도가 가장 높은 것으로 나타났다.

Table 6.

Multi-level unconditional model test results for health status variables

<표 7>은 건강상태를 나타내는 종속변수(EQ-5D, EQ-VAS, BMI)별로 독립변수를 포함한 4수준 분석모형의 결과이다. 신체활동 모형과 마찬가지로 위계 수준만 고려한 무제약 모형과 비교할 때 최종 다수준 모형의 적합도가 상대적으로 높은 것을 알 수 있다. <표 7>의 건강 수준 분석모형에서 L.R. Test는 OLS 회귀분석 모형보다 다수준 분석 모형이 적합함을 보였다. 예를 들어 비만도 모형의 경우 시간(0.05%)과 시간-행정동(0.5%) 수준의 설명력은 크지 않지만 시간-행정동-가구 수준을 연계할 경우 개인의 비만도를 12.1% 설명할 수 있음을 의미한다. 시간-행정동-가구 수준의 ICC 값은 모형4(EQ-5D)의 경우 2.48% 그리고 모형5(EQ-VAS)의 경우 5.51%로 나타났다.

Table 7.

Multi-level analysis for health status variables

개인 및 가구특성을 살펴보면 신체활동 모형과 비슷한 결과를 얻을 수 있었다. 남성은 여성보다 전체적인 건강상태(EQ-5D)와 삶의 질(EQ-VAS)은 높지만, 비만도(BMI)는 더 높은 것으로 나타났다. 이러한 결과는 남성이 여성보다 활동제한이 없고 삶의 질이 높다는 결과를 보고한 성현곤(2011)의 연구와 비슷한 결과로 볼 수 있다. 그리고 나이가 많을수록 활동 제한과 더불어 전반적인 건강 수준과 삶의 질이 낮아지며 비만도가 증가하는 것으로 나타났다.

교육수준의 경우 학력이 높은 사람일수록 전반적으로 건강 수준과 삶의 질이 높으며 비만도는 낮은 것으로 분석되었다. 직업적 특성에서는 비사무직에 비해 사무직과 기타 직업군의 사람들은 자신의 전체적인 건강상태를 낮게 평가하고 있음을 알 수 있다. 결혼한 사람은 그렇지 않은 사람보다 건강 수준과 삶의 질이 높고 더 비만한 것으로 나타났다.

평소 술을 자주 마시는 사람은 전반적인 건강상태에서 술을 마시지 않는 사람과 별 차이를 보이지 않았으며, 한 달에 한번 이하의 음주를 하는 사람의 삶의 질은 음주하지 않는 사람보다 높은 것으로 나타났다. 그러나 한 달에 한 번 이상의 음주는 통계적으로 유의하지 않았다. 비만도의 경우 일주일에 2~3회 음주를 하는 사람의 경우 음주를 하지 않는 사람보다 비만도가 높은 것으로 나타났다. 흡연의 경우 하루에 담배를 20개(1갑) 정도까지 많이 피우는 사람일수록 전반적인 건강상태나 삶의 질이 낮은 것을 볼 수 있으며, 특히 흡연은 삶의 질과 음(-)의 관계로 통계적으로 매우 유의하게 나타났다. 그러나 흡연과 비만도와의 관계는 통계적 유의성을 보이지 않았다.

또 다른 지표인 주관적인 건강인식수준이 높은 사람은 전반적인 건강 수준과 삶의 질이 높은 것으로 나타났으며, 소득수준도 양(+)의 관계를 갖고 있는 것으로 나타났다. 신체활동 특성으로 중강도를 제외한 저강도(보행) 또는 고강도 신체활동 참여시간이 긴 사람일수록 전반적인 건강 수준과 삶의 질이 높은 것으로 나타났다. 이러한 결과는 주관적인 건강인식 수준, 소득, 신체활동 등과 건강 수준이 일정한 연계가 있음을 의미한다.

그러나 음주, 흡연, 주관적 건강인식, 가구소득 등은 모형6의 비만도(BMI) 모형에서는 통계적 유의성이 매우 낮았다. 그리고 고강도 신체활동 참여시간이 긴 사람은 비만도가 높게 나타났다. 이는 고강도 신체활동 참여시간이 비만도를 높이는 것이 아니라, 비만도가 높은 사람이 스스로가 비만임을 인지하고 고강도 신체활동에 참여하고 있기 때문으로 판단된다.

결론적으로 개인 및 가구특성에서 개인의 기본적인 특성인 성별, 나이, 직업 등을 제어한 상태에서 주관적인 건강인식 수준이 높을수록, 소득이 높을수록, 그리고 저강도 신체활동인 보행시간이 길고, 고강도 신체활동 시간이 길수록 전반적인 건강 수준(EQ-5D)과 삶의 질(EQ-VAS)이 높은 것으로 나타났다.

다음으로 물리적 환경특성을 살펴보면 우선, 인구나 고용밀도 변수의 경우 건강 수준을 다룬 세 가지 모형에서 통계적으로 유의미한 관계가 나타나지 않았다. 이는 인구밀도가 삶의 질(EQ-VAS)에 양(+)의 영향을 미치는 것으로 보고한 성현곤(2011)의 연구와 인구밀도가 높을수록 비만도가 낮은 관계를 보인 이경환(2012)의 연구와도 차이가 있는 부분이다. 반면, 인구밀도나 주거밀도가 개인의 건강 수준에 영향이 없다고 보고한 최태규·김흥순(2013)의 연구와 비슷한 결론이다. 이러한 상반된 분석결과는 분석의 단위와 분석자료 그리고 분석 변수와 분석방법론에 따라 다를 수 있으므로 추가적인 연구가 필요할 것으로 판단된다.

주택 유형에서는 단독주택보다 아파트에 사는 사람일수록 비만도가 낮은 것으로 나타났다. 이는 주택 유형의 차이보다는 주택유형에 의한 경제적 특성과 관련이 있을 것으로 판단된다. 토지이용 측면에서는 다가구·다세대·연립주택의 밀도가 높은 지역에 사는 사람은 전반적인 건강 수준(EQ-5D)은 좋지 않으나, 비만도(BMI)는 낮은 것으로 나타났다. 업무시설의 밀도가 높은 지역에 사는 사람은 비만도가 더 높게 나타나는데, 이는 보행활동 참여시간이 짧은 이유와 연계해서 생각해볼 수 있다. 주거, 상업, 업무용도에 대한 엔트로피지수인 토지이용혼합도는 전반적인 건강 수준이나 삶의 질에는 유의미하게 나타나지 않았지만, 비만도(BMI)에 영향을 미치는 요인으로 나타났다. 토지이용혼합도가 높은 지역에 사는 사람일수록 비만도가 낮은 것으로 나타났는데, 이러한 결과는 토지이용혼합도와 비만과의 유의미한 관계를 설명한 성현곤(2011)과 이경환(2012)의 분석결과와 유사한 부분이다.

공원의 경우 근린공원의 면적이 넓은 지역에 사는 사람은 전반적인 건강상태(EQ-5D)가 좋은 것으로 나타났다. 근린공원이 보행과 신체활동 등에 영향을 미치고 이는 개인의 건강 수준에 영향을 미쳤을 가능성이 크다. 더불어 공원의 심미적, 환경적 요인이 개인의 건강상태에도 영향을 미쳤다고 이해할 수 있다. 그러나 어린이공원의 면적이 넓은 지역 거주자의 경우 근린공원과는 반대의 효과가 나타나는데, 이는 다가구·다세대·연립주택과 상당수의 어린이공원이 인접한 점 등을 고려할 때, 다가구·다세대·연립주택이 밀집한 지역 특성과 중첩된 결과로 볼 수 있다. 그러나 이 부분은 추가적인 연구를 통해 구체적인 관계를 분석할 필요가 있다.

마지막으로 시설특성으로 제시한 자전거도로 연장이 긴 지역에 사는 사람일수록 비만도가 낮은 것으로 분석되었다. 이는 자전거도로 연장이 유의미하지 않은 것으로 보고한 김은정·강민규(2011b)의 결과와는 다르나 선행연구에서는 분석단위가 시·군·구로 넓은 영역임을 고려했을 때 나타날 수 있는 차이로 볼 수 있다. 결과를 해석하면, 자전거도로가 비만도에 직접 영향을 미쳤다기보다는 자전거도로가 주거지 근처에 위치할 경우 이용 가능성이 증대되어 자전거도로를 이용하는 사람이 많아지므로 간접적인 영향을 미쳤다고 볼 수 있다. 더불어 자전거도로는 신체활동 모형에서 도출된 결과와 함께 살펴보았을 때, 그 주변에 구축된 여러 환경이 사람들에게 신체활동 장소로써 제 기능을 하였다고 볼 수 있다. 전체적으로 신체활동과 건강 수준 두 가지 모형에서 밀도, 주택유형, 토지이용 유형별 연면적, 공원 특성 등을 구분한 것은 의미가 있음을 확인할 수 있다.