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 변수를 변경하는데 질문이 있습니다. [2] 최00 2019.01.21
1153 SOS file UnuDIo 2019.09.02
1152 클리블랜드 점그래프 그리기 오류 관련 [5] Mefirst~! 2016.02.09
1151 인류 역사상 최고의 천재였다는 남자.jpg 소꿈노리 2018.05.25
1150 ANCOVA에 대해서 여쭤보고 싶습니다. [5] 고든 2015.09.13
1149 교수님 데이터 편집 질문드립니다. [1] file 백지백지 2019.08.02
1148 ggplot에서 y axis limit정하는 방법 여쭙니다. [2] fascinate23 2018.10.26
1147 matchIt 함수 사용시 질문 있습니다. [1] YJK 2019.02.23
1146 time ROC 에 대해 질문드립니다. file YH 2018.04.19
1145 r studio server 가 오늘(3/15) 만료된다고 합니다 [2] 스튜던트 2018.03.15
1144 web-r 에서 NRI, IDI 할수 있는지요,, [4] hawk1227 2016.04.07
1143 PSM 에서 standardized difference 구할 수 있는 방법 여쭙습니다. [1] shiny365 2018.01.29
1142 mytable에서 하는 통계법 질문 [1] 수수 2019.07.11
1141 [책 소개] Mastering RStudio – Develop, Communicate, and Collaborate with R file fermat39 2015.12.18
1140 2 X 2 table에서 categorical variable을 비교한 P 값 [2] namuhanayeyo 2016.05.01
1139 어제 강의 너무 잘 들었습니다 [5] 마술샘 2016.02.29
1138 ggplot 이용시 자꾸 오류가 납니다~ [1] 돌님 2017.10.23
1137 R markdown 관련 질문드립니다. [1] Lloyd 2016.04.06
1136 로지스틱 회귀분석에 관하여 질문드립니다. [2] tympa 2015.12.09
1135 nomogram [5] 횽횽 2018.05.16