polytree {PolyTree}

Causal Polytree Estimation

Description

Estimates directed causal polytree from data, using algorithm developed in Chatterjee and Vidyasagar (2022).

Usage

polytree(x)

Arguments

x

Data matrix, whose rows are i.i.d. data vectors generated from the model.

Value

A directed polytree estimated from the input data, as an igraph object.

References

Sourav Chatterjee and Mathukumalli Vidyasagar (2022). Estimating large causal polytrees from small samples. Available at https://arxiv.org/abs/2209.07028

Examples

p <- 10
n <- 200
x <- matrix(nrow = n, ncol = p)
for (i in 1:n) {
   x[i,1] = rnorm(1)
   for (j in 2:p) {
      x[i,j] = (x[i,j-1] + rnorm(1))/sqrt(2)
   }
}
p <- polytree(x)