Unimodality

From HandWiki
(Redirected from Strongly unimodal)
Short description: Property of having a unique mode or maximum value

In mathematics, unimodality means possessing a unique mode. More generally, unimodality means there is only a single highest value, somehow defined, of some mathematical object.[1]

Unimodal probability distribution

File:Normal distribution pdf.svg

Figure 2. A simple bimodal distribution.
Figure 3. A bimodal distribution. Note that only the largest peak would correspond to a mode in the strict sense of the definition of mode

In statistics, a unimodal probability distribution or unimodal distribution is a probability distribution which has a single peak. The term "mode" in this context refers to any peak of the distribution, not just to the strict definition of mode which is usual in statistics.

If there is a single mode, the distribution function is called "unimodal". If it has more modes it is "bimodal" (2), "trimodal" (3), etc., or in general, "multimodal".[2] Figure 1 illustrates normal distributions, which are unimodal. Other examples of unimodal distributions include Cauchy distribution, Student's t-distribution, chi-squared distribution and exponential distribution. Among discrete distributions, the binomial distribution and Poisson distribution can be seen as unimodal, though for some parameters they can have two adjacent values with the same probability.

Figure 2 and Figure 3 illustrate bimodal distributions.

Other definitions

Other definitions of unimodality in distribution functions also exist.

In continuous distributions, unimodality can be defined through the behavior of the cumulative distribution function (cdf).[3] If the cdf is convex for x < m and concave for x > m, then the distribution is unimodal, m being the mode. Note that under this definition the uniform distribution is unimodal,[4] as well as any other distribution in which the maximum distribution is achieved for a range of values, e.g. trapezoidal distribution. Usually this definition allows for a discontinuity at the mode; usually in a continuous distribution the probability of any single value is zero, while this definition allows for a non-zero probability, or an "atom of probability", at the mode.

Criteria for unimodality can also be defined through the characteristic function of the distribution[3] or through its Laplace–Stieltjes transform.[5]

Another way to define a unimodal discrete distribution is by the occurrence of sign changes in the sequence of differences of the probabilities.[6] A discrete distribution with a probability mass function, [math]\displaystyle{ \{p_n : n = \dots, -1, 0, 1, \dots\} }[/math], is called unimodal if the sequence [math]\displaystyle{ \dots, p_{-2} - p_{-1}, p_{-1} - p_0, p_0 - p_1, p_1 - p_2, \dots }[/math] has exactly one sign change (when zeroes don't count).

Uses and results

One reason for the importance of distribution unimodality is that it allows for several important results. Several inequalities are given below which are only valid for unimodal distributions. Thus, it is important to assess whether or not a given data set comes from a unimodal distribution. Several tests for unimodality are given in the article on multimodal distribution.

Inequalities

Gauss's inequality

A first important result is Gauss's inequality.[7] Gauss's inequality gives an upper bound on the probability that a value lies more than any given distance from its mode. This inequality depends on unimodality.

Vysochanskiï–Petunin inequality

A second is the Vysochanskiï–Petunin inequality,[8] a refinement of the Chebyshev inequality. The Chebyshev inequality guarantees that in any probability distribution, "nearly all" the values are "close to" the mean value. The Vysochanskiï–Petunin inequality refines this to even nearer values, provided that the distribution function is continuous and unimodal. Further results were shown by Sellke and Sellke.[9]

Mode, median and mean

Gauss also showed in 1823 that for a unimodal distribution[10]

[math]\displaystyle{ \sigma \le \omega \le 2 \sigma }[/math]

and

[math]\displaystyle{ |\nu - \mu| \le \sqrt{\frac{3}{4}} \omega , }[/math]

where the median is ν, the mean is μ and ω is the root mean square deviation from the mode.

It can be shown for a unimodal distribution that the median ν and the mean μ lie within (3/5)1/2 ≈ 0.7746 standard deviations of each other.[11] In symbols,

[math]\displaystyle{ \frac{|\nu - \mu|}{\sigma} \le \sqrt{\frac{3}{5}} }[/math]

where | . | is the absolute value.

In 2020, Bernard, Kazzi, and Vanduffel generalized the previous inequality by deriving the maximum distance between the symmetric quantile average [math]\displaystyle{ \frac{ q_\alpha + q_{(1-\alpha)} }{ 2 } }[/math] and the mean,[12]

[math]\displaystyle{ \frac{ \left| \frac{ q_\alpha + q_{(1-\alpha)} }{2} - \mu \right| }{ \sigma } \le \left\{ \begin{array}{cl} \frac{\sqrt[]{\frac{4}{9(1-\alpha)}-1} \text{ } + \text{ } \sqrt[]{\frac{1-\alpha}{1/3+\alpha}}}{2} & \text{for }\alpha \in \left[\frac{5}{6},1\right)\!, \\ \frac{\sqrt[]{\frac{3 \alpha}{4-3\alpha}} \text{ } + \text{ } \sqrt[]{\frac{1-\alpha}{1/3+\alpha}}}{2} & \text{for }\alpha \in \left(\frac{1}{6},\frac{5}{6}\right)\!,\\ \frac{\sqrt[]{\frac{3 \alpha}{4-3\alpha}} \text{ } + \text{ } \sqrt[]{\frac{4}{9 \alpha} -1}}{2} & \text{for }\alpha \in \left(0,\frac{1}{6}\right]\!. \end{array} \right. }[/math]

It is worth noting that the maximum distance is minimized at [math]\displaystyle{ \alpha=0.5 }[/math] (i.e., when the symmetric quantile average is equal to [math]\displaystyle{ q_{0.5} = \nu }[/math]), which indeed motivates the common choice of the median as a robust estimator for the mean. Moreover, when [math]\displaystyle{ \alpha = 0.5 }[/math], the bound is equal to [math]\displaystyle{ \sqrt{3/5} }[/math], which is the maximum distance between the median and the mean of a unimodal distribution.

A similar relation holds between the median and the mode θ: they lie within 31/2 ≈ 1.732 standard deviations of each other:

[math]\displaystyle{ \frac{|\nu - \theta|}{\sigma} \le \sqrt{3}. }[/math]

It can also be shown that the mean and the mode lie within 31/2 of each other:

[math]\displaystyle{ \frac{|\mu - \theta|}{\sigma} \le \sqrt{3}. }[/math]

Skewness and kurtosis

Rohatgi and Szekely claimed that the skewness and kurtosis of a unimodal distribution are related by the inequality:[13]

[math]\displaystyle{ \gamma^2 - \kappa \le \frac{ 6 }{ 5 } = 1.2 }[/math]

where κ is the kurtosis and γ is the skewness. Klaassen, Mokveld, and van Es showed that this only applies in certain settings, such as the set of unimodal distributions where the mode and mean coincide.[14]

They derived a weaker inequality which applies to all unimodal distributions:[14]

[math]\displaystyle{ \gamma^2 - \kappa \le \frac{ 186 }{ 125 } = 1.488 }[/math]

This bound is sharp, as it is reached by the equal-weights mixture of the uniform distribution on [0,1] and the discrete distribution at {0}.

Unimodal function

As the term "modal" applies to data sets and probability distribution, and not in general to functions, the definitions above do not apply. The definition of "unimodal" was extended to functions of real numbers as well.

A common definition is as follows: a function f(x) is a unimodal function if for some value m, it is monotonically increasing for x ≤ m and monotonically decreasing for x ≥ m. In that case, the maximum value of f(x) is f(m) and there are no other local maxima.

Proving unimodality is often hard. One way consists in using the definition of that property, but it turns out to be suitable for simple functions only. A general method based on derivatives exists,[15] but it does not succeed for every function despite its simplicity.

Examples of unimodal functions include quadratic polynomial functions with a negative quadratic coefficient, tent map functions, and more.

The above is sometimes related to as strong unimodality, from the fact that the monotonicity implied is strong monotonicity. A function f(x) is a weakly unimodal function if there exists a value m for which it is weakly monotonically increasing for x ≤ m and weakly monotonically decreasing for x ≥ m. In that case, the maximum value f(m) can be reached for a continuous range of values of x. An example of a weakly unimodal function which is not strongly unimodal is every other row in Pascal's triangle.

Depending on context, unimodal function may also refer to a function that has only one local minimum, rather than maximum.[16] For example, local unimodal sampling, a method for doing numerical optimization, is often demonstrated with such a function. It can be said that a unimodal function under this extension is a function with a single local extremum.

One important property of unimodal functions is that the extremum can be found using search algorithms such as golden section search, ternary search or successive parabolic interpolation.[17]

Other extensions

A function f(x) is "S-unimodal" (often referred to as "S-unimodal map") if its Schwarzian derivative is negative for all [math]\displaystyle{ x \ne c }[/math], where [math]\displaystyle{ c }[/math] is the critical point.[18]

In computational geometry if a function is unimodal it permits the design of efficient algorithms for finding the extrema of the function.[19]

A more general definition, applicable to a function f(X) of a vector variable X is that f is unimodal if there is a one-to-one differentiable mapping X = G(Z) such that f(G(Z)) is convex. Usually one would want G(Z) to be continuously differentiable with nonsingular Jacobian matrix.

Quasiconvex functions and quasiconcave functions extend the concept of unimodality to functions whose arguments belong to higher-dimensional Euclidean spaces.

See also

References

  1. Weisstein, Eric W.. "Unimodal". http://mathworld.wolfram.com/Unimodal.html. 
  2. Weisstein, Eric W.. "Mode". http://mathworld.wolfram.com/Mode.html. 
  3. 3.0 3.1 A.Ya. Khinchin (1938). "On unimodal distributions" (in ru). Trams. Res. Inst. Math. Mech. (University of Tomsk) 2 (2): 1–7. 
  4. Hazewinkel, Michiel, ed. (2001), "Unimodal distribution", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4, https://www.encyclopediaofmath.org/index.php?title=U/u095330 
  5. Vladimirovich Gnedenko and Victor Yu Korolev (1996). Random summation: limit theorems and applications. CRC-Press. ISBN 0-8493-2875-6.  p. 31
  6. Medgyessy, P. (March 1972). "On the unimodality of discrete distributions". Periodica Mathematica Hungarica 2 (1–4): 245–257. doi:10.1007/bf02018665. http://www.akademiai.com/content/j5012306777g764n/. 
  7. Gauss, C. F. (1823). "Theoria Combinationis Observationum Erroribus Minimis Obnoxiae, Pars Prior". Commentationes Societatis Regiae Scientiarum Gottingensis Recentiores 5. 
  8. D. F. Vysochanskij, Y. I. Petunin (1980). "Justification of the 3σ rule for unimodal distributions". Theory of Probability and Mathematical Statistics 21: 25–36. 
  9. Sellke, T.M.; Sellke, S.H. (1997). "Chebyshev inequalities for unimodal distributions". American Statistician (American Statistical Association) 51 (1): 34–40. doi:10.2307/2684690. 
  10. 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
  11. Basu, S.; Dasgupta, A. (1997). "The Mean, Median, and Mode of Unimodal Distributions: A Characterization". Theory of Probability & Its Applications 41 (2): 210–223. doi:10.1137/S0040585X97975447. http://epubs.siam.org/doi/pdf/10.1137/S0040585X97975447. 
  12. Bernard, Carole; Kazzi, Rodrigue; Vanduffel, Steven (2020). "Range Value-at-Risk bounds for unimodal distributions under partial information". Insurance: Mathematics and Economics 94: 9–24. doi:10.1016/j.insmatheco.2020.05.013. 
  13. Rohatgi, Vijay K.; Székely, Gábor J. (1989). "Sharp inequalities between skewness and kurtosis". Statistics & Probability Letters 8 (4): 297–299. doi:10.1016/0167-7152(89)90035-7. 
  14. 14.0 14.1 Klaassen, Chris A.J.; Mokveld, Philip J.; Van Es, Bert (2000). "Squared skewness minus kurtosis bounded by 186/125 for unimodal distributions". Statistics & Probability Letters 50 (2): 131–135. doi:10.1016/S0167-7152(00)00090-0. 
  15. "On the unimodality of METRIC Approximation subject to normally distributed demands.". Method in appendix D, Example in theorem 2 page 5. http://homepage.univie.ac.at/thibaut.barthelemy/METRIC.pdf. 
  16. "Mathematical Programming Glossary.". https://glossary.informs.org/indexVer1.php?page=U.html. 
  17. Demaine, Erik D.; Langerman, Stefan (2005). "Optimizing a 2D Function Satisfying Unimodality Properties". in Brodal, Gerth Stølting; Leonardi, Stefano (in en). Algorithms – ESA 2005. Lecture Notes in Computer Science. 3669. Berlin, Heidelberg: Springer. pp. 887–898. doi:10.1007/11561071_78. ISBN 978-3-540-31951-1. https://link.springer.com/chapter/10.1007/11561071_78. 
  18. See e.g. John Guckenheimer; Stewart Johnson (July 1990). "Distortion of S-Unimodal Maps". Annals of Mathematics. Second Series 132 (1): 71–130. doi:10.2307/1971501. 
  19. Godfried T. Toussaint (June 1984). "Complexity, convexity, and unimodality". International Journal of Computer and Information Sciences 13 (3): 197–217. doi:10.1007/bf00979872.