R語言--傾向性評分matchit

rm(list = ls())
library("MatchIt")
library("tableone")
library('survival')
library('Zelig')
data(lalonde)
head(lalonde)
set.seed(1234)
m.out <- matchit(treat ~ educ + black + hispan, data = lalonde,
                 method = "nearest",ratio=4,caliper=0.001)
m.data <- match.data(m.out)
table(m.data$treat,m.data$educ)
summary(m.out)
plot(m.out,type="jitter")
plot(m.out,type="QQ")
plot(m.out,type="hist")
table1 <- CreateTableOne(vars = c('educ', 'black', 'hispan','married','nodegree'),
                         data = m.data,
                         strata = 'treat')
table1
table(m.data$treat,m.data$married)
d1 <- m.out$match.matrix   #行名稱（實驗組）對應值（配對的id）#