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
1194 PSM 시 범주형 변수가 총 10가지 일 때 어떻게 돌려야 하는지요 [1] 긍정적인마음가짐 2023.05.06
1193 첫번째 행을 라벨로 사용 버튼 클릭시 에러 문의 [3] file 수담이 2023.05.04
1192 안녕하세요? web-R table 그리고 관련하여 문의드립니다. [4] mkw7654 2023.05.02
1191 메타분석 U3 그래프 관련 문의드립니다 secret 화이팅 2023.05.02
1190 교수님 안녕하세요~~ secret 메시12 2023.04.29
1189 연속형 변수인데 범주형으로 인지됩니다. [2] secret shararaku 2023.04.28
1188 통계마당 [2] secret ejo 2023.04.27
1187 moonBook 오류 관련 재문의드립니다. [1] secret holyhyun 2023.04.24
1186 moonBook 사용시 오류메시지로 문의드립니다. [1] holyhyun 2023.04.24
1185 에러 문의 [2] secret 화이팅 2023.04.19
1184 Cox regression 에 대한 문의입니다. [2] file saga 2023.04.18
1183 PSM 문의 [3] secret starwars 2023.04.13
1182 ggPS, ggPSMSummary function [6] 로렌쵸메디치 2023.04.12
1181 Survival 분석에 대해서 문의 드립니다. 정훈신 2023.04.05
1180 안녕하세요. 분석 툴에 대해 문의드립니다. [1] 스매직맨 2023.04.04
1179 드디어 publish 되었네요. [2] 정지윤 2023.04.04
1178 연구동향분석 시 논문 데이터 정제 관련 [5] vtlbtl 2023.03.31
1177 안녕하세요 교수님 에러 문의드립니다 [2] secret 소소 2023.03.26
1176 코딩값이 여러개인 변수에 대한 전체 p값을 알 수 있는 방법이 궁금합니다 쥬쥬즈주 2023.03.22
1175 PSM 매칭 버튼 보이지 않음. secret swpapa 2023.03.19