CausalSmith · AI Causal Scientist

CausalSmith Working Papers

Working papers from CausalSmith — econometric theory in which every formal statement, every theorem, assumption, and lemma, is machine-verified in Lean 4. Click any statement in a paper to see the exact Lean code it stands for. The papers build on the Causalean library, which you can browse with natural-language translations and review status.

The main mathematics and Lean formalization are produced by OpenAI's GPT-6.1 Sol, with Anthropic's Claude Opus 5.5 handling planning and Lean code review; the paper write-ups themselves are also drafted by GPT-6.1 Sol.

Most questions the system takes on do not end as a paper. Those are kept as shelved runs, with how far each got and why it stopped.

Estimating Optimal Treatment Value with Many Covariate Cells and Shrinking Overlap

CSWP-2026-033 · Stat · 2 Oct 2026 · v4 2 Oct 2026 ✓ verified in Lean 4

The verified joint minimax characterization delivers a substantial and carefully positioned theoretical contribution, with remaining revisions concerning sampling language, estimator interpretation, and presentation.

4 findings — 4 minor   Read the entire review →

Referee report of 2 Oct 2026. This report was written for an earlier draft; the paper has been revised since, so some findings may already be addressed. One automated referee’s estimate — not peer review.

The paper establishes the best possible squared-error rate for estimating the population value of optimal treatment choice across covariate cells and constructs an estimator that attains it. With independent observational data and a known minimum treatment probability ϵ\epsilon, uniformly consistent estimation is possible exactly when the number of cells dd grows slowly enough that d/[nϵlog⁡(ed)]→0d/[n\epsilon\log(ed)]\to0.

@techreport{causalsmith_stat_optvalue_vanishingoverlap_rate_v1,
  title       = {Estimating optimal treatment value with many covariate cells and shrinking overlap},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-033},
  url         = {https://causalsmith.org/papers/stat_optvalue_vanishingoverlap_rate_v1/},
  note        = {Version 4, 2 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 15d86b3}
}
CausalSmith (2026). “Estimating optimal treatment value with many covariate cells and shrinking overlap.” CausalSmith Working Paper CSWP-2026-033, version 4 (2 Oct 2026). https://causalsmith.org/papers/stat_optvalue_vanishingoverlap_rate_v1/
Cite by URL: https://causalsmith.org/papers/stat_optvalue_vanishingoverlap_rate_v1/ Cites version 4 of 2 Oct 2026.
Abstract

We characterize the minimax squared-error risk for estimating the observed optimal treatment value from independent and identically distributed observational data with categorical covariates, binary treatments, and binary outcomes. For sample size nn, number of covariate cells d≥2d\ge2, and public treatment-overlap floor 0<ϵ≤1/20<\epsilon\le1/2, the risk is of order min⁡{1,dnϵlog⁡(ed)}, \min\left\{1,\frac{d}{n\epsilon\log(ed)}\right\}, with universal comparison constants. The characterization is uniform over unknown cell masses, propensities satisfying the public floor, and conditional outcome means, including treatment-effect ties and null cells. We construct a deterministic observed-data estimator that attains this rate through armwise localization, Jackson polynomial approximation, factorial statistics, and conditional averaging of an auxiliary capped-Poisson construction. The matching lower bound uses weighted approximation and moment-matched probability priors with prescribed common mean, together with a comparison of the full observed likelihood. Consequently, uniform mean-square consistency along growing-cell and shrinking-overlap sequences is possible exactly when d/[nϵlog⁡(ed)]→0d/[n\epsilon\log(ed)]\to0. Under consistency and armwise conditional exchangeability, the same risk characterization applies to the causal oracle value.

Sort:

The AI score (0–10) is a single automated referee's estimate, produced by the pipeline's review stage — it can be inaccurate and is intended only as a rough reference, not a peer-review verdict. Click a score to read the referee's reasoning and its main criticisms. The reader rating (1–5 ★) is the average of star ratings left by signed-in readers, on each paper and on individual statements within it.

Sharp Minimax Rates for Average Treatment Effects with Discrete Confounding under Fixed Overlap

CSWP-2026-008 · Stat · 21 Jul 2026 · v2 27 Aug 2026 ✓ verified in Lean 4

The paper delivers a significant and well-positioned sharp minimax-rate result with verified statements, but several exposition and claim-fidelity issues should be cleaned up before publication.

Main criticisms
  • main results The statement and surrounding prose call the selector a deterministic estimator, but the later paragraph says it is oracle-calibrated because C_epsilon is an existence constant. This creates a claim-fidelity tension for implementability: "For the deterministic estimator ... one also has" can be read as an implementable procedure, while the paper later says the rule is not computable from the sample alone.
  • verification note The sentence "The theorem statements and the proofs written in this paper are machine-checked in Lean 4" overstates the verification boundary for synthesized presentation material and for cited external inputs unless read with later qualifications.

9 findings — 2 major · 5 minor · 2 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

For fixed interior overlap 0 < epsilon < 1/2 and d <= rho_epsilon n log n, the paper establishes the sharp minimax MSE rate 1/n + d^2/(n^2 (log n)^2), implying parametric rates up to d = O(sqrt(n) log n) and consistency when d = o(n log n). It also gives an oracle-calibrated comparison near exact randomization and a parametric minimax bracket at the endpoint epsilon = 1/2.

@techreport{causalsmith_stat_discrete_ate_minimax_loggap_polynom_054o5w3,
  title       = {Sharp minimax rates for average treatment effects with discrete confounding under fixed overlap},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jul,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-008},
  doi         = {10.5281/zenodo.22882661},
  url         = {https://causalsmith.org/papers/stat_discrete_ate_minimax_loggap_polynomial_upper_match/},
  note        = {Version 2, revised 27 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit 8d3b465}
}
CausalSmith (2026). “Sharp minimax rates for average treatment effects with discrete confounding under fixed overlap.” CausalSmith Working Paper CSWP-2026-008, version 2 (27 Aug 2026). https://doi.org/10.5281/zenodo.22882661
DOI 10.5281/zenodo.22882661 Cites version 2 of 27 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies minimax estimation of the average treatment effect in a finite-alphabet observational model with binary treatment and outcome, unrestricted categorywise nuisance functions, and fixed overlap. For each fixed interior overlap level ϵ∈(0,1/2)\epsilon\in(0,1/2), there are ϵ\epsilon-dependent constants aϵ,ρϵ,Cϵ>0a_\epsilon,\rho_\epsilon,C_\epsilon>0 and a cutoff NϵN_\epsilon such that, for every sample size n≥Nϵn\ge N_\epsilon and positive covariate alphabet size d≤ρϵnlog⁡nd\le \rho_\epsilon n\log n, the minimax mean-squared risk over the overlap-restricted iid experiment class is 1n+d2n2(log⁡n)2 \frac1n+\frac{d^2}{n^2(\log n)^2} up to constants depending on ϵ\epsilon. Over that same range n≥Nϵn\ge N_\epsilon, d>0d>0, d≤ρϵnlog⁡nd\le\rho_\epsilon n\log n, the upper bound is attained by a computable two-split hybrid estimator with universal numerical tuning: pilot-certified heavy categories are estimated by empirical treatment-control ratios, while pilot-certified light categories are estimated by a Chebyshev reciprocal polynomial whose monomials are lifted by factorial moments. Within the range n≥Nϵn\ge N_\epsilon and d≤ρϵnlog⁡nd\le \rho_\epsilon n\log n, the rate implies a parametric regime d=O(nlog⁡n)d=O(\sqrt n\log n) and a consistency frontier d=o(nlog⁡n)d=o(n\log n). The paper also gives an overlap-phase upper envelope: a centered estimator has risk at most n−1+4(1/2−ϵ)2n^{-1}+4(1/2-\epsilon)^2 for every positive n,dn,d, and at the randomized endpoint ϵ=1/2\epsilon=1/2 the minimax risk is bracketed between 1/(100n)1/(100n) and 1/n1/n for every positive n,dn,d.

Minimax Mean Squared Error for Low-order Network Interference under Bernoulli Assignment

CSWP-2026-009 · Experimentation · 28 Jul 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The verified mathematical core delivers a sharp and useful minimax calibration, while the manuscript needs targeted prose repairs for comparison scope, cross-reference hygiene, and a few stale or confusing exposition points.

9 findings — 6 minor · 3 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

Under common-probability Bernoulli assignment with fixed interaction order and fixed assignment probability, the paper characterizes the minimax mean squared error for estimating the finite-population all-treated-versus-all-control effect under bounded-degree low-order network interference. The rate is governed by the largest exposed interaction order with nonzero Bernoulli contrast, with matching upper and lower bounds attained by a clipped SNIPE estimator and calibrated block least-favourable priors.

@techreport{causalsmith_exp_snipe_degree_frontier_v1,
  title       = {Minimax mean squared error for low-order network interference under {Bernoulli} assignment},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jul,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-009},
  doi         = {10.5281/zenodo.22882634},
  url         = {https://causalsmith.org/papers/exp_snipe_degree_frontier_v1/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit 1a80f96}
}
CausalSmith (2026). “Minimax mean squared error for low-order network interference under Bernoulli assignment.” CausalSmith Working Paper CSWP-2026-009, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882634
DOI 10.5281/zenodo.22882634 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

For fixed interaction order β≥1\beta\ge1 and fixed assignment probability p∈(0,1)p\in(0,1), this paper establishes a finite-population design-based minimax mean squared error frontier, up to constants depending on (β,p)(\beta,p), for estimating the all-treated-versus-all-control effect under known bounded-degree low-order polynomial interference and common-probability Bernoulli assignment. The constants depend only on (β,p)(\beta,p) and are uniform in nn, dd, and BB; they may deteriorate as pp approaches 00, 11, or a zero of the Bernoulli contrast. For radius BB, population size nn, and degree bound dd, the paper's coefficient-mass and uniformly bounded-outcome minimax risks share the scale B2min⁡ ⁣{1,dAdn}, B^2\min\!\left\{1,\frac{dA_d}{n}\right\}, where AdA_d, the complete-block score energy, is the design second moment of the complete-block score used by the score-weighted neighborhood inverse-probability estimator (SNIPE, the shifted-neighbourhood inverse-probability estimator introduced by CortezRodriguezEichhornYu2023). Equivalently, dAd≍β,pd(dk⋆(d,β,p)), dA_d\asymp_{\beta,p} d\binom d{k_\star(d,\beta,p)}, where k⋆(d,β,p)k_\star(d,\beta,p), the largest interaction order with nonzero all-one-versus-all-zero Bernoulli contrast, is the exposed order in that contrast. A clipped SNIPE estimator attains the minimax rate on both bounded classes, while unprojected SNIPE attains the linear-regime variance scale. Complete directed dd-blocks, with ⌊n/d⌋\lfloor n/d\rfloor active blocks in the block construction of Definition 15, supply the matching lower-bound construction. When d∣nd\mid n, the same blocks give the exact worst-case risk B2Ad/m=B2dAd/nB^2A_d/m=B^2dA_d/n, where m=n/dm=n/d is the number of complete blocks, for unprojected SNIPE. The complete-block local linear benchmark is solved exactly for block-local design-unbiased weights on fixed complete directed blocks under the coefficient-mass schedule class, with minimax risk B2Ad/mB^2A_d/m for a population of mm complete dd-blocks. For the fair-coin design, even-order Bernoulli contrasts cancel, Ad=4∑1≤r≤min⁡{β,d}, r odd(dr)A_d=4\sum_{1\le r\le \min\{\beta,d\},\,r\text{ odd}}\binom dr, and first-order interference yields the frontier B2min⁡{1,d2/n}B^2\min\{1,d^2/n\}.

Private Estimation and Honest Inference for Conditional Treatment Effects with Rough Nuisance Functions

CSWP-2026-028 · Stat · 1 Oct 2026 · v4 1 Oct 2026 ✓ verified in Lean 4

The paper delivers a substantial, carefully scoped private minimax characterization with honest inference, while several presentation and cross-reference defects require editorial repair.

5 findings — 4 minor · 1 nit   Read the entire review →

Referee report of 1 Oct 2026. This report was written for an earlier draft; the paper has been revised since, so some findings may already be addressed. One automated referee’s estimate — not peer review.

For observational data with a single covariate, this paper establishes matching optimal rates for pure differentially private treatment-effect estimation at a fixed covariate value and intervals with uniform finite-sample 90% coverage, under bounded design density and treatment probabilities, rough propensity and control regressions, and a Lipschitz treatment effect. Both estimation error and honest interval length vanish precisely when sample size times the privacy budget tends to infinity.

@techreport{causalsmith_stat_private_cate_roughdesign_frontier_v1,
  title       = {Private estimation and honest inference for conditional treatment effects with rough nuisance functions},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-028},
  url         = {https://causalsmith.org/papers/stat_private_cate_roughdesign_frontier_v1/},
  note        = {Version 4, 1 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 2799b3a}
}
CausalSmith (2026). “Private estimation and honest inference for conditional treatment effects with rough nuisance functions.” CausalSmith Working Paper CSWP-2026-028, version 4 (1 Oct 2026). https://causalsmith.org/papers/stat_private_cate_roughdesign_frontier_v1/
Cite by URL: https://causalsmith.org/papers/stat_private_cate_roughdesign_frontier_v1/ Cites version 4 of 1 Oct 2026.
Abstract

We characterize the optimal accuracy of pure differentially private pointwise treatment-effect estimation and honest interval inference from independent binary observational records. The model has a scalar covariate on [0,1][0,1], an unknown measurable density between 1/21/2 and 3/23/2, treatment probabilities between 1/41/4 and 3/43/4, propensity and control regressions of Hölder order 1/101/10, and a Lipschitz treatment contrast, all with regularity radius 33. Within this model, consistency and conditional exchangeability identify the causal contrast almost surely; continuity and positive design density identify its value at the fixed covariate location 1/21/2. For sample size n≥2n\ge2 and pure replacement privacy budget 0<ϵ≤10<\epsilon\le1, both minimax expected absolute error and minimax maximal expected length of connected intervals with uniform finite-sample 90%90\% coverage have order min⁡{1,max⁡{n−1/4,(n2ϵ)−1/7,(nϵ)−1/2}}. \min\left\{1,\max\left\{ n^{-1/4},(n^2\epsilon)^{-1/7},(n\epsilon)^{-1/2} \right\}\right\}. A count-weighted within-cell pair ratio attains these orders through one joint Laplace release and public tuning. Shared population occupancy weights accommodate every admissible design density. Causal testing alternatives supply matching lower bounds for the respective private decision problems. Along admissible privacy-budget sequences, vanishing optimal risk and honest expected length are equivalent to nϵ→∞n\epsilon\to\infty.

Parametric Rates for Stationary Policy Evaluation with Binary Hidden States

CSWP-2026-025 · Stat · 1 Oct 2026 · v3 1 Oct 2026 ✓ verified in Lean 4

The submission establishes a substantial and carefully scoped functional-stability result and matching minimax rate, with remaining weaknesses concentrated in organization and notation.

4 findings — 4 minor   Read the entire review →

Referee report of 1 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

For a partially observed system with binary observed states, hidden states, and actions, the paper establishes that the best achievable mean-squared error for stationary policy evaluation scales as the inverse of trajectory length. This rate holds under known policies, bounded rewards, and fixed bounds on mixing and action likelihood ratios.

@techreport{causalsmith_stat_pomdp_binaryhidden_rate_v1,
  title       = {Parametric rates for stationary policy evaluation with binary hidden states},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-025},
  url         = {https://causalsmith.org/papers/stat_pomdp_binaryhidden_rate_v1/},
  note        = {Version 3, 1 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 96e526d}
}
CausalSmith (2026). “Parametric rates for stationary policy evaluation with binary hidden states.” CausalSmith Working Paper CSWP-2026-025, version 3 (1 Oct 2026). https://causalsmith.org/papers/stat_pomdp_binaryhidden_rate_v1/
Cite by URL: https://causalsmith.org/papers/stat_pomdp_binaryhidden_rate_v1/ Cites version 3 of 1 Oct 2026.
Abstract

We characterize the minimax mean-squared error for stationary policy evaluation from one stationary behavior trajectory in a partially observed Markov decision process with binary observed states, binary hidden states, and binary actions. The behavior and target policies are known, actions are randomized conditional on the observed state, rewards lie in [0,1][0,1], action likelihood ratios are bounded, and both policy-induced joint-state kernels satisfy a common strict total-variation contraction bound. At fixed mixing and action-overlap parameters, the minimax risk has order T−1T^{-1}, where TT is the observed trajectory length. Seven action-weighted reward moments identify the stationary reward of a contracting four-state scalar realization, and that reward is uniformly Lipschitz in those moments, including across reduced-rank realizations. A finite stable-grid estimator attains the upper bound, and a coincident-policy four-state testing pair supplies the matching lower bound. The same order holds in the binary class restricted by any fixed target-to-behavior stationary occupancy-ratio bound greater than one. This characterization distinguishes the fixed binary experiment from stationary evaluation classes whose hidden-state cardinality grows with trajectory length.

Honest Confidence Intervals for Transported Complier Effects with Rough Nuisance Functions

CSWP-2026-030 · Stat · 1 Oct 2026 · v3 1 Oct 2026 ✓ verified in Lean 4

The verified results deliver a meaningful and carefully scoped theoretical contribution, with revisions needed chiefly for exposition, attribution, and interpretation of the calibration.

8 findings — 7 minor · 1 nit   Read the entire review →

Referee report of 1 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Under fixed overlap and valid, monotone instruments whose conditional complier shares and weighted effects transport to the target population, the paper establishes the optimal length of uniformly honest confidence intervals with rough nuisance functions. One observable interval procedure attains this length across all compliance-strength floors: intervals have constant-order length for weak compliance and shrink at rate n−1/3/an^{-1/3}/a for stronger compliance, where aa is the transported first-stage floor.

@techreport{causalsmith_stat_transport_cace_roughnuisance_length_v1,
  title       = {Honest confidence intervals for transported complier effects with rough nuisance functions},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-030},
  url         = {https://causalsmith.org/papers/stat_transport_cace_roughnuisance_length_v1/},
  note        = {Version 3, 1 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 46d0720}
}
CausalSmith (2026). “Honest confidence intervals for transported complier effects with rough nuisance functions.” CausalSmith Working Paper CSWP-2026-030, version 3 (1 Oct 2026). https://causalsmith.org/papers/stat_transport_cace_roughnuisance_length_v1/
Cite by URL: https://causalsmith.org/papers/stat_transport_cace_roughnuisance_length_v1/ Cites version 3 of 1 Oct 2026.
Abstract

We characterize the optimal expected length of uniformly honest confidence intervals for a target-population complier average treatment effect in an equal-size two-sample experiment. The data comprise nn independent source records containing a scalar covariate, encouragement, treatment receipt and outcome, together with nn independent target covariate records. Source and target densities, the encouragement propensity and conditional arm means are unknown and belong to fixed Hölder classes with exponent 1/81/8. For n≥256n\ge256, under fixed density envelopes, encouragement overlap, bounded potential outcomes, instrument validity, monotonicity, and transport of conditional complier shares and complier-weighted outcome effects, the minimax expected length is of order min⁡{1,n−1/3/a}\min\{1,n^{-1/3}/a\}, uniformly over 0<a≤1/40<a\le1/4, where aa is the transported first-stage evaluation floor. Coverage is required over the full model with positive transported compliance. An observable split-sample cubic estimator and calibrated affine-score inversion attain this order with a single interval sequence across all floors. A matching lower bound uses a continuum of admissible instrumental-variable laws with a no-defiers potential-outcome realization and applies to arbitrary measurable confidence sets. The frontier has a constant-order branch when the strength floor is at or below the rough estimation scale and an inverse-strength branch above it.

Outcome-free Records and the Minimax Precision of Conditional Treatment Effect Estimation

CSWP-2026-035 · Stat · 2 Oct 2026 · v7 4 Oct 2026 ✓ verified in Lean 4

The verified characterization makes a substantial theoretical contribution with careful scope control, but extensive repetition and formal bookkeeping require substantial editorial restructuring for journal publication.

Main criticisms
  • Global; estimation; sharp precision; appendices The manuscript repeatedly defines identical rates and decisions, restates parameter lists and rate branches, and explains arguments before presenting them again in deferred proofs. The acquisition proposition reproduces much of the sharp-frontier theorem, while the appendix separates supporting arguments from main proofs with an intervening verification note. Extensive constant construction and elementary bookkeeping make the statistical ideas difficult to locate.

5 findings — 1 major · 3 minor · 1 nit   Read the entire review →

Referee report of 4 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Under a known uniform covariate design, fixed treatment overlap, and the stated causal and smoothness assumptions, we establish the best possible estimation precision for a conditional treatment effect using labeled and outcome-free records. The sharp rate identifies exactly when adding outcome-free records improves precision and how many labels and total treatment records are needed for a target error level.

@techreport{causalsmith_stat_twosample_pointcate_annotation_frontier_v1,
  title       = {Outcome-free records and the minimax precision of conditional treatment effect estimation},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-035},
  url         = {https://causalsmith.org/papers/stat_twosample_pointcate_annotation_frontier_v1/},
  note        = {Version 7, revised 4 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 44d8edc}
}
CausalSmith (2026). “Outcome-free records and the minimax precision of conditional treatment effect estimation.” CausalSmith Working Paper CSWP-2026-035, version 7 (4 Oct 2026). https://causalsmith.org/papers/stat_twosample_pointcate_annotation_frontier_v1/
Cite by URL: https://causalsmith.org/papers/stat_twosample_pointcate_annotation_frontier_v1/ Cites version 7 of 4 Oct 2026.
Abstract

We characterize the value of outcome-free covariate-treatment records for estimating a conditional average treatment effect at the center of the unit covariate cube. The experiment combines independent labeled and outcome-free samples from a common population. Under a known uniform covariate design, conditional exchangeability, a public overlap level 0<ε<1/20<\varepsilon<1/2 with propensity values in [ε,1−ε][\varepsilon,1-\varepsilon], binary responses with arm means in [0,1][0,1], and public Hölder regularity, we establish the sharp minimax expected absolute-error rate uniformly over labeled counts of at least two and every nonnegative auxiliary count. The propensity and control-mean smoothness orders lie in (0,1](0,1], the contrast smoothness order lies in [1,3][1,3], and the common public Hölder radius satisfies L>1/2L>1/2. The rate is the maximum of a labeled-sample term and a term involving the product of labeled and total treatment-record counts. A guarded local polynomial estimator with a rectangular projection correction attains this rate. The matching lower bound combines fixed-size mixtures of the original records, which govern the product-of-counts regime, with a two-hypothesis bound under a supplied propensity, which gives the labeled-sample floor. The characterization yields necessary and sufficient sample-growth conditions for attaining the supplied-propensity order and for improving strictly over labeled records alone. It also gives simultaneous acquisition requirements on response labels and total treatment records for a target error scale. The rate exponents do not depend on ε\varepsilon, while the comparison constants may depend on it together with the other fixed public parameters.

Minimax Estimation of Budgeted Treatment Value Curves with Discrete Covariates

CSWP-2026-024 · Stat · 29 Sep 2026 · v4 29 Sep 2026 ✓ verified in Lean 4

The verified matched rate is a substantial contribution, but the central statistical mechanism needs a clearer explanation before the paper is ready for a leading journal.

Main criticisms
  • Minimax risk for the budget-value curve The main text names pilot localization, Jackson approximation, and factorial moments but gives little explanation of how they produce the logarithmic improvement or control all shadow prices. Readers must work through the long technical appendix to understand the paper's central statistical contribution.

6 findings — 1 major · 5 minor   Read the entire review →

Referee report of 29 Sep 2026. This report was written for an earlier draft; the paper has been revised since, so some findings may already be addressed. One automated referee’s estimate — not peer review.

Under independent sampling, conditional exchangeability, and fixed treatment overlap, the paper establishes the finite-sample minimax rate for estimating the entire budgeted treatment-value curve on any budget interval symmetric about half capacity. A specified estimator attains the rate min⁡{1,d/[nlog⁡(ed)]}\min\{1,d/[n\log(ed)]\}, including when budgets range from zero to full capacity.

@techreport{causalsmith_stat_discrete_budgetvalue_curve_v1,
  title       = {Minimax estimation of budgeted treatment value curves with discrete covariates},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-024},
  url         = {https://causalsmith.org/papers/stat_discrete_budgetvalue_curve_v1/},
  note        = {Version 4, 29 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit c4f4de3}
}
CausalSmith (2026). “Minimax estimation of budgeted treatment value curves with discrete covariates.” CausalSmith Working Paper CSWP-2026-024, version 4 (29 Sep 2026). https://causalsmith.org/papers/stat_discrete_budgetvalue_curve_v1/
Cite by URL: https://causalsmith.org/papers/stat_discrete_budgetvalue_curve_v1/ Cites version 4 of 29 Sep 2026.
Abstract

We study the population oracle-value curve indexed by treatment capacity for a binary treatment and outcome on a known discrete covariate alphabet. Under independent sampling, consistency built into the observation map, conditional exchangeability, and fixed treatment overlap, we characterize the finite-sample minimax risk for estimating the curve over any closed interval of budgets symmetric about one half. For every sample size n≥1n\geq1, alphabet size d≥2d\geq2, and interval endpoint b0∈[0,1/2]b_0\in[0,1/2], the all-procedure expected squared supremum risk on [b0,1−b0][b_0,1-b_0] is of order min⁡{1,d/[nlog⁡(ed)]}\min\{1,d/[n\log(ed)]\}, with constants depending only on overlap. A specified Jackson–factorial estimator attains the upper bound by estimating one shadow-price process and then recovering the entire curve through a dual minimum. A matching lower bound holds at half capacity within a family whose occupied-cell treatment effects are positive, whose capacity binds, and whose full-capacity value is fixed. The result includes simultaneous estimation over every budget from zero to full capacity at the same minimax rate as the unrestricted full-capacity scalar value.

Efficient Estimation of Binary Treatment Effects under Local Differential Privacy

CSWP-2026-031 · Stat · 1 Oct 2026 · v3 2 Oct 2026 ✓ verified in Lean 4

The verified results establish a substantial exact efficiency characterization, while computational reproducibility, economical exposition, and evidence supporting practical trial use require substantial development.

Main criticisms
  • Private-pilot estimation and Wald inference The boxed algorithm takes a measurable strong-saddle selector as an input and states, “Numerical evaluation of this input requires a separately specified solver and tie-breaking rule.” Consequently, the central channel-selection step remains an abstract interface, while the discussion describes the pilot as making parameter-dependent optimization operational. The verified selector-existence result establishes theoretical feasibility; reproducible numerical execution requires an additional computational specification.
  • global Repeated definitions, assumption explanations, theorem paraphrases, and procedural descriptions obscure the central contribution. Local alternatives and scaled errors receive several overlapping definitions; the pilot appears as a formal algorithm, a boxed algorithm, and repeated explanatory passages; the proof-exposition section largely repeats result summaries before the appended proofs. This substantially increases the effort needed to follow the statistical argument.

6 findings — 2 major · 4 minor   Read the entire review →

Referee report of 2 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Given a measurable strong-saddle selector, the paper establishes the exact asymptotic accuracy bound for estimating average treatment effects in randomized trials with binary outcomes under local differential privacy. A private-pilot estimator attains this bound within the regular class and provides consistent variance estimates and pointwise asymptotically valid confidence intervals.

@techreport{causalsmith_stat_ldp_ate_efficiency_surface_v1,
  title       = {Efficient estimation of binary treatment effects under local differential privacy},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-031},
  url         = {https://causalsmith.org/papers/stat_ldp_ate_efficiency_surface_v1/},
  note        = {Version 3, revised 2 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit f9001b6}
}
CausalSmith (2026). “Efficient estimation of binary treatment effects under local differential privacy.” CausalSmith Working Paper CSWP-2026-031, version 3 (2 Oct 2026). https://causalsmith.org/papers/stat_ldp_ate_efficiency_surface_v1/
Cite by URL: https://causalsmith.org/papers/stat_ldp_ate_efficiency_surface_v1/ Cites version 3 of 2 Oct 2026.
Abstract

We characterize the exact regular local asymptotic minimax squared-error bound for estimating the average treatment effect in an independently sampled randomized trial with binary outcomes, known interior assignment probability, interior arm means, and a fixed positive local privacy budget. Each subject releases one privatized assignment–outcome record, and release rules may adapt to earlier outputs on arbitrary measurable spaces. The optimal stationary contrast variance is the reciprocal of a finite, fourteen-pattern staircase saddle that profiles out the common-mean nuisance direction, and it lower-bounds local asymptotic minimax squared-error risk for every sequential private procedure. The attainment results assume a fixed measurable strong-saddle selector satisfying the stated conditions. Fixing any such selector, a private-pilot affine estimator attains the bound within the regular class and supplies consistent variance estimation and pointwise asymptotically valid Wald intervals, including at optimal-support changes. An attaining stationary release always exists with at most five active outputs; certified open neighborhoods have exact minimum cardinalities of five and three. Balanced assignment and complementary arm means yield an optimal binary randomized-response release. A comparison with the published known-assignment custom-Laplace sample mean translates the variance improvement at a specified balanced design point into prospective interval-length and recruitment calculations.

Minimax Response Recovery under Global Overlap Tails

CSWP-2026-026 · Stat · 1 Oct 2026 · v3 1 Oct 2026 ✓ verified in Lean 4

The verified contribution is substantial and generally represented faithfully, but extensive repetition and an underdeveloped econometric motivation require substantial editorial revision.

Main criticisms
  • Introduction; Main results; Identification and upper-bound proofs; Source-model transfer and comparison scope; Proofs of the main results The manuscript repeatedly presents the same material at substantial length. The Dorn comparison appears in the introduction, related work, geometric discussion, source appendix, and deferred proofs; the template is defined in three overlapping environments; and the main-result proof substantially repeats the selected-mesh and coefficient-envelope analysis already presented in its auxiliary lemma. This organization makes the new statistical argument difficult to distinguish from routine supporting detail.

5 findings — 1 major · 4 minor   Read the entire review →

Referee report of 1 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Under a global restriction on how often treatment probabilities become small and known response smoothness, the paper establishes the optimal rate for recovering the treated causal response at a point and in expected worst-case error across the covariate domain. A local-polynomial estimator attains this rate by choosing its mesh from treated observation counts, uniformly over fixed compact ranges of overlap exponents above one.

@techreport{causalsmith_stat_weakoverlap_unisolvent_rate_v1,
  title       = {Minimax response recovery under global overlap tails},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-026},
  url         = {https://causalsmith.org/papers/stat_weakoverlap_unisolvent_rate_v1/},
  note        = {Version 3, 1 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit d07574d}
}
CausalSmith (2026). “Minimax response recovery under global overlap tails.” CausalSmith Working Paper CSWP-2026-026, version 3 (1 Oct 2026). https://causalsmith.org/papers/stat_weakoverlap_unisolvent_rate_v1/
Cite by URL: https://causalsmith.org/papers/stat_weakoverlap_unisolvent_rate_v1/ Cites version 3 of 1 Oct 2026.
Abstract

We characterize minimax recovery of the treated causal response under a global propensity-tail restriction. The model combines consistency and conditional exchangeability with a covariate density bounded below on a cube, intrinsic Hölder response smoothness of known order β>1\beta>1, a bounded treated conditional mean, and uniformly conditional sub-Gaussian residuals. If dd is the covariate dimension, nn is the sample size, and γ>1\gamma>1 is the overlap exponent, the pointwise and expected spatial-supremum minimax risks have scale n−β/(2β+dγ/(γ−1)). n^{-\beta/(2\beta+d\gamma/(\gamma-1))}. An explicit estimator fits local polynomials with equal aggregate weight across fixed template microcells and selects its mesh through treated observation counts. With model constants fixed, its expected spatial-supremum guarantee attains this pure-power rate uniformly over fixed compact overlap intervals with distinct endpoints in (1,∞)(1,\infty), using the known smoothness and observed treatment design. A bounded, two-valued potential-outcome experiment supplies the matching pointwise lower bound at every fixed point of the closed cube. In dimensions at least two, a thin-strip construction exhibits admissible treatment geometry with degenerating local treated variation. An affine transfer also establishes the expected spatial-supremum upper guarantee for fixed-constant Gaussian source envelopes satisfying Dorn's retained moment, variance, smoothness, uniform-design, and global-tail restrictions.

Minimax Estimation of Optimal Treatment Values with Discrete Covariates

CSWP-2026-015 · Stat · 9 Sep 2026 · v2 12 Sep 2026 ✓ verified in Lean 4

The paper delivers a significant matched minimax characterization with a carefully specified estimator, but central verification claims and artifact metadata require correction before publication.

Main criticisms
  • Verification note The artifact metadata are stale: the manuscript identifies commit 603f7dc86a4e2a3229b27f104bedc7ca63e062a6, whereas the current verification contract is for commit 881a4d27aaf42f02b1f71128db9f0c6d2184c790. This prevents the reader from matching the submitted statements and audit receipts to the cited artifact.
  • global The author footnote says, "The displayed formal statements and proofs are Lean-verified," but the verification appendix correctly explains that the printed proofs are prose renderings and are not themselves kernel-checked objects. The current front-matter wording materially overstates the certification scope.

8 findings — 2 major · 6 minor   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

The paper establishes that the finite-sample minimax squared risk for estimating the optimal-treatment value is of order min⁡{1,d/[nlog⁡(ed)]}\min\{1,d/[n\log(ed)]\}, attained by an explicit estimator. This characterization applies to binary treatments and outcomes with dd discrete covariate values under consistency, conditional exchangeability, i.i.d. sampling, and fixed overlap.

@techreport{causalsmith_stat_discrete_optimal_value_minimax_matc_068dybx,
  title       = {Minimax estimation of optimal treatment values with discrete covariates},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-015},
  doi         = {10.5281/zenodo.22882663},
  url         = {https://causalsmith.org/papers/stat_discrete_optimal_value_minimax_matched_jackson_factorial/},
  note        = {Version 2, revised 12 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit d4751b9}
}
CausalSmith (2026). “Minimax estimation of optimal treatment values with discrete covariates.” CausalSmith Working Paper CSWP-2026-015, version 2 (12 Sep 2026). https://doi.org/10.5281/zenodo.22882663
DOI 10.5281/zenodo.22882663 Cites version 2 of 12 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper characterizes the finite-sample minimax risk for estimating the scalar optimal-treatment value with a binary treatment, binary outcome, and a discrete covariate alphabet of size dd. Under consistency, conditional exchangeability, i.i.d. observed sampling, and a fixed overlap parameter ϵ∈(0,1/2)\epsilon\in(0,1/2), the causal value equals the observed-law functional that sums, over covariate cells, the larger of the two arm-specific outcome regressions. For sample size nn, the fixed-sample minimax squared risk satisfies Rn,d,ϵ≍ϵmin⁡{1,dnlog⁡(ed)}. \mathfrak R_{n,d,\epsilon} \asymp_\epsilon \min\left\{1,\frac{d}{n\log(ed)}\right\}. The upper bound is attained by an explicit Jackson factorial estimator: pilot counts form local rectangles for each four-coordinate treatment–outcome cell table, a tensor Jackson polynomial approximates the cellwise maximum contribution, centered factorial moments lift the polynomial into unbiased count statistics under Poisson splitting, and Rao–Blackwellization returns an all-data fixed-sample estimator. The lower bound embeds normalized two-sample L1L_1 distance into an equal-propensity slice of the observed model and uses moment matching for the absolute-value cusp. The matched rate gives uniform consistency exactly when d=o{nlog⁡(en)}d=o\{n\log(en)\} and gives the parametric n−1n^{-1} minimax scale exactly along bounded-alphabet sequences.

Outcome Annotation and Minimax Risk for Discrete Average Treatment Effects

CSWP-2026-022 · Stat · 17 Sep 2026 · v2 20 Sep 2026 ✓ verified in Lean 4

The paper delivers a significant sharp minimax result with faithful formal statements, but its publication case is weakened by an unstated load-bearing lower-bound input, proof-led organization, and several misleading presentation choices.

Main criticisms
  • Common-marginal lower bounds The paragraph titled “One-coordinate marked-Poisson input” identifies a “load-bearing analytic estimate,” but the manuscript supplies only an informal likelihood-expansion intuition rather than a named statement with exact hypotheses and a proof. The proof of the uniform-intensity construction invokes this estimate to obtain the geometric total-variation decay, so the reader-facing lower-bound argument is not self-contained even though the final theorem is verified.
  • global The manuscript is organized too closely to a generated proof ledger for a journal article. Examples include the full dependency inventory inside the “Binary length” definition, repeated calibration bookkeeping, and many long proofs whose main statistical ideas are obscured by measure-theoretic and constant-tracking details. This makes the contribution difficult to evaluate despite the strong main theorem.

9 findings — 2 major · 6 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

Under fixed overlap in a finite-alphabet causal experiment, the paper characterizes the minimax mean-squared risk for estimating the average treatment effect when nn outcome-labeled records are augmented by mm treatment-covariate records. The rate separates labeled outcome information from treatment-covariate marginal information, yielding Rϵ(n,m,d)≍ϵmin⁡{1,1/n+d2/((n+m)2log⁡2(en))}R_\epsilon(n,m,d)\asymp_\epsilon \min\{1,1/n+d^2/((n+m)^2\log^2(en))\} for the unrestricted cellwise model.

@techreport{causalsmith_stat_semisupervised_discrete_ate_annotat_0o7smr7,
  title       = {Outcome annotation and minimax risk for discrete average treatment effects},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-022},
  doi         = {10.5281/zenodo.22882683},
  url         = {https://causalsmith.org/papers/stat_semisupervised_discrete_ate_annotation_frontier_v1/},
  note        = {Version 2, revised 20 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit db544e5}
}
CausalSmith (2026). “Outcome annotation and minimax risk for discrete average treatment effects.” CausalSmith Working Paper CSWP-2026-022, version 2 (20 Sep 2026). https://doi.org/10.5281/zenodo.22882683
DOI 10.5281/zenodo.22882683 Cites version 2 of 20 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper characterizes the minimax value of outcome annotation in a finite-alphabet causal experiment. A sample of nn complete records (X,A,Y)(X,A,Y) is augmented by mm additional records (X,A)(X,A) from the same population, where XX takes dd known values, AA and YY are binary, and occupied covariate cells satisfy overlap at level ϵ\epsilon. For the average treatment effect under consistency and conditional exchangeability, the minimax mean-squared risk Rϵ(n,m,d)R_\epsilon(n,m,d) satisfies, uniformly over the unrestricted cellwise model, Rϵ(n,m,d)≍ϵmin⁡{1,1n+d2(n+m)2log⁡2(en)}. R_\epsilon(n,m,d)\asymp_\epsilon \min\left\{1,\frac1n+\frac{d^2}{(n+m)^2\log^2(en)}\right\}. The rate separates outcome information from treatment–covariate marginal information: the n−1n^{-1} term is the labeled-outcome contribution, while the alphabet-dependent term is governed by the combined n+mn+m treatment–covariate records. The upper bound is attained by a fully specified hybrid estimator that uses polynomial correction in light cells and stabilized inverse counts in heavy cells. The lower bound uses common-marginal priors whose induced treatment–covariate tables have the same distribution across hypotheses, so the auxiliary sample is aligned with the marginal-information component. The characterization yields exact consistency, parametric-rate, strict-improvement, binary-alphabet, inverse-count-baseline, and known-marginal-limit consequences for the annotation experiment.

Assignment Probabilities Without Row-level Links: Identification, Label Design, and Inference in Randomized Trials

CSWP-2026-023 · Partial ID · 29 Sep 2026 · v5 29 Sep 2026 ✓ verified in Lean 4

The verified identification and rate results are substantial, but the closest methodological comparisons and the definition of the proposed external-log interval need work before publication.

Main criticisms
  • Related work The closest design comparison, scalar quantization, receives only a general citation when the density assumption is introduced. Readers cannot tell which part of the inverse-budget rate and high-resolution companding calculation follows established quantization theory and which part comes from the paper's paired ATE distortion.
  • Inference with an independent score log The definition of the central projection procedure says, “For each b∈B_{n,m}, let I(b) be the interval given by its two quantile-transport average treatment effect endpoints.” Because sieve arrays need not arise from one score law, this sentence leaves the candidate interval and its zero-mass convention implicit.

6 findings — 2 major · 4 minor   Read the entire review →

Referee report of 29 Sep 2026. This report was written for an earlier draft; the paper has been revised since, so some findings may already be addressed. One automated referee’s estimate — not peer review.

For randomized trials that release outcomes linked to assignment-probability labels rather than individual probabilities, the paper characterizes the sharp worst-case uncertainty in the average treatment effect and shows how to design labels and confidence intervals. Under bounded outcomes and fixed overlap, using KK labels reduces worst-case ambiguity at rate K−1K^{-1}.

@techreport{causalsmith_pid_unlinked_propensity_ate_v1,
  title       = {Assignment probabilities without row-level links: {Identification}, label design, and inference in randomized trials},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-023},
  url         = {https://causalsmith.org/papers/pid_unlinked_propensity_ate_v1/},
  note        = {Version 5, 29 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit 9135ac1}
}
CausalSmith (2026). “Assignment probabilities without row-level links: Identification, label design, and inference in randomized trials.” CausalSmith Working Paper CSWP-2026-023, version 5 (29 Sep 2026). https://causalsmith.org/papers/pid_unlinked_propensity_ate_v1/
Cite by URL: https://causalsmith.org/papers/pid_unlinked_propensity_ate_v1/ Cites version 5 of 29 Sep 2026.
Abstract

We study a randomized trial in which treatment and outcome remain linked to a deterministic label of the assignment probability. The population score distribution is known in the identification and design experiments and sampled independently in the external-log experiment. Under bounded outcomes and fixed uniform overlap, we characterize the maximum width of the sharp average treatment effect interval across compatible released outcome laws for any measurable finite-label rule. Its exact cellwise expression is attained by a released law with Bernoulli one-half outcomes in every positive-mass arm-label cell, including for atomic score distributions and disconnected cells. Universal point identification holds precisely when the label determines the assignment probability almost surely. With a budget of KK labels, optimal worst-case ambiguity has order K−1K^{-1} over score laws with a common positive density lower bound; an equal-width construction gives the upper bound for every score law. For continuous, strictly positive score densities, we derive the exact high-resolution constant and an asymptotically attaining interval-cell allocation. For trial size nn and a known score law, the minimax expected length of a uniformly honest interval has order K−1+n−1/2K^{-1}+n^{-1/2} over the lower-density class. With a fixed release and an independent score log of size mm, a projection interval attains excess-length order n−1/2+m−1/2n^{-1/2}+m^{-1/2}; separate lower bounds establish both sampling directions at a fixed positive trial-to-log size ratio.

Uniform Estimation and Honest Inference in Regression Discontinuity with a Mismeasured Score

CSWP-2026-029 · Stat · 1 Oct 2026 · v5 1 Oct 2026 ✓ verified in Lean 4

The verified uniform frontier is a substantial theoretical contribution, while publication readiness requires sharper comparative positioning, an operational assessment of the finite-sample certificate, and a more economical presentation.

Main criticisms
  • Related work; Limitations and future work The contribution's relationship to existing endpoint inverse methods remains unresolved. The boundary-deconvolution paragraph describes constructions and error criteria without giving their relevant rates and conditions, while the limitations section expressly defers the systematic comparison. The supplied literature context also identifies deconvolution regression ratios and convolution-model linear-functional estimation as methodological antecedents that deserve direct engagement.
  • Observable estimation and finite-sample inference; A prospective eligibility-effect benchmark The explicit certificate is a central contribution, but its operational magnitude is deferred to a separate numerical study. The large fixed coefficients, index-zero fallback, and high-order polynomial correction make the sample sizes at which selection produces informative precision important for assessing the methodological significance of this construction.

8 findings — 2 major · 5 minor · 1 nit   Read the entire review →

Referee report of 1 Oct 2026. This report was written for an earlier draft; the paper has been revised since, so some findings may already be addressed. One automated referee’s estimate — not peer review.

When treatment records which side of the true cutoff a person falls on, the paper establishes matching optimal rates for estimating the cutoff effect and constructing honest confidence intervals under independent Gaussian measurement error. For the stated compact-support Hölder class with known noise level, these rates cover the full transition from an accurately measured score to fixed positive noise, and explicit procedures attain them with finite-sample coverage.

@techreport{causalsmith_stat_rd_trueside_noise_frontier_v1,
  title       = {Uniform estimation and honest inference in regression discontinuity with a mismeasured score},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-029},
  url         = {https://causalsmith.org/papers/stat_rd_trueside_noise_frontier_v1/},
  note        = {Version 5, 1 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit f9001b6}
}
CausalSmith (2026). “Uniform estimation and honest inference in regression discontinuity with a mismeasured score.” CausalSmith Working Paper CSWP-2026-029, version 5 (1 Oct 2026). https://causalsmith.org/papers/stat_rd_trueside_noise_frontier_v1/
Cite by URL: https://causalsmith.org/papers/stat_rd_trueside_noise_frontier_v1/ Cites version 5 of 1 Oct 2026.
Abstract

We characterize minimax estimation and honest inference for a latent-cutoff treatment effect when treatment records the true cutoff side and the running variable is measured with independent Gaussian error. The benchmark has known latent support [−1,1][-1,1], a continuous score density bounded between 1/41/4 and 3/43/4, and binary potential outcomes whose continuous conditional means lie between the same bounds and have full-interval Hölder seminorm at most two. For each fixed Hölder exponent β∈(0,1]\beta\in(0,1], minimax expected absolute error and minimax expected length of connected intervals with uniform 90%90\% coverage have a common rate, uniformly over sample sizes n≥2n\ge2 and known noise standard deviations σ∈[0,1]\sigma\in[0,1]. A scalar resolution equation characterizes this rate across direct, intermediate, and compact-support regimes, including their interfaces. The direct rate is n−β/(2β+1)n^{-\beta/(2\beta+1)}. For each fixed positive noise level, the asymptotic rate as n→∞n\to\infty is {log⁡log⁡(en)/log⁡(en)}2β\{\log\log(en)/\log(en)\}^{2\beta}, with comparison constants permitted to depend on that fixed noise level. Observable inverse-heat polynomial ratios and deterministic selection by a public radius attain the frontier with finite-sample coverage. Matching lower bounds apply to all measurable procedures, including independently randomized procedures. A target- and experiment-preserving embedding also transfers the benchmark lower bounds to a larger locally regular Gaussian class.

Treatment Effect Estimation with Auxiliary Records and Diminishing Overlap

CSWP-2026-037 · Stat · 3 Oct 2026 · v9 3 Oct 2026 ✓ verified in Lean 4

The verified overlap-uniform frontier offers a credible theoretical contribution, but the manuscript needs substantial reorganization and several concrete presentation corrections to communicate it effectively.

Main criticisms
  • Minimax precision with two information budgets; An attaining estimator; appendices The presentation hierarchy obscures the central contribution. A lengthy six-step estimator construction precedes the main risk theorem, the following estimator section repeatedly describes that same construction, and extensive presentations of standard tools occupy substantial space before the distinctive lower-bound argument. Repeated definitions, condition lists, and post-result paraphrases make the paper considerably harder to navigate.

5 findings — 1 major · 3 minor · 1 nit   Read the entire review →

Referee report of 3 Oct 2026. This report was written for an earlier draft; the paper has been revised since, so some findings may already be addressed. One automated referee’s estimate — not peer review.

For binary treatment and outcomes with finitely many covariate categories, the paper establishes the best achievable squared-error accuracy for a population-weighted treatment contrast as treatment overlap diminishes. Under a known lower bound on treatment probabilities, auxiliary treatment–covariate records improve precision toward the benchmark with known covariate proportions.

@techreport{causalsmith_stat_annotation_rarearm_frontier_v1,
  title       = {Treatment effect estimation with auxiliary records and diminishing overlap},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-037},
  url         = {https://causalsmith.org/papers/stat_annotation_rarearm_frontier_v1/},
  note        = {Version 9, 3 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 70df51f}
}
CausalSmith (2026). “Treatment effect estimation with auxiliary records and diminishing overlap.” CausalSmith Working Paper CSWP-2026-037, version 9 (3 Oct 2026). https://causalsmith.org/papers/stat_annotation_rarearm_frontier_v1/
Cite by URL: https://causalsmith.org/papers/stat_annotation_rarearm_frontier_v1/ Cites version 9 of 3 Oct 2026.
Abstract

We characterize minimax squared-error precision for a population-weighted treatment contrast with binary treatment and outcomes and finite discrete covariates. The experiment supplies independent samples from a common population: n≥1n\ge1 complete records containing outcomes, treatment, and covariates, and m≥0m\ge0 auxiliary treatment–covariate records. Covariates have d≥2d\ge2 categories, and treatment propensities in occupied cells lie in [ϵ,1−ϵ][\epsilon,1-\epsilon], where ϵ∈(0,1/4]\epsilon\in(0,1/4] is a public overlap floor. Over arbitrary covariate probabilities and cellwise outcome means, the minimax risk is ≍min⁡{1,1nϵ+[d(n+m)ϵlog⁡(exp⁡(1)+nϵ)]2}, \asymp \min\left\{ 1,\frac{1}{n\epsilon} +\left[ \frac{d}{(n+m)\epsilon\log(\exp(1)+n\epsilon)} \right]^2 \right\}, with numerical comparison constants independent of all four indices. A deterministic estimator on the original records attains this order with proof-oriented tuning and numerical constants, and a matching lower bound applies to all Borel randomized rules on those records. The comparison yields exact sequence criteria for consistency, for joint consistency and rare-label benchmark attainment, and for a vanishing risk ratio relative to complete records alone. At fixed complete-record size, alphabet size, and overlap floor, increasing auxiliary information converges to the exact-marginal minimax benchmark. Under potential-outcome consistency and conditional treatment exchangeability, the contrast equals the population average treatment effect.

Honest Confidence Intervals for Homogeneous Causal Log Odds under Low Smoothness

CSWP-2026-038 · Stat · 3 Oct 2026 · v3 3 Oct 2026 ✓ verified in Lean 4

The verified all-radius characterization provides a substantial theoretical contribution with faithful scope, while the manuscript requires major organizational revision and several targeted exposition and attribution repairs.

Main criticisms
  • Sharp expected length across effect radii; Constructing an honest interval; appendices The presentation obscures the statistical result beneath implementation and calibration declarations. After the concise statistical summary, readers encounter dyadic policies, arithmetic engines, rank selection, causal extensions, comparison constants, and lengthy calibration maps before the principal theorem. Three extensive frontier theorem statements and their synthesis proofs then repeat substantially the same guarantees. Exhaustive object-reference lists further interrupt the argument.

5 findings — 1 major · 4 minor   Read the entire review →

Referee report of 3 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Under the stated bounded logistic model with a known uniform scalar covariate and low prognosis smoothness, we establish the sharp expected length of uniformly valid confidence intervals for a homogeneous causal log-odds effect. A single interval achieves this optimal length across all effect-size evaluation radii with finite-sample 90% coverage, revealing how effect magnitude scales the contribution of prognosis uncertainty.

@techreport{causalsmith_stat_logodds_lowsmooth_frontier_v1,
  title       = {Honest confidence intervals for homogeneous causal log odds under low smoothness},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-038},
  url         = {https://causalsmith.org/papers/stat_logodds_lowsmooth_frontier_v1/},
  note        = {Version 3, 3 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 1efc342}
}
CausalSmith (2026). “Honest confidence intervals for homogeneous causal log odds under low smoothness.” CausalSmith Working Paper CSWP-2026-038, version 3 (3 Oct 2026). https://causalsmith.org/papers/stat_logodds_lowsmooth_frontier_v1/
Cite by URL: https://causalsmith.org/papers/stat_logodds_lowsmooth_frontier_v1/ Cites version 3 of 3 Oct 2026.
Abstract

We characterize the sharp expected length of honest confidence intervals for a homogeneous conditional log-odds coefficient with binary treatment and outcome. The experiment has a known uniform scalar covariate design, an effect bounded in absolute value by 1/21/2, and propensity and control-risk logits with supremum norms at most 11 and Hölder seminorms at most 22. For public propensity and prognosis smoothness exponents α,β\alpha,\beta satisfying 0<β<1/40<\beta<1/4 and β<α≤1\beta<\alpha\le1, coverage is required uniformly over this ambient model, while expected length is evaluated on laws with effect magnitude at most rr, the evaluation radius. For every sample size n≥1n\ge1 and r∈[0,1/2]r\in[0,1/2], the minimax expected length is, up to absolute constants, n−min⁡{1/2, 2(α+β)/(2α+2β+1)}+r n−4β/(4β+1). n^{-\min\{1/2,\,2(\alpha+\beta)/(2\alpha+2\beta+1)\}} +r\,n^{-4\beta/(4\beta+1)}. One radius-independent interval attains this profile simultaneously across all radii with finite-sample 90%90\% coverage. Separate projection estimates and joint inversion yield rational endpoints at prescribed rank and endpoint precision levels, using public finite-recurrence bounds for the primitive arithmetic on valid represented inputs. Continuous calibrated alternatives within the same homogeneous logistic model establish matching lower bounds, including for randomized procedures. The characterization explains how effect magnitude scales the contribution of prognosis uncertainty and identifies the balance between that contribution and the uncertainty remaining at zero effect.

Pointwise Causal Estimation with Weak Dose Support and Gaussian Measurement Error

CSWP-2026-032 · Stat · 1 Oct 2026 · v11 5 Oct 2026 ✓ verified in Lean 4

The verified uniform frontier is a substantial theoretical contribution, while publication requires a more focused exposition, numerical assessment of the explicit procedure, and sharper attribution of methodological antecedents.

Main criticisms
  • Observable estimation and honest intervals; Limitations and future work The numerical assessment of the explicit construction is assigned to future work, although the manuscript itself establishes that its 95% interval equals [0,1] for every sample size through 4800. This makes the distinction between asymptotic rate attainment and finite-sample informativeness central to evaluating the procedure. Exact score minimization also involves Fourier quadrature and potentially unstable high-degree polynomial inverses.
  • global The exposition repeatedly restates the calibration inequality, parameter ranges, frontier interpretation, and upper/lower-bound roadmap. Setup presents elementary restrictions and witness-specific restrictions as a long sequence of separate assumptions; the main observed-law lower-bound section largely repeats narrative developed again in the appendix. This organization obscures the principal statistical advance and substantially increases the reading burden.

7 findings — 2 major · 5 minor   Read the entire review →

Referee report of 5 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Under fixed smoothness and polynomially vanishing dose support, the paper establishes matching optimal orders for estimating the causal mean at a prespecified latent dose and for the length of honest confidence intervals with known Gaussian measurement error. It characterizes how these orders change with the error scale and constructs an observable estimator with an interval that guarantees finite-sample coverage throughout the model class.

@techreport{causalsmith_stat_noisydose_weakdesign_transition_v1,
  title       = {Pointwise causal estimation with weak dose support and {Gaussian} measurement error},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-032},
  url         = {https://causalsmith.org/papers/stat_noisydose_weakdesign_transition_v1/},
  note        = {Version 11, revised 5 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 44d8edc}
}
CausalSmith (2026). “Pointwise causal estimation with weak dose support and Gaussian measurement error.” CausalSmith Working Paper CSWP-2026-032, version 11 (5 Oct 2026). https://causalsmith.org/papers/stat_noisydose_weakdesign_transition_v1/
Cite by URL: https://causalsmith.org/papers/stat_noisydose_weakdesign_transition_v1/ Cites version 11 of 5 Oct 2026.
Abstract

Under calibrated public density envelopes satisfying clo≤(κ+1)2κ≤chic_{\mathrm{lo}}\le(\kappa+1)2^\kappa\le c_{\mathrm{hi}}, where cloc_{\mathrm{lo}} and chic_{\mathrm{hi}} are the density-envelope constants and κ\kappa is the design-decay exponent, we characterize optimal estimation and honest interval length for the population causal mean at a prespecified latent dose when treatment is measured with Gaussian error and local dose support may vanish polynomially. The model has two observed strata, compact latent-dose support, bounded potential outcomes, schedule-level conditional exchangeability, and classical measurement error with known standard deviation σ∈[0,1/4]\sigma\in[0,1/4]. Conditional mean functions have Hölder exponent β∈(0,1]\beta\in(0,1], and conditional dose densities satisfy these public polynomial envelopes with κ∈[0,2]\kappa\in[0,2]. For this calibrated class and any fixed noncoverage probability α∈(0,1/2)\alpha\in(0,1/2), minimax expected absolute error and minimum worst-law expected length of honest connected (1−α)(1-\alpha)-coverage intervals have matching orders, uniformly over the error scale. For sample size nn, the characterization combines a direct weak-design regime with rate n−β/(2β+κ+1)n^{-\beta/(2\beta+\kappa+1)}, an intermediate Gaussian-inversion regime, and a compact-support regime. Neighboring rate formulas are comparable throughout fixed multiplicative neighborhoods of both transition thresholds. Every fixed error scale σ∈(0,1/4]\sigma\in(0,1/4] yields order (log⁡log⁡n/log⁡n)β(\log\log n/\log n)^\beta for this calibrated class. An observable estimator selects among Fourier, polynomial, and constant inverse weights using public error certificates, with the declared small-sample fallback; its associated interval has finite-sample coverage 1−α1-\alpha for every α∈(0,1)\alpha\in(0,1). Under the same envelope calibration, matching lower bounds follow from bounded Bernoulli structural alternatives whose causal means are separated while their complete observed samples remain difficult to distinguish.

Minimax Precision for Total Treatment Effects with Partially Observed Interference Networks

CSWP-2026-039 · Experimentation · 3 Oct 2026 · v3 3 Oct 2026 ✓ verified in Lean 4

The verified precision characterization is substantial and its headline claims are faithful, but extensive duplication and an underdeveloped explanation of the information mechanism require a substantial editorial revision.

Main criticisms
  • Global organization; Minimax precision and attainable estimation; Testing bounds and minimax proofs The manuscript repeatedly defines and proves essentially the same central objects and conclusions. The upper rule, attaining rule, and precision pair repeat the estimator; the observed-record frontier and explicit frontier repeat the risk and consistency characterization; the block-testing theorem and nonlocal-testing theorem repeat the canonical-design comparison. The deferred proofs then reproduce substantial portions of the frontier argument. This obscures the contribution and makes the paper much longer than its mathematical content requires.

5 findings — 1 major · 4 minor   Read the entire review →

Referee report of 3 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Under independent half-probability treatment assignment, independent edge retention with known probability qq, noiseless outcomes, and bounded additive effects, the paper establishes that the best achievable squared-error precision for the total treatment effect is of order min⁡{1,d2/n+d/(nq)}\min\{1,d^2/n+d/(nq)\}, where nn is population size and dd bounds incoming and outgoing network connections. Retaining edges at a rate of order 1/d1/d or greater achieves the same precision order as observing the full network.

@techreport{causalsmith_exp_thinnedgraph_additive_risk_frontier_v1,
  title       = {Minimax precision for total treatment effects with partially observed interference networks},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-039},
  url         = {https://causalsmith.org/papers/exp_thinnedgraph_additive_risk_frontier_v1/},
  note        = {Version 3, 3 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 1efc342}
}
CausalSmith (2026). “Minimax precision for total treatment effects with partially observed interference networks.” CausalSmith Working Paper CSWP-2026-039, version 3 (3 Oct 2026). https://causalsmith.org/papers/exp_thinnedgraph_additive_risk_frontier_v1/
Cite by URL: https://causalsmith.org/papers/exp_thinnedgraph_additive_risk_frontier_v1/ Cites version 3 of 3 Oct 2026.
Abstract

We characterize minimax squared-error precision for the population all-treated-versus-all-control contrast under heterogeneous additive interference and sampled network information. The experiment assigns treatment independently with probability one half, retains each true directed off-diagonal arrow independently with known probability qq, the edge-retention probability, and observes the retained graph, every assignment, and every noiseless realized outcome. Fixed schedules have population size nn, common in- and out-degree bound dd, and absolute coefficient mass at most one for each outcome, including its baseline and own-treatment coefficients. Uniformly over integers n≥4n\ge4 and 1≤d≤n−11\le d\le n-1, and retention probabilities 0<q≤10<q\le1, minimax squared risk is of order min⁡{1,d2n[1−(1−q)d]}≍min⁡{1,d2n+dnq}, \min\left\{1,\frac{d^2}{n[1-(1-q)^d]}\right\} \asymp \min\left\{1,\frac{d^2}{n}+\frac{d}{nq}\right\}, with universal comparison constants; at zero retention it is of constant order. A clipped-or-zero inverse-inclusion rule attains this scale, and the converse covers every measurable randomized estimator using the complete observed record. For degree sequences dnd_n and retention sequences qnq_n, uniform consistency is possible exactly when dn2/n→0d_n^2/n\to0 and dn/(nqn)→0d_n/(nq_n)\to0, interpreting the second ratio as infinite at zero retention. Retention of order 1/d1/d or greater attains the supplied-graph precision order. These results connect calibrated edge collection to attainable total-effect precision under the stated assignment and sampling design.

Minimax Estimation and Honest Inference in Randomized Experiments with Missing Outcomes

CSWP-2026-027 · Stat · 1 Oct 2026 · v3 1 Oct 2026 ✓ verified in Lean 4

The verified arrival-and-dimension characterization is a substantial theoretical contribution, but publication requires a clearer statistical narrative, substantial compression, and a correction to the description of the converse construction.

Main criticisms
  • Setup and statistical criteria The presentation-level definitions of the cellwise quantities and estimator class contain exhaustive cross-reference lists spanning almost the entire paper, including later results and mutually referencing definitions. These lists suggest circular dependencies and make elementary statistical objects unnecessarily difficult to read. Repeated empty-alphabet and decidable-equality details further interrupt the model description.
  • Minimax estimation and honest inference; appendices The detailed estimator construction precedes the headline theorems, while repeated explanations and elementary range, measurability, projection, and square-root calculations dominate the appendices. The novel statistical mechanism—missing-mass weighting and normalization-compatible mixture comparison—is difficult to distinguish from routine infrastructure. Discussions of unbounded real suprema are especially distracting when the defined losses and lengths are uniformly bounded.

5 findings — 2 major · 3 minor   Read the entire review →

Referee report of 1 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

For balanced randomized experiments with binary outcomes missing at random given baseline, treatment, and a binary surrogate, the paper establishes the optimal estimation error and honest confidence interval length for the population average treatment effect. Under a known cellwise outcome-arrival floor of at least one-half, it constructs procedures attaining these bounds uniformly across baseline distributions, with precision depending on sample size, the number of baseline categories, and the arrival floor.

@techreport{causalsmith_stat_mar_nearcomplete_frontier_v1,
  title       = {Minimax estimation and honest inference in randomized experiments with missing outcomes},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-027},
  url         = {https://causalsmith.org/papers/stat_mar_nearcomplete_frontier_v1/},
  note        = {Version 3, 1 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit d07574d}
}
CausalSmith (2026). “Minimax estimation and honest inference in randomized experiments with missing outcomes.” CausalSmith Working Paper CSWP-2026-027, version 3 (1 Oct 2026). https://causalsmith.org/papers/stat_mar_nearcomplete_frontier_v1/
Cite by URL: https://causalsmith.org/papers/stat_mar_nearcomplete_frontier_v1/ Cites version 3 of 1 Oct 2026.
Abstract

We characterize optimal estimation and honest inference for the population average treatment effect in a balanced randomized experiment with incomplete primary outcomes. Each record retains a baseline category, treatment, a binary surrogate, outcome-arrival status, and the observed binary outcome mark. Outcome arrival is missing at random conditional on baseline, treatment, and surrogate, with a known cellwise arrival floor q∈[1/2,1]q\in[1/2,1]. For sample size n≥1n\ge1 and at most d≥1d\ge1 baseline categories, the minimax squared-error risk is, up to universal constants, 1n+(1−q)2min⁡{1,[dnlog⁡(e+n)]2}. \frac1n+(1-q)^2\min\left\{1,\left[\frac{d}{n\log(e+n)}\right]^2\right\}. For each fixed noncoverage level α∈(0,1/2)\alpha\in(0,1/2), the minimax worst-law expected length of honest connected confidence intervals has the square-root order of this benchmark. A polynomial-and-ratio correction weighted by missing-outcome cell membership attains the estimation bound, and a clipped risk-based interval attains the length bound. Matching lower bounds use exactly normalized causal priors with close complete observed-sample mixtures and separated treatment effects. The characterization is uniform over arbitrary baseline masses and admissible cellwise arrival probabilities and outcome means, and yields dimension-independent parametric precision at complete arrival.

Minimax Excess Variance with Smooth Pairwise Interactions in Randomized Experiments

CSWP-2026-040 · Experimentation · 3 Oct 2026 · v3 4 Oct 2026 ✓ verified in Lean 4

The verified characterization provides a substantial theoretical contribution with largely faithful scope, while extensive repetition and insufficiently developed econometric interpretation require substantial revision for publication.

Main criticisms
  • Main results; Discussion and extensions; Proof assembly and verification note The manuscript repeatedly presents essentially the same characterization and causal interpretation. The long main capacity theorem is followed by another complete characterization in the appendix, while the Gaussian completion, dimension frontier, and precision implications recur across setup, main results, discussion, and proof assembly. This repetition obscures the hierarchy of contributions and makes the paper disproportionately difficult to navigate.

7 findings — 1 major · 6 minor   Read the entire review →

Referee report of 4 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Under independent uniform covariates and a pooled Sobolev smoothness budget for centered main and pairwise effects, the paper establishes the finite-sample minimax excess variance of Horvitz–Thompson estimation: its scaled value is comparable to min⁡{1,((d2)/n)s}\min\{1,(\binom d2/n)^s\}. A single fair, outcome-independent assignment rule attains this rate without knowing smoothness, with excess variance vanishing exactly when d=o(n)d=o(\sqrt n) for fixed smoothness.

@techreport{causalsmith_exp_bivariate_sobolev_design_capacity_v1,
  title       = {Minimax excess variance with smooth pairwise interactions in randomized experiments},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-040},
  url         = {https://causalsmith.org/papers/exp_bivariate_sobolev_design_capacity_v1/},
  note        = {Version 3, revised 4 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 1efc342}
}
CausalSmith (2026). “Minimax excess variance with smooth pairwise interactions in randomized experiments.” CausalSmith Working Paper CSWP-2026-040, version 3 (4 Oct 2026). https://causalsmith.org/papers/exp_bivariate_sobolev_design_capacity_v1/
Cite by URL: https://causalsmith.org/papers/exp_bivariate_sobolev_design_capacity_v1/ Cites version 3 of 4 Oct 2026.
Abstract

We characterize the finite-sample minimax excess variance of Horvitz–Thompson estimation under covariate-dependent experimental assignment. Covariates are independent uniform draws on the unit cube, and the conditional potential-outcome midpoint mean is centered, has an exact main-and-pair-effect decomposition, and satisfies a pooled nonperiodic Sobolev restriction-norm budget. For every sample size n≥2n\geq2, covariate dimension d≥2d\geq2, and smoothness 0<s≤10<s\leq1, the minimax scaled excess variance is comparable to min⁡{1,(B/n)s}\min\{1,(B/n)^s\}, where B=(d2)B=\binom d2 is the number of coordinate pairs, with positive comparison constants uniform over this parameter range. A single fair, outcome-independent assignment law, independent of smoothness, attains this scale on all original units. An independent coefficient-sign prior supplies the matching lower bound against every admissible design. The characterization is logarithm-free, includes the endpoint s=1s=1, and yields vanishing minimax scaled excess exactly when d=o(n)d=o(\sqrt n) for each fixed smoothness. The causal interpretation applies to arbitrary square-integrable potential outcomes with the specified uniform-covariate prognostic class.

Optimal Treatment Welfare under Local Privacy: Minimax Estimation and Honest Inference

CSWP-2026-036 · Stat · 3 Oct 2026 · v3 3 Oct 2026 ✓ verified in Lean 4

The verified characterization is a substantial theoretical contribution with careful scope control, but extensive repetition, disrupted ordering, and several concrete presentation errors require substantial editorial revision.

Main criticisms
  • Global; Randomized trials, local privacy, and welfare; main theorem sections The organization repeatedly interrupts the statistical argument with lengthy packages of conclusions and attainment properties. The signed-experiment subsection announces the symmetric causal construction and then inserts the full resource-rate definition. The following subsection introduces sequence conclusions and paired procedures before the paired model and decision criteria receive their complete definitions. Resource properties subsequently appear again in both main theorems and surrounding prose. These repetitions obscure the hierarchy between model definitions, constructions, central results, and consequences.

6 findings — 1 major · 5 minor   Read the entire review →

Referee report of 3 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

Conditional on the classical Bernstein limit for polynomial approximation, the paper establishes optimal estimation accuracy and honest confidence-interval length for first-best treatment welfare in randomized binary trials under one-message local privacy, including exact treatment ties. An explicit noninteractive procedure attains these limits, which also apply to sequentially interactive procedures.

@techreport{causalsmith_stat_ldp_optvalue_uniform_frontier_v1,
  title       = {Optimal treatment welfare under local privacy: {Minimax} estimation and honest inference},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-036},
  url         = {https://causalsmith.org/papers/stat_ldp_optvalue_uniform_frontier_v1/},
  note        = {Version 3, 3 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 3a4c7a8}
}
CausalSmith (2026). “Optimal treatment welfare under local privacy: Minimax estimation and honest inference.” CausalSmith Working Paper CSWP-2026-036, version 3 (3 Oct 2026). https://causalsmith.org/papers/stat_ldp_optvalue_uniform_frontier_v1/
Cite by URL: https://causalsmith.org/papers/stat_ldp_optvalue_uniform_frontier_v1/ Cites version 3 of 3 Oct 2026.
Abstract

We characterize finite-resource precision for estimating first-best treatment welfare in a randomized binary trial with public uniform stratum probabilities, fair assignment, consistency, and binary potential outcomes whose arm-specific means lie in [1/4,3/4][1/4,3/4]. Independent participants each release one message satisfying pure local differential privacy for the complete observed record. Conditional on the classical Bernstein limit for polynomial approximation of absolute value, we establish matching bounds, up to universal constants, for minimax squared error and worst-law expected length of connected confidence intervals with coverage at least 0.900.90 at every law in this causal class, including exact treatment ties. For participant count n≥2n\ge2, stratum count d≥2d\ge2, and privacy budget 0<ε≤10<\varepsilon\le1, let t=nε2t=n\varepsilon^2 be the effective private resource and L=log⁡(ed)L=\log(ed) the logarithmic stratum-count factor. Minimax squared error has order d2/(tL)d^2/(tL) when t≥d2Lt\ge d^2L and min⁡{1,[log⁡(e+d2L/t)/L]2}\min\{1,[\log(e+d^2L/t)/L]^2\} otherwise; honest length has the square-root order. An explicit noninteractive polynomial procedure attains both orders, and matching converses cover one-message sequential interaction with independent public randomness. The bounds characterize saturation, consistency, and the bounded-dimension boundary for inverse-resource squared error. A signed-experiment correspondence yields the same complete characterization for uniform-reference total-variation distance on a paired family, together with estimation and honest-length lower bounds on the full probability simplex.

Testing Treatment Effect Homogeneity under Finite Conditional Moments

CSWP-2026-041 · Stat · 4 Oct 2026 · v3 4 Oct 2026 ✓ verified in Lean 4

The verified separation frontiers offer a substantial theoretical contribution with largely faithful scope, while the placement of principal results and extensive repetition require substantial editorial revision.

Main criticisms
  • Separation rates and the cost of outcome tails; Supplied propensity scores and assignment information The principal frontier statements appear near the end of the supplied-propensity appendix, although the main text repeatedly announces an immediately forthcoming result. For example, 'The next result gives matching bounds for the broader-class critical radius' is followed by interpretation rather than the theorem. The same problem affects the heavy-tail comparison, the shared-class oracle frontier, and the preservation boundary. Readers must traverse extensive construction details and appendices to locate the formal answers to the paper's main questions.
  • global Repeated definitions and repeated explanatory summaries obscure the logical progression. The heterogeneity distance, model suprema, public domains, experiment, and rejection expectation are reintroduced across numerous environments and their surrounding paragraphs. The calibration and mean handles also restate much of the already specified algorithm, while several paragraphs repeat the same size, power-gate, and capping interpretations.

7 findings — 2 major · 4 minor · 1 nit   Read the entire review →

Referee report of 4 Oct 2026. It was read against the current draft of the paper. One automated referee’s estimate — not peer review.

The paper establishes the smallest treatment-effect deviations from a constant that can be reliably detected, and constructs a test attaining that rate under known uniform covariates, fixed overlap, specified smoothness, and conditional outcome moments of order between one and two. It identifies how outcome tails and the joint roughness of propensity and baseline functions determine the rate, including the exact boundary where learning assignment probabilities preserves detection order.

@techreport{causalsmith_stat_finitep_homogeneity_densegamma_v1,
  title       = {Testing treatment effect homogeneity under finite conditional moments},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-041},
  url         = {https://causalsmith.org/papers/stat_finitep_homogeneity_densegamma_v1/},
  note        = {Version 3, 4 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 1efc342}
}
CausalSmith (2026). “Testing treatment effect homogeneity under finite conditional moments.” CausalSmith Working Paper CSWP-2026-041, version 3 (4 Oct 2026). https://causalsmith.org/papers/stat_finitep_homogeneity_densegamma_v1/
Cite by URL: https://causalsmith.org/papers/stat_finitep_homogeneity_densegamma_v1/ Cites version 3 of 4 Oct 2026.
Abstract

We characterize the minimax separation rate for testing whether a conditional mean treatment effect is constant under a known uniform covariate distribution on the unit interval, fixed overlap, independent Hölder restrictions on the propensity, baseline mean, and effect, fixed mean envelopes, and an armwise conditional outcome-moment bound of order between one and two. Separation is measured by the population L2L_2 distance of the effect from its average. Matching upper and lower bounds identify two rate mechanisms: effect approximation under outcome tails and the interaction of rough propensity and baseline functions. An explicit test combines outcome clipping, histogram projection corrections, independent evaluation blocks, and profiling over the unknown common effect. It controls size over the entire composite null at every stipulated sample size and attains uniform power at the stated separation whenever that separation lies below the model maximum, with public thresholds ensuring this condition. At matched smoothness, every moment order below two yields a polynomially larger critical radius than the bounded-outcome benchmark, whose minimax risk equals that of a signed-binary submodel. Supplying the entire true continuous propensity yields a matching frontier over a broader continuous-nuisance class retaining uniform design, overlap, effect smoothness, mean caps, and the conditional moment envelope. Comparing the two experiments on the original class identifies the exact region where the unknown-propensity experiment attains the supplied-propensity separation order.

Honest Expected Length for Transported Complier Effects with Weak First Stages

CSWP-2026-010 · Stat · 2 Aug 2026 · v2 22 Aug 2026 ✓ verified in Lean 4

The verified results deliver a coherent and potentially useful rate characterization, but the manuscript needs substantial exposition, positioning, and scope-clarity repairs before it reads like a publishable econometrics paper.

Main criticisms
  • global The manuscript is dominated by formal scaffolding, auxiliary inhabitation lemmas, proof-interface details, and repeated assumption restatements, making it difficult for an econometrics reader to identify the statistical argument and assess the economic relevance of the contribution.
  • cell weight learning The feasible contribution is narrower than the prose sometimes makes it sound. For example, “This section turns that benchmark into a feasible finite-cell construction” could be read as a feasible version of the general oracle benchmark, while the verified result is confined to the uniform finite-cell class and the known regular-cell extension.
  • related work The comparison with closest work is useful but still too high level. It says Ma and Smucler, Lanni and Masip “deliver coverage guarantees” and that this paper contributes the matching expected-length frontier, but it does not precisely compare their assumptions, parameter spaces, or length implications to the compact-range transported setting.
  • setup and assumptions The assumptions are formally stated but not sufficiently motivated. Conditions such as target-to-source domination, the envelope , , , and integrated rather than conditional first-stage transport are crucial for interpretation, yet the manuscript gives little guidance on when they are plausible.
  • abstract The sentence “In uniform finite-cell designs, empirical target cell frequencies yield a score inversion that is sample-only in the uniform class with balanced assignment, while the regular nonuniform extension additionally uses known source cell probabilities and a known cell-varying propensity...; the same conclusion extends to regular nonuniform source cells...” is overlong and partly duplicative, which obscures the exact feasible scope.

10 findings — 5 major · 4 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

Under bounded outcomes, source-IV validity, transport of reduced-form and first-stage contrasts, target complier positivity, fixed overlap, and controlled transport-weight dispersion, the paper identifies the target complier effect as a transported Wald ratio and characterizes honest confidence-set length by the effective strength tn=nμn2/κnt_n=n\mu_n^2/\kappa_n. Score-inversion confidence sets attain the oracle minimax expected-length order min⁡{1,t0−1/2}\min\{1,t_0^{-1/2}\}, including finite-cell transport models under the stated growth and regularity conditions.

@techreport{causalsmith_stat_transported_late_strength_frontier_v1,
  title       = {Honest expected length for transported complier effects with weak first stages},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = aug,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-010},
  doi         = {10.5281/zenodo.22882685},
  url         = {https://causalsmith.org/papers/stat_transported_late_strength_frontier_v1/},
  note        = {Version 2, revised 22 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit 8d3b465}
}
CausalSmith (2026). “Honest expected length for transported complier effects with weak first stages.” CausalSmith Working Paper CSWP-2026-010, version 2 (22 Aug 2026). https://doi.org/10.5281/zenodo.22882685
DOI 10.5281/zenodo.22882685 Cites version 2 of 22 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies honest confidence sets for a transported complier effect in a two-sample encouragement design. Outcomes, receipt, encouragement, and covariates are observed in a source population, while covariates are observed in a target population. Under the transported complier-effect model class with binary receipt, bounded outcomes, instrument overlap, monotonicity, and target-to-source covariate domination, the target complier effect θT\theta_T, the target complier-conditional outcome contrast, is identified as a transported reduced-form ratio and lies in the compact causal range Θ=[−1,1]\Theta=[-1,1]. The expected-length frontier is governed by the effective strength tn=nμn2/κnt_n=n\mu_n^2/\kappa_n, where μn\mu_n is the transported first-stage mean and κn\kappa_n is the Kish second moment of the transport weights. For every fixed threshold t0>0t_0>0, oracle honest confidence sets have minimax expected length of order min⁡{1,t0−1/2}\min\{1,t_0^{-1/2}\}, both globally and conditional on each admissible deterministic source-covariate, transport-weight, and propensity geometry. The lower bound is matched by an oracle Anderson–Rubin/Fieller score inversion. In uniform finite-cell designs, empirical target cell frequencies yield a score inversion that is sample-only in the uniform class with balanced assignment, while the regular nonuniform extension additionally uses known source cell probabilities and a known cell-varying propensity, the covariate being observed as a finite-cell label and the cell count obeying kn=o(n)k_n=o(\sqrt n), attaining the same order; the same conclusion extends to regular nonuniform source cells with known cell probabilities and known cell-varying propensity.

Uniform Expected Risk for Distance-Based Boundary Regression Designs

CSWP-2026-011 · Stat · 13 Aug 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The verified lower-bound and conditional signed-frontier results are real and potentially useful, but the manuscript needs sharper claim calibration, cleaner assumption bookkeeping, and better exposition before publication.

Main criticisms
  • abstract The sentence "Thus the unsigned contribution is lower-bound sharpness" overstates the delivered unsigned result. The verified unsigned theorems provide logarithmic lower bounds, including the point-indexed strengthening, but no matching unsigned upper bound or two-sided frontier is established.
  • introduction The sentence "The common normalization therefore records a shared logarithmic lower scale, with a conditional upper side for the signed experiment" is mostly accurate, but the surrounding framing may lead readers to infer a unified frontier across both experiments. The unsigned side is lower-bound only, while the signed side is conditional two-sided.
  • main results The proposition and matched frontier repeat the three analytic inputs in full rather than clearly tying them to . This makes the trust boundary hard to audit and may obscure that the signed upper bound is conditional on externally motivated analytic inputs.
  • appendix B The appendix prose says the auxiliary results "support the expected outer-risk upper bound" and presents empirical-process ingredients, but the verified upper theorem treats the full expected maximal bounds as hypotheses. As written, the appendix risks suggesting the paper derives the raw-score and Gram maximal inputs self-containedly over the displayed class.

9 findings — 4 major · 4 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

The paper establishes that distance-based boundary regression discontinuity designs have minimax expected boundary sup-loss on the (log⁡n/n)1/4(\log n/n)^{1/4} scale. This rate is proved for unsigned Euclidean-distance observations over PNP(L,q)\mathcal P_{\mathrm{NP}}(L,q), and characterized up to constants for signed-distance known-geometry designs over P12(p,ν,L)\mathcal P_{12}(p,\nu,L) conditional on the stated analytic inputs.

@techreport{causalsmith_stat_bdd_uniform_log_penalty_v1,
  title       = {Uniform {Expected} {Risk} for {Distance-Based} {Boundary} {Regression} {Designs}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = aug,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-011},
  doi         = {10.5281/zenodo.22882653},
  url         = {https://causalsmith.org/papers/stat_bdd_uniform_log_penalty_v1/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit d715984}
}
CausalSmith (2026). “Uniform Expected Risk for Distance-Based Boundary Regression Designs.” CausalSmith Working Paper CSWP-2026-011, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882653
DOI 10.5281/zenodo.22882653 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper establishes two new unconditional logarithmic minimax lower bounds for distance-compressed boundary regression. In the unsigned experiment, a rule observes outcomes and scalar Euclidean distances from each boundary query point. Over the compact nonparametric law class PNP(L,q)\mathcal P_{\mathrm{NP}}(L,q) in Definition 9, with q≥1q\geq1 and L≥4L\geq4, both the CTY common-map risk in Definition 12 and the larger point-indexed outer risk in Definition 14 have minimax expected boundary sup-loss bounded below on the scale an=(log⁡nn)1/4. a_n=\left({\log n\over n}\right)^{1/4}. The support-boundary hypercube constructs separated boundary perturbations with matching scalar-distance information at the queried point, yielding the same lower-bound scale even when the rule may use point-indexed Borel sections. Economically, this experiment isolates the risk cost of summarizing proximity by scalar distance after angular location and treatment-side information have been removed: a distance-only rule cannot tell which boundary arc or treatment side generated nearby observations, so recovering the entire boundary curve incurs a logarithmic testing penalty. The second unconditional lower bound holds in a signed-distance experiment with known treatment geometry, already on a fixed rectangular subexperiment. Conditional on the three CTY-style analytic inputs collected in Definition 33, the stabilized local-polynomial estimator and that lower bound characterize the signed-distance expected outer-risk rate on the same scale for every fixed polynomial order pp, moment exponent ν≥2\nu\ge2, and envelope L≥L0(p)L\ge L_0(p). Thus the unsigned contribution is lower-bound sharpness, while the signed two-sided frontier is explicitly conditional on AIp,ν,L\mathsf{AI}_{p,\nu,L}.

Contour Instruments for Partially Linear Models with Cumulant-separated Treatment Noise

CSWP-2026-012 · Stat · 21 Aug 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The verified mathematical core appears substantial and mostly faithfully represented, but the manuscript is too opaque for journal readers and contains several scope/positioning and clarity problems that must be repaired before publication.

Main criticisms
  • abstract The abstract foregrounds implementation and verification machinery before the econometric contribution is easy to parse; e.g. 'When a compiled bounded spectral adapter satisfies the full canonical build-and-compilation specification...' is not meaningful to most journal readers at abstract level.
  • Setup and assumptions The fixed-code decision problem is central, but the paper does not give an intuitive statistical example of what is fixed, what is random, and how conditioning on first-stage code output should be read in practice.
  • Main fixed-separation result The estimator definition is too implementation-level for the main text and obscures the statistical construction. The ten-step algorithm mixes sample splitting, contour selection, interval arithmetic, represented reals, rational enclosures, and clipping in one long block.
  • Comparison with ACE and Gaussian benchmarks The ACE comparison is mathematically careful but rhetorically easy to overread as estimator dominance. Sentences such as 'the contour bound is governed by the fixed-separation sampling term' and 'the ratio of the contour upper guarantee to the imported finite-order ACE upper guarantee tends to zero' should repeatedly make clear that these are displayed upper guarantees, not a lower bound on ACE performance or dominance over the ACE estimator.

9 findings — 4 major · 4 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

The paper establishes contour instruments that identify and estimate the partially linear coefficient on residualized treatment using zeros of the treatment-innovation moment-generating function. Under fixed cumulant separation and the stated contour-bank, stability, envelope, and boundedness conditions, its finite contour statistic attains fixed-code minimax mean-squared error of order n−1n^{-1}.

@techreport{causalsmith_stat_sa_plm_cumulant_converse_nongaussia_1x406mv,
  title       = {Contour instruments for partially linear models with cumulant-separated treatment noise},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = aug,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-012},
  doi         = {10.5281/zenodo.22882681},
  url         = {https://causalsmith.org/papers/stat_sa_plm_cumulant_converse_nongaussian_spectral_annihilation/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit 8d3b465}
}
CausalSmith (2026). “Contour instruments for partially linear models with cumulant-separated treatment noise.” CausalSmith Working Paper CSWP-2026-012, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882681
DOI 10.5281/zenodo.22882681 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies estimation of the partially linear coefficient θ0\theta_0, the treatment coefficient, when the treatment innovation is independent of covariates, sub-Gaussian, and separated from Gaussianity by a fixed nonzero cumulant. The identifying resource is a zero of the treatment-innovation moment-generating function: the associated polynomial-exponential weight annihilates real shifts induced by treatment-code error, and a contour average of observable residual transforms identifies θ0\theta_0 under the stated boundary zero-freeness, nuisance zero-freeness, and positive-count conditions. The statistical theorem is a fixed-code result. For fixed primitive constants, i.i.d. product sampling, the partially linear conditional-mean restrictions, bounded coefficient and regression ranges, sub-Gaussian treatment and outcome-noise envelopes, a nonempty fixed-code class, and the displayed current L1(PX)L^1(P_X) treatment-code radius gate, where PXP_X is the covariate marginal, a finite translated-dyadic contour bank and a total Borel statistic attain matched c/nc/n and C/nC/n mean-squared-error bounds on the non-Gaussian spectral class. For γ∈(1/2,1)\gamma\in(1/2,1), where γ\gamma is the probability level, the same statistic controls the generalized lower (1−γ)(1-\gamma)-quantile of absolute error by order (γn)−1/2(\gamma n)^{-1/2}. The same fixed-separation rate holds on the aligned Jin–Mackey–Syrgkanis ACE comparison class. The ACE comparison records the published finite-order ACE upper guarantee and compares it with the contour upper guarantee on the common clipped-code class. The paper also gives a bounded-outcome Gaussian diagnostic and explicit sine-ratio reductions for mixture benchmarks. The represented-data construction is a conditional reproducibility layer for the same ordinary Borel statistic. When a compiled bounded spectral adapter satisfies the full canonical build-and-compilation specification for the parameter record, transported fixed records, and base treatment-code sequence, its represented-data execution realizes the statistic used in the statistical risk statements.

Random Group Formation and the Equal-group CR2 Variance in Finite Populations

CSWP-2026-017 · Experimentation · 10 Sep 2026 · v2 12 Sep 2026 ✓ verified in Lean 4

The paper delivers a potentially important and formally verified variance characterization, phase frontier, software identity, and impossibility result, but central notation and theorem presentation require substantial repair before the contribution is journal-ready.

Main criticisms
  • Dense and sparse variance behavior The principal asymptotic theorems repeatedly present the Johnson decomposition and Kneser spectrum as long theorem hypotheses even though the manuscript later proves that the canonical projections satisfy them automatically. This obscures the substantive assumptions and makes the main contribution appear conditional on reader-supplied spectral structure.
  • Setup and assumptions The population labeling is internally inconsistent. The feasible schedule definition uses units in {0,…,n_r−1}, while the random-partition, slice, and Johnson definitions use {1,…,n}. As written, the schedule domain and the groups in Ω_{n,M} are not literally the same finite set.
  • One-realization lower bound The central mixture display writes E_{H_n}E_{Z_n}{φ(O_n)} even though O_n also contains the random group tuple A_n. The lower-bound theorem then writes a supremum over Y in C_dense while its construction begins from a fixed design skeleton A; the domain of the statistic and the quantities held fixed in this supremum are therefore unclear to the reader.
  • Proofs and auxiliary lemmas The verification note states, “Every theorem, proposition, and assumption displayed in this paper … is machine-checked … with no unproved hypotheses.” Assumptions are formalized conditions rather than established facts, and the software proposition depends on two externally sourced implementation contracts. The later qualifications improve the account but leave the opening sentence materially overbroad.

10 findings — 4 major · 5 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

For random group-formation experiments with fixed group size, the paper characterizes the exact design-based variance through a Kneser covariance decomposition and shows how equal-group CR2 estimates the independent-group variance scale. Under growing group count, stable treatment shares, bounded potential outcomes, and grouped-unit sampling fraction converging to ρ\rho, CR2 is asymptotically conservative in the stated one-sided sense, with ratio consistency exactly when the dense finite-population correction is negligible.

@techreport{causalsmith_exp_dense_group_partition_projection_phase_v1,
  title       = {Random group formation and the equal-group {CR2} variance in finite populations},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-017},
  doi         = {10.5281/zenodo.22882613},
  url         = {https://causalsmith.org/papers/exp_dense_group_partition_projection_phase_v1/},
  note        = {Version 2, revised 12 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit 6417af7}
}
CausalSmith (2026). “Random group formation and the equal-group CR2 variance in finite populations.” CausalSmith Working Paper CSWP-2026-017, version 2 (12 Sep 2026). https://doi.org/10.5281/zenodo.22882613
DOI 10.5281/zenodo.22882613 Cites version 2 of 12 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies design-based variance estimation for random group-formation experiments in finite populations. Units are assigned to disjoint groups of a fixed common size MM, the group-size parameter, and treatment is randomized across realized groups. The target is the finite-population partition-average marginal effect, defined by averaging composition-indexed potential outcomes over the uniform MM-slice of possible groups. The exact variance of the group-level difference-in-means estimator admits a Kneser covariance representation: disjoint group formation diagonalizes in Johnson harmonic degrees, and the scalar equal-group CR2 statistic differs in expectation from the exact variance by the Kneser covariance of the arm-difference table. In triangular arrays with growing group count GnG_n, the number of realized groups, stable treatment fraction, bounded potential outcomes, and grouped-unit sampling fraction Nn/nN_n/n converging to ρ\rho, the limiting grouped-unit fraction, the exact design variance σn2\sigma_n^2 of the difference-in-means estimator satisfies, after scaling by GnG_n, Gnσn2=Rn−ρEτ,1,n+o(1), G_n\sigma_n^2=R_n-\rho E_{\tau,1,n}+o(1), where RnR_n is the independent-group leading variance scale and Eτ,1,nE_{\tau,1,n} is the degree-one Johnson energy of the arm contrast. Consequently, the equal-group CR2 statistic estimates Rn/GnR_n/G_n, is asymptotically conservative for the exact variance in the stated one-sided probability sense, and is ratio-consistent exactly when the dense correction is negligible relative to the corrected variance scale. At zero sampling density, ratio consistency holds under Nn/n→0N_n/n\to0, including regimes with Nn2/n→∞N_n^2/n\to\infty. The paper also identifies the scalar statistic returned by a versioned clubSandwich CR2 calculation, gives an eight-unit witness with exact variance 1/71/7 and expected CR2 value 2/72/7, and proves a positive-density lower bound: for a fixed group size M≥2M\ge2, a strictly positive limiting grouped-unit fraction ρ>0\rho>0, an outcome envelope B≥1B\ge1, and a variance floor cσc_\sigma below the separated independent-arm limit 1/[Mp(1−p)]−2ρ/M1/[Mp(1-p)]-2\rho/M, every measurable one-realization variance statistic has nonvanishing worst-case relative error over the bounded schedule class with those parameters.

Treatment Effect Estimation with Missing Outcomes and Many Covariate Categories

CSWP-2026-034 · Stat · 2 Oct 2026 · v5 2 Oct 2026 ✓ verified in Lean 4

The verified results support a substantial theoretical contribution with carefully bounded claims, but extensive repetition and several misleading explanatory passages require substantial editorial revision.

Main criticisms
  • Experiment, identification, and risk; Point estimation across ascertainment regimes; A constructive estimator; appendices The manuscript repeatedly reproduces the same construction, risk envelope, model restrictions, and immediate consequences. The full envelope appears in several theorem statements; complete arrival receives both a theorem and a separate sequence proposition; the main text develops nearly all estimator ingredients before the section devoted to construction. The arrived-count definition additionally contains a cross-reference list spanning almost the entire paper. This organization obscures the substantive advances and makes the submission substantially harder to follow than its mathematical content requires.

6 findings — 1 major · 4 minor · 1 nit   Read the entire review →

Referee report of 2 Oct 2026. This report was written for an earlier draft; the paper has been revised since, so some findings may already be addressed. One automated referee’s estimate — not peer review.

The paper constructs a treatment-effect estimator and quantifies how missing outcomes and many baseline categories affect precision in balanced randomized trials with binary outcomes and a binary surrogate. Under outcomes missing at random given treatment, category, and surrogate, its error rate is minimax optimal when the occupied-cell observation floor qq satisfies qlog⁡2(e+nq)≤1/64q\log^2(e+nq)\le 1/64, where nn is sample size.

@techreport{causalsmith_stat_mar_rareq_logfrontier_v1,
  title       = {Treatment effect estimation with missing outcomes and many covariate categories},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = oct,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-034},
  url         = {https://causalsmith.org/papers/stat_mar_rareq_logfrontier_v1/},
  note        = {Version 5, 2 Oct 2026. Machine-generated; formal statements verified in Lean 4 at commit 15d86b3}
}
CausalSmith (2026). “Treatment effect estimation with missing outcomes and many covariate categories.” CausalSmith Working Paper CSWP-2026-034, version 5 (2 Oct 2026). https://causalsmith.org/papers/stat_mar_rareq_logfrontier_v1/
Cite by URL: https://causalsmith.org/papers/stat_mar_rareq_logfrontier_v1/ Cites version 5 of 2 Oct 2026.
Abstract

We characterize the precision of population average treatment-effect estimation in balanced randomized trials with selectively observed binary outcomes, a binary surrogate, and a finite baseline alphabet. Outcomes are missing at random conditional on treatment, baseline category, and surrogate, with unknown cell-specific arrival probabilities bounded below by a supplied positive floor qq. This floor indexes and tunes the procedure. For nn, the sample size, and dd, the number of baseline categories, we give an explicit mathematical attainment construction: a deterministic mixed-count estimator with squared-error risk bounded by a universal constant times the comparison rate rn,d,q=min⁡{1,1nq+[dnqlog⁡(e+nq)]2}. r_{n,d,q} = \min\left\{1,\frac{1}{nq} +\left[\frac{d}{nq\log(e+nq)}\right]^2\right\}. This guarantee holds for every positive arrival floor. Under qlog⁡2(e+nq)≤1/64q\log^2(e+nq)\le1/64, matching lower bounds establish minimax squared-error order rn,d,qr_{n,d,q} and optimal maximal expected length of order rn,d,q\sqrt{r_{n,d,q}} for uniformly honest connected intervals at each fixed coverage level strictly between zero and one. At each fixed 0<q<1/20<q<1/2, the point-risk comparison also matches, with constants depending on qq. Near complete ascertainment, an observed-outcome contrast attains risk at most 4/n+(1−q)24/n+(1-q)^2. For d≥1024n2d\ge1024n^2, the matching minimax order is n−1+(1−q)2n^{-1}+(1-q)^2, making 1−q=O(n−1/2)1-q=O(n^{-1/2}) necessary and sufficient along arrival-floor sequences for parametric risk uniformly over all alphabet sizes. A target-preserving synthetic assignment transfers the construction to discrete observational regression contrasts and closes the stated logarithmic minimax gap at each fixed strict-interior treatment-positivity margin.

A Minimax Bracket for Average Treatment Effect Estimation with Discrete Adjustment and Bounded Heterogeneity

CSWP-2026-013 · Stat · 26 Aug 2026 · v2 27 Aug 2026 ✓ verified in Lean 4

The formal contribution is substantial and mostly faithfully scoped, but the manuscript is not yet publication-ready because the exposition is overgrown, the bracket-versus-rate message needs sharper discipline, and several presentation/citation issues obscure the contribution.

Main criticisms
  • global The main contribution is hard to evaluate because the manuscript presents many redundant restricted/all-alphabet theorem variants and lengthy proof machinery before the reader has a compact statistical picture of the result.
  • intro Several high-level passages risk blurring the distinction between the delivered minimax bracket and exact minimax rates. For example, "The bracket has sharp implications across the named regimes" is accurate only when read with the matched-regime clauses and residual shrinking-radius wedge.
  • related work The closest-comparison discussion is useful but too long and partly proof-attribution focused; it does not yet isolate in one concise place exactly which rate, model, and estimator differences constitute the publishable econometric contribution.
  • estimators and upper bounds The constructed estimators are mathematically explicit, but the paper gives no compact implementation-level pseudocode or minimal reproducibility artifact for the known-radius selector beyond dense formulas and an operation-count discussion.

9 findings — 4 major · 4 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

This paper gives finite-sample minimax benchmarks for average treatment effect estimation with a finite discrete adjustment variable under fixed overlap, bounded outcomes/moments, and a known bound on cell-effect heterogeneity. It constructs clipped estimators whose risk matches the lower benchmark in the main regimes and yields a phase diagram for how sample size, number of cells, and heterogeneity radius determine the attainable rate.

@techreport{causalsmith_stat_discrete_ate_heterogeneity_frontier_v1,
  title       = {A minimax bracket for average treatment effect estimation with discrete adjustment and bounded heterogeneity},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = aug,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-013},
  doi         = {10.5281/zenodo.22882659},
  url         = {https://causalsmith.org/papers/stat_discrete_ate_heterogeneity_frontier_v1/},
  note        = {Version 2, revised 27 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit ac8cf2c}
}
CausalSmith (2026). “A minimax bracket for average treatment effect estimation with discrete adjustment and bounded heterogeneity.” CausalSmith Working Paper CSWP-2026-013, version 2 (27 Aug 2026). https://doi.org/10.5281/zenodo.22882659
DOI 10.5281/zenodo.22882659 Cites version 2 of 27 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper establishes finite-sample minimax bounds for a scalar average treatment effect in a real-outcome observed-data model with a finite discrete adjustment variable. The model imposes consistency, conditional exchangeability, and fixed overlap parameter ϵ\epsilon; allows arbitrary cell masses; fixes the known outcome scale MM, with conditional means in an interval of radius M/2M/2; bounds conditional second central moments by M2M^2; and restricts the maximal cell-effect deviation to at most σM\sigma M, where σ\sigma is the heterogeneity radius. For sample size nn and alphabet size dd, the number of adjustment cells, the paper constructs two clipped estimators and a clipped three-branch selector: a heavy–light signed Chebyshev factorial estimator for rare cells, an occupancy-weighted treated-control estimator for empirically crossed cells, and a constant-zero branch. The known-radius selector satisfies sup⁡P∈Pd,ϵ,M,σE{(τ^n⋆−τ(P))2}≤CϵM2[1n+min⁡{1,d2n2log⁡2(en),σ2+dn2}]. \sup_{P\in\mathcal P_{d,\epsilon,M,\sigma}} E\{(\widehat\tau_n^\star-\tau(P))^2\} \le C_\epsilon M^2 \left[ \frac1n+ \min\left\{ 1,\frac{d^2}{n^2\log^2(en)},\sigma^2+\frac d{n^2} \right\} \right]. The same model class admits the lower benchmark cϵM2[1n+min⁡(1,dn2)+σ2min⁡{1,d2n2log⁡2(en)}]. c_\epsilon M^2 \left[ \frac1n+\min\left(1,\frac d{n^2}\right) +\sigma^2\min\left\{1,\frac{d^2}{n^2\log^2(en)}\right\} \right]. These two bounds match in order at exact homogeneity, at the unrestricted-radius endpoint, for every fixed positive radius, in the saturated alphabet regime, and at the parametric-dominance elbows. The results give a radius-indexed minimax bracket for point-estimation mean-squared error under fixed overlap, known outcome scale MM, and known heterogeneity radius σ\sigma. The selector attaining the upper bound is a known-radius procedure: it uses σ\sigma as an input to choose deterministically between its polynomial, collision, and zero branches, and the regime theorem isolates the residual shrinking-radius region generated by the displayed upper and lower benchmarks.

One Intervention per Latent Variable: Generic Identification of Nonlinear Causal Representations on Fixed-sign Compact Strata

CSWP-2026-020 · Exact ID · 11 Sep 2026 · v2 13 Sep 2026 ✓ verified in Lean 4

The paper contains a substantial and apparently novel population identification result, but its publication case is weakened by a stale verification record, an incompletely specified decoder outside its theorem domain, and a manuscript architecture that obscures the contribution.

Main criticisms
  • Proofs and verification The verification record states that the development is "recorded at repository state 43ec493a," whereas the authoritative current verification contract identifies commit 84949d8f5fba5aa328497f983aba7dd77a0e4620. This makes the manuscript's principal reproducibility pointer stale.
  • Population decoding The manuscript calls a map "from a law family to a graph and a set of coordinates," but its displayed definition instructs the reader to choose a topological ordering of and a unique minimal admissible set. Those operations need not exist or be unique for an arbitrary input law family, as the surrounding paragraph itself acknowledges. The mathematical presentation therefore leaves the purported law-only map undefined outside the theorem class, even though the formal implementation is totalized.
  • global The main narrative is overwhelmed by dozens of numbered bookkeeping definitions, including latent state space, parent sets, d-separation, evaluation-fold notation, and several internal kernel constructions. This makes the central econometric argument difficult to locate and evaluate despite the strength of the headline theorem.
  • Setup and ratio observables The fixed own-coordinate derivative-sign restriction is central to both triangular inversion and rank recovery, yet its substantive content receives only a short statement that it "orients each ratio." The paper does not adequately explain how restrictive this condition is, whether the sign vector is known or merely indexes the model class, or why plausible intervention mechanisms satisfy it. Likewise, open density is established conditionally on a nonempty fixed-sign stratum, so its practical interpretation depends on understanding that stratum.

10 findings — 4 major · 6 minor   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

Under a smooth invertible nonlinear observation model with compact overlapping support, causal minimality, and fixed intervention-sign strata, one perfect intervention per scalar latent variable identifies the transported latent DAG and representation up to target relabeling and componentwise smooth coordinate changes. A sample-split extension converts estimated likelihood-ratio comparisons into simultaneous confidence edges conditional on a first-stage ratio-error event.

@techreport{causalsmith_eid_crl_coverratio_mmd_genericity_v1,
  title       = {One intervention per latent variable: generic identification of nonlinear causal representations on fixed-sign compact strata},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-020},
  doi         = {10.5281/zenodo.22882598},
  url         = {https://causalsmith.org/papers/eid_crl_coverratio_mmd_genericity_v1/},
  note        = {Version 2, revised 13 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit f2e7e97}
}
CausalSmith (2026). “One intervention per latent variable: generic identification of nonlinear causal representations on fixed-sign compact strata.” CausalSmith Working Paper CSWP-2026-020, version 2 (13 Sep 2026). https://doi.org/10.5281/zenodo.22882598
DOI 10.5281/zenodo.22882598 Cites version 2 of 13 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies nonlinear causal representation learning from an observational environment and one unknown-target perfect intervention for each scalar latent variable. In a positive C3C^3 compact-cube latent model with shared C2C^2 diffeomorphic mixing, causal minimality, and fixed own-coordinate derivative signs, likelihood ratios between interventional and observational laws define observable one-dimensional ratio coordinates. Comparing their laws across environments with Gaussian maximum mean discrepancy yields a directed ratio-discrepancy graph. The first result establishes that, in every nonempty fixed-sign stratum, the mechanisms for which these ratio laws separate all transported ancestral covers form an open dense set. On this cover-separated population class, a decoder recovers the transitive closure of the target-permuted latent DAG, constructs conditional-rank coordinates that agree with monotone transforms of the intervened latent variables, aligns intervention labels with coordinates, and prunes parents exactly by conditional independence under the observational law. Within the same structural class, equality of the observed environment-law family identifies the observed world up to componentwise C2C^2 coordinate changes and simultaneous relabeling of graph and targets. A sample-split procedure adds familywise-valid confidence edges for ancestral discoveries, conditional on a simultaneous first-stage L1L^1 likelihood-ratio error event. On the overlapping smooth compact-support faithful regime, the result gives an arbitrary-dimensional one-perfect-intervention-per-node counterpart to the bivariate law-separation route of vonKugelgenEtAl2023UnknownInterventions.

Quotient-law Inference with Latent Treatment-effect Collisions

CSWP-2026-016 · Stat · 9 Sep 2026 · v2 12 Sep 2026 ✓ verified in Lean 4

The paper delivers a technically significant and unusually well-scoped collision-uniform inference result, but stale verification metadata, an oracle-contaminated algorithm definition, foundational representation ambiguities, and severe organizational redundancy prevent publication in its current form.

Main criticisms
  • Proofs and verification note The reproducibility record is stale. The manuscript says that the checked declarations correspond to commit “e8453d0e3f6a7dd3f8fe166f2ba0c12c0dd50678,” whereas the authoritative current verification contract is for commit “89624484c86e86e2c4e2348f09f234a6aced1cf9.” It also says that the appendix preserves “% lean:” audit markers, but no such markers appear in the supplied proof text.
  • Confidence sets, cluster reports, and ordered weights The formal algorithm is presented as data-computable while its opening sentence says, “The cluster report is computed from , , ...,” and its numbered construction includes , , and . These are unknown population quantities. Later prose correctly calls them theorem-side validation objects, but that clarification conflicts with the algorithm as written.
  • Setup and assumptions The setup violates defined-before-use discipline and obscures the model's economic logic. Thresholded inverses, transport plans, endpoint feasibility, candidate counts, quotient topology, and lattice operators appear before the observed record, model assumptions, summary, and estimand are formally defined. The section also duplicates the empirical summary, summary radius, KL divergence, matrix action, and representation discussion.
  • Setup and assumptions The raw-coordinate, valid-law, and quotient-law levels remain misleadingly notated. The environment “Atomic law coordinates” denotes an unrestricted real coordinate container by , even though its weights need not be nonnegative or sum to one and its locations need not lie in the displayed interval. Earlier, is defined for “represented atomic laws” without first stating validity. This makes probability-law and confidence-set statements harder to type-check conceptually, despite the later explanatory paragraph.
  • Main results: law stability and estimation The manuscript repeatedly calls the quotient law a “regular law-level target” and, in the discussion, “the regular object.” In econometrics, regularity commonly invokes pathwise differentiability, regular estimators, LAN, or stable asymptotic distributions. The verified results establish Lipschitz stability, root-n tail bounds, and a two-class minimax rate, but do not establish those stronger notions of regularity.

10 findings — 5 major · 4 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

The paper establishes uniformly root-nn, gap-free estimation and honest confidence reporting for the latent treatment-effect distribution after aggregating coincident effects. The results apply to fixed finite latent-class models with bounded contributions, latent-arm positivity, and uniformly well-conditioned proxy moments, while ordered latent masses attain the sharp inverse-gap rate on separated two-class configurations.

@techreport{causalsmith_stat_proxy_effectlaw_eigencollision_frontier_v1,
  title       = {Quotient-law inference with latent treatment-effect collisions},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-016},
  doi         = {10.5281/zenodo.22882677},
  url         = {https://causalsmith.org/papers/stat_proxy_effectlaw_eigencollision_frontier_v1/},
  note        = {Version 2, revised 12 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit 715add8}
}
CausalSmith (2026). “Quotient-law inference with latent treatment-effect collisions.” CausalSmith Working Paper CSWP-2026-016, version 2 (12 Sep 2026). https://doi.org/10.5281/zenodo.22882677
DOI 10.5281/zenodo.22882677 Cites version 2 of 12 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies proxy-based causal inference for latent-class mean treatment effects in a fixed finite-class model. Under the Virk–Mazaheri–Wu proxy moment structure with fixed latent cardinality kk (the number of latent classes), fixed proxy dimensions, bounded observable contributions, supplied latent-arm positivity and proxy-rank margins, the target is the quotient latent-effect law νP\nu_P, the population-weighted distribution of latent-class average treatment effects with masses aggregated at coincident class means. The main population result establishes a gap-free Lipschitz modulus from the observable five-block summary S(P)S(P), the vector of proxy moment summaries, to νP\nu_P in W1W_1, the one-Wasserstein distance. The modulus covers homogeneous, partially colliding, and separated effect configurations within the same uniformly conditioned class. The statistical results turn this modulus into root-nn law estimation and confidence reporting. A nearest-summary repair estimator and a structured-lattice estimator, calibrated by the supplied constants k,dx,dz,L,π0,σ0k,d_x,d_z,L,\pi_0,\sigma_0, are Borel sample maps and achieve uniform root-nn W1W_1 tail bounds over the model class. In the fixed-dimensional unit-cost exact-real model, the lattice estimator has polynomial candidate growth in nn and centers an honest finitely represented exact-real Wasserstein outer confidence set. The associated cluster report merges empirically unresolved support points and reports compatible aggregate mass intervals. On separated gap-local strata, effect-ordered latent masses have clipped inverse-gap risk min⁡{1,(n g)−1}\min\{1,(\sqrt n\,g)^{-1}\}, where gg is the local effect-gap scale. Local and same-class lower bounds from explicit two-class proxy experiments establish root-nn sharpness for quotient-law estimation and the inverse-gap rate for ordered weights on the displayed two-class specialization.

Latent Stationary Overlap and Minimax Off-policy Evaluation in Finite POMDPs

CSWP-2026-021 · Stat · 13 Sep 2026 ✓ verified in Lean 4

The verified all-radius minimax frontier is a substantial contribution, but the manuscript requires major restructuring and several foundational clarifications before its scope and argument are publication-ready.

Main criticisms
  • Setup and estimands The statistical experiment is not coherently specified in the manuscript. The model-class tuple contains “L(O_T),” an observed-data law, while the Markov, sequential-ignorability, and stationary-start restrictions are statements about the latent full trajectory. Later the appendix separately introduces a full law and observation map, but those objects are needed to understand the model class itself.
  • Setup and estimands The presentation-level definition “Behavior and target stationary joint-state laws” defines P_p as a Markov kernel and d_p as its stationary law before b and e have been required to be probability vectors and before contraction has supplied existence and uniqueness. For a raw policy carrier, P_p need not yet be stochastic, so the definition is logically premature as written.
  • Proofs and auxiliary constructions The appendix reproduces lengthy, mechanically structured proof transcripts whose local bookkeeping overwhelms the econometric argument. The central ideas—radius-sensitive bias, window covariance control, terminal-depth filtering, KL calibration, and two-point reduction—are difficult to recover, despite the formal correctness of the statements.

13 findings — 3 major · 8 minor · 2 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

This paper characterizes the minimax risk for stationary off-policy evaluation in finite POMDPs under stationary start, bounded rewards, known observed-state policies, one-step policy overlap, geometric mixing, and latent stationary overlap. It shows that clipped partial-history importance weighting attains the resulting risk surface, including the Hu-Wager partial-history rate for fixed latent overlap constant and the parametric transition on small overlap-distance scales.

@techreport{causalsmith_stat_pomdp_latent_overlap_minimax_v1,
  title       = {Latent stationary overlap and minimax off-policy evaluation in finite {POMDPs}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-021},
  doi         = {10.5281/zenodo.22882673},
  url         = {https://causalsmith.org/papers/stat_pomdp_latent_overlap_minimax_v1/},
  note        = {Version 1, 13 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit f2e7e97}
}
CausalSmith (2026). “Latent stationary overlap and minimax off-policy evaluation in finite POMDPs.” CausalSmith Working Paper CSWP-2026-021, version 1 (13 Sep 2026). https://doi.org/10.5281/zenodo.22882673
DOI 10.5281/zenodo.22882673 Cites version 1 of 13 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper characterizes the observable-data minimax risk for stationary off-policy evaluation in finite partially observed Markov decision processes. The data are one behavior-policy trajectory of length TT, the trajectory horizon, with bounded rewards and binary actions; policies are known functions of the observed state; the joint observed–latent Markov state satisfies stationary start, sequential ignorability, one-step policy overlap L=exp⁡(ζ)L=\exp(\zeta), where ζ\zeta is the log-overlap scale, geometric total-variation contraction α=exp⁡(−1/t0)\alpha=\exp(-1/t_0), where t0t_0 is the mixing scale, and latent stationary overlap de(s)≤Cdb(s)d_e(s)\le C d_b(s), where CC is the stationary-overlap constant and ded_e and dbd_b are the target and behavior stationary joint-state laws. For sufficiently large TT, and cardinality-uniformly over finite observed and latent alphabets, the minimax squared-error risk is comparable, uniformly over C≥1C\ge1, to T−1+T−βqC2(1−β),qC=C−1C,β=22+t0ζ, T^{-1}+T^{-\beta}q_C^{2(1-\beta)}, \qquad q_C=\frac{C-1}{C}, \qquad \beta=\frac{2}{2+t_0\zeta}, where qCq_C is the overlap-distance radius, the normalized distance of the stationary density-ratio bound from unit overlap, and β\beta is the rate exponent set by the mixing and log-overlap scales, up to constants depending only on t0t_0 and ζ\zeta. For every fixed C>1C>1, the surface has the Hu–Wager partial-history exponent T−βT^{-\beta}; the parametric term changes the rate on shrinking overlap-distance scales qC≲T−1/2q_C\lesssim T^{-1/2}. A clipped partial-history importance weighting (PHIW) estimator with radius-calibrated history depth attains this surface when the depth is calibrated from (t0,ζ,C)(t_0,\zeta,C). A signed-depth finite POMDP pair supplies matching hard alternatives while satisfying the same overlap and contraction conditions. A finite controlled insulin-policy demonstration shows that bounded latent stationary occupancy and bounded action overlap can hold together in a stylized treatment model.

Sharp Ordinal Benefit Bounds for Survivor Compliers under Selection and Noncompliance

CSWP-2026-019 · Partial ID · 10 Sep 2026 · v2 12 Sep 2026 ✓ verified in Lean 4

The paper contains a strong verified identification and computation result, but its publication case requires a more self-contained model presentation, reproducible algorithm description, sharper competitor comparison, and clearer positioning of the highly conservative inference result.

Main criticisms
  • Setup and assumptions The central observable object is not defined operationally in the main text. The definition says that ell_i(x) and h_j(x) are "obtained from the corresponding conditional instrument contrasts," but the actual four conditional probabilities and their signs appear only deep in the proof. A reader therefore cannot reconstruct the capacity vector, check its observable implications, or implement the estimator from the definition that every subsequent theorem invokes.
  • Related work The dedicated literature section identifies the relevant strands but gives the closest competing IV ordinal-benefit result, Gabriel–Sachs–Jensen, only a brief descriptive mention. It likewise summarizes Dong–Heiler without comparing their identified objects, observable inputs, sharpness arguments, and inferential conditions theorem by theorem. Because existing work already supplies response-function IV bounds for finite ordinal outcomes and selected-complier distributional bounds, the current prose leaves the incremental novelty insufficiently adjudicated.
  • global The manuscript is organized around a large inventory of theorem-interface objects rather than the economic argument. Separate definitions of each potential variable are followed by a second system definition; threshold aggregates are defined after their first complete definition; cost objects repeat the algorithm; and the guard-event prose is followed by a formal definition saying only "all guard-event indices." This repetition obscures the short conceptual chain from observed IV contrasts to capacities, transport, sharp bounds, and inference.
  • Computation and witness The claimed linear-time algorithm is supported by the verified theorem, but the displayed construction is not reproducible as an algorithmic specification. Phrases such as "construct a maximum subflow" and "allocate the residual mass" omit the ordering, tie rules, pointer updates, stopping conditions, and residual invariants that later carry the complexity and optimality argument.
  • Estimation and guarded inference The finite-sample guard is mathematically valid, but its headline placement gives it similar prominence to the sharp identification theorem even though the paper's own calculation shows that the reported set remains [0,1] at remotely realistic sample sizes in the smallest model configuration. The practically informative inferential result presently delivered is the fixed-law directional limit, while the face-aware multiplier object is diagnostic and the uniform guard is best interpreted as an existence or fallback guarantee.

13 findings — 5 major · 7 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

Under the stated instrumental-variables independence, monotonicity, selection, overlap, and positive survivor-complier mass conditions, the paper gives sharp finite-sample formulas for the identified interval of the survivor-complier probability that treatment strictly improves an ordered outcome. It also constructs endpoint-attaining latent laws and provides plug-in estimation with a guarded confidence interval that uniformly contains the full sharp interval under fixed overlap and positive mass.

@techreport{causalsmith_pid_slate_benefit_partialtransport_v1,
  title       = {Sharp ordinal benefit bounds for survivor compliers under selection and noncompliance},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-019},
  doi         = {10.5281/zenodo.22882645},
  url         = {https://causalsmith.org/papers/pid_slate_benefit_partialtransport_v1/},
  note        = {Version 2, revised 12 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit 8f179a2}
}
CausalSmith (2026). “Sharp ordinal benefit bounds for survivor compliers under selection and noncompliance.” CausalSmith Working Paper CSWP-2026-019, version 2 (12 Sep 2026). https://doi.org/10.5281/zenodo.22882645
DOI 10.5281/zenodo.22882645 Cites version 2 of 12 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper characterizes sharp bounds on the survivor-complier strict-benefit probability in a finite ordered-outcome instrumental-variables model with treatment-induced selection. The target is the probability θ\theta, the survivor-complier benefit probability, that treatment raises the ordered outcome among units whose treatment receipt is shifted by the instrument and whose outcome would be observed under either treatment state. Under conditional instrument independence, consistency and exclusion through received treatment, instrument overlap, treatment monotonicity, covariate-specific weak selection monotonicity, and positive aggregate survivor-complier mass, observed instrument contrasts identify two cellwise selected-complier outcome capacities. Their smaller total is the survivor-complier mass, including cells with zero selected-complier capacity gap. The sharp identified set is the interval obtained by solving a branch-free exact-mass partial-transport problem in each covariate cell and aggregating the resulting threshold-cut values. The lower and upper endpoints have closed finite formulas, and endpoint-attaining full latent laws establish sharpness. The same construction gives sparse endpoint witnesses computable in O(∣X∣K)O(|\mathcal X|K) operations for finite covariate support X\mathcal X and KK ordered outcome levels, with O(∣X∣K2)O(|\mathcal X|K^2) work when full coupling matrices are materialized. A three-level synthetic witness yields the sharp interval [0,0.7][0,0.7]. For i.i.d. samples from the fixed finite-slate structural class with fixed finite support, prescribed screening thresholds ηn\eta_n satisfying ηn→0\eta_n\to0 and n ηn→∞\sqrt n\,\eta_n\to\infty, a common instrument-overlap bound, and a uniform positive lower bound on aggregate survivor-complier mass, a projected screened plug-in estimator has a fixed-law Hadamard directional limit and a deterministic guarded confidence interval contains the full sharp identified interval uniformly.

Graph-Adaptive Bernoulli Design for Bipartite Interference

CSWP-2026-006 · Experimentation · 16 Jul 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The formal results appear internally sound and carefully caveated, but the paper does not yet demonstrate enough econometric value or practical relevance for a leading journal.

Main criticisms
  • abstract; Introduction; Discussion and extensions The paper still asks readers to care about a new experimental design rule without showing that it improves any decision-relevant object in realistic bipartite experiments. The manuscript says the criterion is useful because it is observable, convex, and conservative, but gives no numerical evidence on realized variance, interval width, exposure frequencies, coverage, sensitivity to the positivity floor, or computation time.
  • Introduction and positioning The novelty case remains too narrow for a leading econometrics journal. The manuscript now correctly states that it does not prove a sharp minimax design theorem, but that leaves the central contribution as optimizing a conservative upper bound whose tightness and practical ranking properties are unknown.
  • Discussion and extensions The design checklist is helpful but still not operational enough for applied use. It tells readers to prespecify acceptable floors, tolerances, and exposure-probability minima, but gives no default diagnostics or calibration guidance for when the envelope is too loose or when clustering should dominate.

8 findings — 3 major · 4 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

Under the stated interference, independence, boundedness, positivity, graph-dependence, feasibility, and nondegenerate-variance assumptions, the paper develops a graph-adaptive choice of heterogeneous Bernoulli treatment probabilities for bipartite experiments, yielding a computable variance envelope, an optimal feasible design, and asymptotically conservative Wald intervals for the all-treated-versus-all-control effect.

@techreport{causalsmith_exp_bipartite_minimax_design_v1,
  title       = {{Graph-Adaptive} {Bernoulli} {Design} for {Bipartite} {Interference}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jul,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-006},
  doi         = {10.5281/zenodo.22882608},
  url         = {https://causalsmith.org/papers/exp_bipartite_minimax_design_v1/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit 8d3b465}
}
CausalSmith (2026). “Graph-Adaptive Bernoulli Design for Bipartite Interference.” CausalSmith Working Paper CSWP-2026-006, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882608
DOI 10.5281/zenodo.22882608 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

We study independent, heterogeneous Bernoulli assignment in bipartite experiments, where intervention units are assigned treatment and outcome units may depend on assignments in known intervention neighborhoods. For the finite-population contrast between all-treated and all-control neighborhood outcomes, we analyze an exposure-weighted Hájek estimator under bipartite neighborhood interference. With bounded potential outcomes, we derive a graph-and-design-dependent envelope that upper-bounds the design-based variance scale of its linearization and can be minimized subject to a positivity-constrained expected-treatment budget. Under the full assumptions of Theorem 4 and Theorem 5, the envelope-optimal design supports asymptotic normality and a graph-only Wald scale yields asymptotically conservative coverage. We also give conditions under which heterogeneous probabilities strictly improve the envelope relative to homogeneous assignment and study a separable overlap-based surrogate under an admissible budget and bounded outcome degree. A complementary unbounded-degree construction shows that degree dispersion and comparable surrogate weights alone do not control this approximation. Together, the results provide an outcome-model-free design criterion with convex, inferential, and separation guarantees, and identify approximation questions for surrogate criteria.

Minimax Inference for Threshold Modified Treatment Policies with Continuous Treatments

CSWP-2026-014 · Stat · 7 Sep 2026 · v2 12 Sep 2026 ✓ verified in Lean 4

The paper delivers a potentially important and rigorously verified minimax frontier, but a contradictory sharp-modulus program, stale verification metadata, understated model scope, and absent finite-sample evidence prevent publication in its current form.

Main criticisms
  • Discussion, limitations, and open questions The sharp-length question is internally incompatible with the object it proposes to calibrate. The manuscript explicitly explains that the radius is unweighted, omits retained-mean and atom-mass uncertainty, and has its natural local-regression scale h_n^β at zero threshold; nevertheless the remark requires its length and 2W_{n,B_n} to match r_n for every regime, including δ_n=0 through the q=0 case. At zero and regular/critical thresholds, the missing atom weight generally makes this component radius much larger than the complete clamp frontier, so the displayed matching limits cannot describe the defined interval. The requirement of one constant c across all phase limits is also unsupported and generally implausible when the relative regular and local contributions vary with the phase limit.
  • Proofs and verification note The manuscript states that the checks were run at commit 6f11a039165a75921bbb6499f3d9bf40f3f404f8, while the authoritative current verification contract identifies commit 84c32d60cb53b2f1cdef917448a45b96b08296ad. The reproduction and archive narrative therefore does not identify the source state supporting the submitted statements.
  • Discussion, limitations, and open questions The paper proposes a nonstandard split-sample estimator, total-Gram stabilization rule, fallback branch, and honest interval, yet supplies only a one-cell exponent calculation. The manuscript itself identifies fallback frequency, realized coverage, and expected length as potentially sensitive to constants and moderate sample size. Without numerical evidence, readers cannot assess whether the proposed procedure is usable in the regimes for which it is advocated.
  • abstract and Introduction The prose repeatedly describes the density restriction as “lower-tail” thinning or as holding “near the lower endpoint,” but the stated model imposes c_-a^κ ≤ π_x(a) ≤ c_+a^κ almost everywhere on the entire interval [0,1]. This understates a material global restriction and can lead readers to believe the frontier has been proved under a merely local boundary condition.

10 findings — 4 major · 5 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

For lower-threshold modified treatment policies with continuous treatments, the paper characterizes the minimax risk and honest confidence-interval length when boundary treatment density vanishes polynomially. Under finite strata, Hölder regression smoothness, and density exponent conditions, a stabilized local-polynomial estimator attains the frontier n−1/2+δnκ+1hnβn^{-1/2}+\delta_n^{\kappa+1}h_n^\beta, with matching lower bounds and regime separation.

@techreport{causalsmith_stat_lmtp_threshold_atom_frontier_v1,
  title       = {Minimax inference for threshold modified treatment policies with continuous treatments},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-014},
  doi         = {10.5281/zenodo.22882669},
  url         = {https://causalsmith.org/papers/stat_lmtp_threshold_atom_frontier_v1/},
  note        = {Version 2, revised 12 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit 6f1bf59}
}
CausalSmith (2026). “Minimax inference for threshold modified treatment policies with continuous treatments.” CausalSmith Working Paper CSWP-2026-014, version 2 (12 Sep 2026). https://doi.org/10.5281/zenodo.22882669
DOI 10.5281/zenodo.22882669 Cites version 2 of 12 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

Continuous-treatment modified treatment policies can map a range of natural treatment values to a single assigned dose. This paper studies the exact lower-threshold clamp, which sends treatments below a deterministic threshold δn\delta_n, the policy threshold, to δn\delta_n and leaves larger treatments unchanged. Under a fixed finite-stratum observed-data model with bounded outcomes, Hölder outcome regression with known exponent β\beta and radius LL, and polynomial lower-tail treatment-density thinning c−aκ≤πx(a)≤c+aκc_-a^\kappa\leq\pi_x(a)\leq c_+a^\kappa with known exponent κ\kappa, the post-policy law contains a moving atom. The target is the observed clamp mean θδn(P)\theta_{\delta_n}(P), the sum of the retained mean above the threshold and an atom-regression product at the boundary. The paper establishes matching minimax absolute-risk and uniformly honest expected-length rates rn=n−1/2+δnκ+1hnβ,nhn2β+1(δn+hn)κ≍1, r_n=n^{-1/2}+\delta_n^{\kappa+1}h_n^\beta, \qquad n h_n^{2\beta+1}(\delta_n+h_n)^\kappa\asymp1, attained by a split-sample total-Gram stabilized local-polynomial estimator and a bias-aware interval. The phase diagram identifies the regular, critical, atom-dominated, fixed-threshold, and zero-threshold regimes induced by δn\delta_n. A full-data lift gives the same frontier a causal interpretation under simultaneous consistency, conditional exchangeability, and structural-mean continuity on the threshold range. A continuity-only model yields the separate frontier n−1/2+δnκ+1n^{-1/2}+\delta_n^{\kappa+1}, with an explicit elbow and fixed-positive-threshold behavior.

Second-order Minimax Risk in Multi-arm Binary Randomization

CSWP-2026-018 · Experimentation · 10 Sep 2026 · v2 12 Sep 2026 ✓ verified in Lean 4

The paper delivers a substantial verified fixed-K minimax-rate result and useful finite-dimensional reductions, but its novelty positioning, organization, and presentation of several certificate and attainment claims require major repair before the contribution is journal-ready.

Main criticisms
  • Related work The discussion of the closest competitor is too cursory to establish the paper's incremental contribution. The sentence “Sudijono, Dobriban, and Tchetgen Tchetgen obtains a sharp two-arm binary expansion” omits the competitor's scalar reduction, exact Airy constant, nonlinear attaining rule, and least-favorable-prior characterization. Consequently, readers cannot see that this paper supplies an internal matching-order converse and a fixed-K extension while leaving the sharp multi-arm constant and optimizer open.
  • global The manuscript repeatedly repackages the same results in the exact-reduction theorem, the multi-arm specialization theorem, the universal-rate theorem, the certificate frontier theorem, and the final attainment theorem. Several packages add automatic consequences such as compactness of the cluster set or repeat complete certificate-transfer statements. This makes the logical contribution difficult to identify and gives bookkeeping results the same prominence as the core econometric theorems.
  • Setup and the finite orbit game The “Hull risk-value inclusion” clause is tautological: it proves membership by choosing the same design and estimator. The subsequent “strict treatment-domain separation” merely observes that a real-valued two-arm schedule set is not a subset of a binary K-arm schedule set with a different ambient type. Neither clause provides a risk comparison or substantive extension, yet both appear inside a main theorem alongside the actual orbit and rate results.
  • Three-arm diagnostics and support-specific conclusions The certificate-frontier statement is type-confused: “there exist sequences G_n, P_n, and L_n such that G_n is an exact rational grid-certificate cluster” treats the scalar normalized grid value G_n as the certificate cluster defined earlier as a collection of primal weights, dual multipliers, an estimator, and derived bounds. The proof repeats the same categorical mismatch.

10 findings — 4 major · 5 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

For fixed finite populations with binary outcomes and multi-arm zero-sum contrasts, the paper reduces the labeled design-based minimax problem to a finite orbit game and identifies the first-order risk constant C0(c)/nC_0(c)/n. Under assignment probabilities proportional to ∣ca∣|c_a|, a clipped shrinkage estimator improves the worst-case risk at the second-order n−4/3n^{-4/3} scale, with matching order certified by an internal two-arm information argument.

@techreport{causalsmith_exp_multiarm_secondorder_minimax_frontier_v1,
  title       = {Second-order minimax risk in multi-arm binary randomization},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = sep,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-018},
  doi         = {10.5281/zenodo.22882624},
  url         = {https://causalsmith.org/papers/exp_multiarm_secondorder_minimax_frontier_v1/},
  note        = {Version 2, revised 12 Sep 2026. Machine-generated; formal statements verified in Lean 4 at commit 715add8}
}
CausalSmith (2026). “Second-order minimax risk in multi-arm binary randomization.” CausalSmith Working Paper CSWP-2026-018, version 2 (12 Sep 2026). https://doi.org/10.5281/zenodo.22882624
DOI 10.5281/zenodo.22882624 Cites version 2 of 12 Sep 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies design-based minimax risk for randomized experiments with a fixed number KK of treatment arms, nn labeled units, binary potential outcomes, and a prespecified nonzero zero-sum treatment contrast. The game, which optimizes jointly over the assignment law and an estimator clipped to the estimand's own range against all complete labeled schedules, has the same value as an explicit finite game whose states are the response-type count vectors over the 2K2^K binary response types. For every fixed nonzero zero-sum contrast cc, the first-order minimax constant is C0(c)=Lc2/4C_0(c)=L_c^2/4, where LcL_c is the ℓ1\ell_1 norm of the contrast, and the contrast-weighted allocation qa⋆∝∣ca∣q_a^\star\propto |c_a| attains the finite-sample upper bound C0(c)/nC_0(c)/n with a projected Horvitz–Thompson rule. The same allocation, paired with an explicit clipped-shrinkage estimator, improves on that envelope by a positive multiple of n−4/3n^{-4/3} for all sufficiently large nn, while an embedded two-arm Bayesian information argument gives a matching O(n−4/3)O(n^{-4/3}) converse. Thus 4/34/3 is the universal second-order exponent for every fixed nonzero zero-sum binary contrast. For rational contrasts the paper builds a finite linear program whose primal and dual solutions supply a feasible procedure and a prior — a matching pair of computable upper and lower bounds on the minimax risk, with all data exactly rational — and transfers those brackets to every real contrast by a square-root risk continuity bound. A three-arm example shows that at n=3n=3, for the contrast (1,−1/2,−1/2)(1,-1/2,-1/2) under one fixed allocation, estimators using the full observed arm label and outcome attain a strictly smaller worst-case risk than estimators using only a scalar signed score.

A Lower-Bound Calibration for Joint Margin–Overlap Decay in Offline Policy Learning

CSWP-2026-001 · Stat · 13 Jun 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The core lower-bound calibration appears correct and potentially useful, but the paper is still too narrow and too conditionally framed, with several claim-fidelity and presentation problems that prevent publication as written.

Main criticisms
  • global The paper's delivered contribution is primarily a two-point observed-law lower-bound calibration plus a conditional upper bound for one supplied-nuisance procedure. As written, this is not yet a complete minimax rate characterization and may be too narrow for a leading econometrics journal unless the paper either proves more or more sharply explains why this lower-bound calibration itself changes practice or theory.
  • related literature The comparison to the closest weak-overlap and clipping literature remains too high level. Table 1 says Liu2026 is closest to the feasible upper-bound side, but it does not compare target, overlap condition, clipping schedule, nuisance assumptions, or rate exponent closely enough to establish novelty.
  • Section 4 The upper bound is still too easy to read as an achievability result. The sentence "Theorem \ref{obj:oeq:feasible-upper} is the promised conditional achievability statement" overstates what is delivered because the theorem assumes supplied nuisance estimators, boundedness, fixed folds, oracle-in-class, and primitive localized envelope and offset-envelope conditions.
  • Verification Note The sentence "The theorem statements and proofs described above are machine-checked in Lean 4" is too broad relative to the verification contract: the Le Cam chi-square testing lemma is listed with no Lean declaration attached.
  • Setup and Assumptions The title "Bounded potential outcomes" and the sentence "Potential-outcome notation Y(a) is used for the boundedness assumption below" do not match the verified assumption, which bounds the observed outcome and observed-law treatment-arm conditional means, not potential outcomes.

10 findings — 5 major · 4 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

This paper proves a minimax regret lower-bound calibration for offline policy learning under joint margin-overlap decay. Its conditional analysis of a specified clipped cross-fitted AIPW rule with supplied nuisances matches the lower-bound exponent in the nonbinding nuisance-and-clipping regime and gives the procedure’s nuisance-limited exponent in the binding regime.

@techreport{causalsmith_stat_policy_regret_margin_overlap_v1,
  title       = {A {Lower-Bound} {Calibration} for {Joint} {Margin--Overlap} {Decay} in {Offline} {Policy} {Learning}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jun,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-001},
  doi         = {10.5281/zenodo.22882671},
  url         = {https://causalsmith.org/papers/stat_policy_regret_margin_overlap_v1/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit 8d3b465}
}
CausalSmith (2026). “A Lower-Bound Calibration for Joint Margin–Overlap Decay in Offline Policy Learning.” CausalSmith Working Paper CSWP-2026-001, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882671
DOI 10.5281/zenodo.22882671 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies offline policy learning for deterministic treatment rules when overlap can deteriorate near either propensity boundary in the same region where the treatment contrast is small. The target is observed-law welfare regret, defined directly from the conditional treatment contrast. The law class imposes bounded outcomes, positivity, a margin condition for small treatment contrasts, a zero-effect convention, and a joint overlap-decay restriction that ties weak treatment-arm information to the small-contrast region. Under the margin-window normalization and the auxiliary calibrations used by the two-point construction, the paper proves an observed-law minimax lower bound Mn(α,γ)≥c n−r⋆(α,γ),r⋆(α,γ)=1+α2+α+βα,γ,βα,γ={0,γ=0,αγ/(α+1),γ>0. M_n(\alpha,\gamma)\ge c\,n^{-r_\star(\alpha,\gamma)}, \qquad r_\star(\alpha,\gamma) = \frac{1+\alpha}{2+\alpha+\beta_{\alpha,\gamma}}, \qquad \beta_{\alpha,\gamma} = \begin{cases} 0,&\gamma=0,\\ \alpha\gamma/(\alpha+1),&\gamma>0. \end{cases} The exponent is obtained by balancing margin mass, local contrast size, and the probability of observing the informative treatment arm. The paper also gives a conditional analysis of a specified clipped cross-fitted AIPW empirical welfare rule with supplied nuisance estimates satisfying explicit rate, boundedness, cross-fitting, and localized empirical-process conditions. In the nonbinding nuisance-and-clipping regime, this conditional upper exponent matches the lower-bound exponent up to logarithmic factors; in the binding regime, the analysis gives the rule's nuisance-limited exponent.

Chebyshev Rollout Schedules for Polynomial Extrapolation under Low-Order Interference

CSWP-2026-004 · Experimentation · 2 Jul 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The verified mathematical core is coherent and potentially useful, but the paper is too narrowly delivered relative to its econometric framing and several prose claims and objects around interference, exact rollout risk, and covariance structure need substantial tightening before publication.

Main criticisms
  • global The title and some framing oversell the result as an interference result. The title is "Chebyshev Rollout Schedules for Polynomial Extrapolation under Low-Order Interference," but the formal layer proves results under an imposed polynomial mean curve and diagonal variance envelope, not under a primitive low-order interference model.
  • setup and discussion The exact nested covariance object is not self-contained. The setup says "The specific assignment law... is not further specified here" and that the exact covariance is "not part of the verified formal layer," but Remark oe:q later refers to "the true design covariance matrix of the monotone Bernoulli rollout," which has not been defined.
  • contribution and significance The delivered contribution is narrower than the econometric motivation suggests. The paper solves the diagonal-envelope linear-unbiased polynomial extrapolation design problem, while exact covariance optimality, primitive interference derivation, finite-population rounding effects, and misspecification risk are outside the formal contribution.
  • setup The notation for the law class suppresses important dependence on the schedule, design, and variance scale. Definition P_beta pins E[bar Y_j]=m_P(p_j), so the class changes with p, yet R_exact writes sup_{P in P_beta} while p is being optimized.
  • setup and assumptions The round-mean variance envelope is central but only heuristically justified. The text says it holds under "weak" or "bounded-neighborhood" dependence, but no formal sufficient condition or citation is provided, and bounded outcomes alone are correctly noted to be insufficient.

11 findings — 5 major · 5 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

Under the stated rollout-consistency, polynomial-mean, variance-envelope, and oversampling conditions, Chebyshev-Lobatto measurement schedules control the variance cost of extrapolating from partial rollout q<1q<1 to full adoption at the minimax exponential rate for the envelope problem. For the exact nested rollout covariance, the same schedule is proved rate-feasible.

@techreport{causalsmith_exp_rollout_chebyshev_minimax_tv_envelop_0fx512j,
  title       = {Chebyshev {Rollout} {Schedules} for {Polynomial} {Extrapolation} under {Low-Order} {Interference}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jul,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-004},
  doi         = {10.5281/zenodo.22882630},
  url         = {https://causalsmith.org/papers/exp_rollout_chebyshev_minimax_tv_envelope_rollout_design/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit c3c6187}
}
CausalSmith (2026). “Chebyshev Rollout Schedules for Polynomial Extrapolation under Low-Order Interference.” CausalSmith Working Paper CSWP-2026-004, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882630
DOI 10.5281/zenodo.22882630 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies how to place measurement rounds in a finite-population rollout experiment when the target is the full-adoption contrast but the rollout budget stops at a treated fraction q<1q<1. Under static rollout consistency, a degree-β\beta polynomial restriction on the rollout mean curve, and a common round-mean variance envelope—restrictions that low-order interference motivates but that we impose rather than derive from a microfounded interference model—unbiased linear estimation reduces to a polynomial extrapolation problem from [0,q][0,q] to 11. For a fixed schedule, the worst-case variance over the positive-semidefinite covariance class meeting this diagonal envelope equals the variance scale times the squared ℓ1\ell^1 (total-variation) norm of the polynomial-exact weights. The resulting diagonal-envelope amplification criterion, defined over linear unbiased estimators, attains its minimax value up to multiplicative constants at shifted Chebyshev-Lobatto measurement fractions when q≤qmax⁡<1q\le q_{\max}<1 and the number of rollout intervals satisfies k≥cβk\ge c\beta for some c>1c>1. The minimax amplification is of order ((1+1−q)2q)2β, \left(\frac{(1+\sqrt{1-q})^2}{q}\right)^{2\beta}, up to constants depending only on qmax⁡q_{\max} and the oversampling ratio. For the exact nested rollout covariance problem, the same Chebyshev schedule is rate-feasible through the envelope upper bound, and exact optimality under the true rollout covariance structure is posed as a separate covariance-specific design question.

Generic Separation of Axis-Normalized Latent-Source Representations by Higher-Order Cumulants

CSWP-2026-005 · Exact ID · 12 Jul 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The verified mathematical core appears correct and carefully scoped, but the manuscript as written is too narrow, hard to interpret economically, and insufficiently positioned to merit publication in a leading econometrics journal without substantial revision.

Main criticisms
  • contribution and positioning The paper still does not make a compelling case that the result is important for econometrics rather than primarily algebraic statistics. The introduction says the target is useful when "an institutional model, factor design, or separate model-selection analysis fixes the number of common shock slots," but it gives no concrete econometric design, empirical model, or identification workflow where excluding the opposite fixed-slot cumulant fiber changes an economic conclusion.
  • introduction and related work The novelty comparison remains too qualitative. The paper repeatedly contrasts order 2m+2 or 2m+1 with the order-2m+3 result of ChenPengHuangEtAl2025Direction, but the reader still cannot tell whether the order difference is due to a genuine information gain, the maintained axis normalization, the source-count convention, the definition of confounding, or the change from a test to a complex/real fiber-exclusion target.
  • proofs of the main results The deferred proofs still read partly like implementation notes rather than mathematical exposition. Phrases such as "retained-coordinate parameter spaces," "truncated coordinate chart," "determinant polynomial," and "moment-interior nonvanishing statement" are used as if they were already familiar mathematical lemmas, but the paper does not state the exact external propositions being invoked in conventional theorem/lemma form.

7 findings — 3 major · 3 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

With a fixed number of middle source slots and a maintained axis normalization, population cumulants generically exclude the opposite representation fiber at order 2m+22m+2, and at order 2m+12m+1 over the real feasible region when m≥3m\ge3.

@techreport{causalsmith_eid_lingam_direction_min_order_v1_trunca_0ijfc5p,
  title       = {Generic {Separation} of {Axis-Normalized} {Latent-Source} {Representations} by {Higher-Order} {Cumulants}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jul,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-005},
  doi         = {10.5281/zenodo.22882604},
  url         = {https://causalsmith.org/papers/eid_lingam_direction_min_order_v1_truncated_cumulant_minimality/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit c3c6187}
}
CausalSmith (2026). “Generic Separation of Axis-Normalized Latent-Source Representations by Higher-Order Cumulants.” CausalSmith Working Paper CSWP-2026-005, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882604
DOI 10.5281/zenodo.22882604 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

We study a population identification question for a bivariate latent linear non-Gaussian model in which the number mm of middle source slots and an axis normalization are maintained inputs. The representation-level target asks whether a truncated joint-cumulant vector through order LL can also arise from the opposite axis-normalized arrow convention with the same mm. At L=2m+2L=2m+2, the unordered finite loading-slope support is generically recovered and the opposite cumulant-map fiber is empty, while the same-arrow parametrization retains ambiguity beyond middle-slot relabelling. The closure of the opposite-arrow compatibility locus has codimension exactly one in each arrow image variety. For m≥3m\ge3, generic real opposite-fiber exclusion already holds at order 2m+12m+1; loading-support recovery uses the apolar order 2m+22m+2. There exist Euclidean-open parameter neighborhoods whose intersections with the feasible regions are nonempty, sharpening the population geometry of the maintained representation class.

Forbidden Comparisons in Fixed-Effect Poisson Difference-in-Differences

CSWP-2026-007 · Panel · 16 Jul 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The core deterministic sign-reversal result appears correct and potentially useful, but the manuscript is not yet strong enough for a leading econometrics journal because the empirical bridge, exposition of the central sign argument, and positioning remain underdeveloped.

Main criticisms
  • intro/setup/discussion The empirical bridge remains incomplete. The manuscript says, "This is the population analogue of the coefficient an applied pooled unit-and-time fixed-effect PPML regression is trying to summarize," while later conceding that "No theorem below connects that sample coefficient to beta_N^star(delta) or beta^star(delta), and no sampling inference is developed." This is honest, but for a leading econometrics journal it leaves the main applied relevance under-supported.
  • main results/appendix The proof exposition of the central global sign result is too compressed. The key step says, "At the frontier fit with treatment coefficient set to zero, the nuisance scores vanish and the treatment score is Phi/(T M)," but the manuscript does not give readers a self-contained derivation of this frontier fit or the concavity/score-crossing argument that turns the score sign into the sign of beta^star.
  • discussion The paper establishes that the sign-reversal region is nonempty, but it gives little sense of how economically relevant the region is. The statement that the example is "relevant for otherwise conventional staggered-adoption designs" is stronger than the current evidence, since the manuscript mostly presents one witness and the algebraic frontier.

6 findings — 3 major · 3 minor   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

In staggered-adoption DiD with proportional effects, pooled fixed-effect Poisson can give a negative limiting treatment coefficient even when every treated cohort-time effect is positive because the coefficient is a misspecified projection with signed residual weights. The paper also proves exact recovery under a common proportional effect and characterizes the sign under the stated collapsed-rank and untreated-mean conditions.

@techreport{causalsmith_panel_ppml_forbidden_comparison_v1,
  title       = {Forbidden {Comparisons} in {Fixed-Effect} {Poisson} {Difference-in-Differences}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jul,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-007},
  doi         = {10.5281/zenodo.22882639},
  url         = {https://causalsmith.org/papers/panel_ppml_forbidden_comparison_v1/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit c3c6187}
}
CausalSmith (2026). “Forbidden Comparisons in Fixed-Effect Poisson Difference-in-Differences.” CausalSmith Working Paper CSWP-2026-007, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882639
DOI 10.5281/zenodo.22882639 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies the population treatment coefficient from pooling a staggered-adoption panel in a unit-and-time fixed-effect Poisson pseudo-likelihood. We show that this projection coefficient can be negative even when every cohort-time proportional treatment effect is positive. The failure reflects heterogeneous effects interacting with the fixed-effect PPML score: a four-cohort design with equal shares and flat untreated means provides an explicit sign reversal, while a positive-weight proportional treatment-on-the-treated target remains positive. The result characterizes the interpretation of the deterministic population projection targeted by the pooled criterion.

Exact Randomized Designs for Two-Block Interference Experiments

CSWP-2026-003 · Experimentation · 2 Jul 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The verified mathematics appears coherent and exact, but the submission is too narrow and too weakly connected to econometric design or inference to earn publication in a leading econometrics journal as written.

Main criticisms
  • global The contribution is too narrow for the econometrics framing as written. The paper explicitly says, "we take the objective as given—a covariance design criterion ... and do not derive it from a specific estimand, estimator, variance, or welfare target," but the title, introduction, and discussion frame the work as a contribution to randomized experiments under interference.
  • abstract The abstract presents the norm term as a "Schatten–2 robustness penalty" and says iid assignment is an "asymptotic target as the robustness weight diverges." The verified layer proves facts about a Schatten–2 norm penalty; it does not establish statistical robustness, minimax robustness, robustness to graph misspecification, or robustness of inference.
  • discussion The paper repeatedly draws design-practice implications that exceed the formal scope. For example, "For interference-aware randomized design, the practical implication is that solving the covariance relaxation is not by itself enough" is true only for this two-block, sign-symmetric, block-exchangeable, no-fixed-count criterion, not for interference-aware randomized design generally.
  • positioning The novelty claim around the parity-truncated characterization is under-supported. The discussion says, "We are not aware of a prior statement of this parity-truncated characterization..." but the paper does not engage enough with finite exchangeability, Rademacher moment problems, or low-dimensional correlation-polytope slices to substantiate novelty at journal standard.
  • verification note The verification note lists Lean version, mathlib revision, namespace, and declaration names, but it does not provide an accessible artifact, build command, or archival identifier. For a paper relying on machine checking, this is insufficient reproducibility information.

13 findings — 5 major · 7 minor · 1 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

This paper characterizes finite randomized {±1}\{\pm1\} implementations of a relaxed covariance-based design problem for two equal homophilous blocks, proving exactness in a strict cut region and giving a finite active-set loss formula. Under the stated parity and robustness conditions, odd block size yields regions where the unique relaxed optimum is separated from the implementable design class by strictly positive loss.

@techreport{causalsmith_exp_design_pm1_exactness_boundary_v1_exp_0bg38kf,
  title       = {Exact {Randomized} {Designs} for {Two-Block} {Interference} {Experiments}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jul,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-003},
  doi         = {10.5281/zenodo.22882616},
  url         = {https://causalsmith.org/papers/exp_design_pm1_exactness_boundary_v1_experimentation/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit c3c6187}
}
CausalSmith (2026). “Exact Randomized Designs for Two-Block Interference Experiments.” CausalSmith Working Paper CSWP-2026-003, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882616
DOI 10.5281/zenodo.22882616 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies when a covariance relaxation for interference-aware randomized design is exactly implementable by a finite {±1}\{\pm1\} assignment law. The setting is a design-based finite-population model with two equal homophilous communities, sign-symmetric assignment, and a block-weighted graph. The design objective combines a graph Laplacian term, a Laplacian-pseudoinverse term, a Schatten–2 robustness penalty, and an aggregate balance penalty. A symmetry reduction shows that the relaxed problem is exactly the optimization over a two-parameter block elliptope, while the implementable problem is the same objective restricted to covariances induced by block-exchangeable sign designs. The relaxation is exact in a strict cut region, where the unique relaxed optimum is generated by a two-point randomized design. Independent fair assignment is implementable throughout the model and is a finite-robustness relaxed optimum precisely on the affine locus a+3b=2ma+3b=2m and r=2b(a+b)r=2b(a+b); elsewhere it emerges as the asymptotic target as the robustness weight diverges. For odd community size, parity creates an open parameter region with a uniquely optimal relaxed covariance and strictly positive implementability loss. Finally, an exact finite active-set formula computes the loss for all admissible m,a,b,r,κm,a,b,r,\kappa, with zero loss for even community size.

A Minimax Lower Bound for Interior Dose-Response Estimation

CSWP-2026-002 · Stat · 29 Jun 2026 · v2 25 Aug 2026 ✓ verified in Lean 4

The verified result is a correct but narrow one-sided lower bound, while the paper’s terminology and positioning still invite a stronger minimax/frontier interpretation than the results deliver.

Main criticisms
  • global The paper’s actual contribution is too narrow for the claimed econometric payoff. The verified mathematical content is a two-point lower bound obtained by perturbing the treatment regression plus algebraic comparisons of exponents; it does not establish the minimax rate, does not prove a same-class upper bound, and leaves the main low-smoothness upper frontier open.
  • global The contribution is overstated by repeated use of sharp/minimax/frontier terminology. Examples include the theorem title "Sharp Pointwise Lower Bound", the theorem title "Smooth Covariate Minimax Floor", and the remarks titled "Beta frontier handle" and "Full beta frontier". The formal layer proves only a one-sided lower bound plus algebraic comparisons; it does not prove sharpness, a minimax rate, or a frontier.
  • main results The imported comparator in "Published HOIF Comparator" remains too strong as written for a self-contained econometrics paper. It introduces a bespoke Bonvini-Kennedy localized-regularity class, asserts subset relations, and summarizes their conditions in ways a reader cannot verify from the present definitions alone.
  • setup and observed-data class The model class definition includes causal assumptions and iid sampling as properties of laws P, which blurs the statistical object. A single-observation observed-data law does not itself satisfy iid sampling, and potential-outcome assumptions are not properties of the observed-data law alone without specifying an augmented law.
  • appendix / auxiliary lemmas and proofs The proof prose contains claims not present as independently stated assumptions or visible lemmas in the manuscript, for example the existence of a "bump Holder gate", membership of constructed witnesses in the class, and KL tensorization details. Even if covered by the verification contract, the submitted paper presents these as ordinary mathematical proofs without enough self-contained detail.

12 findings — 5 major · 5 minor · 2 nit   Read the entire review →

The paper may have been revised since this report. One automated referee’s estimate — not peer review.

Under the stated Hölder smoothness, positivity, boundedness, and strict baseline slack assumptions, the minimax MSE for estimating an interior continuous-treatment dose-response partial mean is bounded below at rate n−2α/(2α+1)n^{-2\alpha/(2\alpha+1)}, independently of treatment-density smoothness β\beta.

@techreport{causalsmith_stat_dose_response_minimax_holder_anisot_16btxfy,
  title       = {A {Minimax} {Lower} {Bound} for {Interior} {Dose-Response} {Estimation}},
  author      = {{CausalSmith}},
  year        = {2026},
  month       = jun,
  institution = {CausalSmith Working Papers},
  number      = {CSWP-2026-002},
  doi         = {10.5281/zenodo.22882583},
  url         = {https://causalsmith.org/papers/stat_dose_response_minimax_holder_anisotropic_converse/},
  note        = {Version 2, revised 25 Aug 2026. Machine-generated; formal statements verified in Lean 4 at commit 8d3b465}
}
CausalSmith (2026). “A Minimax Lower Bound for Interior Dose-Response Estimation.” CausalSmith Working Paper CSWP-2026-002, version 2 (25 Aug 2026). https://doi.org/10.5281/zenodo.22882583
DOI 10.5281/zenodo.22882583 Cites version 2 of 25 Aug 2026. One DOI per paper; the record it opens lists every version.
Abstract

This paper studies lower bounds for estimating an interior continuous-treatment dose-response partial mean. The target is θP(t0)=∫μP(t0,x)pX,P(x) dx, \theta_P(t_0)=\int \mu_P(t_0,x)p_{X,P}(x)\,dx, interpretable as a causal dose-response mean under consistency, no unmeasured confounding, and local positivity. Over the same Hölder dose class Pα,β,s(M,c0,ε0,t0)\mathcal P_{\alpha,\beta,s}(M,c_0,\varepsilon_0,t_0), under iid sampling, bounded outcomes, an interior evaluation dose, the stated anisotropic Hölder restrictions, and a strict-slack baseline condition, we prove the same-class minimax lower bound Rn{Pα,β,s(M,c0,ε0,t0),t0}≥c n−2α/(2α+1) R_n\{\mathcal P_{\alpha,\beta,s}(M,c_0,\varepsilon_0,t_0),t_0\} \ge c\,n^{-2\alpha/(2\alpha+1)} for all sufficiently large nn. The lower-bound exponent is the one-dimensional treatment-regression exponent and holds for every fixed β>0\beta>0, although the constant and slack-baseline feasibility may depend on β\beta. We compare this lower floor with the published higher-order influence-function benchmark ρn=n−2α/(2α+1)∨n−2/(1+d/(4s)+1/α). \rho_n= n^{-2\alpha/(2\alpha+1)} \vee n^{-2/(1+d/(4s)+1/\alpha)}. When d≤4sd\le 4s, the lower-bound exponent equals the exponent of this published comparator, whose upper theorem concerns a distinct localized-regularity class. When 4s<d4s<d, the benchmark is governed by the covariate-smoothness term and has a strictly smaller exponent than the same-class lower-floor exponent. The paper establishes a same-class lower bound and an exact algebraic comparison with the external HOIF benchmark; a matching same-class upper analysis would complete the minimax characterization.