발간사 7
Chapter 1. 성과연구 15
개요 16
1.1. 성과연구의 정의 16
1.2. 성과연구에서의 연구설계 17
1.3. 성과연구에서의 바이어스 18
참고문헌 20
Chapter 2. 교란요인 21
개요 22
2.1. 교란요인의 정의 22
2.2. 교란요인 선정 방법 24
2.2.1. 인과성 도식도의 활용 24
2.2.2. 통계적 방법의 활용 26
2.3. 교란요인을 고려한 설계 및 분석 방법 32
2.4. 효과변경인자 33
2.4.1. 효과변경작용 33
2.4.2. 효과변경작용을 확인하는 방법 35
2.4.3. 효과변경인자에 대한 분석방법 38
2.4.4. 효과변경인자를 고려한 연구분석 사례 38
참고문헌 41
Chapter 3. 제한 43
개요 44
3.1. 제한방법 44
3.2. 예제 46
참고문헌 49
Chapter 4. 매칭 50
개요 51
4.1. 개요 51
4.2. 매칭변수의 선정 52
4.3. 매칭 방법 53
4.3.1. 빈도매칭 53
4.3.2. 개별매칭 54
4.4. 매칭 후 공변량 균형 검토 60
4.4.1. 공변량 균형 검토 60
4.4.2. 최적 균형 61
참고문헌 62
Chapter 5. 층화 64
개요 65
5.1. 층화변수의 범주 선정 65
5.2. 단일 요약 추정치 계산 66
5.2.1. 표준화 추정치 계산 66
5.2.2. 통합 추정치 계산 67
참고문헌 76
Chapter 6. 다변수모형 77
개요 78
6.1. 다변수모형의 종류 78
6.2. 선형회귀모형 80
6.2.1. 기본개념 81
6.2.2. 치료효과의 추정 82
6.2.3. 모형구축 82
6.2.4. 모형진단 83
6.2.5. 자료진단 87
6.2.6. 유의성 검정 및 모형적합도 87
6.2.7. 프로그램 88
6.2.8. 예제 89
6.3. 일반화 선형모형 94
6.3.1. 로지스틱 회귀분석 95
6.3.2. 포아송 회귀분석 98
6.3.3. 유의성 검정 및 모형적합도 100
6.3.4. 프로그램 103
6.3.5. 예제 104
6.4. 콕스비례위험모형 109
6.4.1. 기본개념 110
6.4.2. 치료효과추정 113
6.4.3. 가정 및 모형적합도 검토 113
6.4.4. 프로그램 115
6.4.5. 예제 115
6.4.6. 콕스비례위험모형을 응용한 기타방법 118
6.5. 조건부 로지스틱 회귀모형 119
6.5.1. 기본개념 121
6.5.2. 예제 122
6.6. 일반화 추정방정식 123
6.6.1. 기본개념 125
6.6.2. 예제 127
참고문헌 132
Chapter 7. 성향점수 134
개요 135
7.1. 성향점수 개요 135
7.2. 이론적 배경 136
7.2.1. 반사실적 인과성 136
7.2.2. 평균치료효과 138
7.2.3. 성향점수의 정의 및 성질 140
7.3. 성향점수의 추정 142
7.3.1. 성향점수 추정모형에 포함되어야 할 변수선정방법 142
7.3.2. 성향점수 추정방법 144
7.4. 성향점수 적용방법의 선정 146
7.5. 성향점수를 이용한 제한 148
7.6. 성향점수를 이용한 매칭 149
7.6.1. 매칭 방법 150
7.6.2. 균형 확인방법 153
7.6.3. 매칭 후 분석 시 고려할 점 154
7.6.4. 성향점수 매칭 적용사례 156
7.7. 성향점수를 이용한 층화분석 160
7.8. 성향점수를 공변량으로 보정하는 방법 164
7.9. 성향점수 가중치 방법 168
7.10. 성향점수 방법의 비교 174
7.11. 성향점수 방법의 확장 178
7.11.1. 복합표본설계 자료 178
7.11.2. 치료량 및 치료방법이 세 개 이상인 경우 183
7.11.3. 베이지안 방법 187
7.11.4. 시간의존 치료와 시간의존 공변량이 있는 경우 188
7.11.5. 결측치가 있는 경우 189
7.12. 성향점수 적용을 위한 소프트웨어 191
7.13. 예제 192
7.13.1. 연구대상자들의 기저 특성 193
7.13.2. 성향점수의 추정 195
7.13.3. 성향점수를 이용한 매칭 196
7.13.4. 성향점수를 이용한 층화분석 201
7.13.5. 공변량 보정방법 206
7.13.6. 성향점수 가중치 방법 207
7.13.7. 치료효과 추정 209
참고문헌 212
Chapter 8. 질병위험점수 218
개요 219
참고문헌 223
Chapter 9. 주변구조모형 224
개요 225
참고문헌 230
Chapter 10. 매칭 추정량 231
개요 232
10.1. 단순 매칭 추정량 232
10.2. 바이어스 보정 매칭 추정량 235
10.3. 분산추정 236
10.4. 예제 236
참고문헌 238
[부록 A] R-MatchIt 패키지 240
[부록 B] Stata-psmatch2 함수 249
[부록 C] 매칭관련 프로그램 252
[부록 D] 우심도자술 자료 253
판권기 3
표 2-1. 교란요인 선정방법-결과변수(사망과의 관련성 확인) 27
표 2-2. 교란요인 선정방법-노출변수(우심도자술: RHC)와의 관련성 확인 28
표 2-3. 교란요인 선정방법-우도비 검정 29
표 2-4. 교란요인 선정방법-추정치 변화 검토 30
표 2-5. 측정된 교란요인의 통제방법 선정 33
표 2-6. 효과변경작용이 있는 경우 분석 결과표 예시 39
표 2-7. 효과변경인자를 고려한 위험도분석 예시 40
표 3-1. 제한된 코호트에서 스타틴 사용과 1년 사망률과의 관련성 48
표 4-1. 매칭의 장점과 단점 52
표 5-1. 번째 층의 빈도 자료 66
표 5-2. 표준화 방법에 의한 두 군의 치료효과 비교 67
표 5-3. 직접 추정법에서의 통합 추정치 68
표 5-4. 맨틀-핸첼 방법에서의 통합 추정치 69
표 5-5. 층화분석 결과 71
표 5-6. 편두통 완화여부 자료 72
표 5-7. SAS 프로그램 및 결과 73
표 5-8. R 프로그램 및 결과 74
표 5-9. Stata 프로그램 및 결과 75
표 6-1. 결과변수에 따른 다변수 모형 79
표 6-2. 선형회귀모형의 가정 진단 83
표 6-3. 다중공선성의 정도를 나타내는 지표 86
표 6-4. 영향력 큰 관측치 통계량 87
표 6-5. SAS의 회귀분석 프로그램 88
표 6-6. R과 Stata의 회귀분석 프로그램 89
표 6-7. 자료의 형태 90
표 6-8. 모수추정 결과 92
표 6-9. Cook''s D 통계량과 DFFITS 93
표 6-10. 모형적합 결과 94
표 6-11. 일반화 선형모형 94
표 6-12. Modified Park test 프로그램 95
표 6-13. 회귀계수에 따른 확률분포 95
표 6-14. 로지스틱 회귀분석 결과제시 예 96
표 6-15. 포아송 회귀분석 결과제시 예 98
표 6-16. 유의성 검정통계량 100
표 6-17. SAS, R, Stata의 일반화 선형모형 프로그램 103
표 6-18. 자료의 형태 104
표 6-19. 유의성 검정통계량 105
표 6-20. 로지스틱 회귀분석 결과 106
표 6-21. 오즈비 106
표 6-22. 호스머-램쇼 검정통계량 106
표 6-23. 자료의 형태 107
표 6-24. 포아송 회귀분석 결과 108
표 6-25. 음이항 회귀모형 적합도 109
표 6-26. 음이항 회귀분석 결과 109
표 6-27. 콕스비례위험모형 결과제시 110
표 6-28. 비례가정 평가방법 114
표 6-29. SAS, R, Stata의 콕스비례위험모형 프로그램 115
표 6-30. 자료의 형태 116
표 6-31. 콕스비례위험모형의 결과 117
표 6-32. SAS, R, Stata의 층화 콕스비례위험모형 프로그램 119
표 6-33. 조건부 로지스틱 회귀분석 결과제시 예 1 120
표 6-34. S AS, R, Stata의 조건부 로지스틱 회귀모형 프로그램 121
표 6-35. 조 건부 로지스틱 회귀모형 예제 자료 122
표 6-36. 조건부 로지스틱 회귀분석 결과 123
표 6-37. GEE 분석 결과 제시 예 125
표 6-38. 작업상관행렬 126
표 6-39. SAS, R, Stata를 이용한 GEE 프로그램 127
표 6-40. 일반화 추정방정식 예제 자료(횡적자료) 127
표 6-41. 일반화 추정방정식 예제 자료(종적자료) 128
표 6-42. 일반화 추정방정식 분석 결과 129
표 7-1. 치료효과의 종류 139
표 7-2. 성향점수 매칭 방법에 따른 바이어스와 분산사이에 균형 152
표 7-3. 전체 대상자들에 대한 두 군의 특성 비교 157
표 7-4. 매칭된 대상자들에 대한 두 군의 특성 비교 158
표 7-5. 전체대상자들에 대한 치료효과 추정 결과 159
표 7-6. 매칭 대상자들에 대한 치료효과 추정 결과 160
표 7-7. 전체 대상자에서의 스타틴 사용군과 비사용군의 기본특성 163
표 7-8. 각 층에서의 스타틴 사용군과 비사용군의 기본특성 163
표 7-9. 스타틴 치료효과 추정 결과 164
표 7-10. 약제 사용군과 비사용군의 기본특성 167
표 7-11. 약제의 치료효과 추정 결과 167
표 7-12. 전체 성향점수를 이용한 치료효과 추정 결과 171
표 7-13. 성향점수가 0.05보다 큰 자료를 이용한 치료효과 추정 결과 172
표 7-14. 성향점수 적용방법의 특징 175
표 7-15. 분석 방법에 따른 스타틴 사용의 치료효과 176
표 7-16. 대상자의 일반적 특징 181
표 7-17. 보정된 표본가중치 적용 후 대상자의 일반적 특징 182
표 7-18. 우울감이 자살생각에 미치는 영향 182
표 7-19. 골다공증 치료효과 추정 결과 186
표 7-20. 소프트웨어별 프로그램 186
표 7-21. 베이지안 성향점수방법과 성향점수방법의 결과 비교 188
표 7-22. 성향점수 방법에서 결측치를 고려한 분석방법 189
표 7-23. 성향점수 적용법에 따른 소프트웨어별 패키지 및 매크로 191
표 7-24. 소프트웨어별 사용가능 옵션 192
표 7-25. RHS 원시자료 형태 193
표 7-26. 대상자의 특성 194
표 7-27. 성향점수 매칭 표본 197
표 7-28. 성향점수 매칭 후 대상자의 특성 198
표 7-29. 성향점수 층화 후 대상자의 특성 203
표 7-30. 성향점수 가중치 방법 적용 후 대상자의 특성 208
표 7-31. 성향점수 분석방법에 따른 치료효과 추정 210
표 8-1. 질병위험점수(DRS)를 이용한 치료효과 222
표 9-1. 주변구조모형 분석 결과 제시 예 229
표 10-1. 한 개의 공변량을 사용한 잠재적 결과 추정공식 233
표 10-2. 공변량이 한 개인 경우의 단순 매칭 추정 예시 233
표 10-3. 공변량이 한 개인 경우의 단순 매칭 추정 예시 234
표 10-4. 여러 가지 평균치료효과의 추정공식 234
표 10-5. 매칭 추정량을 위한 Stata 프로그램 및 결과 237
그림 2-1. 무작위배정 비교임상시험에서의 치료효과 예시 22
그림 2-2. 교란요인의 예시 23
그림 2-3. 인과성 도식도의 예시 25
그림 2-4. 교란요인의 종류와 통제방법 32
그림 2-5. 효과변경인자를 파악하기 위한 도식적 탐색방법(2차원): 효과변경인자에 따른 효과차이 35
그림 2-6. 효과변경인자를 파악하기 위한 도식적 탐색방법(3차원): 효과변경인자에 따른 효과차이 37
그림 3-1. 모집단 제한 단계 47
그림 3-2. 신환자 정의 단계 예시(한국보건의료연구원, 2011) 48
그림 4-1. 개별매칭 방법의 구분 54
그림 6-1. 다변수모형 분석 절차 80
그림 6-2. 선형회귀모형 결과 제시 81
그림 6-3. 산점도 91
그림 6-4. 잔차산점도 및 정규확률그림 결과 92
그림 6-5. 시간의존형 자료 예시 111
그림 6-6. 시간의존형 자료 예시 117
그림 6-7. 조건부 로지스틱 회귀분석 결과제시 예 2 120
그림 7-1. 성향점수 방법의 적용 단계 136
그림 7-2. 반사실적 조건 137
그림 7-3. 반사실성의 대체 138
그림 7-4. 치료군과 비치료군의 성향점수의 분포 제시 예 146
그림 7-5. 성향점수 분포 148
그림 7-6. 층별 스타틴 사용군과 비사용군의 성향점수 상자그림 162
그림 7-7. 공변량 균형 확인을 위한 분위수 회귀분석 결과 제시 예 166
그림 7-8. t-PA 치료군과 비치료군의 성향점수 분포 170
그림 7-9. 로지스틱 회귀모형에서 치료효과의 바이어스 178
그림 7-10. 캠페인의 노출정도에 따라 마약사용에 대한 생각 184
그림 7-11. 성향점수 분포 195
그림 7-12. 매칭 전후 성향점수 분포(좌: 매칭전, 우: 매칭후) 199
그림 7-13. 매칭 전후 성향점수의 누적 분포 함수(좌: 매칭전, 우: 매칭후) 199
그림 7-14. 매칭 전후 공변량 분포확인-cardiovascular diagnosis 200
그림 7-15. 매칭 전후 공변량 분포확인-neurological diagnosis 200
그림 7-16. 매칭 전후 공변량 분포확인-DNR status on day 1 200
그림 7-17. 매칭 전후 공변량 분포확인-APACHE scor 201
그림 7-18. 매칭 전후 공변량 분포확인-Mean blood pressure 201
그림 7-19. 층화 후 성향점수 분포확인 205
그림 7-20 층화 후 공변량 분포확인(APACHE score) 206
그림 7-21. 공변량 균형확인 207
그림 7-22. 성향점수 매칭 전후 카플란-마이어 생존곡선 210
그림 7-23. 성향점수 가중치 방법에 따른 카플란-마이어 생존곡선 211
그림 7-24. 성향점수 층화 후 카플란-마이어 생존곡선 211
그림 8-1. 질병위험점수 분포 222
그림 9-1. 고정-치료연구의 인과성 도식도 225
그림 9-2. 시간에 따라 변하는 치료와 교란요인의 인과성 도식도(A: 아스피린 사용, L: 심근경색, Y: 심장사망) 227
부록 표 1. summary() 명령어에서 제공하는 척도 241
부록 표 2. match.data() 옵션 244
부록 표 3. matchit 함수 옵션 245
부록 표 4. distance 방법 246
부록 표 5. 매칭 방법에 따른 추가옵션 247
부록 표 6. summary() 옵션 248
부록 표 7. plot() 옵션 248
부록 표 8. 매칭 방법에 다른 psmatch2 명령어 249
부록 표 9. psmatch2 함수 옵션 250
부록 표 10. psmatch2 함수에서 제공하는 변수 251
부록 표 11. R 매칭 관련 패키지 링크 252
부록 표 12. Stata 매칭 관련 패키지 링크 252
부록 표 13. SAS 매칭 관련 매크로 링크 252
부록 표 14. 우심도자술 자료 변수목록 253
부록 그림 1. summary() 결과 중 전체자료와 매칭된 자료에 대한 통계량 242
부록 그림 2. summary() 결과 중 매칭 전후 표본수 242
부록 그림 3. 지터그림 243
부록 그림 4. Q-Q그림과 히스토그램 243
