AP Statistics Curriculum 2007 EDA Center
From Socr
(→Other Measures of Centrality) |
m (→Medimean) |
||
(24 intermediate revisions not shown) | |||
Line 2: | Line 2: | ||
===Measurements of Central Tendency=== | ===Measurements of Central Tendency=== | ||
- | There are three main features of all populations (or data samples) that are always critical in understanding and interpreting their distributions. These characteristics are '''Center''', '''Spread''' and '''Shape'''. The main | + | There are three main features of all populations (or data samples) that are always critical in understanding and interpreting their distributions. These characteristics are '''Center''', '''Spread''' and '''Shape'''. The main measures of centrality are '''mean''', '''median''' and '''mode'''. |
- | Suppose we are interested in the long-jump performance of some students. We can carry an experiment by randomly selecting | + | Suppose we are interested in the long-jump performance of some students. We can carry an experiment by randomly selecting eight male statistics students and ask them to perform the standing long jump. In reality every student participated, but for the ease of calculations below we will focus on these eight students. The long jumps were as follows: |
<center> | <center> | ||
{| class="wikitable" style="text-align:center; width:75%" border="1" | {| class="wikitable" style="text-align:center; width:75%" border="1" | ||
Line 11: | Line 11: | ||
| 74 || 78 || 106 || 80 || 68 || 64 || 60 || 76 | | 74 || 78 || 106 || 80 || 68 || 64 || 60 || 76 | ||
|}</center> | |}</center> | ||
- | |||
===Mean=== | ===Mean=== | ||
- | The '''sample-mean''' is the arithmetic average of a finite sample of numbers. In the long-jump example, the sample-mean is calculated as follows: | + | The '''sample-mean''' is the arithmetic average of a finite sample of numbers. For instance, the mean of the sample <math>x_1, x_2, x_3, \cdots, x_{n-1}, x_n</math> (short hand notation: <math>\{x_{i}\}_{i=1}^n</math>), is given by: |
+ | <center><math>\bar{x}={1\over n}\sum_{i=1}^{n}{x_{i}}.</math> | ||
+ | </center> | ||
+ | |||
+ | In the long-jump example, the sample-mean is calculated as follows: | ||
<center><math>\overline{y} = {1 \over 8} (74+78+106+80+68+64+60+76)=75.75 in.</math></center> | <center><math>\overline{y} = {1 \over 8} (74+78+106+80+68+64+60+76)=75.75 in.</math></center> | ||
Line 37: | Line 40: | ||
The '''modes''' represent the most frequently occurring values (The numbers that appear the most). The term mode is applied both to probability distributions and to collections of experimental data. | The '''modes''' represent the most frequently occurring values (The numbers that appear the most). The term mode is applied both to probability distributions and to collections of experimental data. | ||
- | For instance, for the [[SOCR_012708_ID_Data_HotDogs | Hot dogs data file]], there | + | For instance, for the [[SOCR_012708_ID_Data_HotDogs | Hot dogs data file]], there appears to be three modes for the calorie variable. This is presented by the histogram of the '''Calorie''' content of all hotdogs, shown in the image below. Note the clear separation of the calories into three distinct sub-populations - the highest points in these three sub-populations are the three modes for these data. |
<center>[[Image:SOCR_EBook_Dinov_EDA_012708_Fig3.jpg|500px]]</center> | <center>[[Image:SOCR_EBook_Dinov_EDA_012708_Fig3.jpg|500px]]</center> | ||
Line 46: | Line 49: | ||
If you remove the student with the long jump distance of 106 and recalculate the median and mean, which one is altered less (therefore is more resistant)? Notice that the mean is very sensitive to outliers and atypical observations, and hence less resistant than the median. | If you remove the student with the long jump distance of 106 and recalculate the median and mean, which one is altered less (therefore is more resistant)? Notice that the mean is very sensitive to outliers and atypical observations, and hence less resistant than the median. | ||
- | === | + | ===Resistant Mean-related Measures of Centrality=== |
The following two sample measures of population centrality estimate resemble the calculations of the [[AP_Statistics_Curriculum_2007_Distrib_MeanVar#Expectation_.28Mean.29 | mean]], however they are much more ''resistant'' to change in the presence of outliers. | The following two sample measures of population centrality estimate resemble the calculations of the [[AP_Statistics_Curriculum_2007_Distrib_MeanVar#Expectation_.28Mean.29 | mean]], however they are much more ''resistant'' to change in the presence of outliers. | ||
Line 52: | Line 55: | ||
<math>\bar{y}_{t,k}={1\over n-2k}\sum_{i=k+1}^{n-k}{y_{(i)}}</math>, where <math>k\geq 0</math> is the trim-factor (large k, yield less variant estimates of center), and <math>y_{(i)}</math> are the order statistics (small to large). That is, we remove the smallest and the largest ''k'' observations from the sample, before we compute the arithmetic average. | <math>\bar{y}_{t,k}={1\over n-2k}\sum_{i=k+1}^{n-k}{y_{(i)}}</math>, where <math>k\geq 0</math> is the trim-factor (large k, yield less variant estimates of center), and <math>y_{(i)}</math> are the order statistics (small to large). That is, we remove the smallest and the largest ''k'' observations from the sample, before we compute the arithmetic average. | ||
- | ====''' | + | ===='''Winsorized k-times mean'''==== |
- | The | + | The Winsorized k-times mean is defined similarly by <math>\bar{y}_{w,k}={1\over n}( k\times y_{(k)}+\sum_{i=k+1}^{n-k}{y_{(i)}}+k\times y_{(n-k+1)})</math>, where <math>k\geq 0</math> is the trim-factor and <math>y_{(i)}</math> are the order statistics (small to large). In this case, before we compute the arithmetic average, we replace the ''k'' smallest and the ''k'' largest observations with the k<sup>th</sup> and (n-k)<sup>th</sup> largest observations, respectively. |
+ | |||
+ | ===Other Centrality Measures=== | ||
+ | The ''arithmetic'' mean answers the question, ''if all observations were equal, what would that value (''center'') have to be in order to achieve the same <u>total</u>?'' | ||
+ | : <math>n\times \bar{x}=\sum_{i=1}^n{x_i}</math> | ||
+ | |||
+ | In some situations, there is a need to think of the average in different terms, not in terms of arithmetic average. | ||
====Harmonic Mean==== | ====Harmonic Mean==== | ||
- | + | If we study speeds (velocities) the ''arithmetic'' mean is inappropriate. However the [http://en.wikibooks.org/wiki/Statistics/Summary/Averages/Harmonic_Mean harmonic mean (computed differently)] gives the most intuitive answer to what the "''middle''" is for a process. The harmonic mean answers the question ''if all the observations were equal, what would that value have to be in order to achieve the same sample <u>sum of reciprocals</u>?'' | |
+ | |||
+ | : ''Harmonic mean'': <math>\hat{\hat{x}}= \frac{n}{\frac{1}{x_1} + \frac{1}{x_2} + \frac{1}{x_3} + \ldots + \frac{1}{x_n}}</math> | ||
====Geometric Mean==== | ====Geometric Mean==== | ||
- | + | In contrast, the [http://en.wikibooks.org/wiki/Statistics/Summary/Averages/Geometric_Mean geometric mean] answers the question ''if all the observations were equal, what would that value have to be in order to achieve the same sample <u>product</u>?'' | |
+ | : ''Geometric mean'': <math>\tilde{x}^n={\prod_{i=1}^n x_i}</math> | ||
+ | |||
+ | : Alternatively: <math>\tilde{x}= \exp \left( \frac{1}{n} \sum_{i=1}^n\log(x_i) \right)</math> | ||
+ | |||
+ | ====Medimean==== | ||
+ | We have already seen the standard definitions for [[AP_Statistics_Curriculum_2007_EDA_Center#Mean|mean]] and [[AP_Statistics_Curriculum_2007_EDA_Center#Median|median]]. There is an alternative ways to define these measures of centrality using optimization theory. | ||
+ | : '''Median''' = \( argmin_x \big ( \sum_{s\in S} {|s-x|^1} \big )\) | ||
+ | : '''Mean''' = \( argmin_x \big ( \sum_{s\in S} {|s-x|^2}\big )\) | ||
+ | : where \( argmin_x \) provides the value x which minimizes the distance function in the operator. | ||
+ | |||
+ | The median and mean, as measures of centrality, represent points that are close to all the values randomly generated by the process. Thus, these parameters minimize the total sum of distances between them and each possible process observation, however, the differences between them is related to the differences in the metric used to compute distances between pairs of numbers. | ||
+ | |||
+ | For the mean, the expression above uses the squared '''Euclidean distance''', \(L_2\) norm, which is procedurally different from the arithmetic averaging, or expectation calculation we have seen below, but generates the same result. The median optimization-based expression relies on the '''Manhattan distance''', or \(L_1\) norm. | ||
+ | |||
+ | : The '''Medimean''' is a fusion of median and mean, (which is another word for average), which is defined by: | ||
+ | : '''Medimean(k)''' = \( argmin_x \big ( \sum_{s\in S} {|s-x|^k} \big )\) | ||
+ | |||
+ | The medimean depends on one parameter ''k'', the distance exponent. For ''k=1'', the medimean coincides with the median, and for ''k=2'', it equals the mean. Thus, the medimean is a continuous and natural extension of the median and the mean for higher order exponents. | ||
+ | |||
+ | The medimean may be used as a compromise between the mean and the median. For example, if you are interested in something like a median, but which moves just a little bit in the direction of far-away elements, if there are any. In that case you may want to use a medimean(1.2) in the analysis. | ||
+ | |||
+ | How do we prove the equivalence of ht standard [[AP_Statistics_Curriculum_2007_EDA_Center#Mean|mean]] and [[AP_Statistics_Curriculum_2007_EDA_Center#Median|median]] definitions and the optimization based definitions of these population parameters? | ||
+ | |||
+ | * For the '''Median''' = \( argmin_x \big ( \sum_{s\in S} {|s-x|^1} \big )\), this translates in the continuous variable case to \( argmin_x \big ( \int_{s\in S} {|s-x|^1P(x)dx} \big )\)= | ||
+ | : \( argmin_x \big ( \int_{s<m} {(m-s)P(s)ds} + \int_{s>m} {(m-s)P(s)ds} \big )\)= | ||
+ | : \( argmin_x \big ( m \int_{s<m} {P(s)ds} - \int_{s<m} {xP(s)ds} + m \int_{s>m} {P(s)ds} -\int_{s>m} {sP(s)ds} \big )\)= | ||
+ | : \( argmin_x \big ( m (F(m) -(1-F(m))) + \int_{s>m} {sP(s)ds} -\int_{s<m} {sP(s)ds} \big )\), where the [[AP_Statistics_Curriculum_2007_Distrib_RV#Probability_density.2Fmass_and_.28cumulative.29_distribution_functions|cumulative distribution function]] \( F(m)= \int_{s<m} {P(s)ds} \). | ||
+ | : And this function, \( m (2F(m) -1) + \int_{s>m} {sP(s)ds} -\int_{s<m} {sP(s)ds}\), is minimized for \(m=F^{-1}(0.5)\), which is the standard definition of the [[AP_Statistics_Curriculum_2007_EDA_Center#Median|median]] . | ||
+ | |||
+ | * For the '''Mean''' = \( argmin_x \big ( \sum_{s\in S} {|s-x|^2}\big )\), observe that minimizing \( \sum_{s\in S} {|s-x|^2}\) is similar to minimizing the [[AP_Statistics_Curriculum_2007_Distrib_MeanVar#Variance|variance]] of the process \( \sum_{s\in S} {(s-x)^2P(x)}\), which is identical to minimizing \(E(X^2)-\mu^2\), for \(\mu\). | ||
+ | |||
+ | ===[[EBook_Problems_EDA_Center | Problems]]=== | ||
<hr> | <hr> |
Current revision as of 01:58, 2 May 2013
Contents |
General Advance-Placement (AP) Statistics Curriculum - Central Tendency
Measurements of Central Tendency
There are three main features of all populations (or data samples) that are always critical in understanding and interpreting their distributions. These characteristics are Center, Spread and Shape. The main measures of centrality are mean, median and mode.
Suppose we are interested in the long-jump performance of some students. We can carry an experiment by randomly selecting eight male statistics students and ask them to perform the standing long jump. In reality every student participated, but for the ease of calculations below we will focus on these eight students. The long jumps were as follows:
74 | 78 | 106 | 80 | 68 | 64 | 60 | 76 |
Mean
The sample-mean is the arithmetic average of a finite sample of numbers. For instance, the mean of the sample (short hand notation: ), is given by:
In the long-jump example, the sample-mean is calculated as follows:
Median
The sample-median can be thought of as the point that divides a distribution in half (50/50). The following steps are used to find the sample-median:
- Arrange the data in ascending order
- If the sample size is odd, the median is the middle value of the ordered collection
- If the sample size is even, the median is the average of the middle two values in the ordered collection.
For the long-jump data above we have:
- Ordered data:
60 | 64 | 68 | 74 | 76 | 78 | 80 | 106 |
- .
Mode(s)
The modes represent the most frequently occurring values (The numbers that appear the most). The term mode is applied both to probability distributions and to collections of experimental data.
For instance, for the Hot dogs data file, there appears to be three modes for the calorie variable. This is presented by the histogram of the Calorie content of all hotdogs, shown in the image below. Note the clear separation of the calories into three distinct sub-populations - the highest points in these three sub-populations are the three modes for these data.
Resistance
A statistic is said to be resistant if the value of the statistic is relatively unchanged by changes in a small portion of the data. Referencing the formulas for the median, mean and mode which statistic seems to be more resistant?
If you remove the student with the long jump distance of 106 and recalculate the median and mean, which one is altered less (therefore is more resistant)? Notice that the mean is very sensitive to outliers and atypical observations, and hence less resistant than the median.
Resistant Mean-related Measures of Centrality
The following two sample measures of population centrality estimate resemble the calculations of the mean, however they are much more resistant to change in the presence of outliers.
K-times trimmed mean
, where is the trim-factor (large k, yield less variant estimates of center), and y_{(i)} are the order statistics (small to large). That is, we remove the smallest and the largest k observations from the sample, before we compute the arithmetic average.
Winsorized k-times mean
The Winsorized k-times mean is defined similarly by , where is the trim-factor and y_{(i)} are the order statistics (small to large). In this case, before we compute the arithmetic average, we replace the k smallest and the k largest observations with the k^{th} and (n-k)^{th} largest observations, respectively.
Other Centrality Measures
The arithmetic mean answers the question, if all observations were equal, what would that value (center) have to be in order to achieve the same total?
In some situations, there is a need to think of the average in different terms, not in terms of arithmetic average.
Harmonic Mean
If we study speeds (velocities) the arithmetic mean is inappropriate. However the harmonic mean (computed differently) gives the most intuitive answer to what the "middle" is for a process. The harmonic mean answers the question if all the observations were equal, what would that value have to be in order to achieve the same sample sum of reciprocals?
- Harmonic mean:
Geometric Mean
In contrast, the geometric mean answers the question if all the observations were equal, what would that value have to be in order to achieve the same sample product?
- Geometric mean:
- Alternatively:
Medimean
We have already seen the standard definitions for mean and median. There is an alternative ways to define these measures of centrality using optimization theory.
- Median = \( argmin_x \big ( \sum_{s\in S} {|s-x|^1} \big )\)
- Mean = \( argmin_x \big ( \sum_{s\in S} {|s-x|^2}\big )\)
- where \( argmin_x \) provides the value x which minimizes the distance function in the operator.
The median and mean, as measures of centrality, represent points that are close to all the values randomly generated by the process. Thus, these parameters minimize the total sum of distances between them and each possible process observation, however, the differences between them is related to the differences in the metric used to compute distances between pairs of numbers.
For the mean, the expression above uses the squared Euclidean distance, \(L_2\) norm, which is procedurally different from the arithmetic averaging, or expectation calculation we have seen below, but generates the same result. The median optimization-based expression relies on the Manhattan distance, or \(L_1\) norm.
- The Medimean is a fusion of median and mean, (which is another word for average), which is defined by:
- Medimean(k) = \( argmin_x \big ( \sum_{s\in S} {|s-x|^k} \big )\)
The medimean depends on one parameter k, the distance exponent. For k=1, the medimean coincides with the median, and for k=2, it equals the mean. Thus, the medimean is a continuous and natural extension of the median and the mean for higher order exponents.
The medimean may be used as a compromise between the mean and the median. For example, if you are interested in something like a median, but which moves just a little bit in the direction of far-away elements, if there are any. In that case you may want to use a medimean(1.2) in the analysis.
How do we prove the equivalence of ht standard mean and median definitions and the optimization based definitions of these population parameters?
- For the Median = \( argmin_x \big ( \sum_{s\in S} {|s-x|^1} \big )\), this translates in the continuous variable case to \( argmin_x \big ( \int_{s\in S} {|s-x|^1P(x)dx} \big )\)=
- \( argmin_x \big ( \int_{s<m} {(m-s)P(s)ds} + \int_{s>m} {(m-s)P(s)ds} \big )\)=
- \( argmin_x \big ( m \int_{s<m} {P(s)ds} - \int_{s<m} {xP(s)ds} + m \int_{s>m} {P(s)ds} -\int_{s>m} {sP(s)ds} \big )\)=
- \( argmin_x \big ( m (F(m) -(1-F(m))) + \int_{s>m} {sP(s)ds} -\int_{s<m} {sP(s)ds} \big )\), where the cumulative distribution function \( F(m)= \int_{s<m} {P(s)ds} \).
- And this function, \( m (2F(m) -1) + \int_{s>m} {sP(s)ds} -\int_{s<m} {sP(s)ds}\), is minimized for \(m=F^{-1}(0.5)\), which is the standard definition of the median .
- For the Mean = \( argmin_x \big ( \sum_{s\in S} {|s-x|^2}\big )\), observe that minimizing \( \sum_{s\in S} {|s-x|^2}\) is similar to minimizing the variance of the process \( \sum_{s\in S} {(s-x)^2P(x)}\), which is identical to minimizing \(E(X^2)-\mu^2\), for \(\mu\).
Problems
References
- SOCR Home page: http://www.socr.ucla.edu
Translate this page: