# 安装包
if (!requireNamespace("ggdag", quietly = TRUE)) {
install.packages("ggdag")
}
# 加载包
library(ggdag)
有向无环图
注记
Hiplot 网站
本页面为 Hiplot Directed Acyclic Graphs
插件的源码版本教程,您也可以使用 Hiplot 网站实现无代码绘图,更多信息请查看以下链接:
有向无环图可视化。
环境配置
系统: Cross-platform (Linux/MacOS/Windows)
编程语言: R
依赖包:
ggdag
数据准备
# 加载数据
<- dagify(
tidy_ggdag ~ x + z2 + w2 + w1,
y ~ z1 + w1 + w2,
x ~ w1 + v,
z1 ~ w2 + v,
z2 ~ ~w2, # bidirected path
w1 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
}
可视化
# 有向无环图
<- ggdag(tidy_ggdag) +
p theme_dag()
p
