Summary
Besides computational complexity, the weak spot of constraint-based causal discovery is the assumption of having an independence oracle. Especially with large conditioning sets, independence testing is difficult. This paper aims to address this limitation by restricting the number of conditioning variables and provides a sound framework that allows us to study the guarantees one can obtain for constraint-based structure learning if one restricts the size of the conditioning sets. The authors show that the results can contain bi-directed edges and relate the obtained graph to maximal ancestral graphs. They further provide a sound algorithm to learn such graphs and evaluate it on a few synthetic and semi-synthetic examples.
The theoretical contributions of the paper are interesting and presented clearly to most parts. The main limitation of the work are the points that I mention with respect to the empirical evaluation.
Strengths
- Rigorous definition of assumptions, goal, and solution that can be obtained under the concept of restricting the size of the conditioning set.
- The new formulation of k-closed graphs was introduced well and discussed in the context of previous notions (ancestral graphs). Further, Markov equivalence and faithfulness are adapted for the introduced graphs.
- The authors propose a learning algorithm k-PC to learn the Markov equivalence class of the introduced k-closure of the true DAG.
- The proposed algorithm has been empirically evaluated to PC, LOCI and NOTEARS, and the code has been provided for reproducibility.
Weaknesses
- The performance gain seems to be only evident in small sample regimes.
- The new approach does not clearly outperform LOCI & it would be interesting to see how PC performs when restricting the conditioning set, as well.
- It would be interesting to study a real-world example, or DAGs with non-linear relations. I could imagine that on the latter, using a smaller conditioning set can be a more substantial advantage, as independence testing with large conditioning sets is especially difficult for arbitrary functional relationships. For more details on the problem and maybe a good source to motivate the proposed algorithm please consider Shah and Peters (2020)
Shah, Peters, The Hardness of Conditional Independence Testing and the Generalised Covariance Measure, The Annals of Statistics, 2020
Minor (suggestions):
- For easier comparison, it might be beneficial to compute the area under the curve for the metrics shown in Figure 3.
Questions
- What is the intuition of step 5 in the algorithm? Why do we have to exclude it for Corollary 4.2? Not much intuition is provided regarding this step, and it just occurs in Corollary 4.2.
- I understand that if we just run the PC algorithm and restrict the conditioning set it does not inherit the same guarantees as PC-k, e.g., it would not find any bi-directed edges. What would be the other differences/failure cases of this naive approach?
- Assume one naively restricts the maximum size of the conditioning set in the PC algorithm, which I think is an option in the pcalg package, for example. How does this version of PC perform if we restrict the conditioning size in a similar manner as for the proposed method (e.g., for the experiments shown in Fig. 3 & 4)?
- How does the computational complexity or empirical runtime compare to PC?
Rating
7: Accept: Technically solid paper, with high impact on at least one sub-area, or moderate-to-high impact on more than one areas, with good-to-excellent evaluation, resources, reproducibility, and no unaddressed ethical considerations.
Confidence
4: You are confident in your assessment, but not absolutely certain. It is unlikely, but not impossible, that you did not understand some parts of the submission or that you are unfamiliar with some pieces of related work.
Limitations
Standard limitations occurring in constraint-based causal discovery (a variant of faithfulness, independence oracle).