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
1114 cox 다변량 분석 결측치 처리 관련 질문드립니다 [1] sungmas 2020.06.08
1113 자료 업로드할 때 뜨는 에러메시지 [1] 배깡 2017.02.11
1112 csv파일 업로드 안되는 문제 [2] file 신돌석 2016.04.25
1111 ggkm download [1] yongjuncha 2018.06.13
1110 아래에 올렸었는데요,, 5Y overall survival 구하는 법이요. [2] file 난둥 2016.05.26
1109 웹에서 하는 ggplot2에서.. [3] suein 2018.03.15
1108 로지스틱 회귀분석에서 독립변수가 범주형 번수일때 OR [3] jh100 2018.12.04
1107 단계구분도 문의드립니다. [5] onepeace00 2017.06.27
1106 로지스틱 - stepwise backward regression 외의 방법이 가능한지요? [2] 트루로맨스 2018.12.16
1105 plot.or그래프 그리는 방법 [4] 두치님 2016.07.04
1104 RPubs에 있는 ROC curve analysis by cardiomoon [2] hawk1227 2016.06.02
1103 안녕하세요 구조방정식(moderated mediation)에 대한 질문들이 있습니다. [2] file 통개 2021.08.19
1102 메타분석 결과 관련 질문드립니다. [2] seokzzn 2017.10.18
1101 결측값 처리와 관련하여 [3] 화모기 2015.10.07
1100 카플란 생존분석에서 로그랭크 검사가 안됩니다. nsstp 2019.05.13
1099 Cox 생존분석 다변량 인자 [2] alse 2018.08.28
1098 메타분석 시 사전-사후 상관계수? [1] 후니22 2018.01.03
1097 forest plot 가 다 보이지 않네요~ [8] 헤세드 2016.11.25
1096 mytable관련 문의드립니다 [4] 매직하마 2018.01.02
1095 표 작성시 질문 [1] simsim 2015.06.17