Struct CausalDiscovery 

pub struct CausalDiscovery {
    pub n_vars: usize,
    pub skeleton: Vec<Vec<bool>>,
    pub directed: Vec<Vec<bool>>,
    pub sep_sets: HashMap<(usize, usize), Vec<usize>>,
    pub alpha: f64,
}

Causal discovery via the PC algorithm.

The PC algorithm learns the structure of a DAG from conditional independence tests on observational data. It produces a Completed Partially Directed Acyclic Graph (CPDAG) representing the Markov equivalence class.

Fields

n_vars: usize

Number of variables.

skeleton: Vec<Vec<bool>>

Adjacency matrix of the skeleton (undirected).

directed: Vec<Vec<bool>>

Directed adjacency: directed[i][j] = true means i → j is oriented.

sep_sets: HashMap<(usize, usize), Vec<usize>>

Separation sets: sep_sets[(i,j)] = the conditioning set that d-separates i and j.

alpha: f64

Significance threshold for independence tests.

Implementations

impl CausalDiscovery

pub fn new(n_vars: usize, alpha: f64) -> Self

Create a new PC algorithm runner.

Arguments

• n_vars — number of observed variables
• alpha — significance level for conditional independence tests

pub fn learn_skeleton(&mut self, data: &[Vec<f64>])

Learn the skeleton from a data matrix using partial correlation tests.

Arguments

• data — data matrix, shape [n_obs][n_vars]

pub fn orient_v_structures(&mut self)

Orient v-structures (colliders) in the skeleton.

For each unshielded triple i — k — j (i and j not adjacent), if k is not in sep(i,j), orient i→k←j.

pub fn apply_meek_rules(&mut self)

Apply Meek’s orientation rules to complete the CPDAG.

Rules R1–R3 propagate orientations to avoid new v-structures and cycles.

pub fn run(&mut self, data: &[Vec<f64>])

Run the full PC algorithm: skeleton → v-structures → Meek rules.

Trait Implementations

impl Clone for CausalDiscovery

fn clone(&self) -> CausalDiscovery

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for CausalDiscovery

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.