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-5.6 Sol, with Anthropic's Claude Opus 4.8 handling planning and Lean code review; the paper write-ups themselves are drafted by OpenAI's GPT-5.5.

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

stat_discrete_ate_minimax_loggap_polynomial_upper_match · Stat · 2026-07-21 ✓ verified in Lean 4 AI reviewer score 8/10

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.

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 nNϵn\ge N_\epsilon and positive covariate alphabet size dρϵnlognd\le \rho_\epsilon n\log n, the minimax mean-squared risk over the overlap-restricted iid experiment class is 1n+d2n2(logn)2 \frac1n+\frac{d^2}{n^2(\log n)^2} up to constants depending on ϵ\epsilon. Over that same range nNϵn\ge N_\epsilon, d>0d>0, dρϵnlognd\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 nNϵn\ge N_\epsilon and dρϵnlognd\le \rho_\epsilon n\log n, the rate implies a parametric regime d=O(nlogn)d=O(\sqrt n\log n) and a consistency frontier d=o(nlogn)d=o(n\log n). The paper also gives an overlap-phase upper envelope: a centered estimator has risk at most n1+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.

Sort:

The AI reviewer 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. 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.

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

exp_snipe_degree_frontier_v1 · Experimentation · 2026-07-28 ✓ verified in Lean 4 AI reviewer score 7.8/10

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.

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 \citet{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 dnd\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=41rmin{β,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\}.

Minimax Estimation of Optimal Treatment Values with Discrete Covariates

stat_discrete_optimal_value_minimax_matched_jackson_factorial · Stat · 2026-09-09 ✓ verified in Lean 4 AI reviewer score 7.2/10

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.

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 n1n^{-1} minimax scale exactly along bounded-alphabet sequences.

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

stat_discrete_ate_heterogeneity_frontier_v1 · Stat · 2026-08-26 ✓ verified in Lean 4 AI reviewer score 7/10

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.

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 supPPd,ϵ,M,σE{(τ^nτ(P))2}CϵM2[1n+min{1,d2n2log2(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,d2n2log2(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.

Honest Expected Length for Transported Complier Effects with Weak First Stages

stat_transported_late_strength_frontier_v1 · Stat · 2026-08-02 ✓ verified in Lean 4 AI reviewer score 6.8/10

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,t01/2}\min\{1,t_0^{-1/2}\}, including finite-cell transport models under the stated growth and regularity conditions.

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,t01/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

stat_bdd_uniform_log_penalty_v1 · Stat · 2026-08-13 ✓ verified in Lean 4 AI reviewer score 6.8/10

The paper establishes that distance-based boundary regression discontinuity designs have minimax expected boundary sup-loss on the (logn/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.

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 q1q\geq1 and L4L\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=(lognn)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 LL0(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

stat_sa_plm_cumulant_converse_nongaussian_spectral_annihilation · Stat · 2026-08-21 ✓ verified in Lean 4 AI reviewer score 6.8/10

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 n1n^{-1}.

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

exp_dense_group_partition_projection_phase_v1 · Experimentation · 2026-09-10 ✓ verified in Lean 4 AI reviewer score 6.8/10

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.

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/n0N_n/n\to0, including regimes with Nn2/nN_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 M2M\ge2, a strictly positive limiting grouped-unit fraction ρ>0\rho>0, an outcome envelope B1B\ge1, and a variance floor cσc_\sigma below the separated independent-arm limit 1/[Mp(1p)]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.

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

eid_crl_coverratio_mmd_genericity_v1 · Exact ID · 2026-09-11 ✓ verified in Lean 4 AI reviewer score 6.6/10

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.

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 \citet{vonKugelgenEtAl2023UnknownInterventions}.

Quotient-law Inference with Latent Treatment-effect Collisions

stat_proxy_effectlaw_eigencollision_frontier_v1 · Stat · 2026-09-09 ✓ verified in Lean 4 AI reviewer score 6.4/10

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.

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,(ng)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

stat_pomdp_latent_overlap_minimax_v1 · Stat · 2026-09-13 ✓ verified in Lean 4 AI reviewer score 6.4/10

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.

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 C1C\ge1, to T1+TβqC2(1β),qC=C1C,β=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 qCT1/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

pid_slate_benefit_partialtransport_v1 · Partial ID · 2026-09-10 ✓ verified in Lean 4 AI reviewer score 6.3/10

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.

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(XK)O(|\mathcal X|K) operations for finite covariate support X\mathcal X and KK ordered outcome levels, with O(XK2)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 ηn0\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

exp_bipartite_minimax_design_v1 · Experimentation · 2026-07-16 ✓ verified in Lean 4 AI reviewer score 6.2/10

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.

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

stat_lmtp_threshold_atom_frontier_v1 · Stat · 2026-09-07 ✓ verified in Lean 4 AI reviewer score 6.2/10

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 n1/2+δnκ+1hnβn^{-1/2}+\delta_n^{\kappa+1}h_n^\beta, with matching lower bounds and regime separation.

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 caκπ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=n1/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 n1/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

exp_multiarm_secondorder_minimax_frontier_v1 · Experimentation · 2026-09-10 ✓ verified in Lean 4 AI reviewer score 6.1/10

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 n4/3n^{-4/3} scale, with matching order certified by an internal two-arm information argument.

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 qacaq_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 n4/3n^{-4/3} for all sufficiently large nn, while an embedded two-arm Bayesian information argument gives a matching O(n4/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

stat_policy_regret_margin_overlap_v1 · Stat · 2026-06-13 ✓ verified in Lean 4 AI reviewer score 6/10

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.

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(α,γ)cnr(α,γ),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

exp_rollout_chebyshev_minimax_tv_envelope_rollout_design · Experimentation · 2026-07-02 ✓ verified in Lean 4 AI reviewer score 5.8/10

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.

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 qqmax<1q\le q_{\max}<1 and the number of rollout intervals satisfies kcβk\ge c\beta for some c>1c>1. The minimax amplification is of order ((1+1q)2q)2β, \left(\frac{(1+\sqrt{1-q})^2}{q}\right)^{2\beta}, up to constants depending only on qmaxq_{\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

eid_lingam_direction_min_order_v1_truncated_cumulant_minimality · Exact ID · 2026-07-12 ✓ verified in Lean 4 AI reviewer score 5.8/10

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 m3m\ge3.

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 m3m\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

panel_ppml_forbidden_comparison_v1 · Panel · 2026-07-16 ✓ verified in Lean 4 AI reviewer score 5.5/10

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.

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

exp_design_pm1_exactness_boundary_v1_experimentation · Experimentation · 2026-07-02 ✓ verified in Lean 4 AI reviewer score 4.2/10

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.

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

stat_dose_response_minimax_holder_anisotropic_converse · Stat · 2026-06-29 ✓ verified in Lean 4 AI reviewer score 4/10

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 n2α/(2α+1)n^{-2\alpha/(2\alpha+1)}, independently of treatment-density smoothness β\beta.

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}cn2α/(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=n2α/(2α+1)n2/(1+d/(4s)+1/α). \rho_n= n^{-2\alpha/(2\alpha+1)} \vee n^{-2/(1+d/(4s)+1/\alpha)}. When d4sd\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.