Chebyshev's inequality

From HandWiki
Short description: Bound on probability of a random variable being far from its mean


In probability theory, Chebyshev's inequality (also called the Bienaymé–Chebyshev inequality) provides an upper bound on the probability of deviation of a random variable (with finite variance) from its mean. More specifically, the probability that a random variable deviates from its mean by more than [math]\displaystyle{ k\sigma }[/math] is at most [math]\displaystyle{ 1/k^2 }[/math], where [math]\displaystyle{ k }[/math] is any positive constant.

The rule is often called Chebyshev's theorem, about the range of standard deviations around the mean, in statistics. The inequality has great utility because it can be applied to any probability distribution in which the mean and variance are defined. For example, it can be used to prove the weak law of large numbers.

Its practical usage is similar to the 68–95–99.7 rule, which applies only to normal distributions. Chebyshev's inequality is more general, stating that a minimum of just 75% of values must lie within two standard deviations of the mean and 88.89% within three standard deviations for a broad range of different probability distributions.[1][2]

The term Chebyshev's inequality may also refer to Markov's inequality, especially in the context of analysis. They are closely related, and some authors refer to Markov's inequality as "Chebyshev's First Inequality," and the similar one referred to on this page as "Chebyshev's Second Inequality."

Chebyshev's inequality is tight in the sense that for each chosen positive constant, there exists a random variable such that the inequality is in fact an equality.[3]

History

The theorem is named after Russian mathematician Pafnuty Chebyshev, although it was first formulated by his friend and colleague Irénée-Jules Bienaymé.[4]:98 The theorem was first stated without proof by Bienaymé in 1853[5] and later proved by Chebyshev in 1867.[6] His student Andrey Markov provided another proof in his 1884 Ph.D. thesis.[7]

Statement

Chebyshev's inequality is usually stated for random variables, but can be generalized to a statement about measure spaces.

Probabilistic statement

Let X (integrable) be a random variable with finite non-zero variance σ2 (and thus finite expected value μ).[8] Then for any real number k > 0,

[math]\displaystyle{ \Pr(|X-\mu|\geq k\sigma) \leq \frac{1}{k^2}. }[/math]

Only the case [math]\displaystyle{ k \gt 1 }[/math] is useful. When [math]\displaystyle{ k \leq 1 }[/math] the right-hand side [math]\displaystyle{ \frac{1}{k^2} \geq 1 }[/math] and the inequality is trivial as all probabilities are ≤ 1.

As an example, using [math]\displaystyle{ k = \sqrt{2} }[/math] shows that the probability that values lie outside the interval [math]\displaystyle{ (\mu - \sqrt{2}\sigma, \mu + \sqrt{2}\sigma) }[/math] does not exceed [math]\displaystyle{ \frac{1}{2} }[/math]. Equivalently, it implies that the probability of values lying within the interval (i.e. its "coverage") is at least [math]\displaystyle{ \frac{1}{2} }[/math].

Because it can be applied to completely arbitrary distributions provided they have a known finite mean and variance, the inequality generally gives a poor bound compared to what might be deduced if more aspects are known about the distribution involved.

k Min. % within k standard
deviations of mean
Max. % beyond k standard
deviations from mean
1 0% 100%
2 50% 50%
1.5 55.56% 44.44%
2 75% 25%
22 87.5% 12.5%
3 88.8889% 11.1111%
4 93.75% 6.25%
5 96% 4%
6 97.2222% 2.7778%
7 97.9592% 2.0408%
8 98.4375% 1.5625%
9 98.7654% 1.2346%
10 99% 1%

Measure-theoretic statement

Let (X, Σ, μ) be a measure space, and let f be an extended real-valued measurable function defined on X. Then for any real number t > 0 and 0 < p < ∞,[9]

[math]\displaystyle{ \mu(\{x\in X\,:\,\,|f(x)|\geq t\}) \leq {1\over t^p} \int_{|f| \geq t} |f|^p \, d\mu. }[/math]

More generally, if g is an extended real-valued measurable function, nonnegative and nondecreasing, with [math]\displaystyle{ g(t) \neq 0 }[/math] then:[citation needed]

[math]\displaystyle{ \mu(\{x\in X\,:\,\,f(x)\geq t\}) \leq {1\over g(t)} \int_X g\circ f\, d\mu. }[/math]

The previous statement then follows by defining [math]\displaystyle{ g(x) }[/math] as [math]\displaystyle{ |x|^p }[/math] if [math]\displaystyle{ x\ge t }[/math] and [math]\displaystyle{ 0 }[/math] otherwise.

Example

Suppose we randomly select a journal article from a source with an average of 1000 words per article, with a standard deviation of 200 words. We can then infer that the probability that it has between 600 and 1400 words (i.e. within k = 2 standard deviations of the mean) must be at least 75%, because there is no more than 1k2 = 1/4 chance to be outside that range, by Chebyshev's inequality. But if we additionally know that the distribution is normal, we can say there is a 75% chance the word count is between 770 and 1230 (which is an even tighter bound).

Sharpness of bounds

As shown in the example above, the theorem typically provides rather loose bounds. However, these bounds cannot in general (remaining true for arbitrary distributions) be improved upon. The bounds are sharp for the following example: for any k ≥ 1,

[math]\displaystyle{ X = \begin{cases} -1, & \text{with probability }\frac{1}{2k^2} \\ 0, & \text{with probability }1 - \frac{1}{k^2} \\ 1, & \text{with probability }\frac{1}{2k^2} \end{cases} }[/math]

For this distribution, the mean μ = 0 and the standard deviation σ = 1/k , so

[math]\displaystyle{ \Pr(|X-\mu| \ge k\sigma) = \Pr(|X| \ge 1) = \frac{1}{k^2}. }[/math]

Chebyshev's inequality is an equality for precisely those distributions that are a linear transformation of this example.

Proof

Markov's inequality states that for any real-valued random variable Y and any positive number a, we have Pr(|Y| ≥a) ≤ E(|Y|)/a. One way to prove Chebyshev's inequality is to apply Markov's inequality to the random variable Y = (Xμ)2 with a = ()2:

[math]\displaystyle{ \Pr(|X - \mu| \geq k\sigma) = \Pr((X - \mu)^2 \geq k^2\sigma^2) \leq \frac{\mathbb{E}[(X - \mu)^2]}{k^2\sigma^2} = \frac{\sigma^2}{k^2\sigma^2} = \frac{1}{k^2}. }[/math]

It can also be proved directly using conditional expectation:

[math]\displaystyle{ \begin{align} \sigma^2&=\mathbb{E}[(X-\mu)^2]\\[5pt] &=\mathbb{E}[(X-\mu)^2\mid k\sigma\leq |X-\mu|]\Pr[k\sigma\leq|X-\mu|]+\mathbb{E}[(X-\mu)^2\mid k\sigma\gt |X-\mu|]\Pr[k\sigma\gt |X-\mu|] \\[5pt] &\geq(k\sigma)^2\Pr[k\sigma\leq|X-\mu|]+0\cdot\Pr[k\sigma\gt |X-\mu|] \\[5pt] &=k^2\sigma^2\Pr[k\sigma\leq|X-\mu|] \end{align} }[/math]

Chebyshev's inequality then follows by dividing by k2σ2.

This proof also shows why the bounds are quite loose in typical cases: the conditional expectation on the event where |X − μ| <  is thrown away, and the lower bound of k2σ2 on the event |X − μ| ≥  can be quite poor.

Extensions

Several extensions of Chebyshev's inequality have been developed.

Selberg's inequality

Selberg derived a generalization to arbitrary intervals.[10] Suppose X is a random variable with mean μ and variance σ2. Selberg's inequality states[11] that if [math]\displaystyle{ \beta \geq \alpha \geq 0 }[/math],

[math]\displaystyle{ \Pr( X \in [\mu - \alpha, \mu + \beta] ) \ge \begin{cases}\frac{ \alpha^2 }{\alpha^2 + \sigma^2} &\text{if } \alpha(\beta-\alpha) \geq 2\sigma^2 \\ \frac{4\alpha\beta - 4\sigma^2}{(\alpha + \beta)^2} &\text{if } 2\alpha\beta \geq 2\sigma^2 \geq \alpha(\beta - \alpha) \\ 0 & \sigma^2 \geq \alpha\beta\end{cases} }[/math]

When [math]\displaystyle{ \alpha = \beta }[/math], this reduces to Chebyshev's inequality. These are known to be the best possible bounds.[12]

Finite-dimensional vector

Main page: Multidimensional Chebyshev's inequality

Chebyshev's inequality naturally extends to the multivariate setting, where one has n random variables Xi with mean μi and variance σi2. Then the following inequality holds.

[math]\displaystyle{ \Pr\left(\sum_{i=1}^n (X_i - \mu_i)^2 \ge k^2 \sum_{i=1}^n \sigma_i^2 \right) \le \frac{1}{k^2} }[/math]

This is known as the Birnbaum–Raymond–Zuckerman inequality after the authors who proved it for two dimensions.[13] This result can be rewritten in terms of vectors X = (X1, X2, ...) with mean μ = (μ1, μ2, ...), standard deviation σ = (σ1, σ2, ...), in the Euclidean norm || ⋅ ||.[14]

[math]\displaystyle{ \Pr(\| X - \mu \| \ge k \| \sigma \|) \le \frac{ 1 } { k^2 }. }[/math]

One can also get a similar infinite-dimensional Chebyshev's inequality. A second related inequality has also been derived by Chen.[15] Let n be the dimension of the stochastic vector X and let E(X) be the mean of X. Let S be the covariance matrix and k > 0. Then

[math]\displaystyle{ \Pr \left( ( X - \operatorname{E}(X) )^T S^{-1} (X - \operatorname{E}(X)) \lt k \right) \ge 1 - \frac{n}{k} }[/math]

where YT is the transpose of Y. The inequality can be written in terms of the Mahalanobis distance as

[math]\displaystyle{ \Pr \left( d^2_S(X,\operatorname{E}(X)) \lt k \right) \ge 1 - \frac{n}{k} }[/math]

where the Mahalanobis distance based on S is defined by

[math]\displaystyle{ d_S(x,y) =\sqrt{ (x -y)^T S^{-1} (x -y) } }[/math]

Navarro[16] proved that these bounds are sharp, that is, they are the best possible bounds for that regions when we just know the mean and the covariance matrix of X.

Stellato et al.[17] showed that this multivariate version of the Chebyshev inequality can be easily derived analytically as a special case of Vandenberghe et al.[18] where the bound is computed by solving a semidefinite program (SDP).

Known correlation

If the variables are independent this inequality can be sharpened.[19]

[math]\displaystyle{ \Pr\left (\bigcap_{i = 1}^n \frac{|X_i - \mu_i|}{\sigma_i} \le k_i \right ) \ge \prod_{i=1}^n \left (1 - \frac{1}{k_i^2} \right) }[/math]

Berge derived an inequality for two correlated variables X1, X2.[20] Let ρ be the correlation coefficient between X1 and X2 and let σi2 be the variance of Xi. Then

[math]\displaystyle{ \Pr\left( \bigcap_{ i = 1}^2 \left[ \frac{ | X_i - \mu_i | } { \sigma_i } \lt k \right] \right) \ge 1 - \frac{ 1 + \sqrt{ 1 - \rho^2 } } { k^2 }. }[/math]

This result can be sharpened to having different bounds for the two random variables[21] and having asymmetric bounds, as in Selberg's inequality. [22]

Olkin and Pratt derived an inequality for n correlated variables.[23]

[math]\displaystyle{ \Pr\left(\bigcap_{i = 1 }^n \frac{|X_i - \mu_i|}{\sigma_i} \lt k_i \right) \ge 1 - \frac{1}{n^2} \left(\sqrt{u} + \sqrt{n-1} \sqrt{n \sum_i \frac 1 { k_i^2} - u} \right)^2 }[/math]

where the sum is taken over the n variables and

[math]\displaystyle{ u = \sum_{i=1}^n \frac{1}{ k_i^2} + 2\sum_{i=1}^n \sum_{j\lt i} \frac{\rho_{ij}}{k_i k_j} }[/math]

where ρij is the correlation between Xi and Xj.

Olkin and Pratt's inequality was subsequently generalised by Godwin.[24]

Higher moments

Mitzenmacher and Upfal[25] note that by applying Markov's inequality to the nonnegative variable [math]\displaystyle{ | X - \operatorname{E}(X) |^n }[/math], one can get a family of tail bounds

[math]\displaystyle{ \Pr\left(| X - \operatorname{E}(X) | \ge k \operatorname{E}(|X - \operatorname{E}(X) |^n )^{ \frac{1}{n} }\right) \le \frac{1 } {k^n}, \qquad k \gt 0, n \geq 2. }[/math]

For n = 2 we obtain Chebyshev's inequality. For k ≥ 1, n > 4 and assuming that the nth moment exists, this bound is tighter than Chebyshev's inequality.[citation needed] This strategy, called the method of moments, is often used to prove tail bounds.

Exponential moment

A related inequality sometimes known as the exponential Chebyshev's inequality[26] is the inequality

[math]\displaystyle{ \Pr(X \ge \varepsilon) \le e^{ -t \varepsilon }\operatorname{E}\left (e^{ t X } \right), \qquad t \gt 0. }[/math]

Let K(t) be the cumulant generating function,

[math]\displaystyle{ K( t ) = \log \left(\operatorname{E}\left( e^{ t x } \right) \right). }[/math]

Taking the Legendre–Fenchel transformation[clarification needed] of K(t) and using the exponential Chebyshev's inequality we have

[math]\displaystyle{ -\log( \Pr (X \ge \varepsilon )) \ge \sup_t( t \varepsilon - K( t ) ). }[/math]

This inequality may be used to obtain exponential inequalities for unbounded variables.[27]

Bounded variables

If P(x) has finite support based on the interval [a, b], let M = max(|a|, |b|) where |x| is the absolute value of x. If the mean of P(x) is zero then for all k > 0[28]

[math]\displaystyle{ \frac{\operatorname{E}(|X|^r ) - k^r }{M^r} \le \Pr( | X | \ge k ) \le \frac{\operatorname{E}(| X |^r ) }{ k^r }. }[/math]

The second of these inequalities with r = 2 is the Chebyshev bound. The first provides a lower bound for the value of P(x).

Finite samples

Univariate case

Saw et al extended Chebyshev's inequality to cases where the population mean and variance are not known and may not exist, but the sample mean and sample standard deviation from N samples are to be employed to bound the expected value of a new drawing from the same distribution.[29] The following simpler version of this inequality is given by Kabán.[30]

[math]\displaystyle{ P( | X - m | \ge ks ) \le \frac 1 {N + 1} \left\lfloor \frac {N+1} N \left(\frac{N - 1}{k^2} + 1 \right) \right\rfloor }[/math]

where X is a random variable which we have sampled N times, m is the sample mean, k is a constant and s is the sample standard deviation.

This inequality holds even when the population moments do not exist, and when the sample is only weakly exchangeably distributed; this criterion is met for randomised sampling. A table of values for the Saw–Yang–Mo inequality for finite sample sizes (N < 100) has been determined by Konijn.[31] The table allows the calculation of various confidence intervals for the mean, based on multiples, C, of the standard error of the mean as calculated from the sample. For example, Konijn shows that for N = 59, the 95 percent confidence interval for the mean m is (mCs, m + Cs) where C = 4.447 × 1.006 = 4.47 (this is 2.28 times larger than the value found on the assumption of normality showing the loss on precision resulting from ignorance of the precise nature of the distribution).

An equivalent inequality can be derived in terms of the sample mean instead,[30]

[math]\displaystyle{ P( | X - m | \ge km ) \le \frac{N - 1} N \frac 1 {k^2} \frac{s^2}{m^2} + \frac 1 N. }[/math]

A table of values for the Saw–Yang–Mo inequality for finite sample sizes (N < 100) has been determined by Konijn.[31]

For fixed N and large m the Saw–Yang–Mo inequality is approximately[32]

[math]\displaystyle{ P( | X - m | \ge ks ) \le \frac 1 {N + 1}. }[/math]

Beasley et al have suggested a modification of this inequality[32]

[math]\displaystyle{ P( | X - m | \ge ks ) \le \frac 1 {k^2( N + 1 )}. }[/math]

In empirical testing this modification is conservative but appears to have low statistical power. Its theoretical basis currently remains unexplored.

Dependence on sample size

The bounds these inequalities give on a finite sample are less tight than those the Chebyshev inequality gives for a distribution. To illustrate this let the sample size N = 100 and let k = 3. Chebyshev's inequality states that at most approximately 11.11% of the distribution will lie at least three standard deviations away from the mean. Kabán's version of the inequality for a finite sample states that at most approximately 12.05% of the sample lies outside these limits. The dependence of the confidence intervals on sample size is further illustrated below.

For N = 10, the 95% confidence interval is approximately ±13.5789 standard deviations.

For N = 100 the 95% confidence interval is approximately ±4.9595 standard deviations; the 99% confidence interval is approximately ±140.0 standard deviations.

For N = 500 the 95% confidence interval is approximately ±4.5574 standard deviations; the 99% confidence interval is approximately ±11.1620 standard deviations.

For N = 1000 the 95% and 99% confidence intervals are approximately ±4.5141 and approximately ±10.5330 standard deviations respectively.

The Chebyshev inequality for the distribution gives 95% and 99% confidence intervals of approximately ±4.472 standard deviations and ±10 standard deviations respectively.

Samuelson's inequality

Main page: Samuelson's inequality

Although Chebyshev's inequality is the best possible bound for an arbitrary distribution, this is not necessarily true for finite samples. Samuelson's inequality states that all values of a sample must lie within N − 1 sample standard deviations of the mean.

By comparison, Chebyshev's inequality states that all but a 1/N fraction of the sample will lie within N standard deviations of the mean. Since there are N samples, this means that no samples will lie outside N standard deviations of the mean, which is worse than Samuelson's inequality. However, the benefit of Chebyshev's inequality is that it can be applied more generally to get confidence bounds for ranges of standard deviations that do not depend on the number of samples.

Semivariances

An alternative method of obtaining sharper bounds is through the use of semivariances (partial variances). The upper (σ+2) and lower (σ2) semivariances are defined as

[math]\displaystyle{ \sigma_+^2 = \frac { \sum_{x\gt m} (x - m)^2 } { n - 1 } , }[/math]
[math]\displaystyle{ \sigma_-^2 = \frac { \sum_{x\lt m} (m - x)^2 } { n - 1 }, }[/math]

where m is the arithmetic mean of the sample and n is the number of elements in the sample.

The variance of the sample is the sum of the two semivariances:

[math]\displaystyle{ \sigma^2 = \sigma_+^2 + \sigma_-^2. }[/math]

In terms of the lower semivariance Chebyshev's inequality can be written[33]

[math]\displaystyle{ \Pr(x \le m - a \sigma_-) \le \frac { 1 } { a^2 }. }[/math]

Putting

[math]\displaystyle{ a = \frac{ k \sigma } { \sigma_- }. }[/math]

Chebyshev's inequality can now be written

[math]\displaystyle{ \Pr(x \le m - k \sigma) \le \frac { 1 } { k^2 } \frac { \sigma_-^2 } { \sigma^2 }. }[/math]

A similar result can also be derived for the upper semivariance.

If we put

[math]\displaystyle{ \sigma_u^2 = \max(\sigma_-^2, \sigma_+^2) , }[/math]

Chebyshev's inequality can be written

[math]\displaystyle{ \Pr(| x \le m - k \sigma |) \le \frac 1 {k^2} \frac { \sigma_u^2 } { \sigma^2 } . }[/math]

Because σu2σ2, use of the semivariance sharpens the original inequality.

If the distribution is known to be symmetric, then

[math]\displaystyle{ \sigma_+^2 = \sigma_-^2 = \frac{ 1 } { 2 } \sigma^2 }[/math]

and

[math]\displaystyle{ \Pr(x \le m - k \sigma) \le \frac 1 {2k^2} . }[/math]

This result agrees with that derived using standardised variables.

Note
The inequality with the lower semivariance has been found to be of use in estimating downside risk in finance and agriculture.[33][34][35]

Multivariate case

Stellato et al.[17] simplified the notation and extended the empirical Chebyshev inequality from Saw et al.[29] to the multivariate case. Let [math]\displaystyle{ \xi \in \mathbb{R}^{n_\xi} }[/math] be a random variable and let [math]\displaystyle{ N \in \mathbb{Z}_{\geq n_\xi} }[/math]. We draw [math]\displaystyle{ N+1 }[/math] iid samples of [math]\displaystyle{ \xi }[/math] denoted as [math]\displaystyle{ \xi^{(1)},\dots,\xi^{(N)},\xi^{(N+1)} \in \mathbb{R}^{n_\xi} }[/math]. Based on the first [math]\displaystyle{ N }[/math] samples, we define the empirical mean as [math]\displaystyle{ \mu_N = \frac 1 N \sum_{i=1}^N \xi^{(i)} }[/math] and the unbiased empirical covariance as [math]\displaystyle{ \Sigma_N = \frac 1 N \sum_{i=1}^N (\xi^{(i)} - \mu_{N})(\xi^{(i)} - \mu_N)^\top }[/math]. If [math]\displaystyle{ \Sigma_N }[/math] is nonsingular, then for all [math]\displaystyle{ \lambda \in \mathbb{R}_{\geq 0} }[/math] then

[math]\displaystyle{ \begin{align} & P^{N+1} \left((\xi^{(N+1)} - \mu_N)^\top \Sigma_N^{-1}(\xi^{(N+1)} - \mu_N) \geq \lambda^2\right) \\[8pt] \leq {} & \min\left\{1, \frac 1 {N+1} \left\lfloor \frac{n_\xi(N+1)(N^2 - 1 + N\lambda^2)}{N^2\lambda^2}\right\rfloor\right\}. \end{align} }[/math]

Remarks

In the univariate case, i.e. [math]\displaystyle{ n_\xi = 1 }[/math], this inequality corresponds to the one from Saw et al.[29] Moreover, the right-hand side can be simplified by upper bounding the floor function by its argument

[math]\displaystyle{ P^{N+1}\left((\xi^{(N+1)} - \mu_N)^\top \Sigma_N^{-1}(\xi^{(N+1)} - \mu_N) \geq \lambda^2\right) \leq \min\left\{1, \frac{n_\xi(N^2 - 1 + N\lambda^2)}{N^2\lambda^2}\right\}. }[/math]

As [math]\displaystyle{ N \to \infty }[/math], the right-hand side tends to [math]\displaystyle{ \min \left\{1, \frac{n_\xi}{\lambda^2}\right\} }[/math] which corresponds to the multivariate Chebyshev inequality over ellipsoids shaped according to [math]\displaystyle{ \Sigma }[/math] and centered in [math]\displaystyle{ \mu }[/math].

Sharpened bounds

Chebyshev's inequality is important because of its applicability to any distribution. As a result of its generality it may not (and usually does not) provide as sharp a bound as alternative methods that can be used if the distribution of the random variable is known. To improve the sharpness of the bounds provided by Chebyshev's inequality a number of methods have been developed; for a review see eg.[11][36]

Cantelli's inequality

Cantelli's inequality[37] due to Francesco Paolo Cantelli states that for a real random variable (X) with mean (μ) and variance (σ2)

[math]\displaystyle{ P(X - \mu \ge a) \le \frac{\sigma^2}{ \sigma^2 + a^2 } }[/math]

where a ≥ 0.

This inequality can be used to prove a one tailed variant of Chebyshev's inequality with k > 0[38]

[math]\displaystyle{ \Pr(X - \mu \geq k \sigma) \leq \frac{ 1 }{ 1 + k^2 }. }[/math]

The bound on the one tailed variant is known to be sharp. To see this consider the random variable X that takes the values

[math]\displaystyle{ X = 1 }[/math] with probability [math]\displaystyle{ \frac{ \sigma^2 } { 1 + \sigma^2 } }[/math]
[math]\displaystyle{ X = - \sigma^2 }[/math] with probability [math]\displaystyle{ \frac{ 1 } { 1 + \sigma^2 }. }[/math]

Then E(X) = 0 and E(X2) = σ2 and P(X < 1) = 1 / (1 + σ2).

An application: distance between the mean and the median

The one-sided variant can be used to prove the proposition that for probability distributions having an expected value and a median, the mean and the median can never differ from each other by more than one standard deviation. To express this in symbols let μ, ν, and σ be respectively the mean, the median, and the standard deviation. Then

[math]\displaystyle{ \left | \mu - \nu \right | \leq \sigma. }[/math]

There is no need to assume that the variance is finite because this inequality is trivially true if the variance is infinite.

The proof is as follows. Setting k = 1 in the statement for the one-sided inequality gives:

[math]\displaystyle{ \Pr(X - \mu \geq \sigma) \leq \frac{ 1 }{ 2 } \implies \Pr(X \geq \mu + \sigma) \leq \frac{ 1 }{ 2 }. }[/math]

Changing the sign of X and of μ, we get

[math]\displaystyle{ \Pr(X \leq \mu - \sigma) \leq \frac{ 1 }{ 2 }. }[/math]

As the median is by definition any real number m that satisfies the inequalities

[math]\displaystyle{ \operatorname{P}(X\leq m) \geq \frac{1}{2}\text{ and }\operatorname{P}(X\geq m) \geq \frac{1}{2} }[/math]

this implies that the median lies within one standard deviation of the mean. A proof using Jensen's inequality also exists.

Bhattacharyya's inequality

Bhattacharyya[39] extended Cantelli's inequality using the third and fourth moments of the distribution.

Let [math]\displaystyle{ \mu = 0 }[/math] and [math]\displaystyle{ \sigma^2 }[/math] be the variance. Let [math]\displaystyle{ \gamma = E[X^3] / \sigma^3 }[/math] and [math]\displaystyle{ \kappa = E[X^4]/\sigma^4 }[/math].

If [math]\displaystyle{ k^2 - k \gamma - 1 \gt 0 }[/math] then

[math]\displaystyle{ P(X \gt k\sigma) \le \frac{ \kappa - \gamma^2 - 1 }{ (\kappa - \gamma^2 - 1) (1 + k^2) + (k^2 - k\gamma - 1) }. }[/math]

The necessity of [math]\displaystyle{ k^2 - k \gamma - 1 \gt 0 }[/math] may require [math]\displaystyle{ k }[/math] to be reasonably large.

In the case [math]\displaystyle{ E[X^3]=0 }[/math] this simplifies to

[math]\displaystyle{ P(X \gt k\sigma) \le \frac{\kappa-1}{\kappa \left(k^2+1\right)-2} \quad \text{for } k \gt 1. }[/math]

Since [math]\displaystyle{ \frac{\kappa-1}{\kappa \left(k^2+1\right)-2} = \frac{1}{2}-\frac{\kappa (k-1)}{2 (\kappa-1)}+O\left((k-1)^2\right) }[/math] for [math]\displaystyle{ k }[/math] close to 1, this bound improves slightly over Cantelli's bound [math]\displaystyle{ \frac{1}{2}-\frac{k-1}{2}+O\left((k-1)^2\right) }[/math] as [math]\displaystyle{ \kappa \gt 1 }[/math].

wins a factor 2 over Chebyshev's inequality.

Gauss's inequality

Main page: Gauss's inequality

In 1823 Gauss showed that for a distribution with a unique mode at zero,[40]

[math]\displaystyle{ P( | X | \ge k ) \le \frac{ 4 \operatorname{ E }( X^2 ) } { 9k^2 } \quad\text{if} \quad k^2 \ge \frac{ 4 } { 3 } \operatorname{E} (X^2) , }[/math]
[math]\displaystyle{ P( | X | \ge k ) \le 1 - \frac{ k } { \sqrt{3} \operatorname{ E }( X^2 ) } \quad \text{if} \quad k^2 \le \frac{ 4 } { 3 } \operatorname{ E }( X^2 ). }[/math]

Vysochanskij–Petunin inequality

Main page: Vysochanskij–Petunin inequality

The Vysochanskij–Petunin inequality generalizes Gauss's inequality, which only holds for deviation from the mode of a unimodal distribution, to deviation from the mean, or more generally, any center.[41] If X is a unimodal distribution with mean μ and variance σ2, then the inequality states that

[math]\displaystyle{ P( | X - \mu | \ge k \sigma ) \le \frac{ 4 }{ 9k^2 } \quad \text{if} \quad k \ge \sqrt{8/3} = 1.633. }[/math]
[math]\displaystyle{ P( | X - \mu | \ge k \sigma ) \le \frac{ 4 }{ 3k^2 } - \frac13 \quad \text{if} \quad k \le \sqrt{8/3}. }[/math]

For symmetrical unimodal distributions, the median and the mode are equal, so both the Vysochanskij–Petunin inequality and Gauss's inequality apply to the same center. Further, for symmetrical distributions, one-sided bounds can be obtained by noticing that

[math]\displaystyle{ P( X - \mu \ge k \sigma ) = P( X - \mu \le -k \sigma ) = \frac{1}{2} P( |X - \mu| \ge k \sigma ). }[/math]

The additional fraction of [math]\displaystyle{ 4/9 }[/math] present in these tail bounds lead to better confidence intervals than Chebyshev's inequality. For example, for any symmetrical unimodal distribution, the Vysochanskij–Petunin inequality states that 4/(9 x 3^2) = 4/81 ≈ 4.9% of the distribution lies outside 3 standard deviations of the mode.

Bounds for specific distributions

DasGupta has shown that if the distribution is known to be normal[42]

[math]\displaystyle{ P( | X - \mu | \ge k \sigma ) \le \frac{ 1 }{ 3 k^2 } . }[/math]

From DasGupta's inequality it follows that for a normal distribution at least 95% lies within approximately 2.582 standard deviations of the mean. This is less sharp than the true figure (approximately 1.96 standard deviations of the mean).

  • DasGupta has determined a set of best possible bounds for a normal distribution for this inequality.[42]
  • Steliga and Szynal have extended these bounds to the Pareto distribution.[43]
  • Grechuk et al. developed a general method for deriving the best possible bounds in Chebyshev's inequality for any family of distributions, and any deviation risk measure in place of standard deviation. In particular, they derived Chebyshev inequality for distributions with log-concave densities.[44]

Related inequalities

Several other related inequalities are also known.

Paley–Zygmund inequality

Main page: Paley–Zygmund inequality

The Paley–Zygmund inequality gives a lower bound on tail probabilities, as opposed to Chebyshev's inequality which gives an upper bound.[45] Applying it to the square of a random variable, we get

[math]\displaystyle{ \Pr( | Z | \gt \theta \sqrt{E[Z^2]} ) \ge \frac{ ( 1 - \theta^2 )^2 E[Z^2]^2 }{E[Z^4]}. }[/math]

Haldane's transformation

One use of Chebyshev's inequality in applications is to create confidence intervals for variates with an unknown distribution. Haldane noted,[46] using an equation derived by Kendall,[47] that if a variate (x) has a zero mean, unit variance and both finite skewness (γ) and kurtosis (κ) then the variate can be converted to a normally distributed standard score (z):

[math]\displaystyle{ z = x - \frac{\gamma}{6} (x^2 - 1) + \frac{ x }{ 72 } [ 2 \gamma^2 (4 x^2 - 7) - 3 \kappa (x^2 - 3) ] + \cdots }[/math]

This transformation may be useful as an alternative to Chebyshev's inequality or as an adjunct to it for deriving confidence intervals for variates with unknown distributions.

While this transformation may be useful for moderately skewed and/or kurtotic distributions, it performs poorly when the distribution is markedly skewed and/or kurtotic.

He, Zhang and Zhang's inequality

For any collection of n non-negative independent random variables Xi with expectation 1 [48]

[math]\displaystyle{ \Pr\left ( \frac{\sum_{i=1}^n X_i }{n} - 1 \ge \frac{1}{n} \right) \le \frac{ 7 }{ 8 }. }[/math]

Integral Chebyshev inequality

There is a second (less well known) inequality also named after Chebyshev[49]

If f, g : [a, b] → R are two monotonic functions of the same monotonicity, then

[math]\displaystyle{ \frac{ 1 }{ b - a } \int_a^b \! f(x) g(x) \,dx \ge \left[ \frac{ 1 }{ b - a } \int_a^b \! f(x) \,dx \right] \left[ \frac{ 1 }{ b - a } \int_a^b \! g(x) \,dx \right] . }[/math]

If f and g are of opposite monotonicity, then the above inequality works in the reverse way.

This inequality is related to Jensen's inequality,[50] Kantorovich's inequality,[51] the Hermite–Hadamard inequality[51] and Walter's conjecture.[52]

Other inequalities

There are also a number of other inequalities associated with Chebyshev:

Notes

The Environmental Protection Agency has suggested best practices for the use of Chebyshev's inequality for estimating confidence intervals.[53]

See also

References

  1. Kvanli, Alan H.; Pavur, Robert J.; Keeling, Kellie B. (2006). Concise Managerial Statistics. cEngage Learning. pp. 81–82. ISBN 978-0-324-22388-0. https://books.google.com/books?id=h6CQ1J0gwNgC&pg=PT95. 
  2. Chernick, Michael R. (2011). The Essentials of Biostatistics for Physicians, Nurses, and Clinicians. John Wiley & Sons. pp. 49–50. ISBN 978-0-470-64185-9. https://books.google.com/books?id=JP4azqd8ONEC&pg=PA50. 
  3. "Error Term of Chebyshev inequality?" (in en). https://math.stackexchange.com/a/776424/352034. 
  4. Knuth, Donald (1997). The Art of Computer Programming: Fundamental Algorithms, Volume 1 (3rd ed.). Reading, Massachusetts: Addison–Wesley. ISBN 978-0-201-89683-1. http://www-cs-faculty.stanford.edu/~uno/taocp.html. Retrieved 1 October 2012. 
  5. Bienaymé, I.-J. (1853). "Considérations àl'appui de la découverte de Laplace". Comptes Rendus de l'Académie des Sciences 37: 309–324. 
  6. Tchebichef, P. (1867). "Des valeurs moyennes". Journal de Mathématiques Pures et Appliquées. 2 12: 177–184. 
  7. Markov A. (1884) On certain applications of algebraic continued fractions, Ph.D. thesis, St. Petersburg
  8. Feller, W., 1968. An introduction to probability theory and its applications, vol. 1. p227 (Wiley, New York).
  9. Grafakos, Lukas (2004). Classical and Modern Fourier Analysis. Pearson Education Inc.. p. 5. 
  10. Selberg, Henrik L. (1940). "Zwei Ungleichungen zur Ergänzung des Tchebycheffschen Lemmas" (in de). Skandinavisk Aktuarietidskrift (Scandinavian Actuarial Journal) 1940 (3–4): 121–125. doi:10.1080/03461238.1940.10404804. ISSN 0346-1238. OCLC 610399869. 
  11. 11.0 11.1 Godwin, H. J. (September 1955). "On Generalizations of Tchebychef's Inequality" (in en). Journal of the American Statistical Association 50 (271): 923–945. doi:10.1080/01621459.1955.10501978. ISSN 0162-1459. http://www.tandfonline.com/doi/abs/10.1080/01621459.1955.10501978. 
  12. Conlon, J.; Dulá, J. H.. A geometric derivation and interpretation of Tchebyscheff's Inequality. http://www.people.vcu.edu/~jdula/WORKINGPAPERS/tcheby.pdf. Retrieved 2 October 2012. 
  13. Birnbaum, Z. W.; Raymond, J.; Zuckerman, H. S. (1947). "A Generalization of Tshebyshev's Inequality to Two Dimensions". The Annals of Mathematical Statistics 18 (1): 70–79. doi:10.1214/aoms/1177730493. ISSN 0003-4851. http://projecteuclid.org/DPubS?service=UI&version=1.0&verb=Display&handle=euclid.aoms/1177730493. Retrieved 7 October 2012. 
  14. Ferentinos, K (1982). "On Tchebycheff type inequalities". Trabajos Estadıst Investigacion Oper 33: 125–132. doi:10.1007/BF02888707. 
  15. Xinjia Chen (2007). "A New Generalization of Chebyshev Inequality for Random Vectors". arXiv:0707.0805v2 [math.ST].
  16. Jorge Navarro (2014). "Can the bounds in the multivariate Chebyshev inequality be attained?". Statistics and Probability Letters 91: 1–5. doi:10.1016/j.spl.2014.03.028. 
  17. 17.0 17.1 Stellato, Bartolomeo; Parys, Bart P. G. Van; Goulart, Paul J. (2016-05-31). "Multivariate Chebyshev Inequality with Estimated Mean and Variance". The American Statistician 71 (2): 123–127. doi:10.1080/00031305.2016.1186559. ISSN 0003-1305. 
  18. Vandenberghe, L.; Boyd, S.; Comanor, K. (2007-01-01). "Generalized Chebyshev Bounds via Semidefinite Programming". SIAM Review 49 (1): 52–64. doi:10.1137/S0036144504440543. ISSN 0036-1445. Bibcode2007SIAMR..49...52V. 
  19. Kotz, Samuel; Balakrishnan, N.; Johnson, Norman L. (2000). Continuous Multivariate Distributions, Volume 1, Models and Applications (2nd ed.). Boston [u.a.]: Houghton Mifflin. ISBN 978-0-471-18387-7. http://www.wiley.com/WileyCDA/WileyTitle/productCd-0471183873.html. Retrieved 7 October 2012. 
  20. Berge, P. O. (1938). "A note on a form of Tchebycheff's theorem for two variables". Biometrika 29 (3/4): 405–406. doi:10.2307/2332015. 
  21. Lal D. N. (1955) A note on a form of Tchebycheff's inequality for two or more variables. Sankhya 15(3):317–320
  22. Isii K. (1959) On a method for generalizations of Tchebycheff's inequality. Ann Inst Stat Math 10: 65–88
  23. Olkin, Ingram; Pratt, John W. (1958). "A Multivariate Tchebycheff Inequality". The Annals of Mathematical Statistics 29 (1): 226–234. doi:10.1214/aoms/1177706720. 
  24. Godwin H. J. (1964) Inequalities on distribution functions. New York, Hafner Pub. Co.
  25. Mitzenmacher, Michael; Upfal, Eli (January 2005). Probability and Computing: Randomized Algorithms and Probabilistic Analysis (Repr. ed.). Cambridge [u.a.]: Cambridge Univ. Press. ISBN 978-0-521-83540-4. http://www.cambridge.org/us/knowledge/isbn/item1171566/?site_locale=en_US. Retrieved 6 October 2012. 
  26. Section 2.1
  27. Baranoski, Gladimir V. G.; Rokne, Jon G.; Xu, Guangwu (15 May 2001). "Applying the exponential Chebyshev inequality to the nondeterministic computation of form factors". Journal of Quantitative Spectroscopy and Radiative Transfer 69 (4): 199–200. doi:10.1016/S0022-4073(00)00095-9. Bibcode2001JQSRT..69..447B.  (the references for this article are corrected by Baranoski, Gladimir V. G.; Rokne, Jon G.; Guangwu Xu (15 January 2002). "Corrigendum to: 'Applying the exponential Chebyshev inequality to the nondeterministic computation of form factors'". Journal of Quantitative Spectroscopy and Radiative Transfer 72 (2): 199–200. doi:10.1016/S0022-4073(01)00171-6. Bibcode2002JQSRT..72..199B. )
  28. Dufour (2003) Properties of moments of random variables
  29. 29.0 29.1 29.2 Saw, John G.; Yang, Mark C. K.; Mo, Tse Chin (1984). "Chebyshev Inequality with Estimated Mean and Variance". The American Statistician 38 (2): 130–2. doi:10.2307/2683249. ISSN 0003-1305. 
  30. 30.0 30.1 Kabán, Ata (2012). "Non-parametric detection of meaningless distances in high dimensional data". Statistics and Computing 22 (2): 375–85. doi:10.1007/s11222-011-9229-0. 
  31. 31.0 31.1 Konijn, Hendrik S. (February 1987). "Distribution-Free and Other Prediction Intervals". The American Statistician 41 (1): 11–15. doi:10.2307/2684311. 
  32. 32.0 32.1 Beasley, T. Mark; Page, Grier P.; Brand, Jaap P. L.; Gadbury, Gary L.; Mountz, John D.; Allison, David B. (January 2004). "Chebyshev's inequality for nonparametric testing with small N and α in microarray research". Journal of the Royal Statistical Society. C (Applied Statistics) 53 (1): 95–108. doi:10.1111/j.1467-9876.2004.00428.x. ISSN 1467-9876. 
  33. 33.0 33.1 Berck, Peter; Hihn, Jairus M. (May 1982). "Using the Semivariance to Estimate Safety-First Rules". American Journal of Agricultural Economics 64 (2): 298–300. doi:10.2307/1241139. ISSN 0002-9092. 
  34. Nantell, Timothy J.; Price, Barbara (June 1979). "An Analytical Comparison of Variance and Semivariance Capital Market Theories". The Journal of Financial and Quantitative Analysis 14 (2): 221–42. doi:10.2307/2330500. 
  35. Neave, Edwin H.; Ross, Michael N.; Yang, Jun (2009). "Distinguishing upside potential from downside risk". Management Research News 32 (1): 26–36. doi:10.1108/01409170910922005. ISSN 0140-9174. 
  36. Savage, I. Richard. "Probability inequalities of the Tchebycheff type." Journal of Research of the National Bureau of Standards-B. Mathematics and Mathematical Physics B 65 (1961): 211-222
  37. Cantelli F. (1910) Intorno ad un teorema fondamentale della teoria del rischio. Bolletino dell Associazione degli Attuari Italiani
  38. Grimmett and Stirzaker, problem 7.11.9. Several proofs of this result can be found in Chebyshev's Inequalities by A. G. McDowell.
  39. Bhattacharyya, B. B. (1987). "One-sided chebyshev inequality when the first four moments are known". Communications in Statistics – Theory and Methods 16 (9): 2789–91. doi:10.1080/03610928708829540. ISSN 0361-0926. 
  40. Gauss C. F. Theoria Combinationis Observationum Erroribus Minimis Obnoxiae. Pars Prior. Pars Posterior. Supplementum. Theory of the Combination of Observations Least Subject to Errors. Part One. Part Two. Supplement. 1995. Translated by G. W. Stewart. Classics in Applied Mathematics Series, Society for Industrial and Applied Mathematics, Philadelphia
  41. Pukelsheim, Friedrich (May 1994). "The Three Sigma Rule" (in en). The American Statistician 48 (2): 88–91. doi:10.1080/00031305.1994.10476030. ISSN 0003-1305. http://www.tandfonline.com/doi/abs/10.1080/00031305.1994.10476030. 
  42. 42.0 42.1 DasGupta, A (2000). "Best constants in Chebychev inequalities with various applications". Metrika 5 (1): 185–200. doi:10.1007/s184-000-8316-9. 
  43. Steliga, Katarzyna; Szynal, Dominik (2010). "On Markov-Type Inequalities". International Journal of Pure and Applied Mathematics 58 (2): 137–152. ISSN 1311-8080. http://ijpam.eu/contents/2010-58-2/2/2.pdf. Retrieved 10 October 2012. 
  44. Grechuk, B., Molyboha, A., Zabarankin, M. (2010). Chebyshev Inequalities with Law Invariant Deviation Measures, Probability in the Engineering and Informational Sciences, 24(1), 145-170.
  45. Godwin H. J. (1964) Inequalities on distribution functions. (Chapter 3) New York, Hafner Pub. Co.
  46. Haldane, J. B. (1952). "Simple tests for bimodality and bitangentiality". Annals of Eugenics 16 (4): 359–364. doi:10.1111/j.1469-1809.1951.tb02488.x. PMID 14953132. 
  47. Kendall M. G. (1943) The Advanced Theory of Statistics, 1. London
  48. He, Simai; Zhang, Jiawei; Zhang, Shuzhong (2010). "Bounding probability of small deviation: a fourth moment approach". Mathematics of Operations Research 35 (1): 208–232. doi:10.1287/moor.1090.0438. 
  49. Fink, A. M.; Jodeit, Max Jr. (1984). "On Chebyshev's other inequality". in Tong, Y. L.; Gupta, Shanti S.. Inequalities in Statistics and Probability. Institute of Mathematical Statistics Lecture Notes - Monograph Series. 5. pp. 115–120. doi:10.1214/lnms/1215465637. ISBN 978-0-940600-04-1. http://projecteuclid.org/euclid.lnms/1215465617. Retrieved 7 October 2012. 
  50. Niculescu, Constantin P. (2001). "An extension of Chebyshev's inequality and its connection with Jensen's inequality". Journal of Inequalities and Applications 6 (4): 451–462. doi:10.1155/S1025583401000273. ISSN 1025-5834. http://emis.matem.unam.mx/journals/HOA/JIA/Volume6_4/462.html. Retrieved 6 October 2012. 
  51. 51.0 51.1 Niculescu, Constantin P.; Pečarić, Josip (2010). "The Equivalence of Chebyshev's Inequality to the Hermite–Hadamard Inequality". Mathematical Reports 12 (62): 145–156. ISSN 1582-3067. http://www.csm.ro/reviste/Mathematical_Reports/Pdfs/2010/2/Niculescu.pdf. Retrieved 6 October 2012. 
  52. Malamud, S. M. (15 February 2001). "Some complements to the Jensen and Chebyshev inequalities and a problem of W. Walter". Proceedings of the American Mathematical Society 129 (9): 2671–2678. doi:10.1090/S0002-9939-01-05849-X. ISSN 0002-9939. https://www.ams.org/journals/proc/2001-129-09/S0002-9939-01-05849-X/. Retrieved 7 October 2012. 
  53. Calculating Upper Confidence Limits for Exposure Point Concentrations at hazardous Waste Sites (Report). Office of Emergency and Remedial Response of the U.S. Environmental Protection Agency. December 2002. http://nepis.epa.gov/Exe/ZyNET.exe/P100CYCE.TXT?ZyActionD=ZyDocument&Client=EPA&Index=2000+Thru+2005&Docs=&Query=&Time=&EndTime=&SearchMethod=1&TocRestrict=n&Toc=&TocEntry=&QField=&QFieldYear=&QFieldMonth=&QFieldDay=&IntQFieldOp=0&ExtQFieldOp=0&XmlQuery=&File=D%3A%5Czyfiles%5CIndex%20Data%5C00thru05%5CTxt%5C00000029%5CP100CYCE.txt&User=ANONYMOUS&Password=anonymous&SortMethod=h%7C-&MaximumDocuments=1&FuzzyDegree=0&ImageQuality=r75g8/r75g8/x150y150g16/i425&Display=p%7Cf&DefSeekPage=x&SearchBack=ZyActionL&Back=ZyActionS&BackDesc=Results%20page&MaximumPages=1&ZyEntry=1&SeekPage=x&ZyPURL#. Retrieved 5 August 2016. 

Further reading

  • A. Papoulis (1991), Probability, Random Variables, and Stochastic Processes, 3rd ed. McGraw–Hill. ISBN 0-07-100870-5. pp. 113–114.
  • G. Grimmett and D. Stirzaker (2001), Probability and Random Processes, 3rd ed. Oxford. ISBN 0-19-857222-0. Section 7.3.

External links