components() to partition graphs into weakly connected components or
strongly connected components of the directed part, including cyclic graphs
of class "UNKNOWN".remove_edges() now removes edges directly in Rust, avoiding graph table
reconstruction and redundant session cloning.nodes() and edges().mutate_caugi() no longer forces simple = TRUE on every target class. A
non-simple graph converted to "UNKNOWN", "ADMG" or "AG" keeps
simple = FALSE, so e.g. a bow (A %-->% B, A %<->% B) converts to a
valid ADMG instead of a simple graph with parallel edges that failed on
build() or any query. The converted graph is now built eagerly, so invalid
conversions fail immediately.
m_separated() now handles ADMG collider paths with multiple bidirected
edges, such as a -> b <-> c <-> d given {b, c}. This also fixes
is_valid_adjustment_admg() and all_adjustment_sets_admg(), which use the
same separation routine in the proper backdoor graph.
dag_from_pdag() no longer adds v-structures the input PDAG did not have.
Dor and Tarsi's condition (b) requires every undirected neighbour of the
candidate sink to be adjacent to all its other neighbours; the previous
implementation only required the undirected neighbours to be pairwise
adjacent, so on e.g. A --> B <-- C with B --- D it oriented D --> B
and created two v-structures. The algorithm now lives in Rust, like every
other graph algorithm in the package.
== and != for caugi objects now compare a symmetric edge as an
unordered pair, so caugi(A %---% B) == caugi(B %---% A) is TRUE as it
should be. This affected ---, <-> and o-o, and any symmetric edge added
with register_caugi_edge(); whether a glyph is symmetric is read from the
graph's own registry rather than hard-coded. Asymmetric glyphs remain
direction-sensitive. testthat's expect_equal() now agrees with ==.
Graphs now keep the edge definitions they were built with after
reset_caugi_registry() and re-registration. Previously ==, add_edges(),
mutate_caugi() and other verbs, write_caugi(), caugi_serialize() and
to_graphml() used the current global registry instead, so they could
compare edges wrongly, fail with Invalid edge code, or write the wrong
glyphs. Adding an edge type the graph was not built with is now an error.
enumerate_dags() to enumerate every DAG in the Markov equivalence
class of a PDAG, and count_dags() to return the MEC size without
materializing every DAG (#297).plot() now automatically bends edges around non-incident nodes that they
would otherwise pass straight through, so edges between collinear nodes (e.g.
within a tier) and edges crossing unrelated nodes stay visible. This is
controlled by edge_style$route (default TRUE); disable it with
edge_style = list(route = FALSE), or per edge type/edge via the usual
edge_style overrides.== and != methods for caugi objects so cg1 == cg2 returns a single
logical comparing graph content (class, nodes, edges, simple) rather than
session identity."CPDAG" graph class support across the constructor
(caugi()), coercion (as_caugi()), and class mutation (mutate_caugi()).
Construction validates the full CPDAG invariant (chordal chain components, an
acyclic component DAG, Meek closure, and strong arrow protection).
generate_graph(class = "CPDAG") now returns a graph with
@graph_class = "CPDAG" instead of "MPDAG", the precise label for the
essential graph of a Markov equivalence class. Predicates defined on PDAGs and
MPDAGs (is_pdag(), is_mpdag(), etc.) continue to accept CPDAGs unchanged.@graph_class = "MPDAG" instead of
"PDAG". This affects the result of meek_closure() and
generate_graph(class = "CPDAG"). Predicates and verbs defined on PDAGs
(is_pdag(), mutate_caugi(), etc.) continue to accept MPDAGs unchanged.aid() is now implemented natively in caugi's Rust backend, removing the
external gadjid dependency (and its pinned git revision and vendored
sources). Results are unchanged. The AID implementation
(src/rust/src/graph/aid.rs) is a derivative of gadjid and is licensed
MPL-2.0; the rest of the crate remains MIT. aid() now takes inputs of class
"DAG" or "CPDAG" (previously "DAG" or "PDAG").LICENSE.note
documents the MPL-2.0 component (the AID files derived from gadjid) and the
licenses of the vendored Rust crates, and their copyright holders are recorded
via cph/ctb roles in Authors@R.adjustment_set(type = "backdoor") now returns an inclusion-minimal backdoor
adjustment set, computed in linear time as a minimal d-separator in the proper
backdoor graph, rather than the full set of parents of the exposure.tools/benchmark/. The harness compares caugi to igraph, bnlearn,
dagitty, ggm, pcalg, pgmpy, and Tetrad on a n ∈ {100, 1000, 10000}
grid. The d-separation benchmark now uses a minimal d-separator computed via
minimal_separator() (previously used a backdoor adjustment set, which is
not in general a d-separating set).hd() returning results that depended on the order in which nodes were
declared. The Hamming distance now aligns nodes by name before comparing, so
logically identical graphs always give the same distance (#323).dag_from_pdag() failing with `from`, `edge`, `to` must be equal
length. when a sink had multiple undirected neighbors (#298).--o) edges to
be plotted as undirected edges.adjustment_set(type = "backdoor") returning an invalid (often empty)
set when a parent of the exposure lies on a backdoor path but is not an
ancestor of the outcome (#308). The result is now always a valid backdoor
adjustment set.is_mag() returning incorrect results for some ancestral graphs
(#309).
Adjacency was tested by binary-searching the concatenation of separately
sorted neighbor buckets, which is not globally sorted, so some adjacent
pairs were missed.tests/, examples/, trybuild fixtures) from the
vendored Rust dependencies so no vendored path exceeds 100 characters. This
silences pak's "very long paths" warning and avoids installation failures on
Windows without long-path support (#319).to_dot() and to_mermaid() (and write_dot()/write_mermaid())
silently converting partial --o and o-o edges into plain directed edges,
dropping the circle endpoints (#307).dag_from_pdag() is now named cg, matching the
convention used by the rest of the package. The previous name PDAG
continues to work as an alias but emits a deprecation warning.caugi_layout_circle() and a "circle" method for caugi_layout() that
places nodes evenly along the perimeter of a circle (#108).list_caugi_edges() function to list all available edge types."MPDAG" graph class support across constructor, class
mutation, and class resolution. class = "AUTO" now resolves Meek-closed
PDAGs to "MPDAG".exogenize() is now implemented in Rust for DAGs, which reduces overhead on larger graphs.normalize_latent_structure() is now implemented in Rust for DAGs for faster latent normalization workflows.minimal_d_separator() is renamed to minimal_separator() and now supports
ADMG and AG inputs (previously DAG-only), returning a minimal m-separator.
Implemented via the unified linear-time algorithm of van der Zander &
Liśkiewicz (UAI 2020). The old name minimal_d_separator() remains as a
deprecated alias.m_separated() on ADMGs: moralization now marries every pair in
pa(v) ∪ sp(v), not just pa(v). The old code missed moral edges from
bidirected co-parents and gave false positives (e.g. claimed Z ⊥ Y | X for
Z -> X -> Y, X <-> Y).is_valid_adjustment_admg() and all_adjustment_sets_admg() to verify
the GAC's m-separation condition in the proper backdoor graph rather than
via a per-neighbour decomposition. The old check trivially accepted neighbours
of X that were themselves in Z, so it falsely classified {C} as a valid
adjustment set in the M-bias ADMG C -> X, C <-> X, C -> Y, C <-> Y, X -> Y
(#277).normalize_latent_structure(), which normalizes the latent structure of a
DAG while preserving the marginal model over observed variables.minimal_d_separator(), which computes a minimal d-separator between
sets of nodes in a DAG, with support for mandatory inclusions and restrictions.posteriors() query function, which is the dual of anteriors(). It
returns all nodes reachable by following paths where every edge is either
undirected or directed away from the source node. For DAGs, posteriors()
equals descendants(). For PDAGs and AGs, it includes both descendants and
nodes reachable via undirected edges.ancestors(), anteriors(), descendants(), and posteriors().
This can be set globally with caugi_options() or locally with the
open = TRUE/FALSE argument. The default remains open = TRUE.is_mpdag() query to check whether a PDAG is closed under Meek's
orientation rules (R1-R4), and meek_closure() to orient all implied edges
until Meek closure.caugi_options() now supports nested key drilling: multiple unnamed arguments
traverse nested options (e.g., caugi_options("plot", "tier_style", "fill")).simple, graph_class, nodes, edges) are sourced from the session.@.state, @name_index_map, @built,
and @ptr now warn on access and return NULL.build and state in
caugi() now warn and are ignored.inplace parameter in verb functions (add_edges(), remove_edges(),
set_edges(), add_nodes(), remove_nodes()). All graph modifications now use
copy-on-write semantics for consistency with R conventions. The parameter is
deprecated and ignored with a warning.all.equal and compare_proxy methods for caugi objects to support
graph-content comparison in tests.asp parameter to plot() for controlling the aspect ratio. When asp = 1,
the plot respects equal units on both axes, preserving the layout
coordinates. Works like base R's asp parameter (y/x aspect ratio) (#195).pdag_to_dag() function that generates a random DAG consistent with a
given CPDAG/PDAG structure if possible (#201).plot() to use incorrect layout if node names were not
in the same order as in the graph object (#198).set_edges() so that it correctly replaces symmetric edges in
simple graphs.all in districts() has been deprecated. Use districts()
without arguments to get all districts.-->),
bidirected (<->), and undirected (---) edges while satisfying ancestral
graph constraints. New functions: is_ag(), is_mag().mode argument to neighbors()/neighbours() to filter neighbors by edge
direction or type ("all", "in", "out", "undirected", "bidirected",
"partial"). This is a structural query, and not a semantic query!neighbors() now supports class = "UNKNOWN" graphs, including mode-based
filtering.simulate_data() that enables simulation from DAGs using SEMs. Standard
linear Gaussian SEMs are defaults, but more importantly custom SEMs are
available."AUTO" parameter for class in caugi objects. This automatically
picks the graph class in order DAG, UG, PDAG, ADMG, AG.is_ag() and is_mag() and m-separation
for AGs.exogenize() function that exogenizes variables for any graph type.
Current implementation is written in R, but it is so simple that it might be
preferable over a Rust implementation. This might be changed later.latent_project() function that does latent projection from DAGs to
ADMGs.write_caugi(), read_caugi(), caugi_serialize(), and
caugi_deserialize(). The format is a versioned JSON schema that captures
graph structure, class, and optional metadata (comments and tags).plot() method for visualizing graphs using various layout algorithms.
The plot is rendered using grid graphics and returns a caugi_plot object
that can be customized with node_style, edge_style, and label_style
arguments. The plot() method accepts layouts as strings, functions, or
pre-computed data.frames.caugi_layout() function to compute node coordinates for graph
visualization.caugi_layout_sugiyama(),
caugi_layout_fruchterman_reingold(), caugi_layout_kamada_kawai(),
caugi_layout_bipartite(), and caugi_layout_tiered(). Each function
provides an API for its specific algorithm.to_dot() and write_dot() functions for exporting caugi graphs to DOT
(graphviz) format. The resulting object is a new S7 class, caugi_export,
which has a knit_print() method for rendering DOT graphs in R Markdown and
Quarto documents.to_graphml, write_graphml,
read_graphml, to_mermaid, write_mermaid, and read_mermaid+ and |
for horizontal arrangement, / for vertical stacking. Compositions can be
nested arbitrarily (e.g., (p1 + p2) / p3).caugi_options() function for setting global defaults for plot
appearance, including composition spacing and default styles for nodes,
edges, labels, and titles.caugi_default_options() function to query or reset to package default options.is_caugi() validation calls internally.caugi_layout_tiered() now returns a tier column and orientation attribute in the layout data.frame, allowing plot() to automatically use tier information without requiring the tiers argument to be passed again.index_name_map parameter from internal .cg_state() function.is_cpdag function that returns TRUE on non-complete PDAGs.shd returning positive values for equivalent graphs given in shuffled order.is_admg(), spouses(), districts(),
and m_separated() (generalization of d-separation for graphs with bidirected
edges).is_valid_adjustment_admg() and all_adjustment_sets_admg() implementing
the Generalized Adjustment Criterion.mutate_caugi() function that allows conversion from one graph type to
another.caugi objects.edges_df argument to caugi() for easier construction from
existing data frames containing the columns from, edge, and to.as_adjacency() and as_igraph() to support bidirected edges.as_caugi() documentation to include "ADMG" as a valid class type for
conversion.caugi.org/.CONTRIBUTING.md.README rewrite.lockBinding and unlockBinding in the package to
silence R CMD check notes.mutate_caugi function that allows conversion from one graph type to another.edges_df argument to caugi for easier construction from existing data frames containing the columns from, edge, and to.caugi in a package vignette to use new conversion functionality.CONTRIBUTING.md to github.as_caugi.dplyr and tibble.data.table.data.tables.caugi.Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.