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
1213 Error in f(x, na.rm = TRUE) : 사용되지 않은 인자 (na.rm = TRUE) [2] twogama2 2017.09.18
1212 paired T test 문의드립니다. [1] ulysses11 2021.09.07
1211 [문의] Web RStudio안에서 R 패키지 설치관련 [1] 보석 2016.02.19
1210 Web-R 에서 logistic regression 시 p value 를 볼 수 있는 방법이 있나요? [2] 전공의 2017.08.15
1209 rstudio-server 패키지 설치오류 file 하리보는맛있어 2019.05.14
1208 Cut off value [2] 변변변수 2017.06.19
1207 안녕하세요 psm은 정회원만 이용 가능한가요? [1] Arete 2020.08.01
1206 문의드릴 사항이 있습니다. [7] 정지윤 2021.02.18
1205 mytable로 만든표를 워드파일로 변환 하는 방법 [3] 정암 2017.06.02
1204 안녕하세요 mytable 함수 관련하여 질문드립니다. [1] 지옥소년 2019.07.26
1203 데이타 merge 방법 [5] nipponica 2017.02.18
1202 로지스틱 회귀분석 김밥친구 2018.10.04
1201 범주형변수와 비정규분포 연속형 변수 간 관계를 볼때 이변량 상관관계 분석이 가능한지요? [1] 아라리 2018.05.30
1200 교수님, 다중공선성 질문드립니다. [2] selteen 2017.11.16
1199 2 way RM ANOVA 문의드립니다. file FAdavid 2021.01.19
1198 web-r에서 psm을 수행하였는데 레퍼런스는 어떻게 표시해야 하나요? [1] 스티브 2020.11.24
1197 Rstudio server 서버 오류 관련해서 문의드립니다. [2] file 노상래 2015.07.22
1196 엑셀파일 data 올릴시 maximal upload size exceeded [1] 딸기공주 2020.07.01
1195 정회원용 서버가 접속이 안됩니다. [1] 마늘맨 2018.02.16
1194 생존분석에서 multivariate analysis에 대해 여쭙니다. [3] joon 2016.12.23