REVIEW 5 major objections 5 minor 10 references
Linear-LLM-SCM: Benchmarking LLMs for Coefficient Elicitation in Linear-Gaussian Causal Models
T0 review · 5 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read This paper claims that current LLMs, given only a DAG and variable semantics, cannot reliably estimate linear-Gaussian causal coefficients and quantifies the failure with new scale-invariant metrics.
desk verdict Useful benchmark for a real gap, though the headline claim would be stronger with a baseline and a check that the acceptance gate doesn't reject the ground truth. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The framework's load-bearing mechanism is node-level decomposition with an acceptance gate. The DAG is traversed in topological order; for each node, the LLM receives only the node's direct parents, their units, value ranges, and a phenomenon description, and must return a concrete linear equation with numeric coefficients plus an intercept. An iterative feedback loop (Algorithm 2) computes the implied range of the target from the parents' ranges and the proposal, and rejects proposals whose implied range falls outside the node's stated bounds, re-prompting up to five times. Results are scored by four metrics: M1 raw L2 distance, M2 L2 distance after normalizing each node's coefficient vecto
What would settle it
Construct a synthetic linear-Gaussian SCM on a familiar DAG with coefficients that could not appear in any training corpus (e.g., random nonzero values), elicit coefficients from the same LLMs under the same prompting, and check whether M2/M3 and M4 stay poor. If the models approximate the held-out coefficients well, the reported failure is an artifact of memorization or constraint crafting; if they still fail, the limitation is robust.
Extended reading notes
Core claim
The central claim is that LLMs, when given a pre-specified DAG and variable semantics, cannot reliably parameterize linear-Gaussian SCMs: the node-normalized error metric M2 ranges from roughly 0.45 to 3.25 across model-DAG pairs (far from the perfect score of 0), and the relative-ordering metric M4 falls when spurious edges are added, showing models fail to assign near-zero coefficients to irrelevant causes. The authors additionally report that two of the three tested LLMs remain stochastic even at temperature zero, producing non-deterministic coefficient estimates, while the smallest model is deterministic but less accurate. The claim is specifically about direct, symbolically elicited coe
Load-bearing premise
The benchmark treats the fitted linear-Gaussian coefficients in the source repository as exact ground truth and assumes the hard value ranges supplied (some handcrafted) agree with those fitted effects; if the ranges conflict or if the LLMs have memorized the source papers, the metrics measure agreement with a particular fitted network or retrieval, not causal reasoning.
Editorial extensions
If this is right
- Scale-invariant metrics (M2–M4) should replace or supplement raw L2 distance in future LLM causal benchmarks, since M1 is dominated by variable ranges and misleads cross-DAG comparisons.
- Predicting full causal effect sizes from prior knowledge alone is currently unsafe for healthcare or public-policy decisions; any LLM-elicited coefficients would need data grounding or strong uncertainty handling.
- Adding spurious edges systematically lowers M4 and raises M2/M3, so downstream causal inference should treat DAG structure as uncertain rather than point-estimated.
- Even with temperature set to zero, larger models show non-deterministic estimates, implying reproducibility cannot be assumed for LLM-based parameterization in deployment.
Reading between the lines
- The benchmark never measures training-data leakage directly; if the source papers for the seven DAGs appear in the models' training corpora, the direct RQ1 scores could be inflated by memorization, making actual zero-shot elicitation competence lower than reported.
- At least one DAG's 'hard constraints' were crafted by the authors rather than taken from the literature; rerunning the framework with different plausible ranges would show whether the acceptance gate, not the model's causal knowledge, drives the scores.
- The unit-tweak experiment found that changing micromolar to nanomolar sometimes improved metrics; a clean test would ask the same model for coefficients in both units and check whether the two answers differ by the correct factor, isolating arithmetic scaling from semantic reasoning.
- The framework could be extended to elicit residual variances and confidence intervals, which the authors explicitly left out, turning the benchmark into a full distributional-parameterization test.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces Linear-LLM-SCM, a framework for eliciting linear-Gaussian structural equation coefficients from LLMs given a DAG, variable descriptions, units, and hard value-range constraints. The framework decomposes the DAG into local parent-child prompts, optionally refines proposals via an iterative acceptance gate (Algorithm 2), and compares LLM outputs against ground-truth parameters from the BnRep repository using four metrics M1–M4. Experiments on seven DAGs with three LLMs (Gemini 2.5 Flash, Llama 3.1 8B, Llama 3.3 70B) are reported for direct estimation (RQ1), unit-tweak robustness (RQ2-I), and spurious-edge robustness (RQ2-II). The authors conclude that LLMs show high variability, scale-sensitive M1, moderate normalized distances, and degradation under structural misspecification, and they caution against deploying LLMs as quantitative causal parameterizers in safety-critical domains.
Significance. If the findings hold, the paper provides a useful, open-sourced benchmark for a relatively underexplored task—quantitative coefficient elicitation for continuous SCMs—and makes a reasonable case that scale-invariant metrics (M3/M4) are more meaningful than raw L2 distance. The framework is model- and DAG-agnostic, and the authors are honest about several limitations, including the linear assumption, limited unit-tweak coverage, and the possibility of non-determinism. The external anchoring of ground truth in BnRep avoids the circularity of self-generated derivations. However, the strength of the central claim depends on three unresolved issues: causal provenance of the BnRep parameters, whether the Algorithm 2 hard-constraint gate is compatible with the ground truth, and the absence of chance-level or statistical baselines for interpreting M2–M4.
major comments (5)
- [Section 4.2, Table 2] Ground-truth validity: The BnRep DAGs are a heterogeneous collection; several entries, e.g., expenditure from Tsagris (2022), are learned predictive networks (FEDHC) rather than established causal graphs. Treating every fitted linear-Gaussian parameter as a causal ground truth conflates conditional association with causal effect size. Since the central claim is about LLMs as 'quantitative causal parameterizers', the benchmark should either restrict to DAGs with explicit causal provenance or add a caveat that it evaluates association-coefficient elicitation. A concrete check would be to compare a subset of coefficients against independent domain-expert values.
- [Algorithm 2; Table 2 (VR column)] Compatibility of the acceptance gate with ground truth: For the expenditure DAG, value ranges are 'Crafted' by the authors, not taken from the literature. Algorithm 2 rejects any proposal whose implied range C1 is not contained in the node's hard range C2. If the BnRep ground-truth coefficients themselves imply C1 values outside C2, then the gate would reject the ground-truth parameter vector, and M1–M4 would penalize LLMs for not reproducing an inadmissible parameterization. The paper does not report whether the ground-truth coefficients pass this compatibility check. Without this, a nontrivial portion of the measured 'failure' may reflect the crafted constraints rather than LLM inability.
- [Section 5.2.2, Table 5] Statistical support for degradation under spurious edges: The text states that adversarial conditions generally yield lower M4 values, but the 95% confidence intervals overlap substantially across conditions. For example, Gemini M4 is 7.520±0.200 for O and 7.320±0.271 for S4; Llama 3.3 M4 is 6.560±0.199 for O and 6.680±0.271 for S4. Overlapping intervals do not support a claim of performance degradation. Report paired significance tests, bootstrap difference distributions, or a non-overlap criterion before concluding that spurious edges degrade parameterization.
- [Section 3.4, Table 3] Missing chance-level baseline: M2–M4 have no null reference. M4 counts the number of nodes with exact effect-size ordering, but the expected value under random ordering depends on the number of parents per node; a score of 1.000 on cachexia1 may be far above chance or simply reflect a single multi-parent node with only two parents. Similarly, M2 and M3 values cannot be judged as 'poor' without knowing the distribution under random coefficients. Add a random-coefficient or permutation baseline for each DAG to calibrate the metrics and support the abstract's claim that LLMs are unreliable.
- [Section 4.3, Section 5.1] Data leakage is invoked but not tested: The authors cite Yang et al. (2023) as motivation for adversarial perturbations to mitigate data leakage, yet they never measure whether direct RQ1 scores are inflated by memorization of published BnRep networks. Given that the same networks appear in the literature, LLMs may be recalling coefficients rather than eliciting them. A control condition using permuted variable names/semantics or held-out DAGs would clarify whether the low accuracy reflects a genuine inability to reason about continuous effect sizes.
minor comments (5)
- [Abstract] The phrase 'seven real-world DAGs effect ground truth' is ungrammatical; should be 'seven real-world DAGs with effect ground truth'.
- [Eq. (4)] Index notation in M4 is inconsistent: the surrounding text uses β_{LLM,j,i} and β_{GT,j,i}, but the equation uses β_{LLM,i,j} and β_{GT,i,j}. Unify the subscript order.
- [Algorithm 2] When the loop budget n is reached, the algorithm returns the last proposal even if it failed the validation check. The text does not state whether this is the last accepted proposal or simply the final sampled one; clarify to avoid ambiguity in interpreting rejected samples.
- [Table 2] The 'VR' column uses 'N (Crafted)' for expenditure, but the fact that the ranges are author-crafted is not explained until the discussion of Algorithm 2. Move or annotate this earlier so readers do not mistake the ranges as literature-sourced.
- [Section 5.1] The interpretation of Llama 3.1 8B's 'determinism' would be stronger if the paper reported repeated runs with different seeds or API call configurations, since temperature=0 is not a guarantee of determinism in all serving stacks.
Circularity Check
No significant circularity; benchmark is externally anchored to BnRep ground truth.
full rationale
The paper is an empirical benchmark, not a derivation chain in which a predicted quantity is equivalent to a fitted input. Ground-truth coefficients are taken from the external BnRep repository (Leonelli, 2025), and LLM outputs are compared with those coefficients via the M1–M4 distance metrics; no step defines the target in terms of the inputs or fits a parameter to the data it later 'predicts.' The only author-supplied values are prompt templates and, for the expenditure DAG, hard value ranges marked 'Crafted' in Table 2. These ranges enter Algorithm 2 as a rejection gate, but they are not used to construct the ground truth or to force the conclusion that LLM coefficients deviate from BnRep parameters. Even if the crafted ranges conflict with the BnRep fitted coefficients, that would be a construct-validity threat—the metrics could partly measure constraint compliance rather than causal parameterization quality—but it is not circular, because the comparison target remains external. Similarly, the possibility of LLM memorization of BnRep networks is an unquantified confound, not an equation-level reduction. There are no load-bearing self-citations: the authors do not rely on their own prior results to justify the benchmark's premise. Hence no circular step can be exhibited with the required specificity.
Assumptions & free parameters
free parameters (4)
- Per-variable hard value-range constraints (R) =
e.g., Share ∈ [0,1], Age ∈ [18,100]; for expenditure DAG: 'Crafted' by authors (Table 2)
- Iterative feedback loop budget (n) =
5
- Prompt composition (persona, phenomenon summary, units, unit-effect justification request) =
hand-designed; one example shown in Figure 1
- Sampling configuration =
n=25 runs, temperature=0
assumptions (4)
- domain assumption Each of the seven phenomena is adequately modeled by the given directed acyclic graph with linear-Gaussian structural equations.
- domain assumption The BnRep repository's stored parameters are the correct ground-truth effect sizes.
- domain assumption LLM JSON outputs can be reliably parsed into structural equations.
- standard math The range-propagation calculation in Algorithm 2 (C1 from parent ranges and proposed betas) is a correct and complete check of the possible value range of a node.
Cite this review
Pith. "Pith review of Linear-LLM-SCM: Benchmarking LLMs for Coefficient Elicitation in Linear-Gaussian Causal Models." pith.science (2026). https://pith.science/paper/QOBTN2A4
@misc{pith2026260210282,
author = {Pith},
title = {Pith review of: Linear-LLM-SCM: Benchmarking LLMs for Coefficient Elicitation in Linear-Gaussian Causal Models},
year = {2026},
howpublished = {\url{https://pith.science/paper/QOBTN2A4}},
note = {Machine review of arXiv:2602.10282}
}
read the original abstract
Large language models (LLMs) have shown potential in identifying qualitative causal relations, but their ability to perform quantitative causal reasoning---estimating effect sizes that parametrize functional relationships---remains underexplored in continuous domains. We introduce Linear-LLM-SCM, a plug-and-play framework for evaluating LLMs on Linear Gaussian structural causal model parametrization when a directed acyclic graph (DAG) is given. The framework decomposes a DAG into local parent-child sets and prompts an LLM to produce a regression-style structural equation per node, which is aggregated and compared against available ground-truth parameters. Our experiments with seven real-world DAGs effect ground truth illustrate limitations of LLMs as quantitative causal parameterizers. Across most models, we observe variability in coefficient estimates and sensitivity to structural perturbations. We open-sourced the framework to further encourage the community to work on studies toward the use of LLM for causal effect elicitation in safety-critical domain, e.g., healthcare.
Figures
Reference graph
Works this paper leans on
-
[2]
URLhttps://proceedings.neurips.cc/paper_files/paper/2023/ file/df88b275bef31ac96c85f0c4013734fc-Paper-Conference.pdf. Gemini Team, Google. Gemini 2.5: Pushing the frontier with advanced reasoning, multimodality, long context, and next generation agentic capabilities, 2025. URLhttps://arxiv.org/ abs/2507.06261. William H Greene. Econometric analysis.Preten...
arXiv 2023
-
[4]
Emre Kiciman, Robert Ness, Amit Sharma, and Chenhao Tan
URLhttps://proceedings.neurips.cc/paper_files/paper/2023/ file/631bb9434d718ea309af82566347d607-Paper-Conference.pdf. Emre Kiciman, Robert Ness, Amit Sharma, and Chenhao Tan. Causal reasoning and large language models: Opening a new frontier for causality.Transactions on Machine Learning Research, 2024. ISSN 2835-8856. URLhttps://openreview.net/forum?id=m...
arXiv 2023
-
[8]
Karan Singhal, Shekoofeh Azizi, Tao Tu, S
URLhttps://proceedings.mlr.press/v202/shi23a.html. Karan Singhal, Shekoofeh Azizi, Tao Tu, S. Sara Mahdavi, Jason Wei, Hyung Won Chung, Nathan Scales, Ajay Tanwani, Heather Cole-Lewis, Stephen Pfohl, Perry Payne, Martin Seneviratne, Paul Gamble, Chris Kelly, Abubakr Babiker, Nathanael Sch ¨arli, Aakanksha Chowdhery, Philip Mansfield, Dina Demner-Fushman, ...
arXiv 2023
-
[9]
Unlike these approaches, we assume pre-specified causal DAG structures without observational data
addressing both domains but requiring numerical observational data with specialized atten- tion mechanisms. Unlike these approaches, we assume pre-specified causal DAG structures without observational data. Nafar et al. (2025) similarly benchmarkeffectestimation from pre-specified DAGs but assume discrete domains with Conditional Probability and sampling ...
2025
-
[2020]
doi: 10.1007/s00291-020-00588-8
ISSN 1436-6304. doi: 10.1007/s00291-020-00588-8. URLhttps://doi.org/10. 1007/s00291-020-00588-8. Xiaoyu Liu, Paiheng Xu, Junda Wu, Jiaxin Yuan, Yifan Yang, Yuhang Zhou, Fuxiao Liu, Tian- rui Guan, Haoliang Wang, Tong Yu, Julian McAuley, Wei Ai, and Furong Huang. Large lan- guage models and causal inference in collaboration: A comprehensive survey. In Luis...
-
[2023]
doi: 10.18653/v1/2023.acl-long.302
Association for Computational Linguistics. doi: 10.18653/v1/2023.acl-long.302. URL https://aclanthology.org/2023.acl-long.302/. Roman Eisner, Cynthia Stretch, Thomas Eastman, Jianguo Xia, David Hau, Sambasivarao Dama- raju, Russell Greiner, David S. Wishart, and Vickie E. Baracos. Learning to predict cancer- associated skeletal muscle wasting from 1H-NMR ...
-
[2024]
URLhttps://proceedings.mlr.press/v235/papamarkou24b.html. Judea Pearl. The seven tools of causal inference, with reflections on machine learning.Commun. ACM, 62(3):54–60, February 2019. ISSN 0001-0782. doi: 10.1145/3241036. URLhttps: //doi.org/10.1145/3241036. Jonathan Richens and Tom Everitt. Robust agents learn causal world models. InThe Twelfth In- ter...
doi:10.1145/3241036 2019
-
[2025]
Gaussian
examine LLMs on causality tasks but focus on causal relationship identification rather than quantitative effect estimation in continuous SCMs. A.2.3 SAMPLES OFDAGSUSED IN THEEXPERIMENTS One may easily find DAGs we used in our experiments from references, but to save time of the readers, we attached the two example DAGs, Figure 2 and Figure 3. GM VB A F GC...
Show all 10 references
-
[7928]
David Ha and J ¨urgen Schmidhuber
URLhttp://jmlr.org/papers/v21/18-668.html. David Ha and J ¨urgen Schmidhuber. World models.CoRR, abs/1803.10122, 2018. URLhttp: //arxiv.org/abs/1803.10122. Jilei Hu, Bin Xiong, Zheng Zhang, and Jing Wang. A continuous Bayesian network regression model for estimating seismic li...
2018 arXiv
-
[8440]
URLhttps://www.sciencedirect.com/ science/article/pii/S2405844023110061
doi: 10.1016/j.heliyon.2023.e23798. URLhttps://www.sciencedirect.com/ science/article/pii/S2405844023110061. Freda Shi, Xinyun Chen, Kanishka Misra, Nathan Scales, David Dohan, Ed H. Chi, Nathanael Sch¨arli, and Denny Zhou. Large language models can be easily distracted by irr...
2023 doi
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.