psm 平衡诊断的密度图

Question

我正在尝试制作一系列漂亮的密度图，以显示 psm 分析中变量的平衡。

full_match <- matchit(buyout_flag ~ tsale + sfincs_avg + logprox + tenure + age +
                     percent_ethwhite_origin + percent_poverty_origin + percent_hs_origin + percent_owner_origin + house_medval_origin + percap_origin,
                   method = "full", distance = "glm", link = "probit", data = hcad_floodp, caliper = 0.1)

基本的绘图功能会产生一个简单的系列（超过三个以上的绘图）。

plot(full_match, type = "density")

钴：bal.plot 产生更好的情节，但我只能弄清楚如何一次做一个。 我可以使用 bal.plot 为所有变量制作一系列或网格图吗？

bal.plot(full_match, var.name = "distance", which = "both")

Answer 1

使用bal.plot()一次不能做多个变量，但可以使用patchwork之类的东西将多个图合并为一个。 请参阅下面的示例。

data("lalonde", package = "MatchIt")
library(cobalt); library(MatchIt)
library(patchwork)

m <- matchit(treat ~ age + educ + race + married + 
               nodegree + re74 + re75, data = lalonde,
             method = "full", link = "probit")

bal.plot(m, "age", which = "both") /
  bal.plot(m, "educ", which = "both") /
  bal.plot(m, "re74", which = "both") +
  plot_layout(guides = "collect")

bal.plot(m, "nodegree", which = "both") /
  bal.plot(m, "race", which = "both") /
  bal.plot(m, "married", which = "both") +
  plot_layout(guides = "collect")

^{由reprex 包于 2022-06-06 创建 (v2.0.1)}

由于尺寸的原因，上面的图看起来不太好，但是您可以更改图形设备的尺寸以使其看起来更好。 为了正确收集图例，需要分别绘制分类协变量和连续协变量。 通过一些额外的编程，您可以真正自定义图看起来更漂亮，例如，通过使用ggplot2::labs()从每个图中删除标题，使用ggplot2::guides(fill = "none")删除图例等。