Y : 암 발생유무 (0,1)
Group : 대조군 (0) , 실험군 (1)
공변량 : x1, x2, x3, x4, x5
1. rawdata
n (%)
| Y=0 (213명) | Y=1 (34명) |
Group=0 (195명) | 165 (84.62) | 30 (15.38) |
Group=1 (52명) | 48 (92.31) | 4 (7.69) |
2. matchit(Group~x1+x2+x3+x4+x5, data=d, ratio=1, method="nearest") 함수를 통해 1:1 매칭 한 결과:
n (%)
| Y=0 (51명) | Y=1 (9명) |
Group=0 (30명) | 25 (49.02) | 5 (55.56) |
Group=1 (30명) | 26 (50.98) | 4 (44.44) |
[질문]. 2번에 매칭된 데이터를 이용하여 Y~Group+x1+x2+x3+x4+x5 모델링을 했을때,
분석결과가 나온다고 해도 그 결과를 신뢰해도 될까요? (event n수가 너무 작아서요,,,,,,,)
Comment 2
-
cardiomoon
2018.01.30 15:44
-
wjdEoddl
2018.01.30 17:16
교수님 말씀처럼 full match로 해보려고 아래의 코드를 실행하였더니 다음과 같은 warning message가 나타났습니다.
matchit(Group~x1+x2+x3+x4+x5, data = d, ratio = 1, method="full");
Warning message:
In optmatch::fullmatch(d, ...) :
Without 'data' argument the order of the match is not guaranteed
to be the same as your original data.
제 개인적인 의견을 말씁드립니다.
1. 가능한 모든 데이터를 사용해서 다중로지스틱휘귀분석을 할 것을 권해드립니다.
2. 꼭 PSM을 하시려면 method를 fullmatch로 해보실 것을 권해드립니다.