Many benchmarks for automated causal inference evaluate a system's performance based on a single numerical output, such as an Average Treatment Effect (ATE). This approach conflates two distinct steps in causal analysis: identification, formulating a valid research design under stated assumptions, and estimation, implementing that design numerically on finite data. We introduce CausalReasoningBenchmark, a benchmark of 174 queries across 139 real-world datasets, curated from 73 peer-reviewed research papers and three widely used causal-inference textbooks. For each query, a system must produce (i) a structured identification specification that names the strategy, the treatment, outcome, and control variables, and all design-specific elements, and (ii) a point estimate with a standard error. By scoring these two components separately, our benchmark enables granular diagnosis: it distinguishes failures in causal reasoning from errors in numerical execution. Baseline results with a state-of-the-art LLM show that, while the model correctly identifies the high-level strategy in 81% of cases, full identification-specification correctness drops to only 38%, revealing that the bottleneck lies in the nuanced details of research design rather than in computation. CausalReasoningBenchmark is publicly available on Hugging Face and is designed to foster the development of more robust automated causal-inference systems.
Figures & tables
Benchmark
# Queries
Real data
ID eval
Quant. eval
Design- specific
Designs covered
QRData [ 64 ]
411
✓
Partial
Mixed
CausalBench [ 89 ]
495
Partial
Mixed
CLadder [ 51 ]
10k
✓
Graph-based
Corr2Cause [ 52 ]
413k
✓
Graph-based
CausalGraph2LLM [ 80 ]
700k+
✓
Graph-based
CausalReasoningBenchmark
174
✓
✓
✓
✓
IV, RDD, DiD, CE, RCT
Table 1: Comparison of CausalReasoningBenchmark with related benchmarks. “ID eval” indicates whether identification is evaluated separately from estimation. “Real data” indicates whether the benchmark uses real-world (non-synthetic) datasets. “Quant. eval” means evaluating causal effect estimation from data (e.g., ATE / ATT / LATE / CATE). “Design-specific” indicates whether the benchmark requires specification of design-specific elements (e.g., instruments, running variables).
Source group
#queries
#datasets
Research papers
120
85
Textbook
54
54
Total
174
139
Table 2: CausalReasoningBenchmark queries and datasets by source group.
Identification strategy
#queries
#datasets
Difference-in-Differences
67
37
Regression Discontinuity
44
39
Instrumental Variable
22
22
Conditional Exogeneity
40
40
RCT
1
1
Table 3: CausalReasoningBenchmark composition by identification strategy.
Source group
DiD
RDD
IV
Cond. Exog.
RCT
Research papers
62
39
19
0
0
Textbook
5
5
3
40
1
Table 4: Query counts by source group and identification strategy.
Metric
Value
Identification Metrics
Strategy correct
81.0% (141/174)
Causal quantity correct
74.7% (130/174)
Treatments correct
87.4% (152/174)
Outcomes correct
96.0% (167/174)
Minimal controlling set included
75.9% (132/174)
Table 5: Aggregate evaluation of the GPT-5.3 baseline on all 174 queries. Identification metrics are exact-match or set-based checks against the gold specification; estimation metrics compare the returned effect and uncertainty to the gold solution. Values in brackets denote the interquartile range.
Design
#Queries
Strategy (%)
Full ID (%)
Med. %Err
Difference-in-Differences
67
94.0
59.7
9.7
Regression Discontinuity
44
95.5
11.4
21.0
Instrumental Variable
22
86.4
36.4
50.1
Conditional Exogeneity
40
40.0
30.0
4.0
RCT
1
100.0
100.0
0.0
Overall
174
81.0
37.9
13.5
Table 6: Per-strategy breakdown of the GPT-5.3 baseline. “Strategy” = fraction with correct strategy label; “Full ID” = fraction with fully correct identification specification; “Med. %Err” = median percentage error on the effect estimate.
Error Family
Error Type
Count
Control variables
Missing required controls
42
Included bad / post-treatment controls
10
Estimand
CATE → LATE
28
CLATE → LATE
7
Other estimand errors
9
Strategy
CE → RCT
24
Table 7: Taxonomy of identification errors across all 108 incorrect cases. A single case may contribute to multiple categories. Counts are grouped by error family.
Appendix figures & tables1 asset
Supplementary material from the paper’s appendix.
Appendix
Paper
Design
Title
#queries
[ 59 ]
DiD
Corporate Board Quotas and Gender Equality Policies in the Workplace
2
[ 75 ]
DiD
Deadly Populism: How Local Political Outsiders Drive Duterte’s War on Drugs in the Philippines
2
[ 23 ]
DiD
Does Compliance Pay? Social Standards and Firm-Level Trade
2
[ 40 ]
DiD
Does Direct Democracy Hurt Immigrant Minorities? Evidence from Naturalization Decisions in Switzerland
2
[ 71 ]
DiD
Education or Indoctrination? The Violent Origins of Public School Systems in an Era of State-Building
2
[ 88 ]
DiD
Elite Cleavage and the Rise of Capitalism under Authoritarianism: A Tale of Two Provinces in China
2
Appendix
Table 8: Research papers included in the benchmark. Each row summarizes queries from one paper and design; papers used under multiple designs appear in separate rows.
Large language models (LLMs) increasingly act as integrated data-science agents, combining abstract reasoning with advanced tool use. Yet the relevant benchmark landscape largely divides into symbolic causal reasoning benchmarks without realistic data analysis or data analysis benchmarks without a principled causal data-generating structure. Furthermore, existing causal evaluation datasets are often restricted to curated examples from existing sources, with diversity coming from limited templatized variations rather than from systematic generation of novel synthetic causal structures. We introduce CausalDS, a benchmark for evaluating causal reasoning in agentic data-science workflows. Each benchmark instance is a scene consisting of a sampled structural causal model (SCM) with generated observational data and an accompanying synthetic natural-language story grounded in a realistic domain. We optionally ground the composition of the benchmark components in empirical distributions obtained from real-world datasets, thus retaining empirical structure while reducing the "causal parrot" risk through completely synthetic generation. From each scene, we then derive tasks spanning all three of Pearl's rungs, with typical data-science prediction tasks appearing as Rung 1. Most tasks include a data science coding component, where the model typically needs to use several tools to arrive at the final answer due to the frequent presence of imperfect observations, which are generated by an observation model. Additionally, recognizing when a question admits no warranted answer and abstaining is treated as a first-class scored outcome. The benchmark thus jointly evaluates symbolic causal reasoning, data science, uncertainty quantification, abstention, and tool use/coding.
Andrej Leban, Yuekai Sun
Department of Statistics University of Michigan Ann Arbor, MI, United States
Existing causal-inference benchmarks for LLMs mostly score method descriptions or whether generated code runs, not whether the executed workflow recovers the target causal estimate. CausalVerify studies this verification problem for structured econometric causal-estimation workflows by separating realistic interpretation from verifiable computation. It pairs 259 published economics papers (reconstructed research question, data description, institutional context) with 100 fixed-seed synthetic scenarios that realise CSV datasets for difference-in-differences, event study, instrumental variables, and regression discontinuity designs. Experiment A (real-paper text agreement) scores method-family and direction agreement against four-LLM consensus labels. Experiment B (synthetic execution) runs model-written R code and checks whether the extracted treatment-effect estimate matches a canonical estimator on the same realised dataset; this execution-grounded correctness layer is L2b+, distinct from L2b, which records only whether the code executes. A calibration arm asks whether self-reported confidence separates correct from incorrect workflows. On Experiment B, seven LLMs reach L2b+ pass rates of 10% to 88% at the default 50% tolerance, and 66 of the 426 workflows that execute (15.5%) return a wrong estimate. Execution ranking (L2b) agrees with L2b+ far better than text-direction scoring (L4): Kendall τ=0.81 and Spearman ρ=0.93, versus Kendall τ between −0.20 and 0.10 for L4. Llama-3.3-70B-Instruct shows the same qualitative gap, and reported confidence does not reliably separate correct from incorrect workflows. The claims are confined to standardized single-shot workflows in these four design families under the evaluated R backend and model panel; the benchmark does not measure general causal-inference ability. Code, data, cached outputs, and a datasheet are released.
Yonghong Zhang, Ricardo Correia, Isabel M. Parra +1
Department of Finance Universidad Autónoma de Madrid Madrid, Spain · Spanish National Research Council (CSIC) Madrid, Spain
In graphical causal model, causal discovery aims to construct a causal graph based on numerical data and domain knowledge in plain text. However, the evaluation of causal discovery methods remains a challenge in the area as the progress of domain researches often makes benchmark causal graphs contain mis-aligned knowledge. This problem especially affects the evaluation of large language model (LLM) based causal discovery methods as they are sensitive to the new discoveries in the literature. This work is the first to systematically study the quality of benchmark causal graphs. Specifically, we design a pipeline that automatically retrieves relevant research papers from scientific databases, and prompts LLMs to check the consistency between the benchmark causal graphs and domain research papers. We evaluate 11 popular real-world benchmarks, for which our pipeline in total proceeds 38,081 domain papers. Our results show that popular benchmarks vary significantly in their consistency with domain research, with clear implications for causal discovery research.