---
title: Neyman–Pearson Causal Discovery
url: https://www.emergentmind.com/topics/neyman-pearson-causal-discovery
type: topic
---

# Neyman–Pearson Causal Discovery

Searching arXiv for the cited papers and topic context.
Neyman–Pearson causal discovery denotes a family of causal inference frameworks that import the logic of classical Neyman–Pearson testing into causal graph learning and causal direction estimation. Across recent work, the central idea is to treat causal discovery as an error-controlled decision problem in which one type of mistake is constrained while another is minimized, or in which competing causal explanations are compared through level-\(\alpha\) tests, directional power, and calibrated rejection behavior. In contemporary formulations, this perspective appears in at least three distinct but related forms: asymmetric edge-error control for graph support recovery, calibrated conditional independence testing inside constraint-based discovery with incomplete data, and paired directional testing for bivariate causal direction estimation [2507.21570], [2605.04838], [2605.13550].

## 1. Conceptual definition and scope

Neyman–Pearson causal discovery is most explicitly defined in work on **unequal edge error tolerance**, where the goal is to **minimize false negatives subject to a user-specified upper bound on false positives**. In that formulation, the causal discovery problem is cast as the analogue of classical Neyman–Pearson testing: among all graph estimators whose false positive rate is at most \(\epsilon\), select one with the smallest false negative rate [2507.21570].

A broader usage of the term includes methods that do not claim a classical Neyman–Pearson lemma or uniformly most powerful rule, but nonetheless formulate causal discovery through level-\(\alpha\) hypothesis tests, directional rejection probabilities, and calibrated inferential modules. Two examples are especially relevant. First, **PAIR-CI** treats conditional independence testing as the inferential engine of PC-style causal discovery and emphasizes that a CI oracle must behave like a level-\(\alpha\) hypothesis test under incomplete data; its contribution is described explicitly through a **Neyman–Pearson lens** [2605.04838]. Second, **CDSP** formulates bivariate causal direction detection as a **pair of hypothesis tests**, defines a causal analogue of power and effect size, and bases direction choice on which directional null is more detectable [2605.13550].

This suggests that “Neyman–Pearson causal discovery” is not a single algorithmic doctrine. Rather, it is a technical orientation in which causal discovery is organized around controlled error rates, explicit null and alternative hypotheses, and the asymmetry between different mistakes. A plausible implication is that the unifying concern is not merely graph recovery accuracy, but whether the inferential decisions embedded in a discovery procedure preserve the meaning of a nominal significance threshold.

## 2. Asymmetric edge-error control in graph support recovery

The most direct formulation appears in **“Causal Link Discovery with Unequal Edge Error Tolerance”** [2507.21570]. The paper studies a DAG
\[
\mathcal{G}=(\mathcal{V},\mathcal{E}), \qquad |\mathcal{V}|=d,
\]
with weighted adjacency matrix
\[
\boldsymbol{A}\in\mathbb{R}^{d\times d},
\]
and defines the support matrix \(\boldsymbol{\chi}\) by
\[
\chi_{i,j}=1 \iff \boldsymbol{A}_{i,j}\neq 0 \ \text{or}\ \boldsymbol{A}_{j,i}\neq 0.
\]
The target is therefore **support / causal link discovery**, not orientation at first [2507.21570].

The aggregate edgewise error rates are defined as
\[
\mbox{false positive rate}: \qquad \epsilon^+ = \frac{\mathbb{E}\sum_{i,j}\mathbbm{1}\{\hat{\boldsymbol{\chi}}_{i,j}=1,\ \boldsymbol{\chi}_{i,j}=0\}}{\mathbb{E}\sum_{i,j}\mathbbm{1}\{\boldsymbol{\chi}_{i,j}=0\}},
\]
and
\[
\mbox{false negative rate}: \qquad \epsilon^- = \frac{\mathbb{E}\sum_{i,j}\mathbbm{1}\{\hat{\boldsymbol{\chi}}_{i,j}=0,\ \boldsymbol{\chi}_{i,j}\neq 0\}}{\mathbb{E}\sum_{i,j}\mathbbm{1}\{\boldsymbol{\chi}_{i,j}\neq 0\}}.
\]
These are explicitly described as **aggregated expected edgewise error rates**, normalized by the expected numbers of absent and present edges, rather than family-wise error probabilities or FDR/FDP [2507.21570].

The central optimization problem is
\[
\begin{aligned}
\inf_{\hat{\boldsymbol{\chi}}}\ \epsilon^- \qquad \textrm{s.t.}\qquad \epsilon^+ \le \epsilon,
\end{aligned}
\]
where \(0<\epsilon<1\) [2507.21570]. This is the paper’s formal definition of **Neyman–Pearson causal discovery**. The constrained quantity is the false positive rate, and the optimized quantity is the false negative rate. The motivation is application-dependent asymmetry: in some domains, declaring a nonexistent causal link is more harmful than omitting a real one.

For each pair \((i,j)\), the edgewise detection problem is written as
\[
H_0:\ \boldsymbol{\chi}_{i,j}=0,\qquad \{X_k\}_{1}^{n}\sim P_{i,j},
\]
\[
H_1:\ \boldsymbol{\chi}_{i,j}=1,\qquad \{X_k\}_{1}^{n}\sim Q_{i,j}.
\]
The paper then proves an analogue of the Neyman–Pearson lemma. The optimal detector is a weighted likelihood-ratio threshold rule involving \(\frac{dP_{i,j}}{dQ_{i,j}}\), prior-induced weights \(\frac{w^-_{i,j}}{w^+_{i,j}}\), a global threshold \(\gamma\), and boundary randomization by \(U\sim \mathrm{Unif}[0,1]\), with \(\gamma,\eta\) selected so that \(\epsilon^+=\epsilon\) [2507.21570]. The presence of randomization is explicitly noted as parallel to classical NP testing.

This is a genuine decision-theoretic formulation, not merely an interpretive analogy. It differs from most score-based or constraint-based graph learning procedures because the false-positive tolerance is fixed in advance and the detector is derived to satisfy that tolerance.

## 3. Information-theoretic characterization and the \(\epsilon\)-CUT algorithm

A major theoretical feature of the same framework is its characterization of performance limits through **Rényi divergence** [2507.21570]. For the optimal detector and any \(\lambda\in[0,1]\), the paper gives an upper bound on \(\epsilon^-\) in terms of \((w^-_{i,j})^{1-\lambda}(w^+_{i,j})^\lambda\), the threshold \(\gamma\), and \(D_\lambda(P_{i,j}\|Q_{i,j})\). It also provides a converse lower bound for any detector satisfying \(\epsilon^+\le \epsilon\), again expressed through \(D_\lambda\), \(D'_\lambda\), and \(D''_\lambda\) [2507.21570]. The stated interpretation is that the tradeoff between constrained false positives and achievable false negatives is fundamentally governed by likelihood-ratio interpolation between the edge-absent and edge-present distributions.

Because the exact optimal detector depends on conditional edge distributions that are mixtures over graph structures and edge weights, it is computationally intractable in general. This motivates the algorithm **\(\epsilon\)-CUT**, designed for **linear additive Gaussian noise models** with equal-variance Gaussian noise:
\[
X_k=\boldsymbol{A}X_k+W_k,
\]
where \(W_k\) are i.i.d. zero-mean Gaussian with covariance \(\sigma^2 I\), and \(\sigma^2\) is assumed known [2507.21570].

The method exploits the fact that, after regressing a node on its true parents, the residual sum of squares has a chi-squared distribution. If \(\hat{\alpha}_i\) is the least-squares estimate using the true parent set \(\mathcal{Z}(i)\), then
\[
\hat{\sigma}^{*2}_i=\sum_{k=1}^n \big(X_k(i)-\hat{\alpha}_i^\top Z_k(i)\big)^2
\]
satisfies
\[
\hat{\sigma}^{*2}_i/\sigma^2 \sim \chi^2_{n-p},
\]
conditioned on the parent observations, where \(p=|\mathcal{Z}(i)|\) [2507.21570].

The resulting procedure tests every unordered pair \((i,j)\) by comparing residual energies under candidate parent sets and declares \(\hat{\chi}_{i,j}=0\) if any candidate pair of parent sets yields sufficiently close residual sums of squares. Otherwise it declares \(\hat{\chi}_{i,j}=1\) [2507.21570]. The threshold is selected via a chi-squared cdf equation so that the false-positive tolerance equals \(\epsilon\).

The main finite-sample guarantee is:

- **Theorem 3**: for any sample size \(n\) and any false positive tolerance \(\epsilon\), the false positive rate of \(\epsilon\)-CUT satisfies
  \[
  \epsilon_1^+ \le \epsilon.
  \]

The paper further states a stronger per-edge result: if there is truly no edge between \(i\) and \(j\), then
\[
\mathbb{P}_{\boldsymbol{A}}(\hat{\boldsymbol{\chi}}_{i,j}=1)\le \epsilon.
\]
Thus the method provides true finite-sample false-positive control for support recovery in the specified linear Gaussian SEM regime [2507.21570].

The framework’s scope is correspondingly narrow. The guarantees are for support rather than orientation; direction estimation is deferred to a post-processing stage in which an unrestricted maximum-likelihood DAG estimate is masked by the support estimate:
\[
\hat{\boldsymbol{A}}=\hat{\boldsymbol{A}}_{ML}\circ \hat{\boldsymbol{\chi}}.
\]
The paper also notes exponential complexity in \(d\), with at most
\[
\frac{d(d-1)}{2}4^{d-2}
\]
linear regressions and overall complexity
\[
\mathcal{O}\big(d^3(d-1)4^{d-2}(n+d)\big)
\]
[2507.21570].

## 4. Neyman–Pearson logic in conditional independence testing for causal graphs

A second line of work applies Neyman–Pearson reasoning not to direct edge detection, but to the **local conditional independence tests** that drive constraint-based causal discovery. **“PAIR-CI: Calibrated Conditional Independence Testing for Causal Discovery with Incomplete Data”** develops a CI test specifically for **constraint-based causal discovery with incomplete data** and frames its contribution through a **Neyman–Pearson lens**: the objective is to provide a CI oracle that behaves like a level-\(\alpha\) test when standard incomplete-data workflows do not [2605.04838].

The paper formalizes the scientific null as
\[
H_0:\ Z \perp Y \mid X,
\]
with alternative
\[
H_1:\ Z \not\perp Y \mid X.
\]
Its strongest motivation is **Proposition 1 (Miscalibration of impute-then-test)**, which states that if the imputation model is asymptotically misspecified and imputation error induces spurious dependence under the imputed distribution, then any consistent CI test satisfies
\[
\lim_{n \to \infty} P_{H_0}(T(\hat{D}_n)\ \textup{rejects}) = 1.
\]
The point is not that the CI test is defective on complete data, but that the imputation step changes the effective null [2605.04838]. This is an explicitly Neyman–Pearson-style failure: the nominal significance level no longer refers to the intended scientific null.

PAIR-CI addresses this by integrating multiple imputation into the inferential construction. It generates
\[
\hat{D}^{(1)},\ldots,\hat{D}^{(M)},
\]
then, on each imputed dataset and fold, compares a **full model** for predicting \(Z\) using \((X,Y)\) with a **partial model** using the same imputed conditioning set and a placebo variable \(U\) produced by conditional permutation. The fold-level loss contrast is
\[
\hat{\mu}_{km} = \ell\!\left(\hat{g}^{ZU}_{km}\right) - \ell\!\left(\hat{g}^{ZY}_{km}\right),
\]
and in empirical form
\[
\hat{\mu}_{km} = \frac{1}{n_k} \sum_{i \in F_k} \left[ \ell\!\left(z_i,\, \hat{g}^{ZU}_{km}(\hat{x}_i, u_{ik})\right) - \ell\!\left(z_i,\, \hat{g}^{ZY}_{km}(\hat{x}_i, y_i)\right) \right].
\]
Aggregation gives
\[
\hat{\mu}_m = \frac{1}{K}\sum_{k=1}^K \hat{\mu}_{km}, \qquad \hat{\mu} = \frac{1}{M}\sum_{m=1}^M \hat{\mu}_m.
\]
The conceptual heart is the **paired design**: both models use the same imputed \(\hat X\), so imputation error is shared and is intended to cancel in the loss difference [2605.04838].

The paper distinguishes the **internal null**
\[
Z \perp Y \mid \hat{X}
\]
from the scientific null
\[
Z \perp Y \mid X.
\]
Conditional permutation gives calibration for the internal null, and imputation consistency extends this to the scientific null [2605.04838].

For inference, the combined variance is
\[
T = \bar{W} + (1 + 1/M)B,
\]
where
\[
\bar{W} = M^{-1}\sum_m n^{-1}K^{-1}\sum_k \hat{\sigma}_{k,m}^2, \qquad B = (M-1)^{-1}\sum_m (\hat{\mu}_m - \hat{\mu})^2.
\]
The studentized statistic is
\[
t_{\mathrm{CI}} = \frac{\hat{\mu}}{\sqrt{T}},
\]
compared against a \(t\)-distribution with **Barnard–Rubin adjusted degrees of freedom** \(\nu\), using the one-sided rejection region
\[
t_{\mathrm{CI}} > t_{\alpha,\nu}.
\]
The paper proves:

- **Proposition 2 (Calibration)**:
  \[
  P(t_{\mathrm{CI}} > t_{\alpha,\nu}) \to \alpha.
  \]

- **Proposition 3 (Consistency)**:
  \[
  P(t_{\mathrm{CI}} > t_{\alpha,\nu}) \to 1 \quad \text{as } n\to\infty.
  \]

- **Theorem 4 (Unified inference under cross-validation and multiple imputation)**:
  \[
  \frac{\hat{\mu}-\mu}{\sqrt{T}} \xrightarrow{d} \mathcal N(0,1).
  \]

- **Corollary 5 (PC consistency)**:
  \[
  \hat G \xrightarrow{P} G_{\mathrm{CPDAG}} \quad \text{as } n\to\infty.
  \]

The paper’s empirical emphasis is on null calibration. In standalone CI experiments with about 30% missingness, the calibration table reports average false positive rates
\[
\text{PAIR-CI}: 0.018 \text{ (Complete)},\ 0.021 \text{ (MAR)},\ 0.036 \text{ (MNAR)},
\]
compared with MNAR false positive rates of \(0.350\) for FZ-single, \(0.282\) for FZ-Rubin, \(0.415\) for imputed GCM, and \(0.447\) for imputed KCI [2605.04838]. The abstract summarizes this as existing imputation-based CI tests exhibiting false positive rates of **28–45%** under MNAR, whereas PAIR-CI averages below the nominal **5%** level. When used inside PC, the method reduces structural Hamming distance by **8% on 10-variable nonlinear graphs**, **15% on 30-variable nonlinear graphs**, and up to **44%** on the **56-variable HAILFINDER network** [2605.04838].

This line of work therefore locates Neyman–Pearson causal discovery at the level of the CI oracle: what matters is not only predictive performance, but whether the local test actually retains its level-\(\alpha\) meaning under the scientific null.

## 5. Directional power and effect-size asymmetry in bivariate discovery

A third formulation appears in **“Causal Discovery via Statistical Power (CDSP)”**, which addresses **bivariate observational** causal direction inference [2605.13550]. The paper does not claim classical NP optimality, but explicitly says that it formulates causal direction detection as a **pair of hypothesis tests** and introduces directional analogues of **power**, **effect size**, and **directional evidence**.

Under the standing assumption that the true direction is \(X \to Y\), the data-generating mechanism is
\[
Y = m(X) + \eta,
\]
with \(\eta \perp X\). The reverse-direction working model is
\[
X = l(Y) + \xi,
\]
where \(l\) is the risk minimizer over a model class \(\mathcal L\), and \(\xi := X-l(Y)\) [2605.13550].

The directional nulls are
\[
\begin{aligned}
H_Y^0 &: X \perp \eta,\; m \in \mathcal{M}, \\
H_Y^1 &: \text{otherwise}, \\
H_X^0 &: Y \perp \xi,\; l \in \mathcal{L}, \\
H_X^1 &: \text{otherwise}.
\end{aligned}
\]
Thus each direction is evaluated by whether a valid structural model in that direction is compatible with the data [2605.13550].

The paper defines **Directional Power** as
\[
\mathbb{P}\Big( (\text{fail to reject }H_Y^0)\wedge(\text{reject }H_X^0) \,\Big|\, X\to Y \Big),
\]
which is the joint success probability of retaining the correct-direction null while rejecting the reverse-direction null [2605.13550].

Its central quantitative objects are the directional effect sizes
\[
\Delta_{Y,n} = \frac{\theta_Y\sqrt{n}}{\sigma_Y}, \qquad \Delta_{X,n} = \frac{\theta_X\sqrt{n}}{\sigma_X},
\]
and the standardized critical values
\[
q_{Y,n}^\alpha = \frac{\sqrt{n}\, c_{Y,n}^\alpha}{\sigma_Y}, \qquad q_{X,n}^\alpha = \frac{\sqrt{n}\, c_{X,n}^\alpha}{\sigma_X}.
\]
The key assumption is the **effect-size asymmetry assumption**
\[
\frac{\theta_X}{\sigma_X} > \frac{\theta_Y}{\sigma_Y} \quad \text{and} \quad q_{X,n}^\alpha - q_{Y,n}^\alpha = o(1).
\]
Operationally, CDSP uses the directional detectability indices
\[
I_{Y,n}:=\frac{\theta_Y - c_{Y,n}^\alpha}{\sigma_Y}, \qquad I_{X,n}:=\frac{\theta_X - c_{X,n}^\alpha}{\sigma_X},
\]
and infers \(X\to Y\) when \(\hat I_{X,n} > \hat I_{Y,n}\) [2605.13550].

The main theorem states that under a specified misspecification regime,
\[
\mathbb{P}\!\left((\text{reject } H_X^0)\wedge(\text{fail to reject } H_Y^0)\right)
>
\mathbb{P}\!\left((\text{reject } H_Y^0)\wedge(\text{fail to reject } H_X^0)\right)
\]
if and only if
\[
I_{X,n} > I_{Y,n},
\]
equivalently,
\[
\Delta_{X,n} - q_{X,n}^\alpha > \Delta_{Y,n} - q_{Y,n}^\alpha.
\]
The paper further proves directional consistency:
\[
\mathbb{P}\bigl(\hat{I}_{X,n} > \hat{I}_{Y,n}\bigr) \to 1
\]
when \(I_X>I_Y\) and the plug-in quantities are consistent [2605.13550].

Uncertainty quantification is not given as a p-value or confidence set for the final direction. Instead the paper defines a bootstrap **directional support probability**
\[
\widehat{P}_{\mathrm{CDSP}}
=
\frac{1}{B}\sum_{b=1}^B \mathbf 1\{\hat I^{(b)}_{X,n} - \hat I^{(b)}_{Y,n}>0\},
\]
for the direction \(X\to Y\), with an analogous expression for \(Y\to X\). The paper interprets this as the probability that the CDSP-favored direction remains favored under resampling [2605.13550].

Empirically, the paper studies model misspecification in
\[
Y = \mathrm{sign}(X-a)\,|X-a|^{d}\beta + \eta,
\]
with \(d \in \{1,1.2,1.25,1.3,1.4,1.5,3\}\). Reported simulation accuracies are \(100\%\) for \(d=1,1.2,1.25,1.3\), \(88\%\) for \(d=1.4\), and \(57\%\) for \(d=1.5\), while LiNGAM drops to \(0\%\) by \(d=1.25\) in the listed settings [2605.13550]. On the **Tübingen Cause–Effect Pairs benchmark**, CDSP achieves
\[
\mathrm{TDR}=62/100=62\%, \qquad \mathrm{FDR}=38/100=38\%,
\]
whereas LiNGAM yields \(44\%\) TDR and \(56\%\) FDR, corresponding to an approximately **18 percentage point** reduction in FDR [2605.13550].

CDSP therefore illustrates a distinct Neyman–Pearson-style idea: causal direction can be chosen by comparing the rejection behavior of two composite directional nulls, not by a single likelihood-ratio test between causal directions.

## 6. Related precursors and adjacent testing-based approaches

Earlier and adjacent work shows that the Neyman–Pearson perspective did not emerge in isolation. **“Causal Discovery from Changes”** studies learning from localized mechanism changes across a sequence of environments [1301.2312]. It is not framed as a Neyman–Pearson method, but its pipeline depends on repeated binary decisions of the form “change” versus “no-change,” based on tests of
\[
H_0: P^{j-1}(X)=P^j(X) \quad \text{vs.} \quad H_1: P^{j-1}(X)\neq P^j(X),
\]
with explicit discussion of significance level \(\alpha\), Type I error, and Type II error [1301.2312]. The method then translates those outcomes into descendant and nondescendant constraints under influentiality and local mechanism-change assumptions. This suggests an important precursor pattern: local hypothesis tests can be endowed with causal semantics and assembled into structural decisions, even if no explicit NP optimality result is proved.

A different adjacent approach is **“Causal generalized linear models via Pearson risk invariance”**, which identifies the causal parent set of a target variable \(Y\) using a hybrid of invariance testing and model selection [2407.16786]. The target obeys
\[
Y \mid X_{PA} = x_{PA} \sim EDF\Big(b(f_{PA}(x_{PA})),a(\phi)\Big),
\]
and the causal model is characterized by a **Pearson risk** identity
\[
\mathbb{E}_{X,Y}
\left[\frac{(Y- \dot{b}(f_{PA}(X_{PA})))^2}{\ddot{b}(f_{PA}(X_{PA}))}\right] = a(\phi),
\]
together with expected-likelihood maximization on the true parent set [2407.16786]. Finite-sample implementation uses the Pearson statistic
\[
R^{+}=\sum_{i=1}^n \left[\frac{(y_i- \dot{b}(\widehat{\beta}_{S}^\top\tilde\psi_S(x_{i,S})))^2}{\ddot{b}(\widehat{\beta }_{S}^\top\tilde\psi_S(x_{i,S}))a(\phi)}\right],
\]
with approximate null distribution
\[
R^{+} \overset{H_0}{\sim} \chi^2_{n-\hat q},
\]
and a two-sided acceptance region [2407.16786]. This is not presented as Neyman–Pearson causal discovery, but it shows a closely related architecture: candidate causal models are screened by a calibrated null, then selected by BIC among the accepted models.

These adjacent cases clarify a common misconception. Neyman–Pearson causal discovery is not restricted to one specific graph-learning algorithm, nor does it require likelihood-ratio optimality in the classical simple-vs-simple sense. The recurring feature is the elevation of statistical testing, calibration, and asymmetric decision costs to first-class causal discovery principles.

## 7. Limitations, controversies, and technical boundaries

Several limitations recur across the literature.

The framework in [2507.21570] is the most classically Neyman–Pearson in form, but its tractable algorithm \(\epsilon\)-CUT is specialized to **linear equal-variance Gaussian SEMs** with known \(\sigma^2\), and the finite-sample guarantee concerns **support / skeleton recovery**, not orientation. The optimal detector itself is generally computationally infeasible because the conditional edge distributions are mixtures over graphs and edge weights. The practical method also has exponential complexity in the number of variables [2507.21570].

PAIR-CI [2605.04838] addresses a different problem: preserving level-\(\alpha\) behavior of CI tests under incomplete data. Its scientific-null calibration theorem formally relies on **imputation consistency** (A1), even though the paper argues that the paired construction is empirically robust beyond strict (A1). The authors are explicit that robustness under MNAR is practical rather than fully covered by the theorem. The method is also computationally heavier than parametric CI tests; on HAILFINDER, about 2400 tests take around 100 minutes per replicate. The paper does not provide a multiplicity treatment for the many sequential CI tests inside PC [2605.04838].

CDSP [2605.13550] likewise stops short of classical NP optimality. It provides no uniformly most powerful result, no global finite-sample Type I error control for the final direction decision, and no p-value for the selected causal direction. Its uncertainty measure is a bootstrap support probability rather than a confidence set. The effect-size asymmetry assumption is central but cannot be checked directly from observed data. Severe misspecification can reverse the asymmetry and induce confidently wrong decisions [2605.13550].

A broader controversy concerns what exactly should be controlled in causal discovery. The explicit NP framework in [2507.21570] constrains aggregated edgewise false positives. PAIR-CI emphasizes local test calibration because, in PC-like procedures, false CI rejections propagate through orientation rules and can be more damaging than false non-rejections [2605.04838]. CDSP instead compares directional success probabilities for paired nulls [2605.13550]. These are different inferential targets. This suggests that “error control” in causal discovery is not monolithic: support recovery, CI oracle behavior, and directional selection each instantiate distinct NP-style objects.

A second misconception is that better predictive performance is sufficient for reliable causal discovery. PAIR-CI explicitly argues that **calibration matters at least as much as raw predictive performance**, because miscalibrated CI tests create spurious edges that then propagate through orientation rules [2605.04838]. CDSP similarly shifts attention from raw dependence asymmetry to threshold-adjusted detectability indices [2605.13550]. The common implication is that causal discovery can fail even when its predictive subroutines are individually strong, if their rejection behavior is misaligned with the scientific nulls.

In summary, Neyman–Pearson causal discovery is best understood as a family of testing-centered formulations of causal inference in which causal claims are governed by constrained false-positive tolerance, level-\(\alpha\) local null calibration, or directional power comparisons. Its most direct expression is the constrained optimization
\[
\inf_{\hat{\boldsymbol{\chi}}}\epsilon^- \quad \text{s.t.}\quad \epsilon^+\le\epsilon
\]
for support recovery [2507.21570]. Its broader significance lies in showing that causal discovery is often limited less by lack of predictive flexibility than by the lack of inferential procedures whose nulls, alternatives, and error rates remain meaningful under the data conditions at hand [2605.04838], [2605.13550].

Source: https://www.emergentmind.com/topics/neyman-pearson-causal-discovery