CausalSmith · seminar slides
Uniform Accuracy Along a Boundary
Distance hides direction; recovering an entire boundary curve requires resolving many hidden local changes at once.
slides for Uniform Expected Risk for Distance-Based Boundary Regression Designs
Research Question
A geographic regression discontinuity design assigns treatment by location. Practitioners want accurate treatment effects along the whole boundary.
Keele and Titiunik (2015) explain the geographic design. Cattaneo et al. (2026) develop identification and local polynomial estimation using signed distance.
The usual procedure fits outcomes against signed distance separately on each side and subtracts the fitted intercepts.
| Experiment | Information at each query | Target | |---|---|---| | Unsigned | Outcomes and Euclidean distances | Boundary regression mean | | Signed | Outcomes, signed distances, query, and known geometry | Boundary treatment effect |
Cattaneo et al. (2026) establish an unsigned common-rule lower scale of n−1/4, and signed-distance bias and stochastic control in probability.
Pointwise accuracy does not control the largest boundary error. Probability bounds alone do not control expected loss on unusual samples.
What does simultaneous recovery cost, and how can an estimator control expected loss?
Key Idea
Pair a smooth change in the regression mean with a change in where observations lie around the query. Distance can hide the pair.
an=(nlogn)1/4.
Here n is sample size and an is the risk normalization.
- Unsigned: expected largest boundary error is asymptotically at least a positive multiple of an, including rules chosen separately at each query.
- Signed: the same lower scale holds with treatment sides and fixed known rectangular geometry.
- Signed attainment: a stabilized local polynomial estimator achieves error at most a constant times an, under distance identification, uniform first-order bias, and joint expected sampling bounds.
Hiding direction leaves fourth-power information about each local change. Resolving many changes adds the logarithm.
Example
Fix support X=[−3,3]2, treated region A1=[−1,1]×[0,2], control region A0=X∖A1, and treatment boundary B=∂A1.
Place separated disks of radius w along the middle of the bottom edge. Their centers are xj.
Each disk carries a binary choice ωj: enable a smooth bump of height Δ, or leave it off. The vector ω specifies the observation law Pω.
τPω(xj)=16xj(1)+ωjΔ.
This is the treatment effect at center j; xj(1) is its horizontal coordinate.
The potential outcomes are binary. The control mean is 1/2; the treatment mean adds a horizontal slope and the enabled bumps.
Positive signed distance identifies the treated semicircle. It still hides whether an observation lies to the left or right within that semicircle.
That remaining ambiguity lets us hide each bump.
Unsigned Model
Let Y be an outcome, X a two-dimensional location, and P their unknown law.
The target is μP(x)=EP[Y∣X=x] on the support boundary bd(XP).
At query x, the rule receives
Vn,P(x)=((Yi,∥Xi−x∥2):1≤i≤n).
These are n independent observations expressed as outcomes and distances; direction is discarded.
For integer q≥1 and envelope L≥4, the class PNP(L,q) requires:
- Compact support inside [−L,L]2, with a Lipschitz boundary.
- Continuous location density between L−1 and L, keeping local observations available.
- Regression derivatives bounded through order q−1, with the highest derivative Lipschitz, preventing arbitrarily narrow changes.
- Continuous conditional variance between L−1 and L, keeping noise bounded and nondegenerate.
Rules may differ across queries, but are chosen independently of the unknown law.
Signed Model
Let Y(0) and Y(1) be potential outcomes and T the treatment indicator.
T=1{X∈A1,P},Y=TY(1)+(1−T)Y(0).
Treatment follows the known assignment region; Y is the observed outcome.
The target is τP(x)=μ1,P(x)−μ0,P(x), where μt,P(x)=EP[Y(t)∣X=x] is the mean for arm t.
At a query on treatment boundary BP, use
DP,x±={1(X∈A1,P)−1(X∈A0,P)}∥X−x∥2
Positive distance denotes treatment; negative distance denotes control. Rules receive outcomes and signed distances, together with the query and known geometry.
The independent-sample class P12(p,ν,L) requires square support, bounded positive continuous density, and a compact interior boundary curve of controlled length.
Arm means have smooth extensions with bounded Lipschitz derivatives through integer order p≥0. Variances lie between L−1 and L; conditional 2+ν moments are bounded by L, with ν≥2.
Here L≥4 is the common envelope. Local geometry must also support fitting on both arms.
Local Geometry
A distance slice is the part of a circle at exactly radius s from query x that lies in arm t:
{z∈At,P:∥z−x∥2=s}
Here z is a location. For every boundary query, both arms, and 0<s≤L−1, require 0<∫slicefP(z)dH1(z)<∞: finite positive density-weighted arc length. The measure H1 measures length along the circle.
A neighborhood condition instead accumulates area over all radii up to h. With K□(u)=1{∣u∣≤1},
h21∫At,PK□({(2t−1)∥z−x∥2}/h)dz≥L−1.
For every query, both arms, and 0<h≤L−1, each arm occupies a uniformly positive fraction of the local disk; dz measures planar area.
Finally, the population matrix Ψt,P,x(h) of local polynomial regressor second moments, normalized by h−2, satisfies λmin(Ψt,P,x(h))≥L−1. This prevents unstable inversion.
In the rectangle, circular arcs give positive slice mass. Both arms retain quarter-disk sectors even at corners, providing neighborhood mass and stable designs.
Identification and Approximation
The signed upper bound needs two population assumptions.
Distance identification connects signed-distance conditional means to the treatment effect:
τP(x)=θ1(0+)−θ0(0−).
Here θ1 and θ0 are conditional outcome means on positive and negative signed distances. Their one-sided limits identify the two arm means at query x. We take this identification input from Cattaneo et al. (2026).
Uniform first-order bias controls the difference between the population degree-p fitted intercepts and τP(x).
n→∞limsupBiasRatiop,ν′,L(hn)≤Cb.
The ratio divides the largest population intercept-difference error by hn, taking the supremum over all boundary queries and all laws in P12(p,ν′,L).
- The same Cb, depending only on p,L, covers every ν′≥2.
- The index ν′ varies the moment class: its laws have bounded conditional 2+ν′ moments. The theorem’s ν remains fixed.
- The bound must hold along every positive decreasing deterministic sequence hn→0 with nhn2→∞.
Thus approximation error is uniformly O(hn), across these classes and bandwidth sequences.
Expected Sampling Bounds
The third additional assumption is joint expected sampling bounds.
Let Ψ be the sample counterpart of the population design matrix Ψ.
- GramDev is the largest entrywise error ∣Ψ−Ψ∣, over arms and boundary queries.
- ScoreDev is the largest error in local polynomial-weighted raw-outcome averages, normalized by hn−2.
For some constant Cm>0, require
EP∗[GramDevn,p,P(hn)]≤Cmnhn2log(hn−1),
EP∗[ScoreDevn,p,P(hn)]≤Cm{nhn2log(hn−1)+n(1+ν)/(2+ν)hn2log(hn−1)}.
These bound expected fluctuations over the whole boundary. The second term accommodates outcome tails; EP∗ is outer expectation, the expected-loss convention for boundary suprema.
Both bounds hold eventually, uniformly over P∈P12(p,ν,L), along every positive deterministic hn→0 satisfying log(hn−1)n(1+ν)/(2+ν)hn2→∞.
All three additional inputs are assumptions.
Unsigned Main Result
Let TnPI contain unsigned-distance rules chosen separately at each query, independently of the unknown law.
RnPI=μnPI∈TnPIinfP∈PNP(L,q)supEP∗[x∈bd(XP)supμn,PPI(x)−μP(x)].
The infimum chooses the best rule, the outer supremum chooses the worst law, and the inner supremum selects the largest boundary error.
informal · Theorem T-2 For integer q≥1 and L≥4, unsigned-distance rules chosen separately at each query have minimax expected largest boundary error asymptotically at least a positive multiple of an.
n→∞liminf(lognn)1/4RnPI=n→∞liminfanRnPI≥cPI.
Here cPI>0 depends on q,L.
One common rule at every query inherits this bound. Choosing rules query by query cannot remove the information barrier.
Signed Main Result
Let Tn12,± contain rules using outcomes, signed distances, the query, and known geometry, chosen independently of the unknown law.
Rn12,±(p,ν,L)=τn∈Tn12,±infP∈P12(p,ν,L)supEP∗[x∈BPsup∣τn,P(x)−τP(x)∣].
This is the best worst-law expected largest treatment-effect error.
informal · Theorem T-3 For integer p≥0, ν≥2, and L≥L0(p)≥48, signed-distance minimax expected largest error is asymptotically at least a positive multiple of an, already with fixed rectangular geometry.
Under distance identification, uniform first-order bias, and joint expected sampling bounds,
c≤n→∞liminfanRn12,±(p,ν,L)≤n→∞limsupanRn12,±(p,ν,L)≤C.
The constants satisfy 0<c≤C. The lower bound is unconditional within the law class; the upper bound uses all three additional inputs.
In the rectangle, the loss includes every edge and corner, although the hard changes lie along one edge.
Angular Cancellation
Return to disk j in the rectangle. Let r be distance from its center and ηj(x) the horizontal direction cosine: negative to the left, positive to the right.
Let φr(r/w) be a smooth radial bump, equal to one at the center and zero beyond radius w.
With fixed location density, enabling the bump adds Δφr(r/w) to the pooled treatment mean throughout the disk.
Instead, change density to (1/36){1+tΔ(r)ηj(x)}. A smooth cutoff turns off the tilt near the center; for r≥256Δ, use tΔ(r)=−32Δφr(r/w)/r.
- Left–right reflection preserves radius and treatment side. The average of ηj is zero, preserving signed-distance probabilities.
- The horizontal baseline varies by rηj/16. Multiplying it by the tilt contributes −2Δφr(r/w)ηj2.
- The semicircle average of ηj2 is one-half, exactly canceling the bump.
For binary outcomes, the success probability determines the whole outcome law.
Let Qj,b be the one-observation outcome–signed-distance law conditional on disk j, under choice b. With E={(y,u):0<u<2CΔ} and C=128,
Qj,0(⋅∩Ec)=Qj,1(⋅∩Ec),
Here y is outcome and u signed distance. Only the small central positive-distance window distinguishes the choices.
Fourth-power Information
Set q=p+1. Smoothness permits the packing
M≥cΔ−1/q,w=AΔ1/q,ρ=36πw2.
Here M counts separated disks, w is their radius, ρ is each disk’s probability, and c,A>0 are fixed construction constants.
Width w keeps the bump derivatives within the smoothness envelope.
Inside the uncancelled window, the binary success probabilities differ by at most Δ and remain between 1/4 and 3/4.
Kullback–Leibler divergence, written KL, measures how distinguishable two laws are. The construction gives
KL(Qj,0,Qj,1)≤C0w2Δ4,KL(Qj,1,Qj,0)≤C0w2Δ4.
Here C0>0 is fixed.
- Squared success-probability separation contributes two powers of Δ.
- The remaining window has area proportional to Δ2, contributing two more.
- Conditioning on the disk divides by its area, proportional to w2.
The hidden change therefore carries fourth-power information.
Unsigned Construction
Use the same bump-and-density pairing on half-disks along the bottom edge of support S=[−1/2,1/2]2.
Take width wn=an1/q, signal δn=c1an, and Mn separated centers xj, with fixed c1>0 sufficiently small.
μω(x)=21+8x(1)+j=1∑Mnωjδnφ(wnx−xj),
This regression mean has a horizontal baseline and enabled radial bumps; φ is the smooth bump used above.
Given X=x, use the outcome law
μω(x)N(1,1)+{1−μω(x)}N(0,1).
The notation N(m,1) means a normal law with mean m and variance one.
The entire mixture law depends linearly on its weight μω(x). The same angular cancellation therefore preserves the full outcome–distance law outside the central window.
One-observation divergence is at most CKLδn4, with fixed CKL>0. For n observations, the construction ensures nCKLδn4≤αlogMn, with 0<α<1/8.
How do these local ambiguities force an error somewhere when every query reuses one sample?
Independent Local Blocks
Distance data at different queries share observations. Independent pointwise tests do not establish a uniform lower bound.
For the argument, draw a Poisson number of observations with mean 2n. Splitting observations among disjoint cells produces independent blocks conditional on the binary choices.
Let Zj collect raw observations in cell j, and Z0 collect the background. Each cell law depends only on its own choice; the background law is common.
Decoder j, which tries to recover choice j, receives:
- Its own block compressed to outcomes and distances, denoted Sj=cj(Zj).
- Every other block and the background with raw locations.
For signed rules, the compression retains signed distance. Other raw blocks let decoder j reconstruct their distances from its query.
This enlarged view reproduces every admissible query rule while keeping its own angular information hidden.
The decoders have different raw-data views. We must compare their joint success, not multiply their individual success probabilities.
Coupled Alternatives
For this argument, Qj,b denotes the compressed sample-block law under choice b.
Its overlap is ρj=1−∥Qj,0−Qj,1∥TV, where total variation is the largest probability difference across events.
Construct a pair of raw blocks (Zj0,Zj1) with the correct laws for choices zero and one, whose compressions agree with probability at least ρj.
First couple the compressed laws on their common probability mass. Then draw raw locations conditional on each compressed block. This preserves both raw-block laws.
Draw these pairs independently across cells, and draw one common background Z0.
For any choice vector ω, select
Zj(ω)=Zjωj,1≤j≤M.
This selects the appropriate raw alternative in each cell. Every fixed ω has exactly its original product experiment.
We can therefore compare all choice vectors using one realized array of paired blocks.
Two-cell Pairing
Fix a coupled array with two cells. Suppose cell 1 overlaps: c1(Z10)=c1(Z11)=S1.
For either fixed choice b∈{0,1} in cell 2:
| Choice vector | Decoder 1 receives | Decoder 2 receives | |---|---|---| | (0,b) | (S1,Z2b,Z0) | (c2(Z2b),Z10,Z0) | | (1,b) | (S1,Z2b,Z0) | (c2(Z2b),Z11,Z0) |
Decoder 2 may see different raw information. Decoder 1 sees exactly the same input, while its true choice changes.
Its guess cannot be correct in both alternatives. Consequently, both decoders can be correct in at most one alternative of each pair.
Pair (0,0) with (1,0), and (0,1) with (1,1): at most half the four choice vectors have simultaneous success.
For any number of cells, the same argument applies: after fixing the coupled array, choose one overlapping coordinate and pair every vector with the vector that flips that coordinate.
Whenever an overlap occurs, at least half the choice vectors incur an error somewhere.
Uniform Decoding Failure
The overlap events are independent because the coupled cell pairs were drawn independently.
The chance that every cell avoids overlap is at most ∏j=1M(1−ρj).
On that event, simultaneous success can be as high as one. Otherwise, the pairing argument bounds success by one-half. Averaging over coupled arrays and uniform binary choices gives
Pr{Ωj=Ωj for every j}≤21{1+j=1∏M(1−ρj)}.
Here Ωj is the true choice and Ωj its decoded value.
If each block has divergence at most κlogM, with 0≤κ<1, the standard divergence-to-overlap bound gives ρj≥21exp{−κlogM}.
Many chances to overlap then imply
Pr{at least one decoding error}≥21[1−exp{−M1−κ/2}].
This lower bound stays away from zero as M grows. It allows all the additional raw information in the decoders’ different views.
Rate Calibration
A signed Poisson cell block has expected size 2nρ.
Its conditional one-observation divergence is at most a constant times Δ4/w2. Since ρ=πw2/36, multiplying by expected block size cancels w2.
Block information is therefore at most a constant times nΔ4.
The packing has M≍Δ−1/q. Calibrate the signal by
nΔ4≍logM,henceΔ≍(nlogn)1/4.
These relations choose a signal comparable to an; a sufficiently small constant keeps block divergence within the simultaneous-testing budget.
- An estimate accurate within half the center separation would decode every choice.
- An error somewhere therefore forces expected maximum estimation error at least a positive multiple of an.
- Keeping the first n Poisson observations recovers the original experiment; the probability of a sample shortage is exponentially small.
Smoothness determines how many bumps fit. Angular cancellation determines the fourth root. Simultaneous recovery determines the logarithm.
Local Polynomial Estimator
For signed attainment, fit degree-p polynomials separately on the two arms within distance h.
Use rp(u)=(1,u,…,up)⊤, with u equal to signed distance divided by h. The arm intervals are I1=[0,∞) and I0=(−∞,0).
Winsorization replaces outcomes outside [−B,B] by the nearest endpoint; denote it by ψB and take B(h)=h−1/3.
mt,xB(h)=nh21i=1∑n1{DP,x,i±∈It}K□(DP,x,i±/h)rp(DP,x,i±/h)ψB(h)(Yi).
This vector averages polynomial-weighted winsorized outcomes. The factor h−2 normalizes a two-dimensional neighborhood.
The sample matrix Ψ uses the same weights with regressor outer products. Set βB=Ψ−1mB when its smallest eigenvalue is at least (2L)−1; otherwise use zero.
With e1=(1,0,…,0)⊤ selecting the intercept, estimate
τn,h12,LP(x)=max{−2L,min{2L,e1⊤(β1,xB(h)−β0,xB(h))}}.
This subtracts the fitted intercepts and clips the effect estimate to its bounded target range.
Residual Control
Let βt,P,x∗(h) be the population coefficient fitted using raw outcomes. Its polynomial-weighted raw residual has mean zero.
Write scoret,x,hB for the sample weighted residual vector: winsorized outcomes minus predictions from β∗.
On the guarded branch, coefficient error is the inverse sample design matrix times this residual vector. The eigenvalue guard bounds amplification.
The moment envelope bounds the absolute conditional winsorization remainder by LB−3, for B≥1 and ν≥2.
For centered fluctuations, we establish
EPn∗[t∈{0,1},x∈BPsupscoret,x,hB−EPn[scoret,x,hB]]≤C{nh2log(B/h)+nh2Blog(B/h)}.
Here Pn is the independent sample law and C=C(p,L). This bounds expected residual fluctuations across both arms and the entire boundary.
We apply van der Vaart and Wellner (1996) for maxima of bounded sample averages. Winsorization bounds individual contributions; local windows have probability proportional to h2.
This proved bound supplies the winsorized residual component. The theorem assumes the joint sampling package; the attainment argument uses its Gram bound.
Expected Risk and Tuning
Expected largest error has four contributions:
| Contribution | Control | |---|---| | Population approximation, O(h) | Uniform first-order bias | | Mean tail remainder, O(B−3) | Conditional outcome moments | | Centered residual fluctuations | Proved expected maximal bound | | Unstable sample designs | Expected Gram bound and clipping |
Small Gram deviation activates every inversion guard. On other samples, clipping bounds the largest loss by 4L. Expected Gram control therefore also pays for unusual samples.
Choose hn=an and Bn=an−1/3. The bandwidth balances first-order bias with {log(hn−1)/(nhn2)}1/2, while Bn−3=an.
Every contribution is at most a constant times an:
P∈P12(p,ν,L)supEP∗[x∈BPsupτn,hn12,LP(x)−τP(x)]≤Can.
This bounds the estimator’s expected largest treatment-effect error uniformly over the law class, eventually in n.
informal · Theorem T-4 For integer p≥0, ν≥2, and L≥4, distance identification, uniform first-order bias, and joint expected sampling bounds give the stabilized estimator expected largest boundary error at most Can eventually.
In the rectangle, this covers every edge and corner. With the lower bound, it gives the matched signed rate for L≥L0(p).
Open Questions
- Establish a matching unsigned upper bound.
- Derive all three signed-distance analytic inputs within the same law class.
- Extend the analysis to other metrics, kernels, boundary regularity, and losses suited to inference.
The common lower scale alone does not rank achievable unsigned and signed risks.
Takeaways
- Hidden direction: angular density reweighting cancels a smooth mean bump in the full distance-observation law outside a small central window.
- Uniform difficulty: fourth-power local information and many separated choices produce an=(logn/n)1/4; coupling and pairing force an error somewhere despite the decoders’ different raw-data views.
- Expected-risk attainment: under distance identification, uniform first-order bias, and joint expected sampling bounds, winsorization, guarded inversion, and clipping attain signed error at most a constant times an.
Appendix: Common-map Lower Bound
One common unsigned-distance rule incurs the logarithmic lower scale over the full nonparametric class.
For every integer q≥1 and every L≥4, there exists a constant cCTY=cCTY(q,L)>0 such that, for the common-map minimax risk RnNP in Definition P-5 and an=(logn/n)1/4, n→∞liminf(lognn)1/4RnNP=n→∞liminfanRnNP≥cCTY.
Appendix: Point-indexed Unsigned Lower Bound
Allowing a separate unsigned-distance rule at each query preserves the logarithmic lower scale.
For every q∈N and L∈R satisfying q≥1 and L≥4, there exists a constant cPI=cPI(q,L)>0 such that, with RnPI the point-indexed outer-expectation minimax risk in Definition P-7 and an=(logn/n)1/4, n→∞liminf(lognn)1/4RnPI=n→∞liminfanRnPI≥cPI.
Appendix: Signed-distance Lower Bound
The signed-distance logarithmic lower bound already holds with one fixed known rectangular geometry.
For every nonnegative integer p, there exists a real number L0=L0(p) such that:
- (Envelope threshold.) L0≥48.
- (Moment exponent.) ν≥2.
- (Uniform envelope.) L≥L0.
- (Risk definition.) Rn12,±(p,ν,L) is the signed-distance minimax risk in Definition P-13.
- (Fixed-geometry risk.) Rn12,±,fix(p,ν,L) is the infimum over the same known-geometry point-indexed decision rules of the same worst-case outer-expected interface sup-loss, with the law supremum further restricted to laws in P12(p,ν,L) satisfying the fixed hard-geometry condition.
For every ν and L satisfying these conditions, there exists a constant c=c(p,ν,L)>0 such that, writing an=(logn/n)1/4, n→∞liminfan−1Rn12,±(p,ν,L)=n→∞liminfanRn12,±(p,ν,L), and c≤n→∞liminfanRn12,±(p,ν,L),c≤n→∞liminfanRn12,±,fix(p,ν,L).
Appendix: Conditional Matched Rate
Under the three analytic inputs, signed-distance minimax expected risk is bounded above and below on the an scale.
For every integer p≥0, there exists a constant L0≥48 such that, for every ν≥2 and every L≥L0, the following conditions imply the frontier conclusion below:
- (Distance identification.) For every A1/A2 law P, every map d satisfying the CTY identification assumptions for p and P, and every x∈BP, there are selected source-coherent signed-distance conditional-mean versions θ1,θ0:R→R whose right and left limits at zero exist and identify τP(x)=u↓0limθ1(u)−u↑0limθ0(u).
- (Uniform first-order bias.) There is a constant depending only on p and L such that, uniformly over every ν′≥2 and every positive antitone deterministic bandwidth sequence hn→0 with nhn2→∞, the finite limsup of the normalized uniform first-order bias ratio is bounded by that constant.
- (Expected maximal bounds.) For this p,ν,L, there is a constant Cmax>0 such that, along every positive bandwidth sequence hn→0 with log(hn−1)n(1+ν)/(2+ν)hn2→∞, for all sufficiently large n and every P∈P12(p,ν,L), the expected outer Gram-deviation supremum satisfies EP∗[t∈{0,1}supx∈BPsupj,ksupΨt,x(hn)jk−Ψt,P,x(hn)jk]≤Cmaxnhn2log(hn−1), and the expected outer centered raw-score supremum satisfies EP∗[t∈{0,1}supx∈BPsup∥mt,x(hn)−mt,P,x(hn)∥]≤Cmax(nhn2log(hn−1)+n(1+ν)/(2+ν)hn2log(hn−1)).
Then there exist constants c and C with 0<c≤C such that, for the signed-distance minimax risk Rn12,±(p,ν,L) of Definition P-13 and an=(logn/n)1/4, c≤n→∞liminfanRn12,±(p,ν,L)≤n→∞limsupanRn12,±(p,ν,L)≤C. Moreover, the explicit winsorized, Gram-stabilized degree-p signed-distance local-polynomial estimator with bandwidth hn=an and winsorization level Bn=an−1/3, clipped to [−2L,2L], has uniform outer risk over P∈P12(p,ν,L) satisfying n→∞limsupansupP∈P12(p,ν,L)EP∗[supx∈BPτn,an12,LP(x)−τP(x)]≤C.