MSc-Abschlussarbeit — Universität zu Köln | Betreuer: Prof. Dr. Dominik Wied | 2026
MSc-Abschlussarbeit — Universität zu Köln | Betreuer: Prof. Dr. Dominik Wied | 2026
Klassische Endogenitätstests sind in großen Stichproben uninformativ — selbst vernachlässigbare Verzerrungen werden mit Wahrscheinlichkeit eins entdeckt und führen zu unnötiger Instrumentalvariablenschätzung, die Präzision kostet ohne inhaltlichen Mehrwert.
Ich entwickle einen Test, der nicht fragt „Ist Endogenität exakt null?", sondern „Verändert die Korrektur die geschätzte bedingte Verteilung um mehr als eine inhaltlich bedeutsame Schwelle Δ?" — eine statistisch wie substanziell fundierte Entscheidungsgrundlage.
Semiparametrische Verteilungsregression · Instrumentalvariablen (Control-Function) · Asymptotische Theorie (Einfluss-Funktionen, sequenzielle Linearisierung) · Monte-Carlo-Simulation in R · Empirische Anwendung (Card 1995, n = 3 010)
Der Test ist asymptotisch pivotal — kein Bootstrap erforderlich. Simulationen bestätigen Niveau α = 0,05 und Konsistenz. Empirisch: Endogenität ist praktisch relevant, aber nur für hochgebildete Arbeitnehmer — ein Befund, den ein klassischer Hausman-Test vollständig verdeckt.
Analyseverfahren, die statistische Signifikanz mit praktischer Bedeutsamkeit gleichsetzen, rauschen über das hinaus, was wirklich zählt. Diese Arbeit demonstriert, wie man eine Entscheidungsregel so formuliert, dass sie inhaltlich kalibriert ist — eine Fähigkeit, die überall dort relevant wird, wo Modellaussagen operative oder rechtliche Konsequenzen haben.
Beide Diagramme zeigen dieselbe gemessene Verzerrung d̂ = 0,026. Der Unterschied liegt allein in der Entscheidungsschwelle. (Monte-Carlo-Illustration: n = 1.000, ρ = 0,5.)
Frage: Ist d > 0?
Jede messbare Abweichung d̂ > 0 führt in großen Stichproben zur Ablehnung — auch wenn die Korrektur inhaltlich bedeutungslos ist.
Frage: Ist d > Δ?
d̂ = 0,026 liegt unterhalb der Relevanzschwelle Δ = 0,05: Verzerrung statistisch vorhanden, aber inhaltlich bedeutungslos. IV-Schätzung würde Präzision ohne Mehrwert kosten.
Das Folgende enthält die formale Theorie, Simulationsergebnisse und die empirische Anwendung der Abschlussarbeit. Die Arbeit ist noch in Bearbeitung; theoretisches Rahmenwerk, Asymptotik und Monte-Carlo-Simulation sind abgeschlossen.
We develop a test for the practical relevance of endogeneity in semiparametric distribution regression (DR). Classical tests ask whether the endogenous regressor is exactly exogenous — a question that is uninformative in large samples because even negligible deviations are detected with probability approaching one. We instead test whether correcting for endogeneity changes the estimated conditional distribution by more than an economically meaningful threshold Δ > 0.
The object of interest is the L²-distance d(x,y₂) = ‖F𝑁(·|x,y₂) − F𝑉(·|x,y₂)‖₂ between the ordinary DR estimator and the control-function IV-DR estimator of Wied (2024). The resulting test statistic is asymptotically pivotal: its limiting distribution under the boundary d = Δ is that of a ratio of Brownian-motion functionals, whose quantiles can be obtained by simulation once and for all. The test is consistent and detects local alternatives at rate n⁻½.
A central challenge in applied microeconometrics is deciding whether to use instrumental variables (IV) estimation. IV is consistent under endogeneity but less precise than OLS. The conventional approach relies on Hausman-type tests of exact exogeneity: if the null is rejected, switch to IV. This creates a well-known trap.
This perspective belongs to the equivalence-testing and minimum-detectable-effect literatures (Wellek 2010; Lakens 2017), but has received little attention in distributional models. This paper fills that gap for semiparametric distribution regression.
Distribution regression (DR), introduced by Foresi & Peracchi (1995) and developed by Chernozhukov, Fernández-Val & Melly (2013), models the entire conditional distribution function of an outcome Y given covariates X:
where Λ is a known link function (e.g. probit or logit) and β(y) varies with the threshold y. Each y-binary regression 1{Yₐ ≤ y} = Λ(Xₐ'β(y)) + Uₐ(y) is estimated separately, yielding a flexible nonparametric description of the full conditional distribution. Unlike quantile regression, DR handles mass points and non-smooth distributions naturally, and avoids bandwidth choices.
When the regressor Y₂ is endogenous — correlated with the latent error U(y) — the standard DR estimator is inconsistent for the structural conditional distribution. The paper adopts a triangular simultaneous-equations structure:
Endogeneity enters through the error decomposition U(y) = α1(y)V + α2(y)ε(y), where α1(y) ≠ 0 induces correlation between the structural error and the first-stage residual V. Using instruments Z and the control-function approach of Rivers & Vuong (1988), Wied (2024) proposes a three-step IV-DR estimator that recovers the structural distribution.
Regress Y₂ₐ on (Xₐ, Zₐ) by OLS to obtain estimated first-stage residuals V̂ₐ = Y₂ₐ − Aₐ'Γ̂.
For each threshold y, estimate the augmented binary choice model that includes V̂ₐ as an additional regressor, treating it as the control function.
Integrate out unobserved heterogeneity V by averaging over the empirical distribution of V̂ₐ, obtaining the average structural distribution F𝑉(y|x,y₂).
The practical endogeneity distance at covariate profile (x,y₂) is the L²-norm of the pointwise difference between the naive and IV-corrected conditional distribution functions over the outcome interval I:
Rejecting H₀ means that correcting for endogeneity changes the estimated conditional distribution by more than Δ over I. Failing to reject means endogeneity is practically irrelevant at that covariate profile, regardless of statistical significance.
The threshold Δ must be chosen before observing the data on substantive grounds. Two natural approaches:
The test statistic is the full-sample integrated squared distance:
The self-normalizing term V̂n(x,y₂) is constructed from sequential (partial-sample) versions of both estimators, adapting the self-normalization device of Dette, Möllenhoff & Wied (2025) to the one-sample correlated setting. The quantile q1−α is that of a Brownian-motion ratio W = B(1) / (∫ε1(B(t)/t − B(1))² dt)1/2 — simulatable to arbitrary precision with q0.95 ≈ 1.755 for ε = 0.1.
The theoretical contribution requires combining three distinct bodies of asymptotic theory, because the naive and IV-DR estimators are estimated from the same data — creating a non-trivial correlation structure absent in the independent two-sample case.
Derives the influence functions φNi(y) and φVi(y) for the naive and IV-DR estimators jointly, accounting for their dependence. The IV influence function requires four steps: first-stage OLS expansion, generated-residual Taylor expansion, Z-estimator expansion for the control-function MLE, and the final averaging step.
States a high-level sequential (partial-sample) linearization condition under which √n{δ̂(t,y) − δ(y)} = (√n/nt) Σi=1nt ξi(y) + op(1) uniformly in (t,y) ∈ [ε,1]×I. Primitive sufficient conditions given in Proposition 5.1.
The key observation is that both √n(T̂n−T) and √n V̂n converge to τ times a Brownian-motion functional. The τ cancels in the ratio (T̂n−T)/V̂n, giving the pivotal limit W — free of nuisance parameters without bootstrap.
Under Assumptions 5.1–5.5, the test has asymptotic level α on the null d ≤ Δ, exact asymptotic level α on the boundary d = Δ, is consistent against d > Δ, and detects local alternatives at rate n−1/2 (Corollary 5.8).
Simulations use a triangular model with one exogenous regressor X, one instrument Z, one endogenous regressor Y₂, and a continuous outcome Y. The endogeneity strength is controlled by ρ ∈ [0,1) — when ρ = 0, Y₂ is exogenous; when ρ > 0, Y₂ is endogenous. Distribution regression is implemented using a probit link, evaluated at the covariate profile (x, y₂) = (0, 0). Sample size n = 1,000; nominal level α = 0.05.
Table 1 shows that even at ρ = 0.5 — substantial correlation between structural error and first-stage residual — the estimated distance d̂ = 0.026 remains below the relevance threshold Δ = 0.05, and the test correctly does not reject.
| ρ | Δ | n | Rejection rate | Mean T̂n | Mean d̂ |
|---|---|---|---|---|---|
| 0.0 | 0.05 | 1000 | 0.00 | 2.82 × 10−⁶ | 0.0014 |
| 0.2 | 0.05 | 1000 | 0.00 | 3.95 × 10−⁶ | 0.0057 |
| 0.5 | 0.05 | 1000 | 0.00 | 7.23 × 10−⁴ | 0.0263 |
| B = 100 Monte Carlo replications. The test does not reject when d < Δ, regardless of statistical significance of endogeneity. | |||||
Table 2 shows power as Δ varies. Rejection probability decreases as Δ increases — the test is sensitive to the threshold rather than merely to the presence of nonzero endogeneity. The transition occurs near Δ ≈ 0.027, consistent with the estimated distance d̂ ≈ 0.027.
| ρ | Δ | n | Rejection rate | Mean d̂ |
|---|---|---|---|---|
| 0.5 | 0.005 | 1000 | 0.65 | 0.0284 |
| 0.5 | 0.010 | 1000 | 0.62 | 0.0268 |
| 0.5 | 0.020 | 1000 | 0.21 | 0.0274 |
| 0.5 | 0.030 | 1000 | 0.01 | 0.0269 |
| 0.5 | 0.040 | 1000 | 0.00 | 0.0279 |
| 0.5 | 0.050 | 1000 | 0.00 | 0.0277 |
| Power decreases monotonically with Δ. Transition near Δ ≈ 0.027 confirms the test targets practical relevance, not mere statistical significance. | ||||
When Δ is set close to the true distance d(x,y₂) ≈ 0.027, rejection rates converge to the nominal level α = 0.05 for both n = 1,000 and n = 2,000, confirming the theoretical prediction of Theorem 5.7.
The test is applied to Card's (1995) returns-to-schooling data (n = 3,010), where education is treated as potentially endogenous and proximity to a four-year college (nearc4) serves as the excluded instrument. Controls include experience, experience squared, Black, SMSA, and South indicators.
The standard linear Wu–Hausman test fails to reject exogeneity (p-value = 0.215), despite the OLS estimate (0.0740) and the IV estimate (0.1323) differing substantially in magnitude. This motivates asking a distributional question: even if mean endogeneity is not statistically detectable, does the IV correction shift the distribution of log wages by a practically relevant amount?
The test is applied at three covariate profiles representing low (educ=10, exper=4), middle (educ=13, exper=8), and high (educ=17, exper=15) education/experience workers. The findings are strikingly heterogeneous.
| Profile | Δ | T̂n | d̂ | Critical value | Reject? |
|---|---|---|---|---|---|
| Low (educ=10) | 0.01 | 0.0219 | 0.1479 | 0.0254 | No |
| 0.02 | 0.0219 | 0.1479 | 0.0257 | No | |
| 0.03 | 0.0219 | 0.1479 | 0.0262 | No | |
| 0.05 | 0.0219 | 0.1479 | 0.0278 | No | |
| Middle (educ=13) | 0.01 | 0.0009 | 0.0304 | 0.0129 | No |
| 0.02 | 0.0009 | 0.0304 | 0.0132 | No | |
| 0.03 | 0.0009 | 0.0304 | 0.0137 | No | |
| 0.05 | 0.0009 | 0.0304 | 0.0153 | No | |
| High (educ=17) | 0.01 | 0.0380 | 0.1949 | 0.0353 | Yes |
| 0.02 | 0.0380 | 0.1949 | 0.0356 | Yes | |
| 0.03 | 0.0380 | 0.1949 | 0.0361 | Yes | |
| 0.05 | 0.0380 | 0.1949 | 0.0377 | Yes | |
| Baseline specification. ε = 0.2, central 10–90% log-wage grid with 31 points, q0.95 ≈ 2.69. Rejection rule: T̂n > Δ² + q1−αV̂n. | |||||
Figures 1 and 2 plot the estimated conditional CDFs of log wages for naive DR (solid) and IV/control-function DR (dashed) at the three covariate profiles. The curves are closest for the middle profile; the low and high profiles show visibly larger discrepancies — especially the high profile, where the IV correction shifts the distribution substantially.
Baseline controls (experience, experience², Black, SMSA, South). Solid: naive DR; dashed: IV-DR. The high-education profile (right panel) shows the largest gap.
Full controls (baseline + smsa66 + region indicators). Qualitative pattern persists: middle profile closest, low and high profiles most discrepant. Formal rejection for the high profile disappears with the richer control set.
Integrates the pointwise test over all covariate profiles with weight measure μ, yielding a global relevance test D² = ∫ ∫ [F𝑁(y|x,y₂) − F𝑉(y|x,y₂)]² dy dμ(x,y₂). Natural choices: empirical, policy-relevant, or population measure.
DR estimators are not guaranteed monotone in y in finite samples. Both monotone rearrangement (Chernozhukov et al. 2010) and isotonic regression (Wied 2024) can be applied — the √n-convergence rate and Gaussian limit carry over.
Inverting the test gives a one-sided lower confidence set [Δ̂α, ∞) for the practical endogeneity distance, providing a data-driven summary of endogeneity magnitude at each covariate profile.
The framework extends to Y₂ ∈ ℝm with m > 1 endogenous regressors and ℓ ≥ m instruments. The first stage becomes a system of linear equations; the asymptotic results carry over under multivariate generalizations.