Exact Search
Algorithm Introduction
Search for the optimal graph using Dynamic Programming (DP 1) or A* search 2.
Usage
from causallearn.search.ScoreBased.ExactSearch import bic_exact_search
dag_est, search_stats = bic_exact_search(X, super_graph, search_method,
use_path_extension, use_k_cycle_heuristic,
k, verbose, include_graph, max_parents)
Parameters
X: numpy.ndarray, shape=(n, d). The data to fit the structure too, where each row is a sample and each column corresponds to the associated variable.
super_graph: numpy.ndarray, shape=(d, d). Super-structure to restrict search space (binary matrix). If None, no super-structure is used. Default is None.
search_method: str. Method of exact search ([‘astar’, ‘dp’]). Default is astar.
use_path_extension: bool. Whether to use optimal path extension for order graph. Note that this trick will not affect the correctness of search procedure. Default is True.
use_k_cycle_heuristic: bool. Whether to use k-cycle conflict heuristic for astar. Default is False.
k: int. Parameter used by k-cycle conflict heuristic for astar. Default is 3.
verbose: bool. Whether to log messages related to search procedure.
max_parents: int. The maximum number of parents a node can have. If used, this means using the k-learn procedure. Can drastically speed up algorithms. If None, no max on parents. Default is None.
Returns
dag_est: numpy.ndarray, shape=(d, d). Estimated DAG.
search_stats: dict. Some statistics related to the seach procedure.
- 1
Silander, T., & Myllymäki, P. (2006, July). A simple approach for finding the globally optimal Bayesian network structure. In Proceedings of the Twenty-Second Conference on Uncertainty in Artificial Intelligence (pp. 445-452).
- 2
Yuan, C., & Malone, B. (2013). Learning optimal Bayesian networks: A shortest path perspective. Journal of Artificial Intelligence Research, 48, 23-65.