Skip to menu

matching을 하려고 합니다.

권혜미 2016.12.02 17:45 Views : 427

Rstudio을 사용하여 matching을 진행하려고 했습니다..



## Propensity score estimation (pscore object returned)
dat.nonrandom <- pscore(formula     = Alc ~ Age + Sex_F0M1 + R_BMI 
                        data        = dat,
                        family      = "binomial",
                        name.pscore = "ps")



## Actual Propensity score matching
res.ps.match <- ps.match(object             = dat.nonrandom,
                         control.matched.by = matched.by,
                         who.treated        = 1,                # value 1 is treated
                         name.match.index   = "match.index",
                         ratio              = 1,              
                         caliper            = "logit",
                         x                  = 0.2,
                         givenTmatchingC    = TRUE,
                         bestmatch.first    = TRUE,
                         setseed            = FALSE,
                         combine.output     = TRUE)
res.ps.match





res.ps.match를 시행하여보면

Matching data:
                                        
Number of treated obs.:              260
Number of matched treated obs.:      260
Number of untreated obs.:           1061
Number of matched untreated obs.:    260
Number of total matched obs.:        520
Number of not matched obs.:          801
Number of matching sets:             260
Number of incomplete matching sets:    0

이러한 결과가 나오는데

## Extract matched data
defMatched <-dat[unlist(res.ps.match[c("match.index","control.index")]), ]
## Construct a table
tabMatched <- CreateTableOne(vars = vars, factorVars = catvars, strata = "Alc", data = defMatched, test = TRUE)
## Show table with SMD
print(tabMatched, smd = TRUE, test=TRUE)

하면

                                       0                1                p      test SMD   
  n                                          434               86                            


로 treated 군이 86명으로 나옵니다ㅠ


저는 어디서 잘못된건지 모르겠습니다.ㅠ

위에 결과보면 260씩 매칭이된거 같은데- matched set을 만들때 잘못된건지 모르겠습니다....


저의 작은 지식으로는 이해가 안됩니다.ㅠ

조언부탁드립니다ㅜㅜ 감사합니다.


No. Subject Author Date
1113 ROC 비교 관하여 질문드립니다!(delong's test) 나현둥 2022.08.09
1112 안녕하세요교수님! VIF 관련하여 질문드립니다. file kms9057 2022.08.08
1111 cox 에서 AIC모델 구하는 방법이 있나요? [1] 쥬쥬즈주 2022.08.03
1110 PSM 중 server disconnected [3] joonho 2022.08.01
1109 서로 다른 set 에서의 kappa 값 비교 papahong 2022.07.30
1108 PSM에 대하여 [1] 김습습 2022.07.29
1107 엑셀의 index 기능을 R 에서 구현하려고 합니다 file emilien 2022.07.28
1106 survival 분석 중 PPT 다운로드 오류 secret 곰곰11 2022.07.26
1105 결과 ppt 다운로드 중 오류발생문제 [2] 곰곰11 2022.07.25
1104 Control/experimental Group이 있을 때 유의미한 결과 도출 방법 [3] 민됴 2022.07.25
1103 Cox uni multivaraite [2] zsdfa 2022.07.14
1102 생존분석 문의드립니다 [1] 최홍재 2022.07.13
1101 인용 관련하여 [1] kms9057 2022.07.12
1100 메타 회귀 분석 시 ggplot 출력에 관해 문의드립니다. endostat 2022.07.09
1099 생존곡선에서 색깔 변화 문의드립니다. file 케로스 2022.07.09
1098 후향적 임상시험에서의 검체수 산정 방법 질문 라라보라 2022.07.05
1097 Kaplan Meier survival curve 질문입니다 [2] cetla 2022.07.03
1096 mytalbe( )패키지 사용중 문의드립니다. 하은 2022.07.02
1095 심평원 R studio 내 lubridate 설치 오류 [4] 꾸이뭉 2022.06.27
1094 안녕하세요 r 코드 중에 도저히 해결이 안나는 문제가 있어서 여기에 문의 드립니다 [1] file 슈퍼맨 2022.06.24