CausalSmith · seminar slides
Minimax Inference for Threshold Clamp Policies
For continuous treatments, we characterize the minimax risk and honest confidence-interval length for a lower-threshold clamp when overlap thins polynomially near the boundary.
slides for Minimax Inference for Threshold Modified Treatment Policies with Continuous Treatments
Overview
- The policy raises all treatments below the threshold δn, the policy threshold, up to δn.
- This creates a moving atom at the boundary of the post-policy treatment law.
- We estimate the clamp mean by combining a retained-course mean, a lower-tail atom mass, and a boundary regression value.
- The frontier has two pieces: ordinary sampling noise and atom-weighted boundary learning.
- The same frontier governs point estimation, honest interval length, and the causal lift.
informal · Theorem T-2 Under the finite-stratum Hölder model with the two-sided polynomial density envelope, the minimax absolute-error risk is at most and at least constant multiples of rn.
Motivation
- Modified treatment policies are a common way to define feasible effects for continuous treatments.
- A threshold rule is natural when doses below a minimum practical level are reassigned to that minimum.
- In applications, the threshold may move with sample size as investigators probe closer to the support boundary.
- The inferential difficulty is concentrated at the boundary value created by the clamp.
- The question is: how much precision is possible when the data become sparse near that boundary?
Related Literature
- Díaz et al. (2021), Williams and Díaz (2023), and Hoffman et al. (2024) make modified treatment policies operational for continuous treatments.
- Kennedy et al. (2017) and Bonvini and Kennedy (2026) study continuous-treatment dose-response inference under regular support.
- van der Laan et al. (2022) study stochastic threshold interventions with efficient inference and bands.
- Gaïffas (2005) gives the closest degenerate-design pointwise regression benchmark.
- Low (1997) and Armstrong and Kolesár (2018, 2016) supply the honest interval and modulus perspective.
- Our contribution is the exact deterministic clamp frontier under polynomial thinning and moving thresholds.
Setup
- Observed data are O=(X,A,Y): finite stratum, continuous treatment, bounded outcome.
- The treatment density in stratum x is πx(a).
- The lower-threshold clamp is dδ(a)=max{a,δ}.
- The observed clamp target is the retained mean above the threshold plus the lower-tail atom mass times the boundary regression.
- The first-order object is the atom at δn: its mass shrinks like δnκ+1, where κ is the thinning exponent.
For P∈M as in Definition P-1, x∈X, and 0≤δ≤δˉ, define qx(δ)=∫[0,δ]πx(a)da,νδ(P)=EP[Y1{A>δ}]. The unsmoothed observed-data clamp target is θδ(P)=νδ(P)+x∈X∑pxqx(δ)μxP(δ).
Assumptions
- We work with fixed finite strata and i.i.d. sampling.
- Every stratum has nonvanishing mass.
- The treatment density obeys a two-sided envelope c−aκ≤πx(a)≤c+aκ on all of [0,1]; it is at the lower endpoint that this thins the design.
- The outcome regression is Hölder smooth with exponent β, the smoothness exponent, and radius L.
- These conditions describe the amount of information near the moving threshold.
For every x∈X and for Lebesgue-almost every a∈[0,1], c−aκ≤πx(a)≤c+aκ.
The regression function μxP satisfies the following condition with parameters β and L. For every x∈X:
- (Continuity.) The map a↦μxP(a) is continuous on [0,1].
- (Range.) For every a∈[0,1], μxP(a)∈[0,1].
- (Regression version.) The function μXP(A) realizes the conditional outcome regression given the complete design: EP[Y∣X,A]=μXP(A)P-almost surely.
- (Hölder Taylor remainder.) Let ℓ=ℓ(β) be the local-polynomial order defined in Definition. For all s,t∈[0,1], with derivatives taken intrinsically on [0,1], μxP(t)−j=0∑ℓj!(μxP)(j)(s)(t−s)j≤L∣t−s∣β.
Information Balance
- The threshold regression is learned locally around δn.
- The bandwidth hn, the local window width, balances local-polynomial bias against local sampling information.
- The resulting Hölder frontier rate combines sampling noise and atom-weighted boundary error.
The information-balance bandwidth hn for the threshold sequence δn is hn=inf{h∈(0,1−δˉ]:nh2β+1(δn+h)κ≥1}, with hn=1−δˉ when the displayed set is empty.
The Hölder minimax rate rn is rn=n−1/2+δnκ+1hnβ.
Estimator
- Split the sample into three fixed blocks.
- Use one block for the retained-course mean above δn.
- Use one block for the empirical lower-tail atom masses.
- Use one block for the local-polynomial boundary regression.
- Stabilize by checking the realized total Gram matrix and falling back to a bounded value on singular local designs.
Main Result
informal · Theorem T-2 The total-Gram estimator attains worst-case absolute error at most Crn, and every estimator has worst-case absolute error at least crn, under the stated Hölder clamp model.
Fix J,β,κ,L,c−,c+,pmin,δˉ,α satisfying the regime conditions:
- (Regime.) J≥1, β>0, κ≥0, L>0, c−>0, c−≤κ+1≤c+, pmin>0, pmin≤1/J, δˉ∈(0,1), and α∈(0,1/2).
Then there exist constants 0<c<C and an amplitude a∈(0,1/4], depending only on the displayed regime constants, such that, for every deterministic threshold sequence (δn)n≥1 with δn∈[0,δˉ] for every n and every deterministic sequence Bn=(I0,n,I1,n,I2,n) of three pairwise-disjoint sample blocks with ∣Ij,n∣≥⌊n/4⌋ for j=0,1,2, for all sufficiently large n, with hn the information-balance bandwidth from Definition P-3 and rn=n−1/2+δnκ+1hnβ as in Definition P-4, the observed minimax absolute-error risk Rn⋆ over the clamp model M in Definition P-1 and the worst-case i.i.d. risk over M of the stabilized estimator using Bn satisfy crn≤Rn⋆≤P∈MsupEP[θn,BnTG−θδn(P)]≤Crn. Moreover, for all sufficiently large n, the lower bound is witnessed inside the same clamp model by:
- (Global shift.) Laws P0,P1∈M satisfying Assumption A-1, with the same (X,A)-design distribution and the same px and πx functions, having Bernoulli outcomes, and obeying μxP1(a)=μxP0(a)+n−1/2for every x and a∈[0,1], with well-posed product chi-squared divergence, P1⊗n≪P0⊗n,∫(dP0⊗ndP1⊗n−1)2dP0⊗n<∞, and cn−1/2≤∣θδn(P1)−θδn(P0)∣,χ2(P1⊗n,P0⊗n)≤C.
- (Localized perturbation.) Laws Q0,Q1∈M satisfying Assumption A-1, with the same (X,A)-design distribution and the same px and πx functions, having Bernoulli outcomes, and admitting a continuous bump b:R→R satisfying b(0)=1, b(u)≥0 for every u, b(u)=0 for u∈/[−1,1], and ∣b(u)∣≤1 for u∈[−1,1], with regression shift μxQ1(t)−μxQ0(t)=ahnβb(hnt−δn) for every x and t∈[0,1]. The product chi-squared divergence is well posed, Q1⊗n≪Q0⊗n,∫(dQ0⊗ndQ1⊗n−1)2dQ0⊗n<∞, and cδnκ+1hnβ≤∣θδn(Q1)−θδn(Q0)∣,χ2(Q1⊗n,Q0⊗n)≤C.
Honest Inference
informal · Theorem T-3 The bias-aware interval has uniform coverage at least 1−α, expected length at most Crn, and every uniformly honest interval has expected length at least crn.
There are constants c,C∈(0,∞), with c<C, for which the following holds.
- (Regime constants.) The number of strata J is a positive integer, β>0, κ≥0, L>0, c−>0, c−≤κ+1≤c+, pmin>0, pmin≤1/J, δˉ∈(0,1), and α∈(0,1/2).
- (Threshold path.) The deterministic thresholds satisfy 0≤δn≤δˉ for every n.
- (Sample splits.) For every n, Bn=(I0,n,I1,n,I2,n) is a deterministic three-way split of the sample indices with pairwise disjoint blocks and ∣Ij,n∣≥⌊n/4⌋ for j=0,1,2.
For all sufficiently large n, set hn to be the information-balance bandwidth for δn from Definition P-3, and set rn=n−1/2+δnκ+1hnβ. Then the stabilized interval CIn has worst-case coverage inf{P⊗n{θδn(P)∈CIn(O1,…,On)}:P∈M}≥1−α, and its worst-case expected length is comparable to rn: crn≤sup{EP⊗n[len(CIn)]:P∈M}≤Crn. Moreover, the minimax honest expected length satisfies crn≤Ln⋆, where Ln⋆ is the infimum, over observed-sample confidence procedures with uniform coverage at least 1−α over M, of their worst-case expected length. Equivalently, every observed-sample confidence procedure Cn with inf{P⊗n{θδn(P)∈Cn(O1,…,On)}:P∈M}≥1−α has worst-case expected length at least crn: crn≤sup{EP⊗n[len(Cn)]:P∈M}.
Phase Diagram
- The threshold path determines which term in rn is visible.
- Near zero, the atom is small enough for root-n behavior.
- Past the critical scale, the atom-weighted boundary regression term governs the rate.
- At a fixed positive threshold, the rate becomes the usual boundary nonparametric rate.
- At threshold zero, the target is the ordinary observed mean.
informal · Theorem T-4 The phase diagram separates regular, critical, vanishing atom-dominated, fixed-threshold, and zero-threshold regimes for rn.
Let J be a positive integer and let β>0,κ≥0,L>0,0<c−≤κ+1≤c+,0<pmin≤1/J,0<δˉ<1,0<α<1/2. Define the phase and edge scales by δcrit,n=n−1/{2(βκ+2β+κ+1)},δedge,n=n−1/(2β+κ+1). Then δedge,nδcrit,n→∞. Moreover, for every deterministic threshold sequence (δn)n≥1 with δn∈[0,δˉ] for every n, let hn be the information-balance bandwidth in Definition P-3, and set an=δnκ+1hnβ,rn=n−1/2+δnκ+1hnβ as in Definition P-4. The following conclusions hold:
- Regular thresholds. If δn/δcrit,n→0, then rn≍n−1/2.
- Critical thresholds. For every c0>0, if δn/δcrit,n→c0, then an≍n−1/2 and rn≍n−1/2.
- Vanishing atom-dominated thresholds. If δn/δcrit,n→∞ and δn→0, then rn≍δnκ+1(nδnκ)−β/(2β+1).
- Fixed thresholds. For every δ0∈(0,δˉ], if δn→δ0, then rn≍n−β/(2β+1).
- Zero threshold identities. For every n and every P∈M satisfying Definition P-1, qx(0)=0for every x,θ0(P)=∫YdP,rnδn=0=n−1/2, where the last identity uses the bandwidth hn at threshold 0.
- Zero threshold estimator. For every n, every P∈M satisfying Definition P-1, and every admissible three-way split B=(I0,I1,I2) as in Definition P-5, the total-Gram estimator in Definition P-5, formed at threshold 0 with bandwidth hn, agrees under the product sampling law generated by P with z↦clamp[0,1](∣I0∣1i∈I0∑Yi).
- Zero threshold stabilized guarantees. For every deterministic sequence B∙=(Bn)n≥1 of admissible three-way split blocks, the stabilized zero-threshold procedure has coverage at least 1−α for all sufficiently large n. Its stabilized worst-case risk at threshold 0 and its stabilized worst-case expected length at threshold 0 with noncoverage level α are both asymptotic to n−1/2.
Calibration
- The one-stratum example makes the phase boundary concrete.
- Take β=κ=1 and π(a)=2a.
- The localized Bernoulli alternatives keep likelihood distance bounded while moving the target by the atom-weighted boundary amount.
- The critical threshold scale is where this movement matches n−1/2.
informal · Theorem T-5 In the one-stratum linear-thinning case, hn≍(nδn)−1/3, Δn≍δn2hn, and δn≍n−1/10 is equivalent to Δn≍n−1/2.
Key Idea
- The clamp mean has two statistically different pieces.
- The retained mean behaves like a bounded sample average.
- The lower-tail atom multiplies the regression value at the threshold.
- A naive boundary plug-in inherits degenerate-design instability near sparse support.
- Total-Gram stabilization uses the realized local design only when it has enough curvature.
- The atom weight shrinks the boundary-regression error from hnβ to δnκ+1hnβ.
Proof Sketch
- Upper bounds decompose the estimator into retained-course error, atom-mass error, and boundary-regression error.
- Polynomial thinning gives the local sample size and Gram curvature scale in the threshold window.
- Hölder smoothness gives a deterministic local-polynomial bias bound.
- The lower bound uses two experiments: a global Bernoulli shift for n−1/2, and a localized threshold bump for δnκ+1hnβ.
- Low-style honest-length lower bounds transfer these testing separations to interval length.
Causal Interpretation
- The full-data class adds potential outcomes through a latent-response representation.
- Consistency links observed outcomes to structural responses at the realized treatment.
- Conditional exchangeability identifies the structural response mean within strata.
- Response continuity aligns the structural mean with the observed regression value on the threshold range.
informal · Theorem T-1 Under the declared full-data causal conditions, the causal clamp mean equals the observed clamp target.
informal · Theorem T-7 Every observed law in the Hölder clamp model has a full-data lift, so the observed and causal minimax criteria coincide.
Causal Frontier
informal · Theorem T-6 The same total-Gram estimator, bias-aware interval, minimax risk rate, and honest-length rate rn hold for the causal clamp mean over the full-data class.
Let J∈N, and let β,κ,L,c−,c+,pmin,δˉ,α∈R. Suppose that
- (Regime constants.) The constants satisfy 0<J,β>0,κ≥0,L>0,c−>0,c−≤κ+1≤c+, pmin>0,pmin≤J1,0<δˉ<1,0<α<21.
- (Threshold sequence.) The deterministic thresholds (δn)n∈N satisfy 0≤δn≤δˉ for every n∈N.
- (Split blocks.) For each n∈N, Bn=(I0,n,I1,n,I2,n) is a deterministic three-way split with pairwise disjoint blocks and ⌊4n⌋≤∣I0,n∣,⌊4n⌋≤∣I1,n∣,⌊4n⌋≤∣I2,n∣.
Let hn, the information-balance bandwidth for δn, be as in Definition P-3, and let rn=n−1/2+δnκ+1hnβ be the frontier rate from Definition P-4. Then there exist constants c,C∈R, with 0<c<C, depending only on these regime constants, such that, for every threshold sequence and every deterministic split-block sequence satisfying the preceding conditions, all sufficiently large n satisfy the following conclusions. The total-Gram estimator θn in Definition P-5, formed with split Bn, order ℓ=⌈β⌉−1, bandwidth hn, and threshold δn, is observed-sample measurable. The corresponding bias-aware interval CIn in Definition P-6, formed with the same split, order, bandwidth, and threshold, is an observed-sample measurable interval. For the causal minimax criteria Rn,F⋆ and Ln,F⋆ in Definition P-9, crn≤Rn,F⋆≤PF∈MFsupE(PF)⊗nθn−ψδn(PF)≤Crn, PF∈MFinf(PF)⊗n{ψδn(PF)∈CIn}≥1−α, PF∈MFsupE(PF)⊗nlen(CIn)≤Crn. Moreover, in the extended nonnegative-real order, after embedding the nonnegative real bounds into R≥0∪{∞}, crn≤Ln,F⋆≤Crn.
Continuity-only Comparison
- The continuity-only model keeps the same design and thinning conditions.
- It replaces quantitative Hölder smoothness with qualitative continuity of the regression extension.
- The estimator uses the retained mean, atom masses, and a fixed bounded fallback for the atom regression.
- The resulting frontier is driven by sampling noise plus the total lower-tail mass.
informal · Theorem T-8 In the continuity-only full-data class, the causal clamp mean equals the continuity-only observed clamp target.
informal · Theorem T-9 The continuity-only observed and causal minimax criteria coincide through observed-margin surjectivity.
informal · Theorem T-10 Over the continuity-only class, observed and causal minimax risk and honest-length rates are characterized by sn=n−1/2+δnκ+1, with the stated elbow and fixed-positive-threshold behavior.
Fix J∈N and constants κ,c−,c+,pmin,δˉ,α∈R. Suppose 0<J,0≤κ,0<c−,c−≤κ+1≤c+,0<pmin≤J1,0<δˉ<1,0<α<21. Then there exist constants c,C∈R with 0<c<C such that the following statements hold.
- (Uniform frontier bounds.) For every deterministic threshold sequence (δn)n≥1 with δn∈[0,δˉ] for all n, and for every deterministic sequence B∙=(Bn)n≥1 of admissible three-way split blocks, all sufficiently large n satisfy the following bounds. With sn=n−1/2+δnκ+1, and with Rn,cont⋆, Ln,cont⋆, Rn,cont,F⋆, and Ln,cont,F⋆ as in Definition P-16, csn≤Rn,cont⋆≤P∈McontsupEP⊗nθn,cont−θδn,cont(P)≤Csn, the interval CIn,cont has uniform coverage at least 1−α over Mcont, its worst-case expected length is at most Csn, and csn≤Ln,cont⋆. Moreover, Rn,cont,F⋆=Rn,cont⋆,Ln,cont,F⋆=Ln,cont⋆, and the same estimator and interval satisfy Rn,cont,F⋆≤PF∈McontFsupE(PobsF)⊗nθn,cont−ψδn(PF)≤Csn, with causal coverage at least 1−α uniformly over McontF and causal worst-case expected length at most Csn.
- (Elbow and null-threshold rates.) For the continuity-only rate in Definition P-12, sn(δn=n−1/(2(κ+1)))≍n−1/2,sn(0)≍n−1/2.
- (Fixed positive thresholds.) For every deterministic threshold sequence (δn)n≥1 with δn∈[0,δˉ] for all n, every δ0∈(0,δˉ], and δn→δ0, Rn,cont⋆→0,Ln,cont⋆→0.
Conclusion
- We characterize estimation and honest inference for the exact lower-threshold clamp under a two-sided polynomial density envelope on [0,1].
- The frontier is rn=n−1/2+δnκ+1hnβ.
- The phase diagram explains how moving thresholds shift the problem between root-n, atom-dominated, and fixed-threshold regimes.
- The causal bridge gives the same rates a full-data interpretation under the stated consistency, exchangeability, and continuity conditions.
- The continuity-only analysis separates the role of qualitative continuity from the quantitative Hölder modulus.