Why Second-Order Remainders Matter
An estimating equation can depend on both a target parameter and nuisance functions. The target is the quantity of interest; nuisance functions help construct a moment that identifies it. In practice, those functions must be estimated. The key question is whether their estimation errors distort the target equation at a consequential order.
Neyman orthogonality makes the population estimating equation locally insensitive, to first order, to perturbations of the nuisance functions. It does not make the equation insensitive to every possible error. After the first-order effect vanishes, a second-order remainder remains. Product-rate conditions control that remainder by requiring the nuisance estimation errors to be sufficiently small jointly.
More precisely, for any admissible direction \(h\), orthogonality means \( \left.\frac{d}{dt}M(\theta_0,\eta_0+th)\right|_{t=0}=0 \). The target equation still needs to identify \(\theta_0\): changing \(\theta\) should change the moment in a way that permits the parameter to be distinguished. Orthogonality concerns perturbations of \(\eta\), not perturbations of \(\theta\).
From Orthogonality to a Remainder Bound
Write an estimated nuisance function as \(\widehat{\eta}=\eta_0+\Delta\eta\). A Taylor expansion of the population moment around the truth has the form
where \(D_{\eta}M[\Delta\eta]\) is the first directional derivative, \(D_{\eta}^{2}M[\Delta\eta,\Delta\eta]\) is the second derivative applied twice to the perturbation, and \(r\) is the higher-order Taylor remainder. The first term is zero because the true moment is zero; the first derivative is zero by orthogonality. If the second derivative is bounded in a suitable norm and the higher-order remainder is controlled, the remaining moment error is quadratic in the nuisance error. With multiple nuisance components, it often takes the more useful form of a product of errors.
That qualification matters. Orthogonality alone is a statement about a derivative at one point. It does not, by itself, guarantee a particular quadratic bound far from that point. A product-rate argument needs both orthogonality and enough regularity to control the remainder over the estimation errors under consideration.
Worked Example: An Exact Product Remainder for a Missing-Data Mean
Suppose \(Y\) is an outcome, \(X\) is a covariate, and \(R\) indicates whether \(Y\) is observed. Assume missingness is conditionally independent of the outcome given \(X\), with \(\pi_0(X)=P(R=1\mid X)\), and define \(m_0(X)=E[Y\mid X]\). Assume \(\pi_0(X)\geq c>0\). For the population mean \(\theta_0=E[Y]\), consider the score \(\psi(W;\theta,m,\pi)=m(X)+\frac{R}{\pi(X)}\{Y-m(X)\}-\theta\).
For any candidate functions \(\widehat m,\widehat\pi\), condition on \(X\). Since \(E[R\{Y-\widehat m(X)\}\mid X]=\pi_0(X)\{m_0(X)-\widehat m(X)\}\), the conditional expectation of the score at \(\theta_0\) is
Taking expectations gives the exact remainder \(R=E[(\widehat\pi-\pi_0)(\widehat m-m_0)/\widehat\pi]\). At the truth, a first-order perturbation of either nuisance alone contributes zero: if \(\widehat m=m_0\) or \(\widehat\pi=\pi_0\), the displayed remainder is zero. This illustrates orthogonality and its product structure directly.
If \(\widehat\pi(X)\geq c/2\) and the errors have finite second moments, the Cauchy–Schwarz inequality gives
Thus, if the two \(L_2\) errors are respectively \(O_p(n^{-0.30})\) and \(O_p(n^{-0.23})\), the remainder is \(O_p(n^{-0.53})\), which is \(o_p(n^{-1/2})\). This conclusion uses the stated lower bound on \(\widehat\pi\); without a denominator bound, the product inequality need not hold.
Product Rates in a Partially Linear Model
A second example shows why a product bound can involve errors from distinct nuisance functions. Suppose \(Y=\theta_0D+g_0(X)+U\), with \(E[U\mid D,X]=0\), and let \(m_0(X)=E[D\mid X]\). Consider the score \(\psi(W;\theta,m,g)=(D-m(X))\{Y-\theta D-g(X)\}\). Assume \(E[(D-m_0(X))^2]>0\), so the residualized regressor has nonzero variation.
At the true parameter, the residual in the second factor is \(U\). The conditional mean restriction implies \(E[U\mid X]=0\) and \(E[(D-m_0(X))U\mid X]=0\). Consequently, the true moment is zero. The derivative of the moment with respect to \(\theta\) at the truth is \(-E[(D-m_0(X))D]=-E[(D-m_0(X))^2]\), which is nonzero by assumption.
Worked Example: Bounding the Partially Linear Remainder
Let \(\widehat m=m_0+\Delta_m\) and \(\widehat g=g_0+\Delta_g\). At \(\theta=\theta_0\), the score becomes \((D-m_0(X)-\Delta_m(X))(U-\Delta_g(X))\). Expanding and taking expectations,
The first term is zero by the conditional mean restriction; the second and third terms are zero by conditioning on \(X\). The remaining term is a product. By Cauchy–Schwarz, \(|M(\theta_0,\widehat m,\widehat g)|\leq\|\Delta_m\|_2\|\Delta_g\|_2\). If the two errors are \(O_p(n^{-0.28})\) and \(O_p(n^{-0.24})\), respectively, the bound is \(O_p(n^{-0.52})=o_p(n^{-1/2})\). Both nuisance functions may converge more slowly than \(n^{-1/2}\); it is their product that controls this population remainder.
The calculation also verifies orthogonality in this model: setting either nuisance perturbation to zero makes the remainder zero, so there is no first-order nuisance contribution. The nonzero derivative with respect to \(\theta\) ensures that a small moment error can be translated into a small target-parameter error, with scale determined by \(E[(D-m_0(X))^2]\).
What the Rate Condition Requires
Suppose a target estimator is based on an empirical estimating equation and the target moment changes at a nonzero rate as \(\theta\) changes. For the nuisance-induced population remainder not to affect a root-\(n\) limit, a common requirement is \(R(\widehat{\eta})=o_p(n^{-1/2})\). For a product remainder, sufficient conditions take the form
This condition concerns the product, not a universal requirement that each nuisance function converge at the same rate. If both errors are \(O_p(n^{-1/4})\), their product is only \(O_p(n^{-1/2})\); that bound alone does not establish the stricter \(o_p(n^{-1/2})\) condition. If one error is \(O_p(n^{-0.30})\) and the other is \(O_p(n^{-0.23})\), their product is \(O_p(n^{-0.53})\), which does meet it.
For a quadratic remainder bounded by \(C\|\widehat{\eta}-\eta_0\|_2^2\), it is sufficient that \(\|\widehat{\eta}-\eta_0\|_2=o_p(n^{-1/4})\). More generally, with several nuisance errors, the relevant condition follows from the actual remainder bound: one component may need a faster rate, or one error may compensate for another. Always derive or establish the bound for the score at hand rather than assuming a product condition automatically applies.
Specify the true parameter and nuisance functions, and verify that the moment is zero at the truth.
Differentiate the population moment in admissible nuisance directions and verify that the derivative vanishes.
Expand the moment or use an exact identity, stating the smoothness, integrability, and denominator conditions used.
Combine the nuisance rates in the bound and check whether the remainder is \(o_p(n^{-1/2})\) when root-\(n\) inference is intended.
Common Pitfall: Orthogonality Is Not a Rate Guarantee
A common mistake is to reason, “The score is orthogonal, so nuisance estimation cannot affect the target.” Orthogonality only eliminates the first-order population effect at the true nuisance functions. A second-order remainder remains, and it can still be too large. For example, two errors each of order \(n^{-0.20}\) have a product of order \(n^{-0.40}\), which is not negligible relative to \(n^{-1/2}\).
There is also a distinction between controlling the population remainder and controlling the full empirical estimating equation. The product-rate bound addresses the former. Sampling fluctuation from estimating and evaluating the score, dependence between observations, and the validity of a limiting distribution require their own conditions. In particular, a product-rate condition by itself does not prove asymptotic normality or justify a confidence interval.
Check Your Understanding
Use the definitions and calculations above to answer these questions.
- What derivative must vanish for an estimating equation to be Neyman-orthogonal at the truth?
- In the missing-data example, why is a lower bound on \(\widehat\pi(X)\) needed to obtain the stated \(L_2\) product bound?
- If two nuisance errors are \(O_p(n^{-0.27})\) and \(O_p(n^{-0.25})\), what order does their product have, and is it \(o_p(n^{-1/2})\)?
- Why does orthogonality alone not imply that nuisance estimation has no effect on the target estimator?
- What additional kinds of conditions, beyond a small population remainder, are needed to support root-\(n\) inference?