Chauvenet's criterion

The idea behind Chauvenet's criterion finds a probability band that reasonably contains all n samples of a data set, centred on the mean of a normal distribution. By doing this, any data point from the n samples that lies outside this probability band can be considered an outlier, removed from the data set, and a new mean and standard deviation based on the remaining values and new sample size can be calculated. This identification of the outliers will be achieved by finding the number of standard deviations that correspond to the bounds of the probability band around the mean (${\displaystyle D_{\mathrm {max} }}$) and comparing that value to the absolute value of the difference between the suspected outliers and the mean divided by the sample standard deviation (Eq.1).

${\displaystyle D_{\mathrm {max} }\geq {\frac {|x-{\bar {x}}|}{s_{x}}}}$

(1)

where

• ${\displaystyle D_{\mathrm {max} }}$ is the maximum allowable deviation,
• ${\displaystyle |\cdot |}$ is the absolute value,
• ${\displaystyle x}$ is the value of suspected outlier,
• ${\displaystyle {\bar {x}}}$ is sample mean, and
• ${\displaystyle s_{x}}$ is sample standard deviation.

In order to be considered as including all ${\displaystyle n}$ observations in the sample, the probability band (centered on the mean) must only account for ${\displaystyle n-{\tfrac {1}{2}}}$ samples (if ${\displaystyle n=3}$ then only 2.5 of the samples must be accounted for in the probability band). In reality we cannot have partial samples so ${\displaystyle n-{\tfrac {1}{2}}}$ (2.5 for ${\displaystyle n=3}$) is approximately ${\displaystyle n}$. Anything less than ${\displaystyle n-{\tfrac {1}{2}}}$ is approximately ${\displaystyle n-1}$ (2 if ${\displaystyle n=3}$) and is not valid because we want to find the probability band that contains ${\displaystyle n}$ observations, not ${\displaystyle n-1}$ samples. In short, we are looking for the probability, ${\displaystyle P}$, that is equal to ${\displaystyle n-{\tfrac {1}{2}}}$ out of ${\displaystyle n}$ samples (Eq.2).

${\displaystyle P={\frac {n-{\tfrac {1}{2}}}{n}}=1-{\tfrac {1}{2n}}}$

(2)

where

• ${\displaystyle P}$ is the probability band centered on the sample mean and
• ${\displaystyle n}$ is the sample size.

The quantity ${\displaystyle {\tfrac {1}{2n}}}$ corresponds to the combined probability represented by the two tails of the normal distribution that fall outside of the probability band ${\displaystyle P}$. In order to find the standard deviation level associated with ${\displaystyle P}$, only the probability of one of the tails of the normal distribution needs to be analyzed due to its symmetry (Eq.3).

${\displaystyle P_{z}={\frac {1}{4n}}}$

(3)

where

• ${\displaystyle P_{z}}$ is probability represented by one tail of the normal distribution and
• ${\displaystyle n}$ = sample size.

Eq.1 is analogous to the ${\displaystyle Z}$-score equation (Eq.4).

${\displaystyle Z={\frac {x-\mu }{\sigma }}}$

(4)

where

• ${\displaystyle Z}$ is the ${\displaystyle Z}$-score,
• ${\displaystyle x}$ is the sample value,
• ${\displaystyle \mu =0}$ is the mean of standard normal distribution, and
• ${\displaystyle \sigma =1}$ is the standard deviation of standard normal distribution.

Based on Eq.4, to find the ${\displaystyle D_{\mathrm {max} }}$ (Eq.1) find the z-score corresponding to ${\displaystyle P_{z}}$ in a ${\displaystyle Z}$-score table. ${\displaystyle D_{\mathrm {max} }}$ is equal to the score for ${\displaystyle P_{z}}$. Using this method ${\displaystyle D_{\mathrm {max} }}$ can be determined for any sample size. In Excel, ${\displaystyle D_{\mathrm {max} }}$ can be found with the following formula: =ABS(NORM.S.INV(1/(4n))).

To apply Chauvenet's criterion, first calculate the mean and standard deviation of the observed data. Based on how much the suspect datum differs from the mean, use the normal distribution function (or a table thereof) to determine the probability that a given data point will be at the value of the suspect data point. Multiply this probability by the number of data points taken. If the result is less than 0.5, the suspicious data point may be discarded, i.e., a reading may be rejected if the probability of obtaining the particular deviation from the mean is less than ${\displaystyle {\tfrac {1}{2n}}}$.

For instance, suppose a value is measured experimentally in several trials as 9, 10, 10, 10, 11, and 50, and we want to find out if 50 is an outlier.

First, we find ${\displaystyle P_{z}}$.

${\displaystyle P_{z}=1-{\frac {1}{4n}}=1-{\frac {1}{4\times 6}}=1-{\frac {1}{24}}\approx .9583}$

Then we find ${\displaystyle D_{max}}$ by plugging ${\displaystyle P_{z}}$ into the Quantile Function.

${\displaystyle D_{max}=Q(P_{z})\approx 1.7317}$

Then we find the z-score of 50.

${\displaystyle z={\frac {50-{\bar {x}}}{s_{x}}}={\frac {50-16.67}{16.34}}\approx 2.04}$

From there we see that ${\displaystyle z>D_{max}}$ and can conclude that 50 is an outlier according to Chauvenet's Criterion.

Another method for eliminating spurious data is called Peirce's criterion. It was developed a few years before Chauvenet's criterion was published, and it is a more rigorous approach to the rational deletion of outlier data.[3] Other methods such as Grubbs's test for outliers are mentioned under the listing for Outlier.

Deletion of outlier data is a controversial practice frowned on by many scientists and science instructors; while Chauvenet's criterion provides an objective and quantitative method for data rejection, it does not make the practice more scientifically or methodologically sound, especially in small sets or where a normal distribution cannot be assumed. Rejection of outliers is more acceptable in areas of practice where the underlying model of the process being measured and the usual distribution of measurement error are confidently known.

• Taylor, John R. An Introduction to Error Analysis. 2nd edition. Sausalito, California: University Science Books, 1997. pp 166–8.
• Barnett, Vic and Lewis, Toby. "Outliers in Statistical Data". 3rd edition. Chichester: J.Wiley and Sons, 1994. ISBN 0-471-93094-6.
• Aicha Zerbet, Mikhail Nikulin. A new statistics for detecting outliers in exponential case, Communications in Statistics: Theory and Methods, 2003, v.32, pp. 573–584.