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
1174 Web R studio 초기화 후 문의드립니다. [2] secret 화모기 2016.02.23
1173 [정보] 구조방정식 모형 분석 관련 참고자료 소개 [2] 보석 2016.02.24
1172 그래프가 이상하게 나와서 어떻게 하면 좋을까요? [3] file osr33260216 2016.02.25
1171 데이터 업로드가 되지 않습니다.. [3] jenni 2016.02.28
1170 어제 강의 너무 잘 들었습니다 [5] 마술샘 2016.02.29
1169 lavaan의 공식 구글 토론실입니다 마술샘 2016.02.29
1168 [오류]web-r 구조방정식 [2] 마술샘 2016.02.29
1167 저 이 홈피 체험 1주일 지난것같은데... [3] osr33260216 2016.03.03
1166 이럴 땐 어떤 통계 기법을 쓰나요? (기초) [1] file 송이버섯 2016.03.03
1165 남북한별 인구를 나타내는 막대 그래프를 만들려는데 왜 바가 않나타나죠? [1] file osr33260216 2016.03.04
1164 제가 금시세 표를 만들었는데 왜 업로드 파일은 이상하게 나왔나요? [1] file osr33260216 2016.03.04
1163 Standard error 질문 드립니다. [7] file R까기 2016.03.05
1162 object 'res' not found [3] 민자 2016.03.06
1161 제가 선생님 말씀대로 금시세표를 csv쉼표분리로 변경했는데 이상하게되었습니다. [1] file osr33260216 2016.03.07
1160 인공지능, 인지 컴퓨팅, 머신 러닝, 딥 러닝 관련 IBM 자료 [1] file Ben 2016.03.09
1159 [RStudio Server] 사용 연장 문의 [2] KB 2016.03.09
1158 descriptive 표만들기 질문입니다. [4] 난둥 2016.03.10
1157 그래프 작성 문의드립니다. [4] file 도치맘 2016.03.10
1156 표만들기 관련 질문입니다. [1] tympa 2016.03.10
1155 머신러닝 레볼루션 행사안내 [1] mike 2016.03.10