Skip to menu

matching을 하려고 합니다.

권혜미 2016.12.02 17:45 Views : 428

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
474 대조군1과 치료군2,3,4 의 PSM 에 관해서 질문있습니다. [4] yessul 2021.11.11
473 생존분석에서 missing value가 있는 경우 coxph 모델링은 어떻게 하는지요.. [2] PipeDragoN 2017.01.18
472 ROC curve 에서 AUC 비교에 관해 질문입니다. [2] Seokin 2020.09.15
471 교수님 propensity score관련하여 질문드립니다. [1] tunneled_vision 2016.07.26
470 ROC 커브 분석 관련 질문 드립니다. [1] odd82 2018.06.30
469 메타분석에서의 개별연구의 질검증 [1] 야스오 2021.08.24
468 ROC 커브 분석에 대하여 문의 드립니다. 선택을 했는데도 버튼이 나타나지 않습니다. [1] DOC 2020.10.24
467 Table 1 만들기 질문 [2] jaychoi 2016.03.13
466 안녕하세요~교수님..통계결과 궁금한점이 있어서요... [1] Beauty 2016.08.22
465 PSM 매칭 파일 크기 [1] 스티브 2021.01.06
464 Propensity matching 시 모형적절성 평가에 대해 여쭙니다. [1] 이상화 2016.10.04
463 표만들기 오류 [1] 김밥친구 2019.03.03
462 R 문제가 도저히 안풀려요.. [1] 자하연미필 2021.10.09
» matching을 하려고 합니다. 권혜미 2016.12.02
460 데이타 업로드후 작업이 안됩니다. T.T [1] 유로정 2016.04.22
459 knit PDF 및 ztable에서 표 출력에 대하여 질문드립니다. [5] file zepiros15 2018.12.22
458 생존분석에서 생존표... [2] 파랑파랑하늘 2016.08.02
457 ROC curve 수정 관련해서 문의드립니다. [2] file hamdgogo 2020.05.01
456 성향점수분석에서 caliper 설정 [1] alse 2017.06.20
455 mytable(), mytable2()관련 질문이 드립니다. [1] file 닉네임은별명 2021.08.03