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
1033 web-R에서 ROC 커브 문의드립니다. [1] file comet23 2021.10.06
1032 정회원 후원하였습니다. [4] secret Dr.Ahn 2021.10.03
1031 통계 오류 및 로지스틱 분석 문의드립니다. secret comet23 2021.10.01
1030 회귀분석 문의 드립니다. secret ulysses11 2021.09.26
1029 정회원 서버에서 disconnect 이됩니다. [2] file 귀차니즘킹 2021.09.23
1028 왜 오류인지만 알려주실 수 있나여.. [1] 몌뭔 2021.09.23
1027 multivaraite에서 confidence interval [2] zsdfa 2021.09.22
1026 복합표본분석 문의드립니다 [1] file Patrick 2021.09.20
1025 PSM관련해서 문의 드립니다. [1] secret realizedp 2021.09.13
1024 논문 revision 중 궁금한 점이 있어 질문합니다. [1] file 맹구모션 2021.09.12
1023 Paired T test 및 다중회귀 분석 문의드립니다. [1] secret ulysses11 2021.09.08
1022 회귀분석 문의 드립니다. [1] ulysses11 2021.09.08
1021 paired T test 문의드립니다. [1] ulysses11 2021.09.07
1020 Web R Table 생성에서 결과가 다르게 나옵니다. [1] file 케로스 2021.09.01
1019 로지스틱 회귀분석 - 에러코드 문의드립니다. secret wave 2021.08.31
1018 가설검증- cfa, sem [1] 통개 2021.08.30
1017 로지스틱 회귀분석 - 에러코드 문의드립니다. [1] wave 2021.08.30
1016 Odds ratio 값을 이용한 메타분석에 대해 문의드립니다. [2] endostat 2021.08.29
1015 mytable 연속형 변수 결과값 질문 [1] 1234 2021.08.28
1014 이런 그래프는 어떻게 그려야 하나요? -ggplot으로 가능한가요? file 보라도리 2021.08.27