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
1174 인공지능, 인지 컴퓨팅, 머신 러닝, 딥 러닝 관련 IBM 자료 [1] file Ben 2016.03.09
1173 [책소개] R in 24 Hours fermat39 2015.12.25
1172 ㅠㅠㅠ 초보 질문입니다 [1] 만고 2021.11.26
1171 파일 다운로드 실패 [2] 김밥친구 2017.09.28
1170 웹에서 하는 R통계 - R코드 [2] 김영진 2015.06.07
1169 R스튜디오 설치 [1] BK 2019.04.10
1168 그냥 넋두리.. [2] 스튜던트 2018.03.25
1167 web-R에서 사용되는 명령어를 볼 수 있는 방법이 있을까요? [1] 통계장님 2018.03.15
1166 구조방정식 관련 여쭙습니다. [9] fascinate23 2018.01.22
1165 스피어만 상관계수(Spearman correlation coefficient) [1] ㅎㅎㅎ 2015.10.15
1164 회귀 분석 에러 [1] file 김승욱 2015.03.26
1163 object 'res' not found [3] 민자 2016.03.06
1162 ztable package 사용과 관련하여 도움을 부탁드립니다. [1] IlKwon 2019.07.29
1161 메타분석 후 결과 forest plot이 위쪽과 아래쪽이 다 출력되지 않는 것에 대한 질문입니다. ^^ [2] file hanclan 2018.08.09
1160 median follow up and median overall survival [2] alse 2017.06.12
1159 변수에 가중치 적용한 붓스트랩 방법 또는 피키지 하고지비 2017.05.17
1158 Propensity score matching 결과 standardized difference 에 대한 질문 nskjh 2020.06.15
1157 Loop 사용 관련하여 문의드립니다. [3] 바슈메트 2020.03.09
1156 Trend test에 대한 질문 하가지만 드립니다. [5] file namuhanayeyo 2017.11.30
1155 메타분석에 관한 문의 드립니다. [1] 봄날 2018.03.08