Skip to content

Huber contamination model

An epsilon-contamination neighbourhood mixes a clean law with arbitrary errors. A 5% point mass shows why mean sensitivity differs from median bias.

The Huber contamination model describes an observed distribution as a mostly trusted model mixed with an arbitrary contaminating distribution. For a quant estimating the centre of a return or forecast-error series, it makes a precise distinction: a small probability of bad observations need not mean a small error in the estimated mean.

A neighbourhood of possible distributions

Write the observed law as

Fϵ=(1−ϵ)P+ϵQ,0≤ϵ<1.F_\epsilon=(1-\epsilon)P+\epsilon Q, \qquad 0\le\epsilon<1.
Here PP is the clean reference distribution, QQ is any probability distribution on the same observation space, and ϵ\epsilon is the contamination weight. Holding PP and ϵ\epsilon fixed while allowing QQ to vary gives the contamination neighbourhood
Hϵ(P)={(1−ϵ)P+ϵQ:Q is a probability distribution}.\mathcal H_\epsilon(P) =\{(1-\epsilon)P+\epsilon Q:Q\text{ is a probability distribution}\}.
This is a mixture distribution model. The contaminating observations need not be conspicuous outliers: QQ can overlap PP or introduce a subtle directional shift. The values and component labels are generally unknown. Prasad and coauthors formulate the model with arbitrary QQ, including subtle departures from the assumed law.
Under the usual independent, identically distributed sampling interpretation, each observation comes from QQ with probability ϵ\epsilon. The number of contaminated component draws in nn observations has a binomial distribution. It is not exactly ϵn\epsilon n in each sample. An adversary who inspects a clean sample and then replaces a fixed number of selected observations defines a different corruption mechanism; guarantees require the stated sampling model.
In a location problem, the goal is often the centre of PP, rather than the mean of the observed mixture. UBC's robust-location notes emphasise estimating the central distribution without having to identify the contaminating component.

The same 5%, two very different means

Consider the synthetic population law

FM=0.95 N(0,1)+0.05 δM,F_M=0.95\,\mathcal N(0,1)+0.05\,\delta_M,
where δM\delta_M puts all its probability at the value MM. The clean normal distribution has centre zero and standard deviation one. These are standardised units, not observed market returns.

For this finite-valued contamination, the population mean is

EFM[X]=0.95(0)+0.05M.\mathbb E_{F_M}[X]=0.95(0)+0.05M.
Contaminating value MMPopulation meanPopulation median
2020110.066010.06601
1000100050500.066010.06601
The median calculation uses probability rather than distance. Below either point mass, the mixture's cumulative distribution function is 0.95Φ(x)0.95\Phi(x), where Φ\Phi is the standard-normal CDF. Its unique median therefore solves
0.95Φ(m)=0.5,m=Φ−1 ⁣(0.50.95)≈0.0660118.0.95\Phi(m)=0.5, \qquad m=\Phi^{-1}\!\left(\frac{0.5}{0.95}\right) \approx0.0660118.
Both 2020 and 10001000 exceed this solution, so moving the point mass between those values leaves the median unchanged. It has nevertheless shifted away from the clean centre zero. The population-median bias calculation also appears in the short course on robust statistics, §3.2.
Synthetic population mixtures contain 95% standard-normal mass and 5% point mass at M. Moving M from 20 to 1000 raises the mean from 1 to 50. Both medians are 0.06601, where the CDF 0.95 Phi(x) crosses 0.5; its value at zero is 0.475.
Synthetic population mixtures contain 95% standard-normal mass and 5% point mass at M. Moving M from 20 to 1000 raises the mean from 1 to 50. Both medians are 0.06601, where the CDF 0.95 Phi(x) crosses 0.5; its value at zero is 0.475.
Original synthetic population comparison. The left panel varies the location of the 5% point mass; the right panel magnifies the common CDF near the median. The panels use different horizontal axes.

This is a sensitivity example, not a recommendation to forecast returns with a median. Mean and median answer different questions, and actual return dependence or genuine tail events require their own model. Nor can the mean formula be applied blindly to every allowed QQ: an arbitrary contaminating distribution may have no finite first moment, leaving the mixture mean undefined or nonfinite.

The model is not the loss function

Huber loss is quadratic near zero and linear in the tails. It is an ingredient for constructing an estimator. The contamination model instead specifies which data distributions the estimator must tolerate. Selecting that loss does not determine QQ, reveal which observations are contaminated, or establish a universal error guarantee. The UBC notes, §1.1 separate assumptions on the distribution from the estimating equation.
An influence function studies an estimator functional's derivative under an infinitesimal point-mass perturbation. A bounded influence function limits this first-order local response under the required regularity conditions. It does not mean zero contamination bias, immunity to all corruptions, or a global finite-sample bound at a fixed contamination fraction. The short course, §4.1 defines that derivative and distinguishes it from breakdown behaviour.
For a refresher on population parameters versus sample summaries, Basic Statistics and Regression introduces both mean and median; the full lesson is available with a free signed-in account.

Sources

Last updated
  • robust statistics
  • contamination
  • mixture model
  • location estimation
  • median
← Back to the glossaryBack to top ↑