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
1154 산점도로 표현하고 싶은데요.. [2] file 난둥 2016.03.11
1153 Table 1 만들기 질문 [2] jaychoi 2016.03.13
1152 생존분석에서 설명변수 없이 생존곡선을 그리고 싶습니다. [4] sabqueen 2016.03.16
1151 그래프 그리기에 대한 질문입니다. [2] file 난둥 2016.03.18
1150 PSM와 c-statisticcs (AUC)에 관련해서 질문드립니다. 김정환 2016.03.22
1149 ANOVA, Kruskal-Wallis 검정 및 다중비교 [3] 소오름 2016.03.31
1148 자료전처리 문의드립니다 (두 명목변수를 조합해서 새 변수 생성) [2] 화모기 2016.03.31
1147 이 데이터의 비교를 하고 싶습니다. [3] file 소오름 2016.04.01
1146 메타분석에 대해 질문이 있습니다 [2] 카우보이 2016.04.01
1145 처음사용하는데 데이터 전처리가 어렵습니다. [1] namuhanayeyo 2016.04.02
1144 밑에서 질문 했던 자료입니다. [1] file 소오름 2016.04.03
1143 표만들기가 안되네요 ㅠㅠ [3] 피터팬79 2016.04.03
1142 자료 업뎃해서 다시 올립니다. [3] file 소오름 2016.04.03
1141 table 만들기에서, 나이 순으로 정렬시에 median 값과 범위를 지정할 수는 없는지요? [2] 혈액 2016.04.06
1140 R markdown 관련 질문드립니다. [1] Lloyd 2016.04.06
1139 ANOVA 사후분석은 어떻게 하나요? [2] drmin 2016.04.07
1138 web-r 에서 NRI, IDI 할수 있는지요,, [4] hawk1227 2016.04.07
1137 로지스틱 회귀 분석에서 반응(종속)변수에 결측치 [1] hawk1227 2016.04.07
1136 서버 접속과 문의 [2] file namuhanayeyo 2016.04.09
1135 파일 올려봅니다. [3] file namuhanayeyo 2016.04.11