回帰分析

作者

宋財泫(関西大学)

View slides in full screen

my_dag <- dagify(
  Y ~ D + X1 + X2,
  X1 ~ Z1 + Z2,
  D ~ Z1 + Z3,
  X2 ~ Z2 + Z3,
  Z1 ~ W,
  Z2 ~ W,
  Z3 ~ W,
  exposure = "D",
  outcome  = "Y",
  coords = list(x = c("Y" = 2, 
                      "X1" = 1, "D" = 2, "X2" = 3,
                      "Z1" = 1, "Z2" = 2, "Z3" = 3,
                      "W" = 2),
                y = c("Y" = 4, 
                      "X1" = 3, "D" = 3, "X2" = 3,
                      "Z1" = 2, "Z2" = 2, "Z3" = 2,
                      "W" = 1))
)

ggdag(my_dag) + theme_dag()

my_dag |>
  adjustmentSets()
{ X1, X2 }
{ X1, Z2, Z3 }
{ X2, Z1, Z2 }
{ Z1, Z3 }
my_dag |>
  ggdag_adjustment_set(shadow = TRUE) +
  theme_dag()