← 인과추론
PART 3 · 12 인과추론 · 05

성향점수매칭

Propensity Score Matching · 성향점수, 매칭, 처치효과 추정

성향점수로지스틱 회귀최근접 이웃 · 칼리퍼SMDATT공통지지

01처치받을 확률 (성향점수)참여자비참여자확률이 비슷한 사람끼리 짝짓기

이미 프로그램에 참여한 사람과 참여하지 않은 사람을 비교해야 하는 경우가 많습니다. 무작위 실험에서는 연구자가 처치집단과 통제집단을 직접 무작위로 구성할 수 있지만, 실제 연구에서는 그렇지 못한 경우가 많습니다. 직업훈련 프로그램에 참여한 사람들의 취업률이 더 높다면, 참여자는 원래 취업 의지가 강하거나 학력이 높았을 수도 있습니다. 이러한 차이를 고려하지 않으면 프로그램의 효과를 과대평가하거나 과소평가할 수 있습니다. 성향점수매칭(PSM)은 관찰된 특성을 바탕으로 처치를 받을 확률이 비슷한 대상끼리 짝을 지어 비교하는 방법입니다. 무작위 실험을 대신할 수 있는 완벽한 방법은 아니지만, 관찰자료에서 집단 간 차이를 줄이기 위한 대표적인 접근입니다.

SECTION 01

성향점수란 무엇인가?

성향점수(Propensity Score)는 대상자가 처치를 받을 확률을 의미하며, 처치 이전에 측정한 특성을 바탕으로 추정합니다.

e(Xi) = P(Ti = 1 | Xi)
Ti = 1: 대상 i가 처치를 받음 · Xi: 처치 이전에 관찰된 특성 · e(Xi): 해당 특성을 가진 대상이 처치를 받을 확률

직업훈련 프로그램의 참여 여부를 예측할 때 나이, 학력, 이전 경력, 이전 소득 등이 관련 변수일 수 있습니다. 이 변수들을 이용해 각 대상자의 성향점수를 추정하면, 참여자와 미참여자 중 처치를 받을 가능성이 비슷한 사람들을 찾을 수 있습니다.

참여자비참여자0.00.20.40.60.81.0추정 성향점수 · 모의 구직자 3,000명(참여 962명) · 참여자는 점수가 높은 쪽에, 비참여자는 낮은 쪽에 몰려 있다
그림 1. 참여자와 비참여자의 성향점수 분포

그림 1의 모의 자료처럼 학력이 높고 이전 소득이 낮은 사람이 훈련에 더 많이 참여했다면, 참여자의 성향점수는 높은 쪽에, 비참여자는 낮은 쪽에 몰립니다. 여기서 중요한 점은 성향점수가 처치의 효과를 나타내는 점수가 아니라는 것입니다. 성향점수는 처치를 받을 가능성을 나타내며, 처치가 결과에 미치는 효과는 매칭 이후 결과를 비교해 추정합니다.

SECTION 02

PSM은 어떻게 작동할까?

PSM은 일반적으로 다음과 같은 단계로 진행됩니다.

① 집단 정의참여자 · 미참여자② 성향점수 추정이전 변수로 참여 확률③ 매칭점수가 비슷한 대상끼리④ 균형 확인SMD로 특성 비교⑤ 결과 비교효과와 불확실성 추정
그림 2. PSM의 단계
  1. 처치집단과 통제집단을 정의합니다. 예: 직업훈련 프로그램 참여자와 미참여자
  2. 성향점수를 추정합니다. 학력, 나이, 경력 등 처치 이전 변수를 이용해 참여 확률을 계산합니다.
  3. 성향점수가 비슷한 대상을 매칭합니다. 참여자와 비참여자 중 특성이 비교적 유사한 대상을 짝짓습니다.
  4. 매칭 후 균형을 확인합니다. 두 집단의 관찰된 특성이 충분히 비슷해졌는지 검토합니다.
  5. 결과 차이를 비교합니다. 연구 목적에 맞는 처치효과를 추정하고 불확실성을 평가합니다.

PSM의 핵심은 단순히 점수가 비슷한 사람을 찾는 데 그치지 않습니다. 매칭이 끝난 뒤 실제로 두 집단의 관찰된 특성이 균형을 이루는지 확인해야 합니다.

SECTION 03

성향점수는 어떻게 계산할까?

성향점수는 일반적으로 로지스틱 회귀분석(범주형자료분석 06장)으로 추정하며, 다른 적절한 예측 모형도 사용할 수 있습니다. 직업훈련 참여 여부를 나이와 학력으로 설명한다면 다음과 같은 모형을 설정할 수 있습니다.

logit(e(Xi)) = β0 + β1나이i + β2학력i
logit: 확률을 로그 오즈로 변환하는 함수 · 추정된 값을 확률로 변환하면 성향점수

가령 다음과 같은 결과가 나왔다고 합시다.

대상실제 참여 여부추정 성향점수
A참여0.72
B미참여0.70
C참여0.35
D미참여0.33
E참여0.91
F미참여0.18
0.00.20.40.60.81.0A 0.72B 0.70C 0.35D 0.33E 0.91F 0.18E와 F: 차이가 커서 비교 어려움참여미참여추정 성향점수 (원고의 가상 자료)
그림 3. 성향점수로 짝짓기

그림 3처럼 A와 B, C와 D는 성향점수가 서로 비슷합니다. 반면 E와 F는 점수 차이가 커서 직접 비교하기 어려울 수 있습니다. 실제 분석에서는 점수 차이뿐 아니라 매칭 기준, 표본 손실, 공통지지 영역 및 공변량 균형을 함께 검토해야 합니다.

SECTION 04

매칭 방법에는 어떤 종류가 있을까?

성향점수를 계산한 뒤에는 어떤 기준으로 대상을 짝지을지 결정해야 합니다.

방법원리특징
최근접 이웃 매칭 (Nearest Neighbor)점수가 가장 가까운 대상을 매칭직관적이고 널리 사용됨
칼리퍼 매칭 (Caliper Matching)점수 차이가 일정 기준 이내인 대상만 매칭지나치게 다른 대상의 매칭을 제한
반경 매칭 (Radius Matching)정해진 범위 안의 비교 대상을 활용조건에 맞는 대상이 여러 명일 수 있음
커널 매칭 (Kernel Matching)여러 비교 대상에 가중치를 부여한 쌍만 선택하지 않고 정보를 활용
최근접 이웃참여자가장 가까운 한 명칼리퍼참여자허용 범위 안에서만커널참여자가까울수록 큰 가중치
그림 4. 매칭 방법의 차이

최근접 이웃 매칭은 성향점수가 0.72인 참여자와 0.70인 미참여자처럼 점수 차이가 가장 작은 대상을 선택합니다. 다만 가장 가까운 대상이라도 실제 점수 차이가 클 수 있으므로, 일정한 허용 범위를 두는 칼리퍼를 함께 적용하기도 합니다. 참여자 한 명당 비교 대상 한 명을 매칭하는 1:1 매칭과 여러 명을 매칭하는 1:k 매칭이 있으며, 비교 대상이 많아지면 활용하는 정보가 늘어날 수 있지만 모든 비교 대상의 특성이 충분히 비슷한 것은 아닐 수 있습니다. 매칭 비율은 표본의 분포와 분석 목적을 고려해 결정합니다.

SECTION 05

매칭 후에는 균형을 확인해야 한다

PSM을 수행했다고 해서 두 집단의 특성이 자동으로 같아지는 것은 아닙니다. 매칭 전에는 참여자와 미참여자의 나이, 학력, 경력 등이 크게 달랐더라도 매칭 후에는 이러한 차이가 줄어들어야 합니다. 이를 공변량 균형(Covariate Balance)이라고 하며, 일반적으로 표준화 평균 차이(SMD)를 활용해 확인합니다.

SMD = (X̄T − X̄C) / √[(sT² + sC²) / 2]
X̄T, X̄C: 처치 · 통제집단의 평균 · sT², sC²: 두 집단의 분산
0.1−0.1-0.5+0.5+1.0나이학력(년)경력(년)이전 소득학습 동기 (미측정)매칭 전매칭 후표준화 평균 차이(SMD) · 측정한 변수는 균형을 이뤘지만, 측정하지 않은 학습 동기는 그대로 남는다
그림 5. 매칭 전후의 공변량 균형

SMD의 절댓값이 작을수록 해당 변수의 집단 간 표준화 차이가 작다는 뜻입니다. 실무에서는 |SMD| < 0.1을 균형을 판단하는 참고 기준으로 자주 사용하지만, 절대적인 합격 기준은 아닙니다. 그림 5에서 학력의 SMD는 매칭 전 0.85에서 매칭 후 0.10로, 이전 소득은 -0.45에서 -0.02로 줄었습니다. 그러나 모형에 넣지 않은 학습 동기는 매칭 후에도 0.41로 남아 있습니다. 매칭은 관찰된 특성의 불균형만 줄이며, 미측정 교란이 없다는 사실까지 입증하지는 않습니다.

균형이 충분하지 않다면 성향점수 모형을 재검토하거나 매칭 기준을 바꾸는 방법을 고려할 수 있습니다. 다만 결과가 좋아 보일 때까지 변수를 임의로 바꾸는 방식은 피하고, 연구 설계와 판단 기준을 명확히 해야 합니다.

SECTION 06

PSM으로 어떤 효과를 추정할까?

PSM은 특히 처치집단의 평균 처치효과인 ATT를 추정하는 데 자주 사용됩니다(02장).

ATT = E[Y(1) − Y(0) | T = 1]

모의 자료에서 참여자 962명 가운데 칼리퍼 안에서 짝을 찾은 827명의 결과는 다음과 같습니다.

구분취업률
훈련 참여자 (매칭됨)62%
매칭된 미참여자49%
차이+12.5%p (95% 신뢰구간 7.6~17.0%p)
실제 ATT 11.2%p단순 비교 (전체)+20.1%pPSM (측정한 4개 변수)+12.5%pPSM (+ 학습 동기까지 측정했다면)+11.0%p취업률 차이 (참여자 − 비교 대상) · 모의 자료에서만 알 수 있는 실제 효과와 비교
그림 6. 매칭 방식에 따른 효과 추정

그림 6처럼 매칭 없이 전체 참여자와 미참여자를 비교하면 차이가 20.1%p로 부풀려지지만, 측정한 네 변수로 매칭하면 12.5%p로 실제 ATT(11.2%p)에 가까워집니다. 남은 차이는 측정하지 않은 학습 동기 때문이며, 학습 동기까지 측정해 매칭할 수 있었다면 11.0%p가 됩니다.

따라서 매칭 후 차이를 인과효과로 해석하려면 관찰된 공변량이 충분히 균형을 이루고, 중요한 교란변수가 누락되지 않았다는 등의 가정이 필요합니다. 또한 이 예처럼 매칭으로 일부 참여자(135명)가 제외되었다면 추정 결과는 원래의 전체 참여자와 정확히 같은 집단을 대상으로 하지 않을 수 있으므로, 분석 결과를 어떤 모집단에 적용할 수 있는지 함께 밝혀야 합니다.

SECTION 07

PSM의 주요 가정과 한계

PSM은 관찰자료의 비교 가능성을 높이는 방법이지만, 몇 가지 중요한 가정이 필요합니다.

조건부 교환가능성

처치 이전의 관찰된 변수들을 통제한 뒤에는 처치 여부와 잠재결과 사이에 중요한 교란이 남아 있지 않아야 합니다. 학력과 나이를 통제했더라도 학습 동기가 훈련 참여와 취업 결과에 모두 영향을 주는데 이를 측정하지 않았다면 교란이 남습니다.

공통지지와 중첩성

참여자와 미참여자 사이에 비교 가능한 성향점수 영역이 존재해야 합니다. 참여자는 0.9에, 미참여자는 0.1에 몰려 있다면 비슷한 대상끼리 매칭하기 어렵습니다.

일관성과 처치 정의

같은 처치로 분류된 대상들이 실질적으로 서로 다른 처치를 받는다면 효과 해석이 복잡해집니다. 프로그램의 내용과 강도가 대상마다 크게 다르다면 처치를 어떻게 정의할지 명확히 해야 합니다.

매칭으로 해결되지 않는 문제

PSM은 관찰된 변수에 기반하므로 관찰되지 않은 변수로 인한 교란을 자동으로 제거하지 않습니다. 매칭 과정에서 표본이 줄어들 수 있고, 매칭 방법과 모형 설정에 따라 결과가 달라질 수 있습니다.

공통지지가 충분함01겹치는 영역에서만 비교 가능공통지지가 부족함01비교할 만한 짝이 거의 없음
그림 7. 공통지지와 중첩성

따라서 PSM은 단순히 성향점수를 계산하고 매칭하는 것에서 끝나지 않습니다. 공변량 균형, 공통지지, 표본 제외, 추정 대상 및 불확실성을 함께 평가해야 합니다.

보고 문장 예시

나이, 학력, 경력, 이전 소득으로 추정한 성향점수를 이용해 직업훈련 참여자와 미참여자를 1:1 최근접 이웃 매칭(칼리퍼 0.2 SD)하였다. 참여자 962명 중 827명이 매칭되었으며, 매칭 후 모든 공변량의 |SMD|는 0.1 이하였다. 매칭된 표본에서 참여자의 취업률은 61.5%로 비교집단(49.1%)보다 12.5%p 높았다(95% 신뢰구간 7.6~17.0%p). 측정되지 않은 학습 동기 등의 교란 가능성은 남아 있다.

✓ 핵심 정리

처치받을 확률이 비슷한 사람끼리 짝지어 관찰된 차이를 줄입니다.


다음은 06 이중차분법(DID)입니다. 정책이나 프로그램의 시행 전후에 처치집단과 비교집단이 얼마나 다르게 변화했는지를 비교하는 방법을 살펴봅니다.

KEY TAKEAWAYS

PSM은 처치받을 확률이 비슷한 대상끼리 짝지어 관찰된 특성의 차이를 줄입니다. 모의 자료에서 단순 비교 +20.1%p는 매칭 후 +12.5%p로 실제 효과(+11.2%p)에 가까워졌지만, 측정하지 않은 학습 동기는 매칭으로 해결되지 않았습니다. SMD로 균형을, 성향점수 분포로 공통지지를 반드시 확인합니다.