有向无环图

作者

[编辑] 郑虎;

[审核] .

注记

Hiplot 网站

本页面为 Hiplot Directed Acyclic Graphs 插件的源码版本教程,您也可以使用 Hiplot 网站实现无代码绘图,更多信息请查看以下链接:

https://hiplot.cn/basic/ggdag?lang=zh_cn

有向无环图可视化。

环境配置

  • 系统: Cross-platform (Linux/MacOS/Windows)

  • 编程语言: R

  • 依赖包: ggdag

# 安装包
if (!requireNamespace("ggdag", quietly = TRUE)) {
  install.packages("ggdag")
}

# 加载包
library(ggdag)

数据准备

# 加载数据
tidy_ggdag <- dagify(
  y ~ x + z2 + w2 + w1,
  x ~ z1 + w1 + w2,
  z1 ~ w1 + v,
  z2 ~ w2 + v,
  w1 ~ ~w2, # bidirected path
  exposure = "x",
  outcome = "y") %>%
  tidy_dagitty()

# 查看数据
head(tidy_ggdag)
$data
# A tibble: 13 × 8
   name       x     y direction to      xend  yend circular
   <chr>  <dbl> <dbl> <fct>     <chr>  <dbl> <dbl> <lgl>   
 1 v     -1.13  -6.01 ->        z1     0.234 -5.62 FALSE   
 2 v     -1.13  -6.01 ->        z2    -1.52  -4.65 FALSE   
 3 w1     0.154 -4.63 ->        x      0.580 -4.16 FALSE   
 4 w1     0.154 -4.63 ->        y     -0.470 -3.58 FALSE   
 5 w1     0.154 -4.63 ->        z1     0.234 -5.62 FALSE   
 6 w1     0.154 -4.63 <->       w2    -0.642 -4.17 FALSE   
 7 w2    -0.642 -4.17 ->        x      0.580 -4.16 FALSE   
 8 w2    -0.642 -4.17 ->        y     -0.470 -3.58 FALSE   
 9 w2    -0.642 -4.17 ->        z2    -1.52  -4.65 FALSE   
10 x      0.580 -4.16 ->        y     -0.470 -3.58 FALSE   
11 y     -0.470 -3.58 <NA>      <NA>  NA     NA    FALSE   
12 z1     0.234 -5.62 ->        x      0.580 -4.16 FALSE   
13 z2    -1.52  -4.65 ->        y     -0.470 -3.58 FALSE   

$dag
dag {
v
w1
w2
x [exposure]
y [outcome]
z1
z2
v -> z1
v -> z2
w1 -> x
w1 -> y
w1 -> z1
w1 <-> w2
w2 -> x
w2 -> y
w2 -> z2
x -> y
z1 -> x
z2 -> y
}

可视化

# 有向无环图
p <- ggdag(tidy_ggdag) +
  theme_dag() 

p
图 1: 有向无环图