Slide 73
Slide 73 text
分析の流れ 73
1.傾向スコアの算出
m.out1 <- matchit(t1 ~ x1 + x2 + x3 + x4 + x5 + x6,
data = data11, replace = TRUE, distance = "glm", method = "nearest")
m.data1 <- match.data(m.out1)
2.バランシングの確認
library(cobalt)
love.plot(m.out1, thresholds = .1)
3.傾向スコアを統制した共分散分析(ATTの推定)
model1 <- lm(y3 ~ t1, data = m.data1, weights =
weights)
model2 <- lm(y3 ~ t1 + x1 + x2 + x3 + x4 + x5 +
x6, data = m.data1, weights = weights)
summary(model1)
summary(model2)
※各共変量の群間差が|d|=0.1未満になることが望ましい
※マッチング後の分析においても、共変量を説明変数として投入した方が良い推定になる