[Stata] optmatch2 명령어
통계 연습 2011. 10. 14. 07:38 |성향 점수(propensity score)를 만드는 데는 psmatch2 명령어가 가장 강력하지만, 만든 다음 짝지은 데이터셋을 분리해내려면 optmatch2 명령어가 편하다. 이 명령어는 isvar 명령어와 함께 동작하므로 미리 설치해두어야 한다.
치료군(treatm)을 성향 점수(propen)에 맞춰 1:1로 데이터셋을 추출하기 위한 명령어는 다음과 같다.
치료군과 비치료군에서 성향 점수가 가장 비슷한 27명을 골랐고, 짝지은 쌍에 대해 새로운 아이디(mid)를 만들었다. 이제 치료군과 비치료군에서 원하는 결과 변수에 대해 통계적 검정을 수행하면 된다.
치료군(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
. drop if mid==.
. save propennew.dta
치료군과 비치료군에서 성향 점수가 가장 비슷한 27명을 골랐고, 짝지은 쌍에 대해 새로운 아이디(mid)를 만들었다. 이제 치료군과 비치료군에서 원하는 결과 변수에 대해 통계적 검정을 수행하면 된다.
'통계 연습' 카테고리의 다른 글
| [Stata] batplot 명령어 (0) | 2011.12.26 |
|---|---|
| [Stata] oddsrisk 명령어 (0) | 2011.10.20 |
| [Stata] 다변수 로지스틱 분석에서 범주형 설명 변수의 빈도표 (0) | 2011.10.18 |
| [Stata] 표준화비 신뢰구간 계산 (0) | 2011.10.14 |
| [Stata] 위험 차(risk difference) 계산 (2) | 2011.10.14 |