안녕하세요
저는 최근 R 공부를 시작하고 있습니다.
제가 가지고 있는 데이터로 PSM을 해보고 있는데
에러를 이해하지 못하겠어서 문의글을 작성하게 되었습니다.
PSM 이후 매칭된 환자의 수가 달라진 것이 에러인데요.
예를 들어, 제가 A 약물을 복용하는 환자군을 기준으로 매칭을 한다고 할때.
실제 데이터에서 A약물을 복용하는 환자는 139명이고, 복용하지 않는 환자는 10000명이라고 하고요.
이때, 1:1로 PSM을 한다고 하면 적어도 복용하는 환자는 139명, 복용하지 않는 환자는 130명 비슷하게 나와야하지 않나요?
근데 제가 시행했을 때는 (MATCHIT, nearest) (변수는 age, sex) (기준은 A약물복용력)
A 약물 복용하는 환자가 201명, 복용하지 않는 환자가 200여명정도가 나옵니다.
뭐가 잘못됐는지 아무리 찾아봐도 잘 모르겠어서요
이렇게 나오는 이유가 무엇인지,
어떻게 수정해야 제대로 나오는지 궁금합니다.
Comment 4
-
cardiomoon
2021.03.23 22:48
-
뽀오또
2021.03.24 18:55
네 제가 현재 가지고 있는 데이터를 몇개의 변수를 지우고 올려보겠습니다.
좀 변수가 많습니다.
저의 데이터는 C302라는 변수가 1인 환자를 기준으로(139명) matching 을 하고 싶습니다.
> match_C3 <- matchit(C302~ BX_AGE + A_SEX, method='nearest', data=C3_4, ratio=1)
> matched <- match.data(match_C3)
> dim(matched)
[1] 276 129
이렇게 나옵니다. 1:1 매칭인데 왜 2배 가까이 더 나오는지 모르겠습니다.
원래 데이터에서는 204:201 정도로 나왔었는데 139명보다 더 많이 나오는 이유도 모르겠습니다.....
뭐가 문제일까요...?
-
cardiomoon
2021.03.25 10:25
올려주신 데이터로 점검해보니 대조군과 치료군 모두 138명 씩 매칭됩니다. 위의 코드에서 summary(match_C3) 해보시면 양군간에 매칭된 수치를 확인하실 수 있습니다. 또한 dim(matched) 는 matched 데이터의 행갯수와 열갯수를 보여주는 명령어로 데이터가 모주 276행 및 129개의 열이라는 뜻입니다. table(matched$C302) 해보시면 대조군, 치료군 모두 138 명인 것을 알 수 있습니다.
-
뽀오또
2021.03.25 15:18
아.....넵~!!
제가 명령어 이해가 부족해서 벌어진 일이군요.
정말 감사합니다~!!!
No. | Subject | Author | Date |
---|---|---|---|
1194 | PSM 시 범주형 변수가 총 10가지 일 때 어떻게 돌려야 하는지요 [1] | 긍정적인마음가짐 | 2023.05.06 |
1193 | 첫번째 행을 라벨로 사용 버튼 클릭시 에러 문의 [3] | 수담이 | 2023.05.04 |
1192 | 안녕하세요? web-R table 그리고 관련하여 문의드립니다. [4] | mkw7654 | 2023.05.02 |
1191 | 메타분석 U3 그래프 관련 문의드립니다 | 화이팅 | 2023.05.02 |
1190 | 교수님 안녕하세요~~ | 메시12 | 2023.04.29 |
1189 | 연속형 변수인데 범주형으로 인지됩니다. [2] | shararaku | 2023.04.28 |
1188 | 통계마당 [2] | ejo | 2023.04.27 |
1187 | moonBook 오류 관련 재문의드립니다. [1] | holyhyun | 2023.04.24 |
1186 | moonBook 사용시 오류메시지로 문의드립니다. [1] | holyhyun | 2023.04.24 |
1185 | 에러 문의 [2] | 화이팅 | 2023.04.19 |
1184 | Cox regression 에 대한 문의입니다. [2] | saga | 2023.04.18 |
1183 | PSM 문의 [3] | starwars | 2023.04.13 |
1182 | ggPS, ggPSMSummary function [6] | 로렌쵸메디치 | 2023.04.12 |
1181 | Survival 분석에 대해서 문의 드립니다. | 정훈신 | 2023.04.05 |
1180 | 안녕하세요. 분석 툴에 대해 문의드립니다. [1] | 스매직맨 | 2023.04.04 |
1179 | 드디어 publish 되었네요. [2] | 정지윤 | 2023.04.04 |
1178 | 연구동향분석 시 논문 데이터 정제 관련 [5] | vtlbtl | 2023.03.31 |
1177 | 안녕하세요 교수님 에러 문의드립니다 [2] | 소소 | 2023.03.26 |
1176 | 코딩값이 여러개인 변수에 대한 전체 p값을 알 수 있는 방법이 궁금합니다 | 쥬쥬즈주 | 2023.03.22 |
1175 | PSM 매칭 버튼 보이지 않음. | swpapa | 2023.03.19 |
안녕하세요? 사용하신 데이터와 R코드를 올려주시면 한번 검토해보겠습니다