성향 점수(propensity score)를 만드는 데는 psmatch2 명령어가 가장 강력하지만, 만든 다음 짝지은 데이터셋을 분리해내려면 optmatch2 명령어가 편하다. 이 명령어는 isvar 명령어와 함께 동작하므로 미리 설치해두어야 한다.

치료군(treatm)을 성향 점수(propen)에 맞춰 1:1로 데이터셋을 추출하기 위한 명령어는 다음과 같다.

. optmatch2 treatm propen, minc(1) maxc(1) gen(mid)

----+--- 1 ---+--- 2  ---+--- 3 ---+--- 4 ---+--- 5
...........................

Matched 27 treated subjects and 27 untreated subjects
Sum of distances within matched sets = 0.0265667

. sort mid treatm

. list id treatm propen mid in 1/10

     +--------------------------------+
     |  id   treatm      propen   mid |
     |--------------------------------|
  1. | 251        0   .11095736     1 |
  2. |  55        1   .11095736     1 |
  3. |  89        0   .05489074     2 |
  4. |  86        1   .05489074     2 |
  5. | 136        0   .20299611     3 |
     |--------------------------------|
  6. |  36        1   .20299611     3 |
  7. |  58        0   .19385163     4 |
  8. | 297        1   .19615819     4 |
  9. |  52        0   .08239614     5 |
 10. | 186        1   .08174527     5 |
     +--------------------------------+

. drop if mid==.

. save propennew.dta 

치료군과 비치료군에서 성향 점수가 가장 비슷한 27명을 골랐고, 짝지은 쌍에 대해 새로운 아이디(mid)를 만들었다. 이제 치료군과 비치료군에서 원하는 결과 변수에 대해 통계적 검정을 수행하면 된다.

Posted by cyberdoc
: