Mode (statistics)

For other uses, see Mode.

The mode is the value that appears most often in a set of data. The mode of a discrete probability distribution is the value x at which its probability mass function takes its maximum value. In other words, it is the value that is most likely to be sampled. The mode of a continuous probability distribution is the value x at which its probability density function has its maximum value, so the mode is at the peak.

Like the statistical mean and median, the mode is a way of expressing, in a single number, important information about a random variable or a population. The numerical value of the mode is the same as that of the mean and median in a normal distribution, and it may be very different in highly skewed distributions.

The mode is not necessarily unique, since the probability mass function or probability density function may take the same maximum value at several points x₁, x₂, etc. The most extreme case occurs in uniform distributions, where all values occur equally frequently. When a probability density function has multiple local maxima it is common to refer to all of the local maxima as modes of the distribution. Such a continuous distribution is called multimodal (as opposed to unimodal).

In symmetric unimodal distributions, such as the normal distribution, the mean (if defined), median and mode all coincide. For samples, if it is known that they are drawn from a symmetric distribution, the sample mean can be used as an estimate of the population mode.

Mode of a sample

The mode of a sample is the element that occurs most often in the collection. For example, the mode of the sample [1, 3, 6, 6, 6, 6, 7, 7, 12, 12, 17] is 6. Given the list of data [1, 1, 2, 4, 4] the mode is not unique - the dataset may be said to be bimodal, while a set with more than two modes may be described as multimodal.

For a sample from a continuous distribution, such as [0.935..., 1.211..., 2.430..., 3.668..., 3.874...], the concept is unusable in its raw form, since no two values will be exactly the same, so each value will occur precisely once. In order to estimate the mode, the usual practice is to discretize the data by assigning frequency values to intervals of equal distance, as for making a histogram, effectively replacing the values by the midpoints of the intervals they are assigned to. The mode is then the value where the histogram reaches its peak. For small or middle-sized samples the outcome of this procedure is sensitive to the choice of interval width if chosen too narrow or too wide; typically one should have a sizable fraction of the data concentrated in a relatively small number of intervals (5 to 10), while the fraction of the data falling outside these intervals is also sizable. An alternate approach is kernel density estimation, which essentially blurs point samples to produce a continuous estimate of the probability density function which can provide an estimate of the mode.

The following MATLAB (or Octave) code example computes the mode of a sample:

X = sort(x);
indices   =  find(diff([X; realmax]) > 0); % indices where repeated values change
[modeL,i] =  max (diff([0; indices]));     % longest persistence length of repeated values
mode      =  X(indices(i));

The algorithm requires as a first step to sort the sample in ascending order. It then computes the discrete derivative of the sorted list, and finds the indices where this derivative is positive. Next it computes the discrete derivative of this set of indices, locating the maximum of this derivative of indices, and finally evaluates the sorted sample at the point where that maximum occurs, which corresponds to the last member of the stretch of repeated values.

Comparison of mean, median and mode

Comparison of common averages of values { 1, 2, 2, 3, 4, 7, 9 }
Type	Description	Example	Result
Arithmetic mean	Sum of values of a data set divided by number of values: $\scriptstyle\bar{x} = \frac{1}{n}\sum_{i=1}^n x_i$	(1+2+2+3+4+7+9) / 7	4
Median	Middle value separating the greater and lesser halves of a data set	1, 2, 2, 3, 4, 7, 9	3
Mode	Most frequent value in a data set	1, 2, 2, 3, 4, 7, 9	2

Use

Unlike mean and median, the concept of mode also makes sense for "nominal data" (i.e., not consisting of numerical values in the case of mean, or even of ordered values in the case of median). For example, taking a sample of Korean family names, one might find that "Kim" occurs more often than any other name. Then "Kim" would be the mode of the sample. In any voting system where a plurality determines victory, a single modal value determines the victor, while a multi-modal outcome would require some tie-breaking procedure to take place.

Unlike median, the concept of mode makes sense for any random variable assuming values from a vector space, including the real numbers (a one-dimensional vector space) and the integers (which can be considered embedded in the reals). For example, a distribution of points in the plane will typically have a mean and a mode, but the concept of median does not apply. The median makes sense when there is a linear order on the possible values. Generalizations of the concept of median to higher-dimensional spaces are the geometric median and the centerpoint.

Uniqueness and definedness

For the remainder, the assumption is that we have (a sample of) a real-valued random variable.

For some probability distributions, the expected value may be infinite or undefined, but if defined, it is unique. The mean of a (finite) sample is always defined. The median is the value such that the fractions not exceeding it and not falling below it are each at least 1/2. It is not necessarily unique, but never infinite or totally undefined. For a data sample it is the "halfway" value when the list of values is ordered in increasing value, where usually for a list of even length the numerical average is taken of the two values closest to "halfway". Finally, as said before, the mode is not necessarily unique. Certain pathological distributions (for example, the Cantor distribution) have no defined mode at all. For a finite data sample, the mode is one (or more) of the values in the sample.

Properties

Assuming definedness, and for simplicity uniqueness, the following are some of the most interesting properties.

All three measures have the following property: If the random variable (or each value from the sample) is subjected to the linear or affine transformation which replaces X by aX+b, so are the mean, median and mode.
However, if there is an arbitrary monotonic transformation, only the median follows; for example, if X is replaced by exp(X), the median changes from m to exp(m) but the mean and mode won't.
Except for extremely small samples, the mode is insensitive to "outliers" (such as occasional, rare, false experimental readings). The median is also very robust in the presence of outliers, while the mean is rather sensitive.
In continuous unimodal distributions the median lies, as a rule of thumb, between the mean and the mode, about one third of the way going from mean to mode. In a formula, median ≈ (2 × mean + mode)/3. This rule, due to Karl Pearson, often applies to slightly non-symmetric distributions that resemble a normal distribution, but it is not always true and in general the three statistics can appear in any order.^[1]^[2]
For unimodal distributions, the mode is within $\sqrt{3}$ standard deviations of the mean, and the root mean square deviation about the mode is between the standard deviation and twice the standard deviation.^[3]

Example for a skewed distribution

An example of a skewed distribution is personal wealth: Few people are very rich, but among those some are extremely rich. However, many are rather poor.

Comparison of mean, median and mode of two log-normal distributions with different skewness.

A well-known class of distributions that can be arbitrarily skewed is given by the log-normal distribution. It is obtained by transforming a random variable X having a normal distribution into random variable Y = e^X. Then the logarithm of random variable Y is normally distributed, hence the name.

Taking the mean μ of X to be 0, the median of Y will be 1, independent of the standard deviation σ of X. This is so because X has a symmetric distribution, so its median is also 0. The transformation from X to Y is monotonic, and so we find the median e⁰ = 1 for Y.

When X has standard deviation σ = 0.25, the distribution of Y is weakly skewed. Using formulas for the log-normal distribution, we find:

\begin{array}{rlll} \text{mean} & = e^{\mu + \sigma^2 / 2} & = e^{0 + 0.25^2 / 2} & \approx 1.032 \\ \text{mode} & = e^{\mu - \sigma^2} & = e^{0 - 0.25^2} & \approx 0.939 \\ \text{median} & = e^\mu & = e^0 & = 1 \end{array}

Indeed, the median is about one third on the way from mean to mode.

When X has a larger standard deviation, σ = 1, the distribution of Y is strongly skewed. Now

\begin{array}{rlll} \text{mean} & = e^{\mu + \sigma^2 / 2} & = e^{0 + 1^2 / 2} & \approx 1.649 \\ \text{mode} & = e^{\mu - \sigma^2} & = e^{0 - 1^2} & \approx 0.368 \\ \text{median} & = e^\mu & = e^0 & = 1 \end{array}

Here, Pearson's rule of thumb fails.

Van Zwet condition

Van Zwet derived an inequality which provides sufficient conditions for this inequality to hold.^[4] The inequality

Mode ≤ Median ≤ Mean

holds if

F( Median - x ) + F( Median + x ) ≥ 1

for all x where F() is the cumulative distribution function of the distribution.

Unimodal distributions

It can be shown for a unimodal distribution that the median $\tilde{X}$ and the mean $\bar{X}$ lie within (3/5)^1/2 ≈ 0.7746 standard deviations of each other.^[5] In symbols,

\frac{\left|\tilde{X} - \bar{X}\right|}{\sigma} \le (3/5)^{1/2}

where |.| is the absolute value.

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

\frac{\left|\tilde{X} - \mathrm{mode}\right|}{\sigma} \le 3^{1/2}.

Confidence interval for the mode with a single data point

It is a common but false belief that from a single observation x we can not gain information about the variability in the population and that consequently that ﬁnite length conﬁdence intervals for mean and/or variance are impossible even in principle.

It is possible for an unknown unimodal distribution to estimate a confidence interval for the mode with a sample size of 1.^[6] This was first shown by Abbot and Rosenblatt^[7] and extended by Blachman and Machol.^[8] This confidence interval can be sharpened if the distribution can be assumed to be symmetrical. It is further possible to sharpen this interval if the distribution is normally distributed.

Let the confidence interval be 1 - α. Then the confidence intervals for the general, symmetric and normally distributed variates respectively are

X \pm \left( \frac{ 2 }{ \alpha } - 1 \right) | X - \theta |

X \pm \left( \frac{ 1 }{ \alpha } - 1 \right) | X - \theta |

X \pm \left( \frac{ 0.484 }{ \alpha } - 1 \right) | X - \theta |

where X is the variate, θ is the mode and || is the absolute value.

These estimates are conservative. The confidence intervals for the mode at the 90% level given by these estimators are X ± 19 | X - θ |, X ± 9 | X - θ | and X ± 5.84 | X - θ | for the general, symmetric and normally distributed variates respectively. The 95% confidence interval for a normally distributed variate is given by X ± 10.7 | X - θ |. It may be worth noting that the mean and the mode coincide if the variates are normally distributed.

The 95% bound for a normally distributed variate has been improved and is now known to be X ± 9.68 | X - θ |^[9] The bound for a 99% confidence interval is X ± 48.39 | X - θ'|

Note

Machol has shown that given a known density symmetrical about 0 that given a single sample value (x) that the 90% confidence intervals of population mean are^[10]

x \pm 5 | x - \nu |

where ν is the population median.

If the precise form of the distribution is not known but it is known to be symmetrical about zero then we have

P( X - k | X - a | \le \mu \le X + k | X - a |) \ge 1 - \frac{ 1 }{ 1 + k }

where X is the variate, μ is the population mean and a and k are arbitrary real numbers.

It is also possible to estimate a confidence interval for the standard deviation from a single observation if the distribution is symmetrical about 0.^[11] For a normal distribution the with an unknown variance and a single data point (X) the 90%, 95% and 99% confidence intervals for the standard deviation are [ 0, 8|X| ], [ 0, 17|X| ] and [ 0, 70|X| ]. These intervals may be shorted if the mean is known to be bounded by a multiple of the standard deviation.

If the distribution is known to be normal then it is possible to estimate a confidence interval for the mean and variance from a simple value.^[12] The 90% confidence intervals are

X - 23.3 | X | \le \mu \le X + 23.3 | X |

\sigma \le 10 | X |

The confidence intervals can be estimated for any chosen range.

This method is not limited to the normal distribution but can be used with any known distribution.

Statistical tests

These estimators have been used to create hypothesis tests for simple samples from normal or symmetrical unimodal distributions.^[13] Let the distribution have an assumed mean ( μ₀ ). The null hypothesis is that the assumed mean of the distribution lies within the confidence interval of the sample mean ( m ). The null hypothesis is accepted if

\mu_0 < \frac{ x + m }{ 2 } \pm k | x - m |

where x is the value of the sample and k is a constant. The null hypothesis is rejected if

\mu_0 > \frac{ x + m }{ 2 } \pm k | x - m |

The value of k depends on the choice of confidence interval and the nature of the assumed distribution.

If the distribution is assumed or is known to be normal then the values of k for the 50%, 66.6%, 75%, 80%, 90%, 95% and 99% confidence intervals are 0.50, 1.26, 1.80, 2.31, 4.79, 9.66 and 48.39 respectively.

If the distribution is assumed or known to be unimodal and symmetrical but not normal then the values of k for the 50%, 66.6%, 75%, 80%, 90%, 95% and 99% confidence intervals are 0.50, 1.87, 2.91, 3.94, 8.97, 18.99, 99.00 respectively.

To see how this test works we assume or know a priori that the population from which the sample is drawn has a mean of μ₀ and that the population has a symmetrical unimodal distribution - a class that includes the normal distribution. We wish to know if the mean estimated from the sample is representative of the population at a pre chosen level of confidence.

Assume that the distribution is normal and let the confidence interval be 95%. Then k = 9.66.

Assuming that the sample is representative of the population, the sample mean ( m ) will then lie within the range determined from the formula:

\mu_0 < \frac{ x + m }{ 2 } \pm 9.66 | x - m |

If subsequent sampling shows that the sample mean lies outside these parameters the sample mean is to be considered to differ significantly from the population mean.

History

The term mode originates with Karl Pearson in 1895.^[14]

References

↑ "Relationship between the mean, median, mode, and standard deviation in a unimodal distribution".
↑ Hippel, Paul T. von (2005). "Mean, Median, and Skew: Correcting a Textbook Rule". J. of Statistics Education 13 (2).
↑ Bottomley, H. (2004). "Maximum distance between the mode and the mean of a unimodal distribution" (PDF). Unpublished preprint.
↑ van Zwet WR (1979) "Mean, median, mode II", Statistica Neerlandica, 33 (1) 1–5
↑ Basu, Sanjib, and Anirban DasGupta. "The mean, median, and mode of unimodal distributions: a characterization." Theory of Probability & Its Applications 41.2 (1997): 210-223.
↑ Edelman, D. (1990). "A confidence interval for the center of an unknown unimodal distribution based on a sample of size 1". The American Statistician 44 (4): 285–287. doi:10.1080/00031305.1990.10475740.
↑ Abbot, J. H.; Rosenblatt, J. (1963). "Two stage estimation with one observation on the first stage". Annals of the Institute of Statistical Mathematics 14 (1): 229–235. doi:10.1007/BF02868644.
↑ Blachman, N. M.; Machol, R. (1987). "Confidence intervals based on one or more observations". IEEE Transactions on Information Theory 33 (3): 373–382. doi:10.1109/TIT.1987.1057306.
↑ Wall, M. M.; Boen, J.; Tweedie, R. (2001). "An effective confidence interval for the mean With samples of size one and two". The American Statistician 55 (2): 102–105. doi:10.1198/000313001750358400.
↑ Machol R (1964) IEEE Trans Info Theor
↑ Rodríguez, C. C. (1996). "Confidence Intervals from one Observation". Maximum Entropy and Bayesian Methods. Fundamental Theories of Physics 70. pp. 175–182. arXiv:bayes-an/9504001. doi:10.1007/978-94-009-0107-0_19. ISBN 978-94-010-6534-4.
↑ Rosenblatt, J. (1966). "Confidence interval for standard deviation from a single observation". Technometrics 8 (2): 367–368. doi:10.1080/00401706.1966.10490358.
↑ Valentine, H. T.; Furnival, G. M.; Gregoire, T. G. (1991). "Confidence Intervals from Single Observations in Forest Research". Forest Science 37 (1): 370–373.
↑ Pearson, Karl (1895). "Contributions to the Mathematical Theory of Evolution. II. Skew Variation in Homogeneous Material", Philosophical Transactions of the Royal Society of London, Ser. A, 186, 343-414

External links

Hazewinkel, Michiel, ed. (2001), "Mode", Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
A Guide to Understanding & Calculating the Mode
Weisstein, Eric W., "Mode", MathWorld.
Mean, Median and Mode short beginner video from Khan Academy

Statistics

Descriptive statistics

Continuous data

Location	Mean arithmetic geometric harmonic Median Mode

Dispersion	Range Standard deviation Coefficient of variation Percentile Interquartile range

Shape	Variance Skewness Kurtosis Moments L-moments

Count data

Index of dispersion

Summary tables

Dependence

Statistical graphics

Data collection

Study design	Effect size Standard error Statistical power Sample size determination

Survey methodology	Sampling stratified cluster Opinion poll Questionnaire

Controlled experiments	Design control optimal Controlled trial Randomized Random assignment Replication Blocking Factorial experiment

Uncontrolled studies	Observational study Natural experiment Quasi-experiment

Statistical inference

Statistical theory

Frequentist inference

Confidence interval Testing hypotheses Power

Unbiased estimators	Mean unbiased minimum-variance Median unbiased

Biased estimators	Maximum likelihood Method of moments Minimum distance Density estimation

Parametric tests	Likelihood-ratio Wald Score

Specific tests

Z (normal) Student's t-test F Shapiro–Wilk Kolmogorov–Smirnov

Goodness of fit	Chi-squared G Sample source (Anderson–Darling) Sample normality (Shapiro–Wilk) Skewness / kurtosis normality (Jarque-Bera) Model comparison (Likelihood-ratio) Model quality (Akaike criterion)

Signed-rank	1-sample (Wilcoxon) 2-sample (Mann–Whitney U) 1-way anova (Kruskal–Wallis)

Bayesian inference

Correlation	Pearson product–moment Partial correlation Confounding variable Coefficient of determination

Regression analysis	Errors and residuals Regression model validation Mixed effects models Simultaneous equations models Multivariate adaptive regression splines (MARS)

Linear regression	Simple linear regression Ordinary least squares General linear model Bayesian regression

Non-standard predictors	Nonlinear regression Nonparametric Semiparametric Isotonic Robust Heteroscedasticity Homoscedasticity

Generalized linear model	Exponential families Logistic (Bernoulli) / Binomial / Poisson regressions

Partition of variance	Analysis of variance (ANOVA, anova) Analysis of covariance Multivariate ANOVA Degrees of freedom

Categorical / Multivariate / Time-series / Survival analysis

Categorical

Multivariate

Time-series

General	Decomposition Trend Stationarity Seasonal adjustment Exponential smoothing Cointegration Structural break Granger causality

Specific tests	Dickey–Fuller Johansen Q-statistic (Ljung–Box) Durbin–Watson Breusch–Godfrey

Time domain	Autocorrelation (ACF) partial (PACF) Cross-correlation (XCF) ARMA model ARIMA model (Box–Jenkins) Autoregressive conditional heteroskedasticity (ARCH) Vector autoregression (VAR)

Frequency domain	Spectral density estimation Fourier analysis Wavelet

Survival

Survival function	Kaplan–Meier estimator (product limit) Proportional hazards models Accelerated failure time (AFT) model First hitting time

Hazard function	Nelson–Aalen estimator

Test	Log-rank test

Applications

Biostatistics	Bioinformatics Clinical trials / studies Epidemiology Medical statistics

Engineering statistics	Chemometrics Methods engineering Probabilistic design Process / quality control Reliability System identification

Social statistics	Actuarial science Census Crime statistics Demography Econometrics National accounts Official statistics Population Psychometrics

Spatial statistics	Cartography Environmental statistics Geographic information system Geostatistics Kriging

Category
Portal
Commons
WikiProject

This article is issued from Wikipedia - version of the Monday, January 11, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.