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
1154 Multiple regression 기본 acs data로 돌려도 오류가 납니다. [3] 케로스 2023.01.22
1153 webrPSM 설치 문제 [1] hawk1227 2023.01.03
1152 정회원용서버 ExploPlots 오류 [1] 버블 2023.01.03
1151 생존분석에서 세군이상 시 사후 검정이 가능할지 여쭙습니다 [1] 정망고 2023.01.03
1150 gaze에서 show.stats Arguments가 작동하지 않습니다. [1] 김찬식 2022.12.31
1149 서버 연결이 자꾸 끊깁니다. [1] secret wave 2022.12.30
1148 PSM, full matching, cox proportional hazard model [4] file 주니 2022.12.28
1147 비밀번호 찾기 오류. [2] ezra006 2022.12.24
1146 그래프 만들던 중 오류 질문드립니다.. [1] file 헝힝헝 2022.12.14
1145 multipleROC 패키지 질문 [2] 강우성 2022.12.14
1144 웹에서 하는 R통계 (무료버전) 업로드 가능한 파일 용량이 궁금합니다. [2] 성열 2022.12.10
1143 데이터프레임 자료를 ggplot2를 이용해 꺾은선그래프 그리는 방법 [1] 다다다다빈 2022.12.06
1142 서버에서 table 만들기 에러 [1] 버블 2022.12.04
1141 mytable 문의드립니다. [1] file seyi0918 2022.11.30
1140 교수님 안녕하세요 moonbook 패키지 mycsv 오류가 나서 문의드립니다. [1] holyhyun 2022.11.15
1139 autoReg(fit, uni=TRUE)시, multivariable에 들어갈 변수들을 manual로 선택할 수 있을까요? [6] 김찬식 2022.11.14
1138 processR 패키지에서 Measurment model 정의가 가능한지요. chrik 2022.11.13
1137 Incidence rate per person-years 및 rate ratio에 대한 95% CI 산출 방법 문의 [1] file 꾸이뭉 2022.11.11
1136 서열로짓 시 평균한계효과 요다 2022.10.30
1135 서버문제로 ppt 다운로드가 안된다고 뜹니다. [5] liverim 2022.10.26