공간자기상관 (Spatial Autocorrelation)

한 줄 정의: 지리적으로 가까운 지점들의 값이 서로 얼마나 비슷하거나 다른 패턴을 보이는지를 수치화한 지표로, 대표적으로 모란의 I(Moran's I) 통계량이 쓰입니다.

쉽게 풀면

전국 시군구별 집값 데이터를 지도에 색칠한다고 해봅시다. 보통 부자 동네 옆에는 부자 동네가, 서민 동네 옆에는 서민 동네가 모여 있는 경향이 있습니다. 즉 "가까운 곳끼리 값이 비슷하다"는 공간적 패턴이 존재하는 것입니다. 이런 패턴을 그냥 우연으로 볼지, 통계적으로 유의한 군집(클러스터)으로 볼지를 판단하는 것이 공간자기상관 분석입니다. 모란의 I 값이 양수이고 크면 "비슷한 값끼리 뭉쳐 있다"는 뜻이고, 음수이면 "높은 값과 낮은 값이 체스판처럼 번갈아 나타난다"는 뜻이며, 0에 가까우면 공간적 패턴 없이 무작위로 흩어져 있다는 뜻입니다.

왜 중요한가

공간 데이터를 다루는 지리학, 역학, 도시계획, 생태학 등에서는 관측치가 지도상 위치를 갖기 때문에 "이웃한 지점은 서로 영향을 주고받는다"는 가정이 자연스럽게 성립합니다. 공간자기상관을 확인하지 않고 분석하면 실제로는 존재하는 지역 군집 효과를 놓치거나, 반대로 무작위 변동을 의미 있는 패턴으로 오인할 위험이 있습니다. 그래서 질병 확산 지도, 부동산 가격 모형, 환경오염 분포 연구 등 공간 데이터를 다루는 거의 모든 논문에서 분석 초반에 공간자기상관 검정을 거칩니다.

논문에서는 이렇게 쓰입니다

"시군구 단위 대기오염 농도에 대해 모란의 I 통계량을 산출한 결과 유의한 양의 공간자기상관이 확인되어 (I = 0.42, p < .001), 인접 지역 간 오염 수준이 군집을 이루는 것으로 나타났다."

이 문장은 대기오염 수치가 지도상에서 무작위로 흩어진 것이 아니라, 오염이 심한 지역끼리 인접해 모여 있는 공간적 패턴을 통계적으로 확인했다는 뜻입니다.

"읍면동별 전염병 신규 확진자 수에 대한 국소 모란의 I(local Moran's I) 분석 결과, 도심 인접 지역에서 고고(High-High) 군집이 유의하게 나타나 확산이 특정 권역에 집중되고 있음을 시사했다."

이 문장은 전체 지역이 아니라 개별 지역 단위에서 어디가 고위험 군집인지를 국소 지표로 세분화해 밝혔다는 뜻입니다.

"토지가격 회귀모형의 잔차에 대해 공간자기상관 검정을 수행한 결과 유의한 양의 상관이 남아 있어, 일반 최소제곱 대신 공간시차모형(spatial lag model)을 적용해 재추정하였다."

이 문장은 회귀분석 이후 잔차에까지 남아있는 공간적 패턴을 점검하고, 이를 반영해 더 적합한 모형으로 바꾸었다는 뜻입니다.

조금 더 깊게 보면

모란의 I는 전역(global) 지표로 전체 지도의 평균적인 군집 경향만 보여주는 반면, 국소 모란의 I(Local Indicators of Spatial Association, LISA)는 개별 지역이 고고(High-High), 저저(Low-Low), 고저(High-Low) 등 어떤 유형의 군집에 속하는지를 개별적으로 알려줍니다. 또한 공간자기상관 계산에는 어떤 지역을 "이웃"으로 볼지를 정하는 공간가중행렬(spatial weight matrix)이 필요하며, 인접 경계 기준으로 정의할 수도 있고 거리 기준으로 정의할 수도 있어 이 선택이 결과에 영향을 줄 수 있습니다.

주의할 점

공간자기상관이 존재하는 데이터에 일반적인 회귀분석을 그대로 적용하면 관측치들이 독립이라는 가정이 깨져 표준오차가 과소추정될 수 있습니다. 이는 시간축에서 관측치 간 상관을 다루는 자기상관 문제와 개념적으로 비슷하며, 공간회귀모형 등으로 보정해야 합니다.