DAG objectUsing the node information contained in the DAG object this function plots the corresponding DAG in a quick and convenient way. Some options to customize the plot are available, but it may be advisable to use other packages made explicitly to visualize DAGs instead if those do not meet the users needs.
# S3 method for DAG
plot(x, layout="nicely", node_size=0.2,
node_names=NULL, node_color="black",
node_fill="red", node_linewidth=0.5,
node_linetype="solid", node_alpha=1,
node_text_color="black", node_text_alpha=1,
node_text_size=8, node_text_family="sans",
node_text_fontface="bold", arrow_color="black",
arrow_linetype="solid", arrow_linewidth=1,
arrow_alpha=1, arrow_head_size=0.3,
arrow_head_unit="cm", arrow_type="closed",
arrow_node_dist=0.03, gg_theme=ggplot2::theme_void(),
include_td_nodes=TRUE, mark_td_nodes=TRUE,
...)Returns a standard ggplot2 object.
A DAG object created using the empty_dag function with nodes added to it using the + syntax. See empty_dag or node for more details.
A single character string specifying the layout of the plot. This internally calls the layout_ function of the igraph package, which offers a great variety of ways to layout the nodes of a graph. Defaults to "nicely". Some other options are: "as_star", "as_tree", "in_circle", "on_sphere", "randomly" and many more. For more details see ?layout_.
Either a single positive number or a numeric vector with one entry per node in the DAG, specifying the radius of the circles used to draw the nodes. If a single number is supplied, all nodes will be the same size (default).
A character vector with one entry for each node in the DAG specifying names that should be used for in the nodes or NULL (default). If NULL, the node names that were set during the creation of the DAG object will be used as names.
A single character string specifying the color of the outline of the node circles.
A single character string specifying the color with which the nodes are filled. Ignored if time-varying nodes are present and both include_td_nodes and mark_td_nodes are set to TRUE.
A single number specifying the width of the outline of the node circles.
A single character string specifying the linetype of the outline of the node circles.
A single number between 0 and 1 specifying the transparency level of the nodes.
A single character string specifying the color of the text inside the node circles.
A single number between 0 and 1 specifying the transparency level of the text inside the node circles.
A single number specifying the size of the text inside of the node circles.
A single character string specifying the family of the text inside the node circles.
A single character string specifying the fontface of the text inside the node circles.
A single character string specifying the color of the arrows between the nodes.
A single character string specifying the linetype of the arrows.
A single number specifying the width of the arrows.
A single number between 0 and 1 specifying the transparency level of the arrows.
A single number specifying the size of the arrow heads. The unit for this size parameter can be changed using the arrow_head_unit argument.
A single character string specifying the unit of the arrow_head_size argument.
Either "open" or "closed", which controls the type of head the arrows should have. See ?arrow.
A single positive number specifying the distance between nodes and the arrows. By setting this to values greater than 0 the arrows will not touch the node circles, leaving a bit of space instead.
A ggplot2 theme. By default this is set to theme_void, to get rid off everything but the plotted nodes (e.g. everything about the axis and the background). Might be useful to change this to something else when searching for good parameters of the number arguments of this function.
Whether to include time-varying nodes added to the dag using the node_td function or not. If one node is both specified as a time-fixed and time-varying node, it's parents in both calls will be pooled and it will be considered a time-varying node if this argument is TRUE. It will, however, also show up if it's argument is FALSE. In this case however, only the parents of that node in the standard node call will be considered.
Whether to distinguish time-varying and time-fixed nodes by fill color. If TRUE, the color will be set automatically using the standard ggplot2 palette, ignoring the color specified in node_fill. Ignored if include_td_nodes=FALSE or if there are no time-varying variables.
Further arguments passed to the layout function specified by the argument of the same name.
Robin Denz
This function uses the igraph package to find a suitable layout for the plot and then uses the ggplot2 package in conjunction with the geom_circle function of the ggforce package to plot the directed acyclic graph defined by a DAG object. Since it returns a ggplot object, the user may use any standard ggplot2 syntax to augment the plot or to save it using the ggsave function.
Note that there are multiple great packages specifically designed to plot directed acyclic graphs, such as the igraph package. This function is not meant to be a competitor to those packages. The functionality offered here is rather limited. It is designed to produce decent plots for small DAGs which are easy to create. If this function is not enough to create an adequate plot, users can use the dag2matrix function to obtain an adjacency matrix from the DAG object and directly use this matrix and the igraph package (or similar ones) to get much better plots.
If the DAG supplied to this function contains time-varying variables, the resulting plot may contain cycles or even bi-directional arrows, depending on the DAG. The reason for that is, that the time-dimension is not shown in the plot. Note also that even though, technically, every time-varying node has itself as a parent, no arrows showing this dependence will be added to the plot.
empty_dag, node, node_td
library(simDAG)
# 2 root nodes, 1 child node
dag <- empty_dag() +
node("age", type="rnorm", mean=50, sd=4) +
node("sex", type="rbernoulli", p=0.5) +
node("smoking", type="binomial", parents=c("sex", "age"), betas=c(1.1, 0.4),
intercept=-2)
if (requireNamespace("ggplot2") & requireNamespace("ggforce")) {
library(ggplot2)
library(igraph)
library(ggforce)
plot(dag)
# get plot using the igraph package instead
g1 <- as.igraph(dag)
plot(g1)
# plot with a time-varying node
dag <- dag +
node_td("lottery", type="time_to_event", parents=c("age", "smoking"))
plot(dag)
}
Run the code above in your browser using DataLab