Tutorials › Econometrics (Graduate) › Inference and Diagnostics from Estimated Influence Functions

Semiparametric Efficiency · Tutorial 804 of 952

Inference and Diagnostics from Estimated Influence Functions

Use estimated influence functions to quantify uncertainty and assess whether standard asymptotic inference is credible.

Advanced 12 min read

What You'll Learn

  • Connect an asymptotically linear estimator to the variance of its influence function.
  • Calculate an estimated standard error and Wald confidence interval from observation-level influence values.
  • State the conditions that justify normal approximation and variance estimation.
  • Diagnose whether a small number of influence values dominates estimated uncertainty.
  • Distinguish a warning sign in a finite sample from evidence that a required moment condition fails.

From an Influence Function to Uncertainty

An estimate alone does not show how precisely a target has been learned. In semiparametric analysis, the influence function links an estimator’s sampling variation to observation-level contributions. Once those contributions can be estimated, their empirical spread provides a standard error and supports confidence intervals and tests. The calculations are direct; their justification depends on the estimator having the right asymptotic behavior and on the contributions having sufficiently well-behaved moments.

Let \(O_1,\ldots,O_n\) be independent, identically distributed observations from a distribution \(P\), and let \(\theta(P)\) be the target parameter. Write \(\phi_P(O)\) for an influence function for this target. It has mean zero under \(P\). The influence functions developed for particular targets, such as means or missing-data means, provide examples of this general object.

Definition: An estimator \(\hat\theta\) is asymptotically linear with influence function \(\phi_P\) if
$$ \sqrt{n}\{\hat\theta-\theta(P)\} = \frac{1}{\sqrt{n}}\sum_{i=1}^{n}\phi_P(O_i)+o_P(1). $$
Here \(o_P(1)\) denotes a term that converges to zero in probability. The influence function has \(E_P\{\phi_P(O)\}=0\); for standard normal inference it must also have finite, positive variance \(V=E_P\{\phi_P(O)^2\}\).

The expansion says that, to first order, the estimation error is an average of independent contributions. If \(V\) is finite and positive, the central limit theorem gives \(\sqrt{n}\{\hat\theta-\theta(P)\}/\sqrt{V}\) an asymptotic standard normal distribution. Thus the approximate variance of \(\hat\theta\) is \(V/n\), and its standard error is \(\sqrt{V/n}\). In practice, \(P\), and therefore \(\phi_P\), is unknown. We use estimated influence values \(\hat\phi_i\), one for each observation, instead.

Estimating the Variance and Constructing Inference

The empirical variance of the estimated influence values estimates the variance of a single contribution. Centering is important: even when the population influence function has mean zero, estimated values need not sum to exactly zero in a finite sample. Define \(\bar{\hat\phi}=n^{-1}\sum_{i=1}^{n}\hat\phi_i\).

Variance estimate and standard error: Use the sample variance of the estimated influence values and divide by the sample size to estimate the variance of \(\hat\theta\):
$$ \hat V = \frac{1}{n-1}\sum_{i=1}^{n} (\hat\phi_i-\bar{\hat\phi})^2, \qquad \widehat{SE}(\hat\theta) = \sqrt{\frac{\hat V}{n}}. $$

Under the assumptions below, an approximate two-sided \(100(1-\alpha)\%\) Wald confidence interval is \(\hat\theta\pm z_{1-\alpha/2}\widehat{SE}(\hat\theta)\), where \(z_{1-\alpha/2}\) is the corresponding standard normal quantile. For a null value \(\theta_0\), the Wald statistic is \((\hat\theta-\theta_0)/\widehat{SE}(\hat\theta)\); its two-sided approximate p-value is \(2\{1-\Phi(|T|)\}\), where \(\Phi\) is the standard normal cumulative distribution function. These are large-sample approximations, not exact finite-sample guarantees.

Conditions for this procedure: The observations are independent and identically distributed; the estimator has the stated asymptotically linear expansion; \(0<V=E_P\{\phi_P(O)^2\}<\infty\); and the estimated influence values consistently reproduce the empirical variance of the true values. One sufficient condition for the last requirement is
$$ \frac{1}{n}\sum_{i=1}^{n} \{\hat\phi_i-\phi_P(O_i)\}^2=o_P(1). $$
In particular applications, asymptotic linearity and this approximation condition depend on how the estimator and any nuisance functions are constructed.

The empirical approximation condition is a useful way to express what is needed of estimated contributions; it does not by itself establish the estimator’s asymptotic linearity. For example, estimating nuisance functions may require additional conditions on their errors and on the estimator’s remainder. When nuisance predictions are obtained by cross-fitting, each observation’s contribution is evaluated using nuisance estimates trained without that observation. Whatever the construction, the standard error is justified only when the full estimator meets the asymptotic conditions, not merely because a software routine returns influence values.

Worked Example: A Wald Interval From Estimated Influence Values

Suppose a hypothetical analysis has \(n=100\) independent observations and estimates a mean outcome in days as \(\hat\theta=2.40\) days. Its estimated influence values have \(\bar{\hat\phi}=0\) days and \(\sum_{i=1}^{100}(\hat\phi_i-\bar{\hat\phi})^2=144\) days squared. Assume the estimator’s asymptotic linearity and the required influence-value approximation are credible, and use a 95% Wald interval.

The estimated variance of one influence value is \[ \hat V=\frac{144}{100-1}=1.4545\text{ days}^2. \] Therefore, \[ \widehat{SE}(\hat\theta) =\sqrt{\frac{1.4545}{100}} =0.1206\text{ days}. \] Using \(z_{0.975}=1.96\), the margin of error is \(1.96(0.1206)=0.2364\) days. The interval is \[ 2.40\pm0.2364=(2.164,\ 2.636)\text{ days}. \]

Under the stated large-sample assumptions, this procedure has approximately 95% coverage across repeated samples. It does not mean there is a 95% probability that the fixed target lies in this particular interval.

Justification for the Plug-in Standard Error

The variance calculation follows from the asymptotic linear representation and consistency of the empirical variance. The reasoning also clarifies why checking the influence values matters.

1
Apply the central limit theorem.
With independent observations, a mean-zero influence function, and finite positive variance \(V\), the normalized sum \(n^{-1/2}\sum_i\phi_P(O_i)\) converges in distribution to \(N(0,V)\). The asymptotic linear expansion transfers this limit to \(\sqrt{n}\{\hat\theta-\theta(P)\}\).
2
Estimate the variance of one contribution.
Because \(E_P\{\phi_P(O)^2\}<\infty\), the law of large numbers applies to the squared influence values. The empirical variance of the true values converges in probability to \(V\).
3
Replace true values by estimated values.
The empirical mean squared difference between \(\hat\phi_i\) and \(\phi_P(O_i)\) converges to zero. By the Cauchy–Schwarz inequality, the difference between their empirical second moments also converges to zero, provided the empirical second moments of the true values remain bounded in probability. Their empirical means converge to zero as well. Hence \(\hat V\) converges in probability to \(V\).
4
Studentize the estimator.
Since \(\hat V\) converges to positive \(V\), the continuous mapping theorem gives \(\sqrt{\hat V}\) converging to \(\sqrt V\). Slutsky’s theorem then yields \(\{\hat\theta-\theta(P)\}/\widehat{SE}(\hat\theta)\) converging in distribution to \(N(0,1)\), which justifies the stated Wald interval and test.

Diagnosing Extreme Contributions and Unstable Variances

An influence value is a first-order contribution to estimation error, not a direct measure of an observation’s substantive importance or a guarantee that deleting it would produce a particular change. Nevertheless, the squared values determine the estimated variance. A small number of large absolute values can therefore dominate the standard error. Inspecting the empirical distribution and the contributions to the sum of squares can reveal instability that a single reported standard error hides.

Worked Example: One Observation Dominates the Estimated Variance

Consider eight hypothetical estimated influence values, in units of the target: \[ -1,\ -1,\ -1,\ -1,\ -1,\ -1,\ -1,\ 7. \] Their sum is zero, so \(\bar{\hat\phi}=0\). Their sum of squares is \(7(1^2)+7^2=7+49=56\). Thus \[ \hat V=\frac{56}{8-1}=8, \qquad \widehat{SE}(\hat\theta)=\sqrt{\frac{8}{8}}=1 \] in the target’s units.

The single value 7 accounts for \(49/56=0.875\), or 87.5%, of the total squared contribution. This is a clear warning that the estimated variance depends heavily on one observation. It does not establish that the value is erroneous or that the standard error is invalid. A careful follow-up checks data quality, the observation’s covariates and fitted nuisance quantities, and whether the influence-function formula is being evaluated as intended. Reporting a sensitivity calculation may help describe dependence on that observation, but deleting it or capping its influence changes the procedure and does not automatically restore valid inference.

Variance instability can also arise when the estimated influence values themselves are sensitive to fitted nuisance functions, or when a few observations have extreme weights or residuals. Useful diagnostics include plotting the values, reporting their range and selected quantiles, and calculating the fraction of the sum of squares contributed by the largest absolute values. Such summaries are diagnostics, not formal tests of the model assumptions. A finite-sample plot cannot establish that a population moment exists or does not exist.

Moment Conditions and Limits of Standard Inference

A finite second moment is central to the standard procedure. The central limit theorem used here requires a finite variance for the influence function. The empirical variance also needs to estimate a finite population variance. If the influence function has infinite variance, the usual square-root-\(n\) normal approximation and its plug-in Wald interval are not justified by this argument. A large observed influence value is a reason to investigate, but the absence of an extreme value in a modest sample cannot verify the required moment condition.

Worked Example: A Mean With an Infinite-Variance Influence Function

Suppose, hypothetically, that an outcome \(X\) has a Pareto distribution with lower bound 1 and tail index \(1.5\), so its density is \(f(x)=1.5x^{-2.5}\) for \(x\geq1\). The target is its mean. For a mean, the influence function is \(\phi_P(X)=X-\mu\), where \(\mu=E_P(X)\).

The mean exists: \[ E_P(X)=\int_1^\infty x(1.5x^{-2.5})\,dx =1.5\int_1^\infty x^{-1.5}\,dx=3. \] But the second raw moment diverges: \[ E_P(X^2) =1.5\int_1^\infty x^{-0.5}\,dx=\infty. \] Consequently \(E_P\{(X-3)^2\}=\infty\), so the influence function does not have finite variance. Even if a finite hypothetical sample yields a calculable sample variance, that number does not estimate a finite population variance under this distribution. The usual standard error and Wald interval from the finite-variance theory are therefore unsupported; merely increasing the reported number of decimal places would not address the problem.

Common Pitfall and Practical Interpretation

A common mistake is to assume that calculating an empirical variance proves standard asymptotic inference is valid. The calculation only produces a number. Its interpretation depends on the estimator’s asymptotic linearity, the relevant moment conditions, and whether estimated influence values adequately approximate the true contributions. If these conditions are doubtful, the interval should not be presented as automatically reliable; explain the concern and investigate the estimator and data-generating features that could cause it.

Keep the units clear throughout. Influence values have the units of the target, their variance has squared units, and the standard error and confidence interval have the target’s units. An interval describes uncertainty under the stated sampling and asymptotic assumptions; it does not remove bias from a misspecified target, guarantee causal interpretation, or establish that the assumptions hold.

Key takeaway: Estimated influence values provide a principled standard error when the estimator is asymptotically linear, the influence function has finite positive variance, and the estimated contributions reproduce its empirical spread. Inspecting extreme contributions and moment plausibility helps assess the calculation, but diagnostics cannot substitute for the required assumptions.

Check Your Understanding

Use the assumptions as well as the arithmetic when answering.

  1. Why is the empirical variance of the influence values divided by \(n\) to estimate the variance of \(\hat\theta\)?
  2. An estimated influence function has sample mean different from zero. Why should the variance calculation center the values?
  3. If one observation contributes 70% of the sum of squared influence values, what does this suggest, and what does it not establish?
  4. Which moment condition supports the standard normal inference described here, and what can go wrong if it fails?
  5. Why does calculating a finite sample standard error not prove that the population influence function has finite variance?
Clinical Biostats