Statistics using excel_succinctly

2
By
Charles Zaiontz
Foreword by Daniel Jebaraj

3
Copyright © 2015 by Syncfusion Inc.
2501 Aerial Center Parkway
Suite 200
Morrisville, NC 27560
USA
All rights reserved.
mportant licensing information. Please read.
This book is available for free download from www.syncfusion.com on completion of a registration form.
If you obtained this book from any other source, please register and download a free copy from
www.syncfusion.com.
This book is licensed for reading only if obtained from www.syncfusion.com.
This book is licensed strictly for personal or educational use.
Redistribution in any form is prohibited.
The authors and copyright holders provide absolutely no warranty for any information provided.
The authors and copyright holders shall not be liable for any claim, damages, or any other liability arising
from, out of, or in connection with the information in this book.
Please do not use this book if the listed terms are unacceptable.
Use shall constitute acceptance of the terms listed.
SYNCFUSION, SUCCINCTLY, DELIVER INNOVATION WITH EASE, ESSENTIAL, and .NET ESSENTIALS are the
registered trademarks of Syncfusion, Inc.
Technical Reviewer: Vittoria Ardino
Copy Editor: Suzanne Kattau
Acquisitions Coordinator: Hillary Bowling, marketing coordinator, Syncfusion, Inc.
Proofreader: Darren West, content producer, Syncfusion, Inc.
I

4
Table of Contents
The Story behind the Succinctly Series of Books..................................................................................9
Chapter 1 Introduction..............................................................................................................................11
What is Statistics?..................................................................................................................................11
Why do Statistics in Excel?....................................................................................................................11
Where is Statistics Used?......................................................................................................................12
What Will This Book Cover? ..................................................................................................................12
Chapter 2 Excel Environment ..................................................................................................................13
The Basics .............................................................................................................................................13
User Interface ........................................................................................................................................13
Array Formulas ......................................................................................................................................15
Array Functions......................................................................................................................................16
Data Analysis Tools ...............................................................................................................................16
Chapter 3 Descriptive Statistics ..............................................................................................................20
Basic Terminology .................................................................................................................................20
Measures of Central Tendency..............................................................................................................21
Mean ....................................................................................................................................................21
Median .................................................................................................................................................22
Mode ....................................................................................................................................................23
Geometric mean ..................................................................................................................................23
Harmonic Mean....................................................................................................................................24
Measures of Variability...........................................................................................................................24
Variance...............................................................................................................................................25
Standard Deviation ..............................................................................................................................27
Squared Deviation ...............................................................................................................................27

5
Average Absolute Deviation ................................................................................................................27
Median Absolute Deviation ..................................................................................................................28
Range ..................................................................................................................................................28
Interquartile Range ..............................................................................................................................28
Ranking..................................................................................................................................................29
MIN, MAX, SMALL, LARGE ................................................................................................................30
RANK ...................................................................................................................................................30
RANK.AVG ..........................................................................................................................................31
PERCENTILE ......................................................................................................................................31
PERCENTILE.EXC..............................................................................................................................32
PERCENTRANK and PERCENTRANK.EXC......................................................................................33
QUARTILE and QUARTILE.EXC ........................................................................................................33
Shape of the Distribution........................................................................................................................34
Symmetry and Skewness ....................................................................................................................34
Kurtosis ................................................................................................................................................35
Graphic Illustrations .............................................................................................................................35
Descriptive Statistics Data Analysis Tool...............................................................................................36
Graphical Representations of Data........................................................................................................38
Frequency Tables ................................................................................................................................38
Histogram.............................................................................................................................................39
Box Plots..............................................................................................................................................39
Outliers...................................................................................................................................................43
Missing Data ..........................................................................................................................................43
Chapter 4 Distributions ............................................................................................................................45
Discrete Distributions.............................................................................................................................45
Continuous Distributions........................................................................................................................46
Excel Distribution Functions ..................................................................................................................47

6
Chapter 5 Normal Distribution.................................................................................................................51
Normal Distribution ................................................................................................................................51
Basic Concepts....................................................................................................................................51
Excel Functions....................................................................................................................................52
Example...............................................................................................................................................52
Standard Normal Distribution...............................................................................................................53
Lognormal Distribution.........................................................................................................................54
Sample Distributions..............................................................................................................................54
Basic Concepts....................................................................................................................................54
Central Limit Theorem .........................................................................................................................55
Hypothesis Testing ..............................................................................................................................55
Chapter 6 Binomial Distribution ..............................................................................................................59
Basic Concepts ......................................................................................................................................59
Excel Functions......................................................................................................................................60
Hypothesis Testing ................................................................................................................................61
Chapter 7 Student’s t Distribution...........................................................................................................62
Basic Concepts ......................................................................................................................................62
Excel Functions......................................................................................................................................63
Hypothesis Testing ................................................................................................................................64
One-Sample Testing............................................................................................................................64
Testing Two Independent Samples .....................................................................................................66
Testing Paired Samples.......................................................................................................................70
Using the T.TEST Function .................................................................................................................73
Confidence Intervals ..............................................................................................................................74
Chapter 8 Chi-square and F Distribution................................................................................................75
Chi-square Distribution ..........................................................................................................................75

7
Basic Concepts....................................................................................................................................75
Excel Functions....................................................................................................................................75
Contingency Tables .............................................................................................................................76
Independence Testing .........................................................................................................................78
F Distribution..........................................................................................................................................79
Basic Concepts....................................................................................................................................79
Excel Functions....................................................................................................................................79
Example...............................................................................................................................................80
Excel Test Function .............................................................................................................................81
Chapter 9 Analysis of Variance ...............................................................................................................82
Introduction ............................................................................................................................................82
One-way ANOVA ...................................................................................................................................82
One-way ANOVA Example..................................................................................................................82
Basic Concepts....................................................................................................................................83
Analysis................................................................................................................................................85
Follow-up Analysis...............................................................................................................................86
Levene’s Test.......................................................................................................................................87
Factorial ANOVA....................................................................................................................................88
Example...............................................................................................................................................88
Basic Concepts....................................................................................................................................89
Analysis................................................................................................................................................91
ANOVA with Repeated Measures .........................................................................................................93
Basic Concepts....................................................................................................................................93
Example...............................................................................................................................................94
Analysis................................................................................................................................................94
Assumptions ........................................................................................................................................96

8
Follow-up Analysis...............................................................................................................................97
Chapter 10 Correlation and Covariance .................................................................................................99
Basic Definitions ....................................................................................................................................99
Scatter Diagram .....................................................................................................................................99
Excel Functions....................................................................................................................................102
Excel Functions with Missing Data ......................................................................................................102
Data Analysis Tools .............................................................................................................................106
Chapter 11 Linear Regression ...............................................................................................................108
Linear Regression Analysis .................................................................................................................108
Regression Line.................................................................................................................................108
Residuals ...........................................................................................................................................110
Model Fit ............................................................................................................................................111
Multiple Regression Analysis...............................................................................................................113

9
The Story behind the Succinctly Series
of Books
Daniel Jebaraj, Vice President
Syncfusion, Inc.
taying on the cutting edge
As many of you may know, Syncfusion is a provider of software components for the
Microsoft platform. This puts us in the exciting but challenging position of always
being on the cutting edge.
Whenever platforms or tools are shipping out of Microsoft, which seems to be about every other
week these days, we have to educate ourselves, quickly.
Information is plentiful but harder to digest
In reality, this translates into a lot of book orders, blog searches, and Twitter scans.
While more information is becoming available on the Internet and more and more books are
being published, even on topics that are relatively new, one aspect that continues to inhibit us is
the inability to find concise technology overview books.
We are usually faced with two options: read several 500+ page books or scour the web for
relevant blog posts and other articles. Just as everyone else who has a job to do and customers
to serve, we find this quite frustrating.
The Succinctly series
This frustration translated into a deep desire to produce a series of concise technical books that
would be targeted at developers working on the Microsoft platform.
We firmly believe, given the background knowledge such developers have, that most topics can
be translated into books that are between 50 and 100 pages.
This is exactly what we resolved to accomplish with the Succinctly series. Isn’t everything
wonderful born out of a deep desire to change things for the better?
S

10
The best authors, the best content
Each author was carefully chosen from a pool of talented experts who shared our vision. The
book you now hold in your hands, and the others available in this series, are a result of the
authors’ tireless work. You will find original content that is guaranteed to get you up and running
in about the time it takes to drink a few cups of coffee.
Free forever
Syncfusion will be working to produce books on several topics. The books will always be free.
Any updates we publish will also be free.
Free? What is the catch?
There is no catch here. Syncfusion has a vested interest in this effort.
As a component vendor, our unique claim has always been that we offer deeper and broader
frameworks than anyone else on the market. Developer education greatly helps us market and
sell against competing vendors who promise to “enable AJAX support with one click” or “turn the
moon to cheese!”
Let us know what you think
If you have any topics of interest, thoughts or feedback, please feel free to send them to us at
succinctly-series@syncfusion.com.
We sincerely hope you enjoy reading this book and that it helps you better understand the topic
of study. Thank you for reading.
Please follow us on Twitter and “Like” us on Facebook to help us spread the
word about the Succinctly series!

11
Chapter 1 Introduction
What is Statistics?
Statistics is a field of study that has two principal objectives:
 Collecting, organizing, interpreting, and presenting data (Descriptive Statistics)
 Using mathematical techniques, especially from probability theory, to make inferences
and/or predictions based on a sample of experimentally observed data (Inferential
Statistics)
Descriptive Statistics involves calculating the mean, median, variance, standard deviation,
and other properties of the data, and presenting this information in ways that make the data
more meaningful such as histograms, box plots, etc.
Inferential Statistics involves analyzing data and inferring characteristics of a general
population based on the same properties of a sample taken from the population. This is what
gives the field of statistics its power since, with a relatively small amount of data, we are able to
make significant assertions even though such inferences are not 100 percent certain but, rather,
are probabilistic in nature.
Why do Statistics in Excel?
There are a number of commonly used, powerful tools for carrying out statistical analyses. The
most popular of these are Statistical Package for the Social Sciences (SPSS), Statistical
Analysis System (SAS), Stata, Minitab and R. Many people choose to use Excel as their
principal analysis tool or as a complement to one of these tools for any of the following reasons:
 It is widely available and so many people already know how to use it
 It is not necessary to incur the cost of yet another tool (as some of the popular tools are
quite expensive)
 It is not necessary to learn new methods of manipulating data and drawing graphs
 It provides numerous built-in statistical functions and data analysis tools
 It is much easier to see what is going on since, unlike the more commonly used
statistical analysis tools, very little is hidden from the user
 It provides the user with a lot of control and flexibility
This makes Excel an ideal tool for quick analyses and even some serious systematic analyses.
However, it has two major shortcomings:
 Many people are not very familiar with the statistical capabilities that are built into Excel
 Excel was not designed to be a comprehensive statistical package and so many
commonly used statistical tests are not provided

12
We will address the first of these shortcomings in this book. We will also give you some
techniques for extending the built-in statistical capabilities included in Excel, which partially
addresses its second shortcoming.
Where Excel provides the statistical analysis capabilities that you need (and, fortunately, many
of the most commonly used tests are included in Excel), it is a great tool to use for the reasons
previously mentioned. Where it does not, there are a number of software packages which
extend Excel’s built-in capabilities. Once such software package is one that I have developed
called Real Statistics. You can download it for free here.
Where is Statistics Used?
Statistics plays a central role in research in the social sciences, pure sciences, and medicine. A
simplified view of experimental research is as follows:
 You make some observations about the world, and then create a theory consisting of a
hypothesis and possible alternative hypotheses that try to explain the observations you
have made
 You then test your theory by conducting experiments. Such experiments include
collecting data, analyzing the results, and coming to some conclusions about how well
your theory holds up
 You iterate this process, observing more about the world and improving your theory
Statistics also plays a major role in decision making for business and government, including
marketing, strategic planning, manufacturing, and finance.
Statistics is a discipline which is concerned with the collection and analysis of data based on a
probabilistic approach. Theories about a general population are tested on a smaller sample and
conclusions are made about how well properties of the sample extend to the population at
large.
What Will This Book Cover?
This book aims to show the reader how to perform statistical analyses using Excel. In particular,
it describes:
 Basic statistical tests and analyses
 Methods for carrying out these tests and analyses using Excel
 Numerous examples
We don’t have space to cover all of the aspects of statistics in Excel but, as you will see, we will
try to give you a good idea of how to do important statistical analyses in Excel.

13
Chapter 2 Excel Environment
The Basics
For the most part, we will assume that the reader is already familiar with Excel. We don’t have
the space here to provide a review of Excel. Instead, we will touch on a few of the capabilities
that are not as commonly used but that are essential for doing statistical analysis in Excel. We
will focus on the most recent versions of Excel in the Windows environment, namely, Excel
2007, 2010, and 2013—although most of what is written will also apply to earlier versions of
Excel in Windows as well as Excel for the Macintosh (especially Excel 2011).
User Interface
Excel works with files called workbooks. Each workbook contains one or more spreadsheets
called worksheets. Each worksheet consists of cells that are organized in a rectangular grid.
The rows of a worksheet are labeled with a number, and the columns are labeled with a letter or
series of letters. The first row is labeled 1, the next row is labeled 2, and so on. The first column
is labeled A, the next column is labeled B, and so on. The column after Z is labeled AA and then
the next columns are labeled AB, AC, and so on. The column after AZ is BA, and then the next
columns are labeled BB, BC, and so on. The column after ZZ is labeled AAA, and then the next
columns are labeled AAB, AAC, and so on.
Each worksheet in Excel 2007/2010/2013 can contain up to 1,048,578 rows and 16,384
columns (i.e., column A through XFD). Each cell can contain a number, text, truth value or a
formula. Cells have both a value and an address (e.g., the cell with address B5 lays in the fifth
row and second column).
A cell range consists of all of the cells in a rectangular area on the worksheet (e.g., the range
B6:C10 consists of the 10 cells in the rectangle whose opposite corners are the cells B6 and
C10. The formula =SUM(B6:C10) has a value which is the sum of the values of all the cells in
the range B6:C10). Formulas can use either absolute addressing, relative addressing, or a
combination of both.
The rectangular grid is only one part of the Excel user interface (UI). The entire UI is as depicted
in Figure 1:

14
Figure 1: Excel user interface
This is the layout that is used in Excel 2007. The layout used in Excel 2010 and Excel 2013 is
almost identical.
The Ribbon Tabs are the top-level menu items. In Figure 1, these are Home, Insert, Page
Layout, Formulas, etc. To access most capabilities in Excel, you click one of these ribbon tabs.
For each tab, a different ribbon will be displayed.
Each ribbon consists of a collection of Excel capabilities that are organized into groups that
correspond to some ribbon tab. For example, the Home ribbon (displayed in Figure1) is
organized into the Clipboard, Font, Alignment, Number, etc. groups. Each group consists of
one or more icons that correspond to some capabilities in Excel. For example, to center the
content of a cell in a worksheet, click that cell and then click the center icon in the Alignment
group on the Home ribbon.
We use the following abbreviation for this sequence of steps: Home > Alignment|Center.
There are also shortcuts for some icons (e.g., to center the contents of a cell, you can click that
cell and then enter Ctrl-E).

15
Array Formulas
Most worksheet formulas used in Excel return a single value that is assigned to the cell which
contains the formula. Excel also allows you to define a formula that assigns values to a range of
cells at the same time. These are called array formulas and they will be used quite often in the
rest of this book.
Example: Calculate the revenues for each item in the worksheet of Figure 2:
Figure 2: Array formulas
Given that the revenue for each item is unit price times quantity, we can enter the formula
=B4*C4 in cell D4 and then copy this formula into cells D5, D6, and D7 (e.g., by clicking cell D4,
pressing Ctrl-C, highlighting the range D5:D7, and pressing Ctrl-V, or by highlighting the range
D4:D7 and pressing Ctrl-D).
Another way to do this is via an array formula using the following steps:
1. Highlight the range D4:D7
2. Enter the array formula =B4:B7*C4:C7
3. Press Ctrl-Shft-Enter
It is essential to press Ctrl-Shft-Enter (step three) and not simply Enter for an ordinary formula.
Note that the formula that appears in the formula bar is {=B4:B7*C4:C7}. The curly brackets
indicate that this is an array formula. If the range B4:B7 is given the name UnitPrice and C4:C7
is given the name Quantity, then the array formula can be entered as =UnitPrice*Quantity (step
two).
The array formula appears in all four cells in the range D4:D7. To make changes to the formula,
you must edit the entire range, not just one, two, or three of these cells. Similarly, you can’t copy
or delete a part of the range but, rather, you must copy or delete the entire range. If you attempt
to modify a part of the range, you will receive an error message. If you get stuck and see a
series of error messages, you just need to press Esc to recover.
You can erase a range that contains an array formula by highlighting the entire range and
pressing Delete. You can write over the array function, replacing it by a value or another
formula. The important thing is to use the entire range and not a part of the range.

16
Note, too, that you can also use array formulas such as {=SUM(B4:B7*C4:C7)}. This returns the
value, which is the sum of the revenues of the four types of equipment. Even though this
formula returns a single value, and so may be placed in a single cell such as D8, it must be
entered as an array formula (since the formula contains an embedded array formula). This
means that you need to type =SUM(B4:B7*C4:C7) and then press Ctrl-Shft-Enter. If you forget
to press Ctrl-Shft-Enter and only press Enter, you will get an error message.
Array Functions
A few of Excel’s built-in functions are array functions where the output of the function is an array
(i.e., a range). These functions are managed as previously described for array formulas.
Example: Change the data range in columns A and B of Figure 3 into an equivalent row range:
Figure 3: Array function
This can be accomplished by means of Excel’s TRANSPOSE array function using the following
steps:
1. Highlight the output range D3:I4
2. Enter the array formula =TRANSPOSE(A3:B8)
3. Press Ctrl-Shft-Enter
Note that the output range (step one) must be of the right size. In this case, the input range is
six rows by two columns, and so the output range must be two rows by six columns. As for array
formulas, the formula bar contains the array formula enclosed in curly brackets. Once again, it is
important to press Ctrl-Shft-Enter.
Data Analysis Tools
Excel provides a number of data analysis tools which are accessible via Data > Analysis|Data
Analysis.

17
If this option is not visible, you may need to first install Excel’s analysis toolpack. This is done by
selecting Office Button > Excel Options > Add-Ins in Excel 2007 or File > Help|Options >
Add-Ins in Excel 2010/2013, and then clicking Go at the bottom of the window. Next, you select
the Analysis ToolPak option in the dialog box that appears, and click OK. You will then be able
to access the data analysis tools.
After selecting Data > Analysis|Data Analysis, you will be presented with the dialog box in
Figure 4:
Figure 4: Data analysis dialog box
You can now select any one of the following options useful for statistical analysis:
 ANOVA: Single Factor
 ANOVA: Two-Factor with Repetition
 ANOVA: Two-Factor without Repetition
 Correlation
 Covariance
 Descriptive Statistics
 Exponential Smoothing
 F-Test: Two Sample for Variances
 Histogram
 Random Number Generation
 Rank and Percentile
 Regression
 Sampling
 t-Test: Paired Two Sample for Means
 t-Test: Two-Sample Assuming Equal Variance
 t-Test: Two-Sample Assuming Unequal Variance
 z-Test: Two-Sample for Means
Each of these options represents a data analysis tool that will be described in this book. Now
suppose, for example, that you choose ANOVA: Single Factor. You will now be presented with
the dialog box shown in Figure 5:

18
Figure 5: Dialog box for ANOVA: Single Factor
The Input Range consists of the Excel range within which the data elements to be analyzed are
stored. Suppose, for example, this data consists of a 4×8 array representing four treats as in
Figure 6:
Figure 6: Sample input range
In this case, you insert the range B2:E9 in the Input Range field (of the dialog box of Figure 5)
and select Columns. If you had assigned a name (e.g., Study1) to the range B2:E9, then you
could put this name in the Input Range field instead of B2:E9.
Alternatively, you could insert B1:E9 in the Input Range field and select the Labels in First
Row check box in the dialog box to indicate that you have included the column headings in the
data range. Note that the participant numbers (in column A) are not used.
If the data were arranged where the treatments are listed by row instead of column, then you
would select Rows and you could select the Labels in First Column check box.

19
The Alpha value (which will be described later in the book) is set to 0.05 by default, although
you have the option to change this to 0.01 or some other value.
You can now select New Worksheet Ply from the Output options (and leave the data field
blank). In this case, a new worksheet will be created (in the tab prior to the current one) and the
ANOVA report will be placed in this worksheet starting at cell A1. You can then copy the results
to the current worksheet (or anywhere else you like).
Alternatively, you can select Output Range or New Workbook, and put the report in some
specific output range that you choose or in a new workbook.

20
Chapter 3 Descriptive Statistics
Basic Terminology
Before exploring Descriptive Statistics, we define the following statistical concepts which will
be used throughout the rest of the book.
Sample: A subset of the data from the population which we analyze in order to learn about the
population. A major objective in the field of statistics is to make inferences about a population
based on properties of the sample.
Random sample: A sample in which each member of the population has an equal chance of
being included and in which the selection of one member is independent from the selection of all
other members.
Random variable: A variable which represents values from a random sample. We will use
letters at the end of the alphabet, especially 𝑥, 𝑦, and 𝑧 as random variables.
Independent random variable: A random variable that is chosen and then measured or
manipulated by the researcher in order to study some observed behavior.
Dependent random variable: A random variable whose value depends on the value of one or
more independent variables.
Discrete random variable: A random variable that takes a discrete set of values. Such random
variables generally take a finite set of values (e.g., heads or tails, brands of cars, etc.) but they
can also include random variables that take a countable set of values (e.g., 0, 1, 2, 3, etc.).
Continuous random variable: A random variable that takes an infinite number of values even
in a finite interval (e.g., the height of a building, a temperature, etc.).
Statistic: A quantity which is calculated from a sample and is used to estimate a corresponding
characteristic (i.e., parameter) about the population from which the sample is drawn.
In the rest of this chapter, we define some statistics which are commonly used to characterize
data. In particular, we define metrics of central tendency (e.g., mean and median), variability
(e.g., variance and standard deviation), symmetry (i.e., skewness) and peakedness (i.e.,
kurtosis).
We also provide some important ways of graphically describing data and distributions including
histograms and box plots.

21
Measures of Central Tendency
We consider a random variable 𝑥 and a data set S= {𝑥1, … , 𝑥 𝑛} of size 𝑛 which contains possible
values of 𝑥. The data can represent either the population being studied or a sample drawn from
the population.
We seek a single measure (i.e., a statistic) which somehow represents the center of the entire
data set S. The commonly used measures of central tendency are the mean, median, and
mode. Besides the normally studied mean (also called the arithmetic mean), we also touch
briefly on two other types of mean: the geometric mean and the harmonic mean.
These measures are summarized in Table 1 where R is an Excel range which contains the data
elements in the sample or population S:
Statistic Excel 2007 Excel 2010/2013
Arithmetic Mean AVERAGE(R) AVERAGE(R)
Median MEDIAN(R) MEDIAN(R)
Mode MODE(R)
MODE.SNGL(R),
MODE.MULT(R)
Geometric Mean GEOMEAN(R) GEOMEAN(R)
Harmonic Mean HARMEAN(R) HARMEAN(R)
Table 1: Measures of central tendancy
Mean
We begin with the mean, the most commonly used measure of central tendency. The mean
(also called the arithmetic mean) of the data set S is defined as follows:
The mean is calculated in Excel using the worksheet function AVERAGE.

22
Example: The mean of S = {5, 2, -1, 3, 7, 5, 0, 2} is (2 + 5 – 1 + 3 + 7 + 5 + 0 + 2) / 8 = 2.875.
We achieve the same result by using the formula =AVERAGE(C3:C10) in Figure 7:
Figure 7: Examples of central tendency in Excel
When the data set S is a population, the Greek symbol µ is used for the mean. When S is a
sample, then the symbol 𝑥̅ is used.
Median
If you arrange the data in S in increasing order, the middle value is the median. When S has an
even number of elements, there are two such values. The average of these two values is the
median.
The median is calculated in Excel using the worksheet function MEDIAN.
Example: The median of S = {5, 2, -1, 3, 7, 5, 0} is 3 since 3 is the middle value (i.e., the 4th of
7 values) in -1, 0, 2, 3, 5, 5, 7. We achieve the same result by using the formula
=MEDIAN(B3:B10) in Figure 7.
Note that each of the functions in Figure 7 ignores any non-numeric values including blanks.
Thus, the value obtained from =MEDIAN(B3:B10) is the same as that for =MEDIAN(B3:B9).
The median of S = {5, 2, -1, 3, 7, 5, 0, 2} is 2.5 since 2.5 is the average of the two middle values
2 and 3 of -1, 0, 2, 2, 3, 5, 5, 7.This is the same result as =MEDIAN(C3:C10) in Figure 7.

23
Mode
The mode of the data set S is the value of the data element that occurs most often.
Example: The mode of S = {5, 2, -1, 3, 7, 5, 0} is 5 since 5 occurs twice, more than any other
data element. This is the result we obtain from the formula =MODE(B3:B10) in Figure 7. When
there is only one mode, as in this example, we say that S is unimodal.
If S = {5, 2, -1, 3, 7, 5, 0, 2}, the mode of 𝑆 consists of both 2 and 5 since they each occur twice,
more than any other data element. When there are two modes, as in this case, we say
that S is bimodal.
The mode is calculated in Excel by the formula MODE. If range R contains unimodal data, then
MODE(R) returns this unique mode. But when R contains data with more than one mode,
MODE(R) returns the first of these modes. For our example, this is 5 (since 5 occurs before 2,
the other mode, in the data set). Thus, MODE(C3:C10) = 5.
If all of the values occur only once, then MODE returns an error value. This is the case for S =
{5, 2, -1, 3, 7, 4, 0, 6}. Thus, MODE(D3:D10) returns the error value #N/A.
Excel 2010/2013 provide an array function MODE.MULT which is useful for multimodal data by
returning a vertical list of modes. When we highlight C19:C20, then enter the array formula
=MODE.MULT(C3: C10), and then press Ctrl-Alt-Enter, we see that both modes are displayed.
Excel 2010/2013 also provides the function MODE.SNGL which is equivalent to MODE.
Geometric mean
The geometric mean of the data set S is defined as follows:
This statistic is commonly used to provide a measure of average rate of growth as described in
the next example.
The geometric mean is calculated in Excel using the worksheet function GEOMEAN.
Example: Suppose the sales of a certain product grow five percent in the first two years and 10
percent in the next two years. What is the average rate of growth over the four years?
If sales in year 1 are $1, then sales at the end of the four years are (1 + .05)(1 + .05)(1 + .1)(1 +
.1) = 1.334. The annual growth rate 𝑟 is that amount such that (1 + 𝑟)4
.= 1.334. Thus, 𝑟 is equal
to . 3341 4⁄
.= .0747.
The same annual growth rate of 7.47 percent can be obtained in Excel using the formula
GEOMEAN(H7:H10) – 1 = .0747.

24
Harmonic Mean
The harmonic mean of the data set S is calculated as follows:
This statistic can be used to calculate an average speed as described in the next example.
The harmonic mean is calculated in Excel using the worksheet function HARMEAN.
Example: If you go to your destination at 50mph and return at 70mph, what is your average rate
of speed?
Assuming the distance to your destination is 𝑑, the time it takes to reach your destination is
𝑑/50 hours and the time it takes to return is 𝑑/70 hours, for a total of 𝑑/50 + 𝑑/70 hours. Since
the distance for the whole trip is 2 𝑑, your average speed for the whole trip is as follows:
This is equivalent to the harmonic mean of 50 and 70, and so can be calculated in Excel as
HARMEAN(50,70) which is HARMEAN(G7:G8) from Figure 7.
Measures of Variability
Let’s now consider a random variable 𝑥 and a data set S = {𝑥1, … , 𝑥 𝑛} of size 𝑛 which contains
possible values of 𝑥. The data set can represent either the population being studied or a sample
drawn from the population.
The mean is the statistic used most often to characterize the center of the data in S. We now
consider the following commonly used measures of variability of the data around the mean,
namely, the standard deviation, variance, squared deviation, and average absolute
deviation.
In addition, let’s also explore three other measures of variability that are not linked to the mean,
namely, the median absolute deviation, range, and interquartile range.
Of these statistics, the variance and standard deviation are the most commonly employed.
These measures are summarized in Table 2 where R is an Excel range which contains the data
elements in the sample or population S:

25
Statistic Excel 2007 Excel 2010/2013 Symbol
Population Variance VARP(R) VAR.P(R) 𝜎2
Sample Variance VAR(R) VAR.S(R) 𝑠2
Population Standard
Deviation
STDEVP(R) STDEV.P(R) 𝜎
Sample Standard
Deviation
STDEV(R) STDEV.S(R) 𝑠
Squared Deviation DEVSQ(R) DEVSQ(R) 𝑆𝑆
Average Absolute
Deviation
AVEDEV(R) AVEDEV(R) 𝐴𝐴𝐷
Median Absolute Deviation See below See below 𝑀𝐴𝐷
Range MAX(R)-MIN(R) MAX(R)-MIN(R)
Interquartile Range
=QUARTILE(R, 3)
– QUARTILE(R, 1)
See below 𝐼𝑄𝑅
Table 2: Measures of variability
Variance
The variance is a measure of the dispersion of the data around the mean. Where S represents
a population, the population variance (symbol 𝝈 𝟐
) is calculated from the population mean µ as
follows:
Where S represents a sample, the sample variance (symbol 𝒔 𝟐
) is calculated from the sample
mean 𝑥̅ as follows:

26
The sample variance is calculated in Excel using the worksheet function VAR.S. The population
variance is calculated in Excel using the function VAR.P. In versions of Excel prior to Excel
2010, these functions are called VAR and VARP.
Example: If S = {2, 5, -1, 3, 4, 5, 0, 2} represents a population, then the variance = 4.25. This is
calculated as follows:
First, the mean = (2+5-1+3+4+5+0+2)/8 = 2.5, and so the squared deviation SS = (2–2.5)2
+
(5–2.5)2
+ (-1–2.5)2
+ (3–2.5)2
+ (4–2.5)2
+ (5–2.5)2
+ (0–2.5)2
+ (2–2.5)2
= 34. Thus, the
variance = SS/n = 34/8 = 4.25.
If instead S represents a sample, then the mean is still 2.5 but the variance = SS/(n–1) = 34/7 =
4.86. These can be calculated in Excel by the formulas VAR.P(B3;B10) and VAR.S(B3:B10) as
shown in Figure 8:
Figure 8: Examples of measures of variability in Excel

27
Standard Deviation
The standard deviation is the square root of the variance. Thus, the population and sample
standard deviations are calculated respectively as follows:
The sample standard deviation is calculated in Excel using the worksheet function STDEV.S.
The population variance is calculated in Excel using the function STDEV.P. In versions of Excel
prior to Excel 2010, these functions are called STDEV and STDEVP.
Example: If S = {2, 5, -1, 3, 4, 5, 0, 2} is a population, then the standard deviation = the square
root of the population variance, that is: √4.25 = 2.06. If S is a sample, then the sample standard
deviation = square root of the sample variance = √4.86 = 2.20.
These are the results of the formulas STDEV.P(B3:B10) and STDEV.S(B3:B10) as shown in
Figure 8.
Squared Deviation
The squared deviation (symbol SS for sum of squares) is most often used in ANOVA and
related tests. It is calculated as follows:
The squared deviation is calculated in Excel using the worksheet function DEVSQ.
Example: If S = {2, 5, -1, 3, 4, 5, 0, 2}, then the squared deviation = 34. This is the same as the
result of the formula DEVSQ(B3:B10) as shown in Figure 8.
Average Absolute Deviation
The average absolute deviation (AAD) of data set S is calculated as follows:
The average absolute deviation is calculated in Excel using the worksheet function AVEDEV.

28
Example: If S = {2, 5, -1, 3, 4, 5, 0, 2}, then the average absolute deviation = 1.75. This is the
same as the result of the formula AVEDEV(B3:B10) as shown in Figure 8.
Median Absolute Deviation
The median absolute deviation (MAD) of data set S is calculated as follows, where 𝑥̃ =
median of the data elements in S:
If R is a range which contains the data elements in S, then the MAD of S can be calculated in
Excel by the array formula:
=MEDIAN(ABS(R-MEDIAN(R)))
Even though the value is presented in a single cell, it is essential that you press Ctrl-Shft-
Enter to obtain the array value, otherwise the result won’t come out correctly.
Example: If S = {2, 5, -1, 3, 4, 5, 0, 2}, then the median absolute deviation = 2 since S = {-1, 0,
2, 2, 3, 4, 5, 5} and so the median of S is (2+3)/2 = 2.5. Thus, MAD = the median of {3.5, 2.5,
0.5, 0.5, 0.5, 1.5, 2.5, 2.5} = {0.5, 0.5, 0.5, 1.5, 2.5, 2.5, 2.5, 3.5}, that is, (1.5+2.5)/2 = 2.
This metric is less affected by extremes in the data (i.e., the tails) because the data in the tails
have less influence on the calculation of the median than they do on the mean.
Range
The range of a data set S is a crude measure of variability and consists simply of the difference
between the largest and smallest values in S.
If R is a range which contains the data elements in S, then the range of S can be calculated in
Excel by the formula:
=MAX(R) – MIN(R)
Example: If S = {2, 5, -1, 3, 4, 5, 0, 2}, then the range = 5 – (-1) = 6
Interquartile Range
The interquartile range (IQR) of a data set S is calculated as the 75th
percentile of S minus the
25th
percentile. The IQR provides a rough approximation of the variability near the center of the
data in S.
If R is a range which contains the data elements in S, then the IQR of S can be calculated in
Excel by the formula:

29
=QUARTILE(R, 3) – QUARTILE(R, 1)
Example: If S = {2, 5, -1, 3, 4, 5, 0, 2}, then the IQR = 4.25 – 1.5 = 2.75
As we will see shortly, in Excel 2010/2011/2013 there are two ways to calculate the quartiles
using QUARTILE.INC, which is equivalent to QUARTILE and QUARTILE.EXC. Using this
second approach, we calculate the IQR for S to be 4.75 – 0.5 = 4.25.
The variance, standard deviation, average absolute deviation, and median absolute deviation
measure both the variability near the center and the variability in the tails of the distribution of
the data. The average absolute deviation and median absolute deviation do not give undue
weight to the tails. On the other hand, the range only uses the two most extreme points and the
IQR only uses the middle portion of the data.
Ranking
Table 3 summarizes the various ranking functions that are available in all versions of Excel for a
data set R. We will describe each of these functions in more detail in the rest of the section, plus
we will describe additional ranking functions that are only available in Excel 2010/2013.
Excel Function Definition Notes
MAX(R)
The largest value in R
(maximum)
=LARGE(R, 1)
MIN(R)
The smallest value in R
(minimum)
=SMALL(R, 1)
LARGE(R, n)
The nth largest element
in R
LARGE(R, COUNT(R)) =
MIN(R)
SMALL(R, n)
The nth smallest element
in R
SMALL(R, COUNT(R)) =
MAX(R)
RANK(c, R, d)
The rank of element c in
R
If d = 0 (or omitted), then the
ranking is in decreasing order;
otherwise, it is in increasing
order
PERCENTILE(R, p)
Element in R at the pth
percentile
0 ≤ p ≤ 1. If PERCENTILE(R,p)
= c, then p% of the data

30
Excel Function Definition Notes
elements are less than c
PERCENTRANK(R, c)
Percentage of elements
in R below c
If PERCENTRANK(R,c) = p,
then PERCENTILE(R, p) = c
QUARTILE(R, n)
Element in R at the nth
quartile
n = 0, 1, 2, 3, 4
Table 3: Ranking functions
MIN, MAX, SMALL, LARGE
MIN(R) = the smallest value in R and MAX(R) = the largest value in R.
SMALL(R, n) = the nth smallest value in R and LARGE(R, n) = the nth largest value in R.
Here, n can take on any value from 1 to the number of elements in R, that is, COUNT(R).
Example: For range R with data elements {4, 0, -1, 7, 5}:
MIN(R) = -1, MAX(R) = 7.
LARGE(R, 1) = 7, LARGE(R, 2) = 5, LARGE(R, 5) = -1
SMALL(R, 1) = -1, SMALL(R, 2) = 0, LARGE(R, 5) = 7
RANK
All versions of Excel contain the ranking function RANK. For versions of Excel starting with
Excel 2010, there is also the equivalent function RANK.EQ.
RANK(c, R, d) = the rank of data element c in R. If d = 0 (or is omitted), then the ranking is in
decreasing order, that is, a rank of 1 represents the largest data element in R. If d ≠ 0, then the
ranking is in increasing order, and so a rank of 1 represents the smallest element in R.
RANK(7, R) = RANK(7, R, 0) = 1
.
RANK(7, R, 1) = 5

31
RANK(0, R) = RANK(0, R, 0) = 4
RANK(0, R, 1) = 2
RANK.AVG
Excel’s RANK and RANK.EQ functions do not take care of ties. For example, if the range R
contains the values {1, 5, 5, 0, 8}, then RANK(5, R) = 2 because 5 is the second-highest ranking
element in range R. But 5 is also the third-highest ranking element in the range, and so for
many applications, it is useful to consider the ranking to be 2.5, namely, the average of 2 and 3.
Versions of Excel starting with Excel 2010 address this issue by providing a new
function RANK.AVG which takes the same arguments as RANK but returns the average of
equal ranks as previously described.
For versions of Excel prior to Excel 2010, you will need to use the following formula to take care
of ties in a similar fashion:
= RANK(c, R) + (COUNTIF(R, c) – 1) / 2
Example: Using the RANK.AVG function, find the ranks of the data in range E17:E23 of Figure
9:
Figure 9: Average ranking
The result is shown in column F of Figure 9. For example, the average rank of 8 (cell E21 or
E22) is 1.5 as shown in cell F21 (or F22) and calculated using the formula
=RANK.AVG(E21,E17:E23). If instead you want the ranking in the reverse order (where the
lowest value gets rank 1), then the results are shown in column G. This time using the formula
=RANK.AVG(E21,E17:E23,1), we see that the rank of 8 is 6.5 as shown in cell G21.
PERCENTILE
For any percentage p (that is, where 0 ≤ p ≤ 1 or equivalently 0% ≤ p ≤
100%), PERCENTILE(R, p) = the element at the pth percentile. This means that if
PERCENTILE(R, p) = c, then p% of the data elements in R are less than c.

32
If p = k/(n–1) for some integer value k = 0, 1, 2, ... n–1 where n = COUNT(R), then
PERCENTILE(R, p) = SMALL(R, k+1) = the k+1th element in R. If p is not a multiple of 1/(n-1),
then the PERCENTILE function performs a linear interpolation as described in the following
example.
Example: For range R with data elements {4, 0, -1, 7, 5}, R’s five data elements divide the
range into four intervals of size 25%, that is, 1/(5-1) = .25. Thus:
PERCENTILE(R, 0) = -1 (the smallest element in R)
PERCENTILE(R, .25) = 0 (the second smallest element in R)
PERCENTILE(R, .5) = 4 (the third smallest element in R)
PERCENTILE(R, .75) = 5 (the fourth smallest element in R)
PERCENTILE(R, 1) = 7 (the fifth smallest element in R)
For other values of p, we need to interpolate. For example:
PERCENTILE(R, .8) = 5 + (7 – 5) * (0.8 – 0.75) / 0.25 = 5.4
PERCENTILE(R, .303) = 0 + (4 – 0) * (0.303 – 0.25) / 0.25 = .85
Of course, Excel's PERCENTILE function calculates all these values automatically without you
having to figure things out.
PERCENTILE.EXC
Starting with Excel 2010, Microsoft introduced an alternative version of the percentile function
named PERCENTILE.EXC.
If n = COUNT(R), then for any integer k with 1 ≤ k ≤ n
PERCENTILE.EXC(R, k/(n+1)) = SMALL(R, k), that is, the kth smallest element in range R.
For 0 < p < 1, if p is not a multiple of 1/(n+1), then PERCENTILE.EXC(R, 𝑝) is calculated by
taking a linear interpolation between the corresponding values in R. For p < 1/(n+1) or p >
n/(n+1), no interpolation is possible and so PERCENTILE.EXC(R, p) returns an error.
Excel 2010/2013 also provides the new function PERCENTILE.INC which is equivalent to
PERCENTILE.
Example: Find the 0 – 100 percentiles in increments of 10% for the data in Figure 10 using both
PERCENTILE (or PERCENTILE.INC) and PERCENTILE.EXC:
Figure 10: Data for percentile calculations

33
The result is shown in Figure 11. For example, the score at the 60th
percentile is 58 (cell P10)
using the formula =PERCENTILE(B3:M3,O10), while it is 59 (cell S10) using the formula
=PERCENTILE.EXC(B3:M3,R10).
Figure 11: PERCENTILE vs. PERCENTILE.EXC
PERCENTRANK and PERCENTRANK.EXC
PERCENTRANK (or PERCENTRANK.INC) and PERCENTRANK.EXC are the inverses of
PERCENTILE and PERCENTILE.EXC. Thus, PERCENTRANK(R, c) = the value p such that
PERCENTILE(R, p) = c. Similarly, PERCENTRANK.EXC(R, c) = the value p such that
PERCENTILE.EXC(R, p) = c.
Example: Referring to Figures 10 and 11, we have:
PERCENTRANK(B3:M3,54) = .4
PERCENTRANK.EXC(B3:M3,S12) = .8
QUARTILE and QUARTILE.EXC
For any integer, n = 0, 1, 2, 3 or 4, QUARTILE(R, n) = PERCENTILE(R, n/4). If c is not an
integer but 0 ≤ c ≤ 4, then QUARTILE(R, c) = QUARTILE(R, INT(c)). Thus:
QUARTILE(R, 0) = PERCENTILE(R, 0) = MIN(R)
QUARTILE(R, 1) = PERCENTILE(R, .25)
QUARTILE(R, 2) = PERCENTILE(R, .5) = MEDIAN(R)
QUARTILE(R, 3) = PERCENTILE(R, .75)
QUARTILE(R, 4) = PERCENTILE(R, 1) = MAX(R)

34
QUARTILE(R, 0) = PERCENTILE(R, 0) = -1
QUARTILE(R, 1) = PERCENTILE(R, .25) = 0
QUARTILE(R, 4) = PERCENTILE(R, 1) = 7
Starting with Excel 2010, Microsoft introduced an alternative version of the quartile function
named QUARTILE.EXC. This function is defined so that QUARTILE.EXC(R, n)
= PERCENTILE.EXC(R, n/4). Excel 2010/2013 also provides the new function QUARTILE.INC
which is equivalent to QUARTILE.
Shape of the Distribution
Symmetry and Skewness
Skewness as a measure of symmetry. If the skewness of S is 0, then the distribution of the data
in S is perfectly symmetric. If the skewness is negative, then the graph is skewed to the left. If
the skewness is positive, then the graph is skewed to the right (see Figures 12 and 14 for
examples).
Excel calculates the skewness of sample S using the formula:
where 𝑥̅ is the mean and 𝑠 is the standard deviation of S. To avoid division by zero, this formula
requires that 𝑛 > 2.
When a distribution is symmetric, the mean = median; when the distribution is positively
skewed, mean > median, and when the distribution is negatively skewed, mean < median.
Excel provides the SKEW function as a way to calculate the skewness of S. For example, if R is
a range in Excel containing the data elements in S, then SKEW(R) = the skewness of S.
There is also a population version of the skewness given by the formula:
This version has been implemented in Excel 2013 using the function SKEW.P.

35
For versions of Excel prior to Excel 2013, you can use the formula:
SKEW(R)*(n–2)/SQRT(n(n–1))
instead of SKEW.P(R) where R contains the data in S = {𝑥1, … , 𝑥 𝑛} and n = COUNT(R).
Kurtosis
Kurtosis as a measure of peakedness or flatness. Positive kurtosis indicates a relatively
peaked distribution. Negative kurtosis indicates a relatively flat distribution.
Kurtosis is calculated in Excel as follows:
where 𝑥̅ is the mean and 𝑠 is the standard deviation of S. To avoid division by zero, this formula
requires that 𝑛 > 3.
Excel provides the KURT function as a way to calculate the kurtosis of S. For example, if R is a
range in Excel containing the data elements in S, then KURT(R) = the kurtosis of S.
Graphic Illustrations
We now look at an example of these concepts using the chi-square distribution:
Figure 12: Example of skewness and kurtosis
Figure 12 contains the graph of two chi-square distributions. We will study the chi-square
distribution later but, for now, note the following values for the kurtosis and skewness:
df = 5 df = 10
Kurtosis 5.4 4.2
Skewness 1.264911 0.894427

36
Figure 13: Comparison of skewness and kurtosis
The red curve ( 𝑑𝑓 = 10) is flatter than the blue curve ( 𝑑𝑓 = 5) which is reflected in the fact that
the kurtosis value of the blue curve is lower.
Both curves are asymmetric and skewed to the right (that is, the fat part of the curve is on the
left). This is consistent with the fact that the skewness for both is positive. But the blue curve is
more skewed to the right, which is consistent with the fact that the skewness of the blue curve is
larger.
Example: Calculate the skewness and kurtosis for the data in range A4:A16 of Figure 14:
Figure 14: Calculation of skewness and kurtosis in Excel
Descriptive Statistics Data Analysis Tool
Excel provides the Descriptive Statistics data analysis tool which produces a summary of the
key descriptive statistics for a data set.
Example: Produce a table of the most common descriptive statistics for the scores in column A
of Figure 15:

37
Figure 15: Output from descriptive statistics data analysis tool
The output from the tool is shown on the right side of Figure 15. To use the tool, select Data >
Analysis|Data Analysis and select Descriptive Statistics. A dialog box appears as in Figure
16:
Figure 16: Descriptive statistics dialog box
Now click Input Range and highlight the scores in column A (i.e., cells A3:A14). If you include
the heading as is done here, check Labels in first row. Since we want the output to start in cell
C3, click Output Range and insert C3 (or click cell C3). Finally, click Summary statistics and
press OK.

38
Note that, had we checked the Kth Largest check box, the output would have also contained
the value for LARGE(A4:A14,k) where k is the number we insert in the box to the right of the
label Kth Largest. Similarly, checking the Kth Smallest check box outputs SMALL(A4:A14,k).
The Confidence Interval for Mean option generates a confidence interval using the t
distribution as explained in Chapter 7.
Graphical Representations of Data
Frequency Tables
When you have a lot of data, it can be convenient to put the data in bins, usually of equal size,
and then create a graph of the number of data elements in each bin. Excel provides the
FREQUENCY(R1, R2) array function for doing this, where R1 = the input array and R2 = the bin
array.
To use the FREQUENCY array function, enter the data into the worksheet and then enter a bin
array. The bin array defines the intervals that make up the bins. For example, if the bin array =
10, 20, 30, then there are four bins, namely, data with values 𝑥 ≤ 10, data with values 𝑥 where
10 < 𝑥 ≤ 20, data with values 𝑥 where 20 < 𝑥 ≤ 30, and, finally, data with values 𝑥 > 30. The
FREQUENCY function simply returns an array consisting of the number of data elements in
each of the bins.
Example: Create a frequency table for the 22 data elements in the range A4:B14 of Figure 17
based on the bin array D4:D7 (the text "over 80" in cell D8 is not part of the bin array):
Figure 17: Example of the FREQUENCY function
To produce the output, highlight the range E4:E8 (i.e., a column range with one more cell than
the number of bins) and enter the formula:
=FREQUENCY(A4:B11,D4:D7)

39
Since this is an array formula, you must press Ctrl-Shft-Enter. Excel now inserts frequency
values in the highlighted range. Here, E4 contains the number of data elements in the input
range with value in the first bin (i.e., data elements whose value is ≤ 20). Similarly, E5 contains
the number of data elements in the input range with value in the second bin (i.e., data elements
whose value is > 20 and ≤ 40). The final output cell (E8) contains the number of data elements
in the input range with value > the value of the final bin (i.e., > 80 for this example).
Histogram
A histogram is a graphical representation of the output from the FREQUENCY function
previously described. You can use Excel’s chart tool to graph the output or, alternatively, you
can use the Histogram data analysis tool to directly accomplish this.
To use Excel’s Histogram data analysis tool, you must first establish a bin array (as for the
FREQUENCY function previously described) and then select the Histogram data analysis tool.
In the dialog box that is displayed, you then specify the input data (Input Range) and bin array
(Bin Range). You have the option to include the labels for these ranges (in which case you click
Labels).
For the data in Figure 17, the Input Range is A4:B14 and the Bin Range is D4:D7 (with the
Labels check box left unchecked). The output is displayed in Figure 18:
Figure 18: Example of the histogram data analysis
You should exercise caution whenever you create histograms in order to present the data in a
clear and accurate way. For most purposes, it is important that the intervals be equal in size
(except for an unbounded first and/or last interval). Otherwise, a distorted picture of the data
may be presented. To avoid this problem, equally-spaced intervals should generally be used.
Box Plots
Another way to characterize data is via a box plot. Specifically, a box plot provides a pictorial
representation of the following statistics: maximum, 75th-percentile, median (50th-percentile),
25th-percentile, and minimum.
As we shall see later, box plots are especially useful when comparing samples and testing
whether or not data is symmetric.
Bin Frequency
20 5
40 11
60 2
80 3
More 1
0
5
10
15
20 40 60 80 More
Frequency
Bin
Histogram
Frequency

40
Example: Create box plots for the three brands shown in range A3:C13 of Figure 19 using
Excel’s charting capabilities:
Figure 19: Box plot data
Select the range containing the data including the headings (A3:C13). Now, create the table in
the range E3:H8. The output in column F corresponds to the raw data from column A. Column G
corresponds to column B and column H corresponds to column C. In fact, once you construct
the formulas for the range F4:F8, you can fill in the rest of the table by highlighting the range
F4:H8 and pressing Ctrl-R.
The formulas for the cells in the range F4:F8 are as follows:
Cell Content
F4 =MIN(A4:A13)
F5 =QUARTILE(A4:A13,1)-F4
F6 =MEDIAN(A4:A13)-QUARTILE(A4:A13,1)
F7 =QUARTILE(A4:A13,3)-MEDIAN(A4:A13)
F8 =MAX(A4:A13)-QUARTILE(A4:A13,3)
Table 4: Formulas in the box plot table

41
Note that, alternatively, you can replace the QUARTILE function in Table 4 by the
QUARTILE.EXC function.
Once you have constructed the table, you can create the corresponding box plot as follows:
1. Select the data range E3:H7. Notice that the headings are included on the range but not
the last row
2. Select Insert > Charts|Column > Stacked Column
3. Select Design > Data|Switch Row/Column if necessary so that the 𝑥 axis represents
the brands
4. Select the lowest data series in the chart (i.e., Min) and set fill to No fill (and, if
necessary, set border color to No Line) to remove the lowest boxes. This is done by
right-clicking on any of the three Min data series boxes in the chart and selecting
Format Data Series. In the resulting dialog box, select Fill|No fill
5. Repeat the previous steps for the lowest visible data series (i.e., Q1-Min). That is, right-
click on the Q1-Min data series and select Format Data Series… > Fill|No fill.
Alternatively, right-click on the Q1-Min data series and press Ctl-Y.
6. With the Q1-Min data series still selected, select Layout > Analysis|Error Bars > More
Error Bar Options. In the resulting dialog box (Vertical Error Bars menu), click Minus
and Percentage, and insert a percentage error of 100%
7. Click the Q3-Med data series (the uppermost one) and select Layout > Analysis|Error
Bars > More Error Bar Options. In the resulting dialog box (Vertical Error Bars menu),
click Plus and Custom, and then click Specify Value. Now specify the range F8:H8,
that is, the last row of the table you previously created, in the dialog box that appears (in
the Positive Error Values field).
8. Remove the legend by selecting Layout > Labels|Legend > None
The resulting box plot is:
Figure 20: Box plot
0
200
400
600
800
1000
1200
1400
1600
1800
2000
Brand A Brand B Brand C

42
For each sample, the box plot consists of a rectangular box with one line extending upward and
another line extending downward (usually called whiskers). The box itself is divided into two
parts. In particular, the meaning of each element in the box plot is described in Table 5:
Element Meaning
Top of upper whisker Maximum value of the sample
Top of box 75% percentile of the sample
Line through the box Median of the sample
Bottom of the box 25% percentile of the sample
Bottom of the lower whisker Minimum of the sample
Table 5: Box plot elements
From the box plot (see Figure 20), we can see that the scores for Brand C tend to be higher
than for the other brands, and those for Brand B tend to be lower. We also see that the
distribution of Brand A is pretty symmetric, at least in the range between the 1st
and 3rd
quartiles,
although there is some asymmetry for higher values (or there is potentially an outlier). Brands B
and C look less symmetric. Because of the long upper whisker (especially with respect to the
box), Brand B may have an outlier.
The approach described in this section works perfectly for non-negative data. When a data set
has a negative value, you need to add -MIN(R) to all the data elements where R is the data
range containing the data. This will shift the y-axis up by –MIN(R). For example, if R ranges
from -10 to 20, then the range described in the chart will range from 0 to 30.
We can also convert the box plot to a horizontal representation of the data (as in Figure 21) by
clicking the chart and selecting Insert > Charts|Bar > Stacked Bar:

43
Figure 21: Horizontal box plot
Outliers
One problem that we face when analyzing data is the presence of outliers or data elements
that are much bigger or much smaller than the other data elements.
For example, the mean of the sample {2, 3, 4, 5, 6} is 4 while the mean of {2, 3, 4, 5, 60} is 14.4.
The appearance of the 60 completely distorts the mean in the second sample. Some statistics,
such as the median, are more resistant to such outliers. In fact, the median for both samples is
4. For this example, it is obvious that 60 is a potential outlier.
One approach for dealing with outliers is to throw away data that is either too big or too small.
Excel provides the TRIMMEAN function for dealing with this issue:
TRIMMEAN(R, p) – calculates the mean of the data in the range R after first throwing away
p% of the data, half from the top and half from the bottom. If R contains n data elements
and k = the largest whole number ≤ np/2, then the k largest items and the k smallest items
are removed before calculating the mean.
Example: Suppose R = {5, 4, 3, 20, 1, 4, 6, 4, 5, 6, 7, 1, 3, 7, 2}. Then TRIMMEAN(R, 0.2)
works as follows. Since R has 15 elements, k = INT(15 * .2 / 2) = 1. Thus, the largest element
(20) and the smallest element (1) are removed from R to get R′ = {5, 4, 3, 4, 6, 4, 5, 6, 7, 1, 3, 7,
2}. TRIMMEAN now returns the mean of this range, namely, 4.385 instead of the mean of R
which is 5.2.
Missing Data
Another problem we face when collecting data is that some data may be missing. For example,
when conducting a survey with 10 questions, perhaps some of the people who take the survey
don’t answer all 10 questions.
0 500 1000 1500 2000
Brand A
Brand B
Brand C

44
The principal way Excel deals with missing data is to ignore the missing data elements. For
example, in the various functions in this chapter (AVERAGE, VAR.S, RANK, etc.), any blank or
non-numeric cells are simply ignored.

45
Chapter 4 Distributions
Discrete Distributions
The (probability) frequency function 𝑓, also called the probability density function
(abbreviated pdf), of a discrete random variable 𝑥 is defined so that for any value 𝑡 in the
domain of the random variable:
𝑓(𝑡) = 𝑃(𝑥 = 𝑡)
where 𝑃(𝑥 = 𝑡) = the probability that 𝑥 assumes the value 𝑡.
The corresponding (cumulative) distribution function 𝐹(𝑥) is defined by
𝐹(𝑡) = ∑ 𝑓(𝑢)
𝑢≤𝑡
for any value 𝑡 in the domain of the random variable 𝑥.
For any discrete random variable defined over the domain S with frequency function 𝑓 and
distribution function 𝐹:
0 ≤ 𝑓(𝑡) ≤ 1 for all 𝑡 in 𝑆
∑ 𝑓(𝑡)
𝑡∈𝑆
= 1
𝐹(𝑡) = 𝑃(𝑥 ≤ 𝑡)
𝑃(𝑡1 < 𝑥 ≤ 𝑡2) = 𝐹(𝑡2) − 𝐹(𝑡1)
A frequency function can be expressed as a table or a bar chart as described in the following
example.
Example: Find the distribution function for the frequency function given in columns A and B in
Figure 22. Also, show the graph of the frequency and distribution functions:

46
Figure 22: Table defining the frequency and distribution functions
Given the frequency function 𝑓(𝑥) defined in the range B4:B11 of Figure 22, we can calculate
the distribution function 𝐹(𝑥) in the range C4:C11 by putting the formula =B4 in cell C4 and the
formula =B5+C4 in cell C5, and then copying this formula into cells C6 to C11 (for example, by
highlighting the range C5:C11 and pressing Ctrl-D).
From the table in Figure 22, we can create the charts in Figure 23:
Figure 23: Charts of frequency and distribution functions
Continuous Distributions
While for a discrete random variable 𝑥, the probability that 𝑥 assumes a value between 𝑎 and 𝑏
(exclusive) is given by ∑ 𝑓(𝑥)𝑎<𝑥<𝑏 , the frequency function 𝑓 of a continuous random variable
can assume an infinite number of values (even in a finite interval), and so we can’t simply sum
up the values in the ordinary way. For continuous variables, the equivalent formulation is that
the probability that 𝑥 assumes a value between 𝑎 and 𝑏 is given by 𝐴𝑟𝑒𝑎 𝑎<𝑥<𝑏 𝑓(𝑥), that is, the
area under the graph of 𝑦 = 𝑓(𝑥) bounded by the x-axis and the lines 𝑥 = 𝑎 and 𝑥 = 𝑏.
0.00
0.05
0.10
0.15
0.20
0.25
0.30
1 2 3 4 5 6 7 8
FrequencyFunction f(x)
0.00
0.20
0.40
0.60
0.80
1.00
1.20
1 2 3 4 5 6 7 8
DistributionFunction F(x)

47
Figure 24: Area under the curve 𝑦 = 𝑓(𝑥)
For a continuous random variable 𝑥, 𝑓 is a frequency function, also called the probability
density function (pdf), provided:
𝑓 is the frequency function, more commonly called the probability density function, for a
particular random variable, provided the area of the region indicated in Figure 24 represents the
probability that 𝑥 assumes a value between 𝑎 and 𝑏, inclusively. 𝑓 only takes non-negative
values and the area between the curve if 𝑦 = 𝑓(𝑥) and the x-axis is 1.
The corresponding (cumulative) distribution function 𝐹 is defined by:
𝐹(𝑥) = Area 𝑡<𝑥 𝑓(𝑡)
For any continuous random variable with distribution function 𝐹:
𝐹(𝑏) = 𝑃(𝑥 < 𝑏)
𝐹(𝑏) − 𝐹(𝑎) = 𝑃(𝑎 < 𝑥 < 𝑏)
Note that the probability that 𝑥 takes any particular value 𝑎 is not 𝑓(𝑎). In fact, for any specific
value 𝑎, the probability that 𝑥 takes the value 𝑎 is considered to be 0.
Essentially, the area under a curve is a way of summing when dealing with an infinite range of
values in a continuum. For those of you familiar with calculus, Area 𝑎<𝑥<𝑏 𝑓(𝑥) = ∫ 𝑓(𝑥)𝑑𝑥
𝑏
𝑎
.
Excel Distribution Functions
The following table provides a summary of the various distribution functions provided by Excel.
The entries labelled PDF/CDF specify the probability density function (where the last argument
for the function is FALSE), as well as the (left-tailed) cumulative distribution function (where the
last argument of the function is TRUE):
y = f(x)
x-axis
x = a x = b

48
Distribution PDF/CDF Inverse Right/Two-tailed Test
Beta BETA.DIST BETA.INV
Binomial BINOM.DIST BINOM.INV
Chi-square CHISQ.DIST CHISQ.INV
CHISQ.DIST.RT
CHISQ.INV.RT
CHISQ.TEST
Exponential EXPON.DIST
F F.DIST F.INV
F.DIST.RT
F.INV.RT
F.TEST
Gamma GAMMA.DIST GAMMA.INV
Hyper-
geometric
HYPGEOM.DIST
Lognormal LOGNORM.DIST LOGNORM.INV
Negative
Binomial
NEGBINOM.DIST
Normal NORM.DIST NORM.INV
Poisson POISSON.DIST
Standard
Normal
NORM.S.DIST NORM.S.INV Z.TEST
Student’s t T.DIST T.INV
T.DIST.RT
T.DIST.2T
T.INV.2T
T.TEST

49
Distribution PDF/CDF Inverse Right/Two-tailed Test
Weibull WEIBULL.DIST
Table 6: Excel distribution functions
The functions listed in Table 6 are not available for versions of Excel prior to Excel 2010. For
earlier releases of Excel, a similar collection of functions is available as summarized in Table 7:
Table 7: Distribution functions for versions of Excel prior to Excel 2010
Distribution CDF PDF Inverse Test
Beta BETADIST BETAINV
Binomial BINOMDIST BINOMDIST CRITBINOM
Chi Square CHIDIST CHIINV CHITEST
Exponential EXPONDIST EXPONDIST
F FDIST FINV FTEST
Gamma GAMMADIST GAMMAINV
Hypergeometric HYPGEOMDIST
Lognormal LOGNORMDIST LOGINV
Negative
Binomial
NEGBINOMDIST
Normal NORMDIST NORMDIST NORMINV
Poisson POISSON POISSON

50
Distribution CDF PDF Inverse Test
Standard
Normal
NORMSDIST NORMSINV ZTEST
Student’s T TDIST TINV TTEST
Weibull WEIBULL

51
Chapter 5 Normal Distribution
Normal Distribution
Basic Concepts
The probability density function of the normal distribution is defined by:
𝑓(𝑥) =
1
𝜎√2𝜋
𝑒
−
1
2
(
𝑥−𝜇
𝜎
)
2
where e is the constant 2.7183..., and π is the constant 3.1415…
The normal distribution is completely determined by the parameters μ and σ. It turns out that μ
is the mean of the normal distribution and σ is the standard deviation.
The graph of the normal distribution has the shape of a bell curve as shown in Figure 25:
Figure 25: Normal curve
Some characteristics of the normal distribution are as follows:
 Mean = median = mode = μ
 Standard deviation = σ
 Skewness = kurtosis = 0
The function is symmetric about the mean with inflection points (i.e., the points where the curve
changes from concave up to concave down, or from concave down to concave up) at 𝑥 = 𝜇 ± 𝜎.
Normal Curve
μ μ+σ μ+2σμ-3σ μ-2σ μ+3σμ-σ
68.3%
95.4%
99.7%

52
As can be seen from Figure 25, the area under the curve in the interval 𝜇 – 𝜎 < 𝑥 < 𝜇 + 𝜎 is
approximately 68.26 percent of the total area under the curve. The area under the curve in the
interval 𝜇 – 2𝜎 < 𝑥 < 𝜇 + 2𝜎 is approximately 95.44 percent of the total area under the curve,
and the area under the curve in the interval 𝜇 – 3𝜎 < 𝑥 < 𝜇 + 3𝜎 is approximately 99.74 percent
of the area under the curve.
As we shall see, the normal distribution occurs frequently and is very useful in statistics.
Excel Functions
Excel provides the following functions regarding the normal distribution:
NORM.DIST(x, μ, σ, cum) where cum takes the values TRUE and FALSE
NORM.DIST(x, μ, σ, FALSE) = probability density function value at x for the normal
distribution with mean μ and standard deviation σ
NORM.DIST(x, μ, σ, TRUE) = cumulative probability distribution function value at x for the
normal distribution with mean μ and standard deviation σ
NORM.INV(p, μ, σ) is the inverse of NORM.DIST(x, μ, σ, TRUE), that is,
NORM.INV(p, μ, σ) = the value x such that NORM.DIST(x, μ, σ, TRUE) = p
These functions are not available in versions of Excel prior to Excel 2010. Instead, these
versions of Excel provide NORMDIST which is equivalent to NORM.DIST and NORMINV, which
is equivalent to NORM.INV.
Example
Create a graph of the distribution of IQ scores using the Stanford-Binet Intelligence Scale.
This distribution is known to have a normal distribution with mean 100 and standard deviation
16. To create the graph, we first must create a table with the values of the probability density
function 𝑓(𝑥) for 𝑥 = 50, 51, …, 150. This table begins as shown on the left side of Figure 26:

53
Figure 26: Distribution of IQ scores
We only show the first 20 entries but, in total, there are 101 entries occupying the range
A4:B104. The values in column B are calculated using the NORM.DIST function. For example,
the formula in cell B4 is:
=NORM.DIST(A4,100,16,FALSE)
We can draw the chart on the right side of Figure 26 by highlighting the range B4:B104 and then
selecting Insert > Charts|Line. The display of the chart can be improved as described
elsewhere in this book.
Standard Normal Distribution
The standard normal distribution is a normal distribution with a mean of 0 and a standard
deviation of 1.
To convert a random variable 𝑥 with normal distribution with mean μ and standard deviation σ to
standard normal form, you use the following linear transformation:
𝑧 =
𝑥 − 𝜇
𝜎

54
The resulting random variable is called a z-score. Excel provides the following function for
calculating the value of 𝑧 from 𝑥, μ and σ:
STANDARDIZE(𝑥, 𝜇, 𝜎) =
𝑥−𝜇
𝜎
In addition, Excel provides the following functions regarding the standard normal distribution:
NORM.S.DIST(z, cum) = NORM.DIST(z,0,1,cum)
NORM.S.INV(p) = NORM.INV(p,0,1)
versions of Excel provide NORMSDIST(z) which is equivalent to NORM.S.DIST(z, TRUE), and
NORMSINV(p) which is equivalent to NORM.S.INV(p). For the pdf, the formula equivalent to
NORM.S.DIST(z, FALSE) is NORMDIST(z, 0, 1, FALSE).
Lognormal Distribution
A random variable 𝑥 is lognormally distributed provided the natural log of 𝑥, ln 𝑥, is normally
distributed.
Note that the lognormal distribution is not symmetric but is skewed to the right. If you have data
that is skewed to the right that fits the lognormal distribution, you may be able to access various
tests described elsewhere in this book that require data to be normally distributed.
Excel provides the following two functions:
LOGNORM.DIST(x, μ, σ, cum) = the lognormal cumulative distribution function with mean μ
and standard deviation σ at x if cum = TRUE, that is, NORMDIST(LN(x), μ, σ, cum), and the
probability density function of the lognormal distribution if cum = FALSE.
LOGNORM.INV(p, μ, σ) = EXP(NORM.INV(p, μ, σ)
versions of Excel provide LOGNORMDIST(x, μ, σ) which is equivalent to LOGNORM.DIST(x, μ,
σ, TRUE), and LOGINV(p, μ, σ) which is equivalent to LOGNORM.INV(p, μ, σ).
Sample Distributions
Basic Concepts
Consider a random sample 𝑥1, 𝑥2, … , 𝑥 𝑛 from a population. As previously described, the mean of
the sample (called the sample mean) is defined as:

55
𝑥̅ =
1
𝑛
∑ 𝑥𝑖
𝑛
𝑖=1
𝑥̅ can be considered to be a number representing the mean of the actual sample taken but it
can also be considered as a random variable representing the mean of any sample of size n
from the population.
If 𝑥1, 𝑥2, … , 𝑥 𝑛 is a random sample from a population with mean μ and standard deviation σ, then
it turns out that the mean of the random variable 𝑥̅ is μ and the standard deviation of 𝑥̅ is 𝜎 √ 𝑛⁄ .
When the population is normal, we have the following stronger result, namely, that the sample
mean 𝑥̅ is normally distributed with mean μ and standard deviation 𝜎 √ 𝑛⁄ .
Central Limit Theorem
Provided n is large enough (generally about 30 or more), even for samples not coming from a
population with a normal distribution, the sample mean 𝑥̅ is approximately normally distributed
with mean μ and standard deviation 𝜎 √ 𝑛⁄ .
This result is called the Central Limit Theorem and is one of the most important and useful
results in statistics.
Hypothesis Testing
We will now briefly introduce the concept of hypothesis testing using the following example.
Example: A company selling electric motors claims that the average life for its motors is at least
120 months. One of its clients wanted to verify this claim by testing a random sample of 48
motors as described in range A3:F10 of Figure 27. Is the company’s claim correct?
Figure 27: One-sample testing of the mean

56
We first note that the average life of the motors in the sample is 125.9375, which is more than
120 months. But is this a statistically significant result or merely due to chance? Keep in mind
that, even if the company’s claim is false, any given sample may have a mean larger than 120
months.
Now, let’s assume that the population of motors is distributed with mean μ and standard
deviation σ. To perform the analysis, we first define the null hypothesis as follows:
H0: μ < 120
The complement of the null hypothesis (called the alternative hypothesis) is therefore:
H1: μ ≥ 120
Usually in an analysis, we are actually testing the validity of the alternative hypothesis by testing
whether or not to reject the null hypothesis. When performing such analyses, there is some
chance that we will reach the wrong conclusion. There are two types of errors:
 Type I – H0 is rejected even though it is true (false positive)
 Type II – H0 is not rejected even though it is false (false negative)
The acceptable level of a Type I error (called the significance level) is designated by alpha (α)
while the acceptable level of a Type II error is designated beta (β).
Typically, a significance level of α = .05 is used (although, sometimes other levels such as α =
.01 may be employed). This means that we are willing to tolerate up to five percent of type I
errors, that is, we are willing to accept the fact that, in 1 out of every 20 samples, we reject the
null hypothesis even though it is true.
For our analysis, we use the Central Limit Theorem. Since the sample size is sufficiently large,
we can assume that the sample means have a normal distribution, with mean μ and standard
deviation 𝜎 √ 𝑛⁄ , as depicted in Figure 28:
Figure 28: Right-tailed critical region
Keep in mind that this curve represents the distribution of all possible values for the sample
mean. We need to perform our test based on the observed sample mean of 125.9375.
α
critical
value

57
The next step in the analysis is to determine the probability that the observed sample mean has
a value of 125.9375 or higher, under the assumption that the null hypothesis is true. If this
probability is low (i.e., less than α), then we will have grounds for rejecting the null hypothesis.
Otherwise, we would not reject the null hypothesis (which is also called retaining the null
hypothesis).
Retaining the null hypothesis is not the same as positing that it is true. It merely means that
based on the one observed sample, we don’t have sufficient evidence for rejecting it.
The p-value (that is, the probability value) is the value p of the statistic used to test the null
hypothesis. If p < α, then we reject the null hypothesis. For our problem, the p-value is the
probability that the observed sample mean has a value of 125.9375 or higher, under the
assumption that the null hypothesis is true. This is the same as testing whether the p-value falls
in the shaded area of Figure 28 (i.e., the critical region) in what is called a one-tailed test.
The p-value for our example is calculated in cell I8 of Figure 28. Since p-value = .042998 < .05
= α, we see that the observed value does fall in the critical region, and so we reject the null
hypothesis. Therefore, we can show that we have sufficient grounds for accepting the
company’s claim that their motors last at least 120 months on average.
Note that the population standard deviation σ is unknown. For the analysis shown in Figure 28,
we made the assumption that the sample standard deviation is a good estimate of the
population standard deviation. For large enough samples, this is a reasonable assumption,
although it introduces additional error. In Chapter 6, we show how to perform the same type of
analysis using the t distribution without having to make this assumption.
To solve this problem, we could also use Excel’s Z.TEST function (called ZTEST for versions of
Excel prior to Excel 2010). See cell I10 in Figure 27.
Z.TEST(R, 𝜇0, 𝜎) = 1 – NORM.DIST(𝑥̅, 𝜇0, 𝜎/√ 𝑛,TRUE) where 𝑥̅ = AVERAGE(R) = the
sample mean of the data in range R and 𝑛 = COUNT(R) = sample size. The third parameter
is optional. When it is omitted, the value of the sample standard deviation of R is used
instead; that is, Z.TEST(R, 𝜇0) = Z.TEST(R, 𝜇0, STDEV.S(R)).
We make one final observation. If we had tried to test the null hypothesis H0: μ = 120, then the
alternative hypothesis would become H0: μ ≠ 120. In this case, we would have two critical
regions where we would reject the null hypothesis as shown in Figure 29:
Figure 29: Two-tailed hypothesis testing
α/2α/2
critical
value
critical
value

58
This is called a two-tailed test, which is more commonly used than the one-tailed test.

59
Chapter 6 Binomial Distribution
Basic Concepts
Suppose an experiment has the following characteristics:
 The experiment consists of n independent trials, each with two mutually exclusive
outcomes (success and failure)
 For each trial, the probability of success is p (and so the probability of failure is 1 – p)
Each such trial is called a Bernoulli trial. Let 𝑥 be the discrete random variable whose value is
the number of successes in 𝑛 trials. Then the probability distribution function for 𝑥 is called the
binomial distribution which is defined as follows:
𝑓(𝑥) = 𝐶(𝑛, 𝑥)𝑝 𝑥(1 − 𝑝) 𝑛−𝑥
where 𝐶(𝑛, 𝑥) =
𝑛!
𝑥!(𝑛−𝑥)!
and n factorial is 𝑛! = 𝑛(𝑛 − 1)(𝑛 − 2) ⋯ 3 ∙ 2 ∙ 1.
𝐶(𝑛, 𝑥) can be calculated in Excel by using the function COMBIN(n,x).
Figure 30 shows a graph of the probability density function for the binomial distribution with n =
10 and p = .25:
Figure 30: Binomial distribution
The mean of this distribution is np and the variance is np(1 – p ).
0.00
0.05
0.10
0.15
0.20
0.25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Binomial Distribution (n = 20, p = .25)

60
Excel Functions
Excel provides the following functions regarding the binomial distribution:
BINOM.DIST(x, n, p, cum) where n = the number of trials, p = the probability of success for
each trial, and cum takes the value TRUE or FALSE
BINOM.DIST(x, n, p, FALSE) = probability density function value at x for the binomial
distribution, that is, the probability that there are x successes in n trials where the probability
of success on any trial is p
BINOM.DIST(x, n, p, TRUE) = cumulative probability distribution value at x for the binomial
distribution, that is, the probability that there are, at most, x successes in n trials where the
probability of success on any trial is p
BINOM.INV(n, p, 1 – α) = the critical value; that is, the smallest value of x such that
BINOM.DIST(x, n, p, TRUE) ≥ 1 – α
For versions of Excel prior to Excel 2010, the following functions are used instead: BINOMDIST
which is equivalent to BINOM.DIST, and CRITBINOM which is equivalent to BINOM.INV.
Excel 2013 introduces the following new function (where x ≤ y ≤ n):
BINOM.DIST.RANGE(n, p, x, y) = the probability there are between x and y successes
(inclusive) in n trials where the probability of success on any trial is p
Thus,
BINOM.DIST.RANGE(n, p, 0, y) = BINOM.DIST(y, n, p, TRUE)
and for x > 0
BINOM.DIST.RANGE(n, p, x, y) = BINOM.DIST(y, n, p, TRUE)
– BINOM.DIST(x – 1, n, p, TRUE)
The y parameter may be omitted, in which case we have:
BINOM.DIST.RANGE(n, p, x) = BINOM.DIST(x, n, p, FALSE)
Example: What is the probability that, if you throw a die 10 times, it will come up a 6 four times?
We can model this problem using the binomial distribution with n = 10 and p = 1/6 as follows:
𝑓(4) = 𝐶(10, 4) (
1
6
)
4
(1 −
1
6
)
10−4
= 0.054266
Alternatively, the problem can be solved using the Excel formula:
BINOM.DIST(4, 10, 1/6, FALSE) = 0.054266

61
The probability that the die will come up a 6 at least four times can be calculated by the Excel
formula:
1 – BINOM.DIST(3, 10, 1/6, TRUE) = 0.06728
Hypothesis Testing
Example: Suppose you have a die and you suspect that it is biased towards the number 3. So
you run an experiment in which you throw the die 10 times and count that the number 3 comes
up four times. Determine whether the die is biased.
The population random variable 𝑥 = the number of times 3 occurs in 10 trials has a binomial
distribution. Let 𝜋 be the population parameter corresponding to the probability of success on
any trial. We define the following null and alternative hypotheses:
H0: 𝜋 ≤ 1/6; that is, the die is not biased towards the number 3
H1: 𝜋 > 1/6
Now, setting α = .05, we have:
𝑃(𝑥 ≤ 4) = BINOM.DIST(4, 10, 1/6, TRUE) = 0.984538 > 0.95 = 1 – α.
And so we reject the null hypothesis with 95 percent level of confidence.
Example: We suspect that a coin is biased towards heads. When we toss the coin nine times,
how many heads need to come up before we are confident that the coin is biased towards
heads?
We use the following null and alternative hypotheses:
H0: 𝜋 ≤ .5
H1: 𝜋 > .5
Using a confidence level of 95 percent (i.e., α = .05), we calculate:
BINOM.INV(n, p, 1 – α) = BINOM.INV(9, .5, .95) = 7
And so 7 is the critical value. If seven or more heads come up, then we are 95 percent confident
that the coin is biased towards heads and so we can reject the null hypothesis.
Note that BINOM.DIST(6, 9, .5, TRUE) = .9102 < .95 while BINOM.DIST(7, 9, .5, TRUE) =
.9804 ≥ .95.

62
Chapter 7 Student’s t Distribution
Basic Concepts
The one sample hypothesis test using the normal distribution as described in Chapter 5 is fine
when one knows the standard deviation of the population distribution and the population is
either normally distributed or the sample is sufficiently large that the Central Limit Theorem
applies.
The problem with this approach is that the standard deviation of the population is generally not
known. One way of addressing this is to use the standard deviation s of the sample instead of
the standard deviation σ of the population, employing the t distribution.
The (Student’s) t distribution with k degrees of freedom has probability density function
given by the formula:
𝑓(𝑥) =
Γ (
𝑘 + 1
2 )
√𝜋𝑘Γ(
𝑘
2)
(1 +
𝑥2
𝑘
)
−(𝑘+1) 2⁄
where Γ(𝑦) is the gamma function.
The overall shape of the probability density function of the t distribution resembles the bell
shape of a normally distribution with mean 0 and variance 1, except that it is a bit lower and
wider. As the number of degrees of freedom grows, the t distribution approaches the standard
normal distribution and, in fact, the approximation is quite close for k ≥ 30:
Figure 31: Chart of the t distribution by degrees of freedom
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
0.45
-4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0
Student's t Distribution
t(df = 1)
t(df = 5)
t(df = 10)
t(df = 30)
std nomal

63
If 𝑥 has a normal distribution with mean μ and standard deviation σ, then for samples of size n
the random variable:
𝑡 =
𝑥̅ − 𝜇
𝑠 √ 𝑛⁄
has a t distribution with n – 1 degrees of freedom. The same is true even when 𝑥 doesn’t have a
normal distribution, provided that n is sufficiently large (by the Central Limit Theorem).
This test statistic is the same as 𝑧 =
𝑥̅− 𝜇
𝜎 √ 𝑛⁄
from the Central Limit Theorem with the population
standard deviation σ replaced by the sample standard deviation s. What makes this useful is
that usually the standard deviation of the population is unknown while the standard deviation of
the sample is known.
Excel Functions
Excel provides the following functions regarding the t distribution:
T.DIST(x, df, cum) = the probability density function value at x for the t distribution with df
degrees of freedom when cum = FALSE and the corresponding cumulative distribution
function at x when cum = TRUE
T.INV(p, df) = the value x such that T.DIST(x, df, TRUE) = p, that is, inverse of T.DIST(x, df,
TRUE)
T.DIST.RT(x, df) = the right tail at x of the t distribution with df degrees of freedom
T.DIST.2T(x, df) = the sum of the right tail of the t distribution with df degrees of freedom at
x plus the left tail at -x, where x ≥ 0 (the function yields an error value when x < 0
T.INV.2T(p, df) = the value xsuch that T.DIST.2T(x, df) = p, that is, inverse of T.DIST.2T(x,
df)
These functions were not available prior to Excel 2010. The following table lists how these
functions are used as well as their equivalents using the functions TDIST and TINV, which were
used prior to Excel 2010:
Tail Distribution Inverse
Right tail T.DIST.RT(x, df)
= TDIST(x, df, 1) if x ≥ 0
= 1-TDIST(-x, df, 1) if x < 0
T.INV(1-α, df) = TINV(2*α, df)

64
Left tail T.DIST(x, df, TRUE)
= 1-TDIST(x, df, 1) if x ≥ 0
= TDIST(-x, df, 1) if x < 0
T.INV(α, df) = -TINV(2*α, df)
Two-tail T.DIST.2T(x, df) = TDIST(x, df, 2) T.INV.2T(α, df) = TINV(α, df)
Table 8: Excel t distribution functions
Hypothesis Testing
One-Sample Testing
As explained previously, the t distribution provides a good way to perform one sample test on
the mean when the population variance is not known, provided the population is normal or the
sample is sufficiently large so that the Central Limit Theorem applies.
It turns out that the t distribution provides good results even when the population is not normal,
and even when the sample is small, provided the sample data is reasonably symmetrically
distributed about the sample mean. This can be determined by graphing the data. The following
are indications of symmetry:
 The box plot is relatively symmetrical; that is, the median is in the center of the box and
the whiskers extend equally in each direction
 The histogram looks symmetrical
 The mean is approximately equal to the median
 The coefficient of skewness is relatively small
Example: A weight reduction program claims to be effective in treating obesity. To test this
claim, 12 people were put on the program and the number of pounds of weight gained/lossed
was recorded for each person after two years. This is shown in columns A and B of Figure 32.
Can we conclude that the program is effective?

65
Figure 32: One sample t test
A negative value in column B indicates that the subject gained weight. We judge the program to
be effective if there is some weight loss at the 95 percent significance level. We choose to
conduct a two-tailed test since there is risk that the program might actual result in weight gain
rather than loss. Thus, our null hypothesis is:
H0: μ = 0; that is, the program is not effective.
From the box plot in Figure 33, we see that the data is quite symmetric and so we use the t test
even though the sample is small:
Figure 33: Box plot of the sample data
-15
-10
-5
0
5
10
15
20
25
30
Weight Loss

66
Column E of Figure 33 contains all the formulas required to carry out the t test. Since Excel only
displays the values of these formulas, we show each of the formulas (in text format) in column G
so that you can see how the calculations are performed.
Thus, where R is B4:B15, we see that 𝑛 = COUNT(R) = 12, 𝑥̅ = AVERAGE(R) = 5.5, s =
STDEV(R) = 11.42, and standard error = 𝑠 √ 𝑛⁄ = 3.30. From this we see that:
𝑡 =
𝑥̅− 𝜇
𝑠 √ 𝑛⁄
=
5.5−0
3.3
= 1.668 with 𝑑𝑓 = 11 degrees of freedom.
Since p-value = T.DIST.2T(t, df) = T.DIST.2T(1.668, 11) = .123 > .05 = α, the null hypothesis is
not rejected. This means there is an 12.3 percent probability of achieving a value for t this high,
assuming that the null hypothesis is true and, since 12.3 percent > 5 percent, we can’t reject the
null hypothesis.
The same conclusion is reached since:
t-crit = T.INV.2T(α, df) = T.INV.2T(.05, 11) = 2.20 > 1.668 = t-obs.
Here we used T.DIST.2T and T.INV.2T since we are employing a two-tailed test.
Testing Two Independent Samples
Assuming equal population variances
Our goal is to determine whether or not the means of two populations are equal given two
independent samples, one from each population. We start by assuming that population
variances are equal even if unknown.
Suppose that 𝑥̅ and 𝑦̅ are the sample means of two independent samples of size 𝑛 𝑥 and 𝑛 𝑦
respectively. Also suppose that 𝑠 𝑥 and 𝑠 𝑦 are the sample standard deviations of two sets of
data. If 𝑥 and 𝑦 are normal or if 𝑛 𝑥 and 𝑛 𝑦 are sufficiently large for the Central Limit Theorem to
hold, then the statistic:
𝑡 =
(𝑥̅ − 𝑦̅) − (𝜇 𝑥 − 𝜇 𝑦)
𝑠√
1
𝑛 𝑥
+
1
𝑛 𝑦
has a t distribution with 𝑛 𝑥 + 𝑛 𝑦 − 2 degrees of freedom where 𝑠2
is the pooled variance defined
by:
𝑠2
=
(𝑛 𝑥 − 1)𝑠 𝑥
2
+ (𝑛 𝑦 − 1)𝑠 𝑦
2
(𝑛 𝑥 − 1) + (𝑛 𝑦 − 1)

67
We can use the statistic t to test the null hypothesis that the two population means are equal or
equivalently that 𝜇 𝑥 − 𝜇 𝑦 = 0. It turns out that, even when 𝑥 and 𝑦 are not normal and 𝑛 𝑥 and 𝑛 𝑦
are not particularly large, as long as the two samples are reasonably symmetric, then this test
yields pretty good results.
Example: A food company wants to determine whether or not their new formula for peanut
butter is significantly tastier than their old formula. They chose two random samples of 10
people and ask the people in the first sample to taste the peanut butter using the existing
formula, and they ask the people in the second sample to taste the peanut butter with the new
formula. They then ask all 20 people to fill out a questionnaire rating the tastiness of the peanut
butter they tasted.
Based on the data in Figure 34, determine whether there is a significant difference between the
two types of peanut butter.
As we can see from Figure 34, the mean score for the new formula is 15 while the mean score
for the old formula is 11. But is this difference just due to random effects or is there a statistically
significant difference?
We also note that the variances for the two samples are 13.3 and 18.8. It turns out that these
two values are close enough to satisfy the equal variance assumption.
Figure 34: Sample data
Finally, if we draw the box plots for the two data sets, we see that both are relatively symmetric
(that is, for each sample the colored areas in each box are approximately equal in size, and the
upper and lower whiskers are fairly equal in length).

68
Thus, the assumptions for the t test are met. To carry out the test, we use Excel’s t-Test: Two
Sample Assuming Equal Variances data analysis tool which, as usual, can be accessed by
selecting Data > Analysis|Data Analysis. The result is shown in Figure 35:
t-Test: Two-Sample Assuming Equal Variances
New Old
Mean 15 11.1
Variance 13.33333 18.76667
Observations 10 10
Pooled Variance 16.05
Hypothesized Mean Difference 0
Df 18
t Stat 2.176768
P(T<=t) one-tail 0.021526
t Critical one-tail 1.734064
P(T<=t) two-tail 0.043053
t Critical two-tail 2.100922
Figure 35: Two-sample t test assuming equal variance
Assuming that we are performing a two-tailed test, we see that the p-value = 0.043 which is less
than the usual alpha value of .05. Thus, we can conclude that there is a significant difference
between the mean scores for the two peanut butter formulations.
Assuming unequal population variances
In the previous section, we assumed that the two population variances were equal (or at least
not very unequal). When the assumption of equal population variances is not met (or when we
don’t have enough evidence to draw this conclusion), we can modify the approach used in the
previous section by employing a slightly different version of the t test.
Suppose that 𝑥̅ and 𝑦̅ are the sample means of two independent samples of size 𝑛 𝑥 and 𝑛 𝑦
respectively. Also suppose that 𝑠 𝑥 and 𝑠 𝑦 are the sample standard deviations of the two sets of
data. If 𝑥 and 𝑦 are normally distributed or if 𝑛 𝑥 and 𝑛 𝑦 are sufficiently large for the Central Limit
Theorem to hold, then the statistic:
𝑡 =
(𝑥̅ − 𝑦̅) − (𝜇 𝑥 − 𝜇 𝑦)
√
𝑠 𝑥
2
𝑛 𝑥
+
𝑠 𝑦
2
𝑛 𝑦
has a t distribution with the following degrees of freedom.

69
𝑑𝑓 =
(
𝑠 𝑥
2
𝑛 𝑥
+
𝑠 𝑦
2
𝑛 𝑦
)
2
(
𝑠 𝑥
2
𝑛 𝑥
)
2
𝑛 𝑥 − 1 +
(
𝑠 𝑦
2
𝑛 𝑦
)
2
𝑛 𝑦 − 1
Example: Repeat the previous analysis using the data in Figure 36:
Figure 36: Revised sample data
This time, we see that there is a big difference between the variances (118.3 vs. 18.8).
Generally, even if one variance is up to four times the other, the equal variance assumption will
give good results. This rule of thumb is clearly violated in this example, and so we need to use
the t test with unequal population variances, namely, Excel’s t-Test: Two Sample Assuming
Unequal Variances data analysis tool.
The output is shown in Figure 37:
t-Test: Two-Sample Assuming Unequal Variances
New Old
Mean 18.9 11.1
Variance 118.3222 18.76667
Observations 10 10
df 12
t Stat 2.106655

70
Figure 37: Two-sample t test assuming unequal variance
Note that the degrees of freedom have been reduced from 18 in Example 1 to 12 in this
example (actually 11.78 although Excel rounds this up to 12). Since for the two-tail test, p-value
= 0.056869 > .05 = α, we cannot reject the null hypothesis that there is no significant difference
between the two formulations.
If we had used the equal variance test, we would have gotten the results shown in Figure 6.
This would have led us to a different conclusion since, for the two-tail test, p-value = 0.049441 <
.05 = α, and therefore we would have concluded that there was a significant difference between
the two formulations.
t-Test: Two-Sample Assuming Equal Variances
New Old
Mean 18.9 11.1
Variance 118.3222 18.76667
Observations 10 10
Pooled Variance 68.54444
df 18
t Stat 2.106655
Figure 38: Same problem assuming equal variance
Testing Paired Samples
In paired sample hypothesis testing, a sample from the population is chosen and two
measurements for each element in the sample are taken. Each set of measurements is
considered a sample. Unlike the hypothesis testing studied so far, the two samples are not
independent of one another.
For example, if you want to determine whether drinking a glass of wine or drinking a glass of
beer has the same or different impact on memory, one approach is to take a sample of, say, 40
people. You then have half of them drink a glass of wine and the other half drink a glass of beer.
You then give each of the 40 people a memory test and compare the results. This is the
approach with independent samples.

71
Another approach is to take a sample of 20 people and have each person drink a glass of wine
and take a memory test. You then have the same people drink a glass of beer and again take a
memory test. Finally, you compare the results. This is the approach used with paired samples.
The advantage of this second approach is that the sample can be smaller. Also, since the
sampled subjects are the same for beer and wine, there is less of a chance that some external
factor (or confounding variable) will influence the result. The problem with this approach is that
it is possible that the results of the second memory test will be lower simply because the person
has imbibed more alcohol. This can be corrected by sufficiently separating the tests, for
example, by conducting the test with beer a day after the test with wine.
It is also possible that the order in which people take the tests influences the result (for example,
the subjects learn something on the first test that helps them on the second test or perhaps
taking the test the second time introduces a degree of boredom that lowers the score). One way
to address these order effects is to have half the people drink wine on Day 1 and beer on Day
2, while for the other half, the order is reversed (this is called counterbalancing).
Obviously, not all experiments can use the paired sample design. For example, if you are
testing differences between men and women, then independent samples will be necessary.
Example: A clinic provides a program to help their clients lose weight and asks a consumer
agency to investigate the effectiveness of the program. The agency takes a sample of 15
people, weighs each person in the sample before the program begins, and then weighs each
person three months later to produce the results in the table in Figure 39. Determine whether
the program is effective.

72
The program is effective if there is a significant reduction in the weight of the people in the
program. Let 𝑥 = the reduction in weight three months after the program starts. The null
hypothesis is:
H0: 𝜇 = 0; that is, any differences in weight is due to chance (two-tailed test).
Essentially, the test to use is the one sample test as previously described, but on the difference
in the weight before the program minus the weight after the program. Since all the values in
column D of Figure 39 are positive, we see that all the subjects in the sample lost weight. But is
the reduction significant?
We could create the analysis (as in Figure 32) using column D of Figure 39 as the single
sample. Instead, we will use Excel’s t-Test: Paired Two Sample for Means data analysis tool.
We begin by selecting Data > Analysis|Data Analysis.
When the dialog box shown in Figure 40 appears, you need to fill in the values as shown in the
figure and click OK:
Figure 40: T-Test: Paired two-sample for means dialog box
The output of the analysis is displayed in Figure 41:
t-Test: Paired Two Sample for Means
Before After
Mean 207.9333 197
Variance 815.781 595
Observations 15 15
Pearson Correlation 0.98372

73
df 14
t Stat 6.6897
P(T<=t) one-tail 5.14E-06
P(T<=t) two-tail 1.03E-05
Figure 41: Paired sample t-test analysis
Since the p-value for the two-tailed test is 1.03E-05, we can conclude that there is a significant
difference between the mean weights. Since the mean afterward is lower than the mean weight
beforehand, we can conclude that there is a significant reduction in weight. Thus, it is fair to say
that the program is effective (at least for the three-month period).
Using the T.TEST Function
In addition to the approaches shown in the last few sections, Excel provides the T.TEST
function to perform the three different types of two sample t tests.
T.TEST(R1, R2, tails, type) = the t-test value for the difference between the means of two
samples R1 and R2 where tails = 1 (one-tailed) or 2 (two-tailed) and type takes the values:
1. The samples have paired values from the same population
2. The samples are from populations with the same variance
3. The samples are from populations with different variances
These three types correspond to the Excel data analysis tools:
1. t-Test: Paired Two Sample for Mean
2. t-Test: Two-Sample Assuming Equal Variance
3. t-Test: Two-Sample Assuming Unequal Variance
The T.TEST function is not available for versions of Excel prior to Excel 2010. For earlier
versions of Excel, the TTEST function is available, which provides identical functionality.
For example, we can repeat the analysis from Figure 35 using the T.TEST function. As we can
see from Figure 42, the p-value for the two-tailed test for two independent samples from
populations with the same variances is 0.043053 (cell E40). This is the same as the result
previously found:
Figure 42: T.TEST for two independent samples

74
Figure 42 also shows the value of the test if we didn’t assume equal variances (cell 39). As you
can see, there is very little difference between the two results.
Confidence Intervals
Until now, we have made what are called point estimates of the mean or differences between
means. We now define the confidence interval which provides a range of possible values
instead of just a single measurement. We can illustrate how this works using the t distribution:
The 1 – α confidence interval for the population mean is:
𝑜𝑏𝑠𝑒𝑟𝑣𝑒𝑑 𝑚𝑒𝑎𝑛 ± 𝑡 𝑐𝑟𝑖𝑡 ∙ 𝑠𝑡𝑑 𝑒𝑟𝑟𝑜𝑟
Example: Calculate the 95 percent confidence interval for the one sample t test example (see
Figure 32):
𝑚𝑒𝑎𝑛 𝑜𝑏𝑠 ± 𝑡 𝑐𝑟𝑖𝑡 ∙ 𝑠𝑡𝑑 𝑒𝑟𝑟 = 5.5 ± 2.201 ∙ 3.297 = 5.5 ± 7.257
This yields a 95 percent confidence interval of (-1.757, 12.757) for the population mean. Since
the interval contains the hypothetical mean of 0, we are justified once again in our decision not
to reject the null hypothesis.
Versions of Excel starting with Excel 2010 provide the following function to calculate the
confidence interval for the t distribution:
CONFIDENCE.T(α, s, n) = k such that ( 𝑥̅ – k, 𝑥̅ + k) is the 1 – α confidence interval for the
population mean; that is, CONFIDENCE.T(α, s, n) = 𝑡 𝑐𝑟𝑖𝑡 ∙ std error, where n = sample size
and s = sample standard deviation.
Thus, for the last example, we calculate CONFIDENCE.T(.05,11.422,12) = 7.257 which yields a
95 percent confidence interval of 5.5 ± 7.257, as before.

75
Chapter 8 Chi-square and F Distribution
Chi-square Distribution
Basic Concepts
The chi-square distribution with k degrees of freedom has the probability density function:
𝑓(𝑥) =
𝑥
𝑘
2
−1
𝑒−
𝑥
2
2 𝑘 2⁄ Γ (
𝑘
2
)
where Γ is the gamma function and k does not have to be an integer and can be any positive
real number. Note that chi-square is commonly written as 𝜒2
.
The following are the graphs of the probability density function with degrees of freedom df= 5
and 10. As df grows larger, the fat part of the curve shifts to the right and becomes more like the
graph of a normal distribution:
Figure 43: Chart of the chi-square distribution
The key relationship between the chi-square and normal distribution is that, if 𝑥 has a standard
normal distribution, then 𝑥2
has a chi-square distribution with one degree of freedom.
Excel Functions
Excel provides the following functions regarding the chi-square distribution:
0.00
0.02
0.04
0.06
0.08
0.10
0.12
0.14
0.16
0.18
0 2 4 6 8 10 12 14 16 18 20 22 24
Chi-square
df = 5
df =10

76
CHISQ.DIST(x, df, cum) = the probability density function value at x for the chi-square
distribution with df degrees of freedom when cum = FALSE and the corresponding
cumulative distribution function at x when cum = TRUE
CHISQ.INV(p, df) = the value x such that CHISQ.DIST(x, df, TRUE) = p, that is, inverse of
CHISQ.DIST(x, df, TRUE)
CHISQ.DIST.RT(x, df) = the right-tail at x of the chi-square distribution with df degrees of
freedom
CHISQ.INV.RT(p, df) = the value 𝑥 such that CHISQ.DIST.RT(x, df) = p, that is, inverse of
CHISQ.DIST.RT(x, df)
functions are used, as well as their equivalents using the functions CHIDIST and CHIINV which
were used prior to Excel 2010:
Right-tail
CHISQ.DIST.RT(x, df) =
CHIDIST(x, df)
CHISQ.INV.RT(α, df) =
CHIINV(α, df)
Left-tail
CHISQ.DIST(x, df, TRUE)
=1-CHIDIST(x, df)
CHISQ.INV(α, df) =
CHIINV(1-α, df)
Table 9: Formulas for the chi-square distribution in Excel
Contingency Tables
The chi-square distribution is commonly used to determine whether or not two sets of data are
independent of each other. Such data are organized in what are called contingency tables as
described in the following example.
Example: A survey is conducted of 175 young adults whose parents are classified either as
wealthy, middle class, or poor to determine their highest level of schooling (e.g., graduated from
university, graduated from high school, or neither). The results are summarized in the
contingency table in Figure 44 (Observed Values). Based on the data collected, is a person’s
level of schooling independent of their parents’ wealth?

77
Figure 44: Contingency Table (observed values)
We set the null hypothesis to be:
H0: Highest level of schooling attained is independent of parents’ wealth.
To use the chi-square test, we need to calculate the expected values that correspond to the
observed values in the table in Figure 44. To accomplish this, we use the fact that if two events
are independent, then the probability that both events occur is equal to the product of the
probabilities that each will occur.
We also assume that the proportions for the sample are good estimates of the probabilities of
the expected values.
We now show how to construct the table of expected values (see Figure 45). We know from
Figure 44 that 45 of the 175 people in the sample are from wealthy families. So the probability
that someone in the sample is from a wealthy family is 45/175 = 25.7 percent. Similarly, the
probability that someone in the sample graduated from university is 68/175 = 38.9 percent. But
based on the null hypothesis, being from a wealthy family is independent of graduating from
university and so, the expected probability of both is simply the product of the two events or
25.7 percent ∙ 38.9 percent = 10.0 percent. Thus, based on the null hypothesis, we expect that
10.0 percent of 175 = 17.5 people are from a wealthy family and have graduated from
university:
Figure 45: Expected values
Thus, for example, cell H6 can be calculated by the formula =H$9*$K6/$K$9. If you highlight the
range H6:J8 and press Ctrl-D and then Ctrl-R you will fill in the rest of the table.

78
Alternatively, you can fill in the table of expected values by highlighting the range H6:J8 and
entering the array formula:
=MMULT(K6:K8,H9:J9)/K9
Independence Testing
The Pearson’s chi-square test statistic is defined as:
∑
(𝑜𝑏𝑠𝑖 − 𝑒𝑥𝑝𝑖)2
𝑒𝑥𝑝𝑖
𝑘
𝑖=1
For sufficiently large samples, the Pearson’s chi-square test statistic has approximately a chi-
square distribution with (row count – 1) (column count – 1) degrees of freedom.
The test requires that the following conditions be met:
 Random sample: Data must come from a random sampling of the population
 Independence: The observations must be independent of each other. This means chi-
square cannot be used to test correlated data (e.g., matched pairs)
 Cell size: All of the expected frequency values in the contingency table are at least 5,
although with large tables, occasional values lower than 5 will probably still yield good
results. With a 2×2 contingency table, it is better to have even larger frequency values.
In Figure 46, we show how to use the chi-square test to carry out the analysis:
Figure 46: Chi-square test
Since p-value = 0.003273 < .05 = α, we reject the null hypothesis and conclude that a person’s
level of schooling is not independent of their parents’ wealth.
Excel provides the following function which carries out the chi-square test for independence for
contingency tables:

79
CHISQ.TEST(R1, R2) = CHISQ.DIST.RT(x, df) where R1 = the array of observed data, R2
= the array of expected value, x is the chi-square statistic, and df = (row count – 1) (column
count – 1).
For our example, CHISQ.TEST(B6:D8,H6:J8) = 0.003273.
The CHISQ.TEST function is only available in versions of Excel starting with Excel 2010. For
previous versions of Excel, the equivalent function CHITEST can be used.
F Distribution
Basic Concepts
If 𝑥1 and 𝑥2 have chi-square distributions with 𝑑𝑓1 and 𝑑𝑓2 degrees of freedom respectively, then
the quotient F = 𝑥1 /𝑥2 has F distribution with 𝑑𝑓1, 𝑑𝑓2 degrees of freedom.
The F distribution is commonly used to compare two variances. It is useful in ANOVA,
regression, and other commonly used tests. In particular, if we draw two independent samples
of size 𝑛1 and 𝑛2 respectively from two different normal populations with the same variance,
then the quotient of the two sample variances has F distribution with 𝑛1 − 1, 𝑛2 − 1 degrees of
freedom.
Excel Functions
Excel provides the following functions regarding the F distribution:
F.DIST(x, df, cum) = the probability density function value at x for the F distribution with df
degrees of freedom when cum = FALSE and the corresponding cumulative distribution
function at x when cum = TRUE
F.INV(p, df) = the value x such that F.DIST(x, df, TRUE) = p, that is, inverse of F.DIST(x, df,
TRUE)
F.DIST.RT(x, df) = the right-tail at x of the F distribution with df degrees of freedom
F.INV.RT(p, df) = the value 𝑥 such that F.DIST.RT(x, df) = p, that is, inverse of
F.DIST.RT(x, df)
functions are used as well as their equivalents using the functions FDIST and FINV that were
available prior to Excel 2010:

80
Right-tail
F.DIST.RT(x, df1, df2) =
FDIST(x, df1, df2)
F.INV.RT(α, df1, df2) =
FINV(α, df1, df2)
Left-tail
F.DIST(x, df1, df2, TRUE) =
1 – FDIST(x, df1, df2)
F.INV(α, df1, df2) =
FINV(1-α, df1, df2)
Table 10: F distribution formulas in Excel
Hypothesis Testing
As we previously observed, if we draw two independent samples of size 𝑛1 and 𝑛2 respectively
from two different normal populations with the same variance, then the quotient of the two
sample variances has F distribution with 𝑛1 − 1, 𝑛2 − 1 degrees of freedom.
This fact can be used to test whether or not the variances of two populations are equal. In order
to deal exclusively with the right-tail of the distribution, when taking ratios of sample variances
we should put the larger variance in the numerator of the quotient.
In order to use this test, the following must hold:
 Both populations are normally distributed
 Both samples are drawn independently from each other.
 Within each sample, the observations are sampled randomly and independently of each
other
Example
A company is comparing two methods for producing high-precision bolts and wants to choose
the method with the least variability. It has taken a sample of the lengths of the bolts using both
methods as given in the table on the left side of Figure 47:

81
Figure 47: F test to compare two variances
To carry out the test, we use Excel’s F-Test Two Sample For Variances data analysis tool
which, as usual, can be accessed by selecting Data > Analysis|Data Analysis. The result is
shown on the right side of Figure 47.
Since p-value = .271364 > .05 = alpha, we cannot reject the null hypothesis. So we must
conclude that there is no significant difference between the two variances.
Excel Test Function
Excel provides the following test statistic function:
F.TEST(R1, R2) = two-tailed F-test comparing the variances of the samples in ranges R1
and R2 = the two-tailed probability that the variance of the data in ranges R1 and R2 are
not significantly different.
Note that, unlike the F-Test data analysis tool or F.DIST, F.DIST.RT, F.INV and F.INV.RT
functions, which are one-tailed, the F.TEST function is two-tailed. For the previous example:
F.TEST(A4:A18,B4:B18) = .5427
which is double the value obtained for the one-tailed p-value in cell E11 of Figure 47.

82
Chapter 9 Analysis of Variance
Introduction
Essentially, analysis of variance (ANOVA) is an extension of two-sample hypothesis testing for
comparing means (when variances are unknown) to more than two samples.
We start with the one-factor case. We will define the concept of factor later but, for now, we
simply view this type of analysis as an extension of the two independent sample t test with equal
population variances.
One-way ANOVA
One-way ANOVA Example
A food company wants to determine whether or not any of their three new formulas for peanut
butter is significantly tastier than their old formula. They create four random samples of 10
people each (one for each type of peanut butter) and ask the people in each sample to taste the
peanut butter for that sample. They then ask all 40 people to fill out a questionnaire rating the
tastiness of the peanut butter they tried.
Based on the data in Figure 48, determine whether there is a significant difference between the
four types of peanut butter:
The null hypothesis for this example is that any difference between the four types of peanut
butter is due to chance, that is:

83
H0: 𝜇1 = 𝜇2 = 𝜇3 = 𝜇4
Basic Concepts
Before we proceed with the analysis for this example, we should review a few basic concepts.
Suppose we have 𝑘 samples, which we will call groups (or treatments). These are the
columns in our analysis (corresponding to the four types of peanut butter in the example). We
will use the index 𝑗 for these. Each group consists of a sample of size 𝑛𝑗. The sample elements
are the rows in the analysis. We will use the index 𝑖 for these.
Suppose the 𝑗th group sample is {𝑥1𝑗, … , 𝑥 𝑛 𝑗 𝑗} and so the total sample consists of all the
elements {𝑥𝑖𝑗: 1 ≤ 𝑖 ≤ 𝑛𝑗, 1 ≤ 𝑗 ≤ 𝑘}. We will use the abbreviation 𝑥̅𝑗 for the mean of the 𝑗th
group sample (called the group mean) and 𝑥̅ for the mean of the total sample (called the total
or grand mean).
Let the sum of squares for the 𝑗th group be 𝑆𝑆𝑗 = ∑ (𝑥𝑖𝑗 − 𝑥̅𝑗)
2
𝑖 . We now define the following
terms:
𝑆𝑆 𝑇 = ∑ ∑(𝑥𝑖𝑗 − 𝑥̅)
2
𝑖𝑗
𝑆𝑆 𝐵 = ∑ 𝑛𝑗(𝑥̅𝑗 − 𝑥̅)
2
𝑗
𝑆𝑆 𝑊 = ∑ 𝑆𝑆𝑗
𝑗
= ∑ ∑(𝑥𝑖𝑗 − 𝑥̅𝑗)
2
𝑖𝑗
𝑆𝑆 𝑇 is the sum of squares for the total sample, that is, the sum of the squared deviations from
the grand mean. 𝑆𝑆 𝑊 is the sum of squares within the groups, that is, the sum of the squared
means across all groups. 𝑆𝑆 𝐵 is the sum of the squares between group sample means, that is,
the weighted sum of the squared deviations of the group means from the grand mean.
We also define the following degrees of freedom, where 𝑛 = ∑ 𝑛𝑗
𝑘
𝑗=1 :
𝑑𝑓𝑇 = 𝑛 − 1 𝑑𝑓𝐵 = 𝑘 − 1 𝑑𝑓 𝑊 = ∑ (𝑛𝑗 − 1)
𝑘
𝑗=1
= 𝑛 − 𝑘
Finally, we define the mean square, as 𝑀𝑆 = 𝑆𝑆/𝑑𝑓, and so:
𝑀𝑆 𝑇 = 𝑆𝑆 𝑇 / 𝑑𝑓𝑇 𝑀𝑆 𝐵 = 𝑆𝑆 𝐵 / 𝑑𝑓𝐵 𝑀𝑆 𝑊 = 𝑆𝑆 𝑊 / 𝑑𝑓 𝑊
We summarize these terms in the following table:

84
𝒅𝒇 𝑺𝑺 𝑴𝑺
𝑇 𝑛 – 1 ∑ ∑(𝑥𝑖𝑗 − 𝑥̅)
2
𝑖𝑗
𝑆𝑆 𝑇 𝑑𝑓𝑇⁄
𝐵 𝑘 – 1 ∑ 𝑛𝑗(𝑥̅𝑗 − 𝑥̅)
2
𝑗
𝑆𝑆 𝐵 𝑑𝑓𝐵⁄
𝑊 𝑛 – 𝑘 ∑ ∑(𝑥𝑖𝑗 − 𝑥̅𝑗)
2
𝑖𝑗
𝑆𝑆 𝑊 𝑑𝑓 𝑊⁄
Table 11: Summary of ANOVA terms
Clearly 𝑀𝑆 𝑇 is the variance for the total sample, 𝑀𝑆 𝑊 is the sum of the group sample
variances, and 𝑀𝑆 𝐵 is the variance for the “between sample”, that is, the variance of
{𝑛1 𝑥̅1, … , 𝑛 𝑘 𝑥̅ 𝑘}.
It is also not hard to show that:
𝑆𝑆 𝑇 = 𝑆𝑆 𝑊 + 𝑆𝑆 𝐵
𝑑𝑓𝑇 = 𝑑𝑓 𝑊 + 𝑑𝑓𝐵
It turns out that, if the null hypothesis is true, then 𝑀𝑆 𝑊 and 𝑀𝑆 𝐵 are both measures of the same
error. Thus, the null hypothesis becomes equivalent to the hypothesis that the population
versions of these statistics are equal, that is:
𝜎 𝐵 = 𝜎 𝑊
We can, therefore, use the F-test described in Chapter 8 to determine whether or not to reject
the null hypothesis. This means that if the 𝑥𝑖𝑗 are independently and normally distributed, and all
the 𝜇 𝑗 are equal (null hypothesis), and all the 𝜎𝑗
2
are equal (homogeneity of variances), then the
test statistic:
𝐹 =
𝑀𝑆 𝐵
𝑀𝑆 𝑊
has an F distribution with 𝑑𝑓𝐵, 𝑑𝑓 𝑊 degrees of freedom.

85
Analysis
To carry out the analysis for the example, we will use the ANOVA: Single Factor data analysis
tool. To access this tool, as usual, press Data > Analysis|Data Analysis and fill in the dialog
box that appears as in Figure 49:
Figure 49: Dialog box for ANOVA: Single-factor data analysis tool
The output is as shown in Figure 50:
Figure 50: ANOVA: Single-factor data analysis

86
All the fields in Figure 50 are calculated as previously described. We can see that the test
statistics is F = 3.206928. Since p-value = F.DIST.RT(3.206928, 3, 36) = 0.034463 < .05 = α,
we reject the null hypothesis that there is no significant difference between the evaluations of
the four different types of peanut butter.
Follow-up Analysis
Although we now know that there is a significant difference between the four types of peanut
butter, we still don’t know where the differences lay.
It appears from Figure 50 that New 2 has a higher rating than the other types of peanut butter.
We can use a t test to determine whether or not there is a significant difference between the
rating for New 2 and New 1, the next highest-rated peanut butter. The result is shown in Figure
51:
New 1 New 2
Mean 13.1 16.6
Variance 21.21111 7.155556
Observations 10 10
Df 14
t Stat -2.07809
Figure 51: Follow-up analysis: New 1 vs. New 2
From Figure 51, we see that there is no significant difference between the two types of peanut
butter based on a two-tailed test (0.056577 > .05). If we instead compare the New 2 with the old
formula, we get the results shown in Figure 52:
Old New 2
Mean 11.1 16.6
Variance 18.76667 7.155556
Observations 10 10

87
Df 15
t Stat -3.41607
Figure 52: Follow-up analysis: Old vs. New 2
This time, we see that there is a significant difference between the two types of peanut butter
using a two-tailed t-test (.003829 < .05).
The problem with this approach is that doing multiple tests incurs higher amounts of what is
called experimentwise error. Remember that when we use a significance level of α = .05, we
accept that five percent of the time we will get a type I error. If we perform three such tests, then
we will essentially increase our overall type I error to 1 – (1 – .05)3
= .14. This means that 14
percent of the time, we will have a type I error which is higher than we would like.
The general approach for addressing this issue is to either reduce α using Bonferroni’s
correction (for example, for three tests we use α/3 = .05/3 = .0167) or use a different type of
test (for example, Tukey’s HSD or REGWQ) which is beyond the scope of this book.
Levene’s Test
As previously mentioned, the ANOVA test requires that the group variances be equal. There is a
lot of leeway here and, even when the variance of one group is four times another, the analysis
will be pretty good. Another option for testing homogeneity of group variances is by using
Levene’s test.
For Levene’s test, the residuals 𝑒𝑖𝑗 of the group means from the cell means are calculated as
follows:
𝑒𝑖𝑗 = 𝑥𝑖𝑗 − 𝑥̅𝑗
An ANOVA is then conducted on the absolute value of the residuals. If the group variances are
equal, then the average size of the residual should be the same across all groups.
There are three versions of the test: using the mean (as described), using the median, or using
the trimmed mean.
Example: Use Levene’s test (residuals from the median) to determine whether the four groups
in the ANOVA example have significantly different population variances.
We begin by calculating the medians for each group (range B14:E14 of Figure 53). For
example, cell B14 contains the formula =MEDIAN(B4:B13).

88
Next, we create the table (in range B19:E29) with the absolute residuals from the median. We
do this by entering the formula =ABS(B4-B$14) in cell B20, highlighting the range B20:E29 and
then pressing Ctrl-R followed by Ctrl-D.
Finally, we use the ANOVA: Single Factor data analysis tool as described earlier in this
chapter on the Input Range B19:E29. The output is shown on the right side of Figure 53:
Figure 53: Levene’s test
Since p-value = .25752 > .05 = α, we cannot reject the null hypothesis. So we can conclude that
there is no significant difference between the four group variances. Thus, the ANOVA test
previously conducted satisfies the homogenity of variances assumption.
Factorial ANOVA
Example
A new fertilizer has been developed to increase the yield on crops. The makers of the fertilizer
want to better understand which of the three formulations (or blends) of this fertilizer are most
effective for wheat, corn, soy beans, and rice (the crops). They test each of the three blends on
five samples of each of the four types of crops. The crop yields for the five samples of each of
the 12 combinations are as shown in Figure 54.

89
Determine whether or not there is a difference between crop yields based on these three types
of fertilizer.
Figure 54: Sample data for Two-Factor ANOVA
Before we proceed with the analysis for this example, we should review a few basic concepts.
Basic Concepts
We now extend the one-factor ANOVA model previously described to a model with more than
one factor.
A factor is an independent variable. A 𝑘 factor ANOVA addresses 𝑘 factors. The model for the
last example contains two factors: crops and blends.
A level is some aspect of a factor; these are what we called groups or treatments in the one
factor analysis previously described. The blend factor contains three levels and the crop factor
contains four levels.
Let’s look at the Two-Factor model in more detail.
Suppose we have two factors, A and B, where factor A has 𝑟 levels and factor B has 𝑐 levels.
We organize the levels for factor A as rows and the levels for factor B as columns. We use the
index 𝑖 for the rows (i.e., factor A) and the index 𝑗 for the columns (i.e., factor B).
This results in an 𝑟 × 𝑐 table whose entries are {𝑋𝑖𝑗: 1 ≤ 𝑖 ≤ 𝑟, 1 ≤ 𝑗 ≤ 𝑐} where 𝑋𝑖𝑗 is a sample
for level 𝑖 of factor A and level 𝑗 of factor B. Here, 𝑋𝑖𝑗 = {𝑥𝑖𝑗𝑘: 1 ≤ 𝑘 ≤ 𝑛𝑖𝑗}. We further assume
that the 𝑛𝑖𝑗 are all equal of size 𝑚 and use 𝑘 as an index to the sample entries.
Fertilizer Wheat Corn Soy Rice
Blend X 123 128 166 151
156 150 178 125
112 174 187 117
100 116 153 155
168 109 195 158
Blend Y 135 175 140 167
130 132 145 183
176 120 159 142
120 187 131 167
155 184 126 168
Blend Z 156 186 185 175
180 138 206 173
147 178 188 154
146 176 165 191
193 190 188 169
Crop

90
We use terms such as 𝑥̅𝑖 (or 𝑥̅𝑖.) as an abbreviation for the mean of {𝑥𝑖𝑗𝑘: 1 ≤ 𝑗 ≤ 𝑟, 1 ≤ 𝑘 ≤ 𝑚}.
We use terms such as 𝑥̅𝑗 (or 𝑥̅.𝑗) as an abbreviation for the mean of {𝑥𝑖𝑗𝑘: 1 ≤ 𝑖 ≤ 𝑐, 1 ≤ 𝑘 ≤ 𝑚}.
We also expand the definitions of sum of squares, mean square, and degrees of freedom for the
one-factor model described earlier to the Two-Factor model as shown in Table 12:
SS Df MS
𝑆𝑆 𝑇 = ∑ ∑ ∑(𝑥𝑖𝑗𝑘 − 𝑥̅)
2
𝑖𝑗𝑘
𝑑𝑓𝑇 = 𝑛 − 1 𝑀𝑆 𝑇 = 𝑆𝑆 𝑇 𝑑𝑓𝑇⁄
𝑆𝑆𝐴 = 𝑚𝑐 ∑(𝑥̅𝑖 − 𝑥̅)2
𝑖
𝑑𝑓𝐴 = 𝑟 − 1 𝑀𝑆𝐴 = 𝑆𝑆𝐴 𝑑𝑓𝐴⁄
𝑆𝑆 𝐵 = 𝑚𝑟 ∑(𝑥̅𝑗 − 𝑥̅)
2
𝑗
𝑑𝑓𝐵 = 𝑐 − 1 𝑀𝑆 𝐵 = 𝑆𝑆 𝐵 𝑑𝑓𝐵⁄
𝑆𝑆𝐴𝐵 = 𝑚 ∑ ∑(𝑥̅𝑖𝑗 − 𝑥̅𝑖 − 𝑥̅𝑗 + 𝑥̅)
2
𝑖𝑗
𝑑𝑓𝐴𝐵 = (𝑟 − 1)(𝑐 − 1) 𝑀𝑆𝐴𝐵 = 𝑆𝑆𝐴𝐵 𝑑𝑓𝐴𝐵⁄
𝑆𝑆 𝑊 = ∑ ∑ ∑(𝑥𝑖𝑗𝑘 − 𝑥̅𝑖𝑗)
2
𝑖𝑗𝑘
𝑑𝑓 𝑊 = 𝑛 − 𝑟𝑐 𝑀𝑆 𝑊 = 𝑆𝑆 𝑊 𝑑𝑓 𝑊⁄
Table 12: Two-Factor ANOVA terminology
Note that, in addition to the A and B factors, the model also contains the intersection between
these factors (labeled AB in the table) which consists of the 𝑚 elements in each of the 𝑟 × 𝑐
combinations of factor levels.
It is not hard to show that:
𝑆𝑆 𝑇 = 𝑆𝑆𝐴 + 𝑆𝑆 𝐵 + 𝑆𝑆𝐴𝐵 + 𝑆𝑆 𝐸
𝑑𝑓𝑇 = 𝑑𝑓𝐴 + 𝑑𝑓𝐵 + 𝑑𝑓𝐴𝐵 + 𝑑𝑓𝐸
It turns out that, if the null hypothesis is true, then 𝑀𝑆 𝑊 and 𝑀𝑆 𝐵 are both measures of the same
error. Thus, the null hypothesis is equivalent to the hypothesis that the population versions of
these statistics are equal, that is:
𝜎 𝐵 = 𝜎 𝑊

91
We can, therefore, use the F-test described in Chapter 8 to determine whether or not to reject
the null hypothesis, that is, if the 𝑥𝑖𝑗 are independently and normally distributed and all 𝜇 𝑗 are
equal (the null hypothesis), and all the 𝜎𝑗
2
are equal (homogeneity of variances), then the test
statistic:
𝐹 =
𝑀𝑆 𝐵
𝑀𝑆 𝑊
has an F distribution with 𝑑𝑓𝐵, 𝑑𝑓 𝑊 degrees of freedom.
Analysis
To carry out the analysis for the last example, we will use the ANOVA: Two Factor with
Replication data analysis tool. To access this tool, as usual, select Data > Analysis|Data
Analysis and fill in the dialog box that appears as in Figure 55:
Figure 55: Dialog box for ANOVA: Two-Factor with Replication analysis tool
After pressing OK, the output shown in Figures 56 and 57 appear:
ANOVA: Two-Factor With Replication
SUMMARY Wheat Corn Soy Rice Total
Blend X
Count 5 5 5 5 20
Sum 659 677 879 706 2921
Average 131.8 135.4 175.8 141.2 146.05
Variance 844.2 707.8 278.7 354.2 782.3658
Blend Y
Count 5 5 5 5 20
Sum 716 798 701 827 3042

92
Average 143.2 159.6 140.2 165.4 152.1
Variance 498.7 978.3 165.7 217.3 511.0421
Blend Z
Count 5 5 5 5 20
Sum 822 868 932 862 3484
Average 164.4 173.6 186.4 172.4 174.2
Variance 443.3 428.8 212.3 175.8 330.6947
Total
Count 15 15 15 15
Sum 2197 2343 2512 2395
Average 146.4667 156.2 167.4667 159.6667
Variance 705.8381 871.0286 605.981 404.9524
Figure 56: Two-Factor ANOVA descriptive statistics
ANOVA
Source of Variation SS df MS F P-value F crit
Sample 8782.9 2 4391.45 9.933347 0.000245 3.190727
Columns 3411.65 3 1137.217 2.572355 0.064944 2.798061
Interaction 6225.9 6 1037.65 2.347138 0.045555 2.294601
Within 21220.4 48 442.0917
Total 39640.85 59
Figure 57: Two-Factor ANOVA
We can now draw some conclusions from the ANOVA table in Figure 57. Since the p-value
(crops) = .0649 > .05 = α, we can’t reject the Factor B null hypothesis. So we can conclude (with
95 percent confidence) that there are no significant differences between the effectiveness of the
fertilizer for the different crops.
Since the p-value (blends) = .00025 < .05 = α, we reject the Factor A null hypothesis and
conclude that there is a statistical difference between the blends.
We also see that the p-value (interactions) = .0456 < .05 = α, and so we conclude that there are
significant differences in the interaction between crop and blend. We can look more carefully at
the interactions by plotting the mean interactions between the levels of the two factors (see
Figure 58). Lines that are roughly parallel indicate the lack of interaction, while lines that are not
roughly parallel indicate interaction.

93
From the first chart, we can see that Brand X has quite a different pattern from the other brands
(especially regarding soy). Although less dramatic, Brand Y is also different from Brand Z
(especially since the line for Brand Y is trending up towards soy but trending down towards
Rice, exactly the opposite of Brand Z):
Figure 58: Interaction plots
Follow-up tests are carried out as for one-way ANOVA, although there are more possibilities.
ANOVA with Repeated Measures
Basic Concepts
ANOVA with repeated measures is an extension of the t test with paired samples test to more
than two variables. As for the paired samples t test, the key difference between this type of
ANOVA and the tests we have considered until now is that the variables under consideration
are not independent.
Some examples of this type analysis are:
 A study is made of 30 subjects, each of whom is asked to take a test on driving a car, a
boat, and an airplane
 A study is made of 30 rats, each of whom is given training once a day for 10 days with
their scores recorded each day
 A study is made of 30 married couples and the husband’s IQ is compared with his wife’s
The important characteristic of each of these examples is that the treatments are not
independent of each other. The most common of these analyses is to compare the results of
some treatment given to the same participant over a period of time (as in the second example).
120
130
140
150
160
170
180
190
Wheat Corn Soy Rice
Blend X
Blend Y
Blend Z
120
130
140
150
160
170
180
190
Blend X Blend Y Blend Z
Wheat
Corn
Soy
Rice

94
Example
A program has been developed to reduce the levels of anxiety for new mothers. In order to
determine whether or not the program is successful, a sample of 15 women was selected and
their level of anxiety was measured (low scores indicate higher levels of anxiety) before the
program. Their level of anxiety was also measured one, two, and three weeks after the
beginning of the program. Based on the data in Figure 59, determine whether the program is
effective in reducing anxiety.
Analysis
To perform ANOVA with repeated measures, we use Excel’s ANOVA: Two Factor without
Replication data analysis tool.
To access this tool, press Data > Analysis|Data Analysis and fill in the dialog box that appears
as in Figure 60:

95
Figure 60: Dialog box for ANOVA: Two-Factor without replication
After pressing OK, the output shown in Figures 61 and 62 appears:
ANOVA: Two-Factor Without Replication
SUMMARY Count Sum Average Variance
1 4 86 21.5 15
2 4 80 20 50
3 4 96 24 12.66667
4 4 80 20 84.66667
5 4 48 12 20.66667
6 4 49 12.25 20.91667
7 4 80 20 16.66667
8 4 97 24.25 14.91667
9 4 78 19.5 27.66667
10 4 90 22.5 28.33333
11 4 78 19.5 28.33333
12 4 107 26.75 6.25
13 4 52 13 23.33333
14 4 55 13.75 40.91667
15 4 29 7.25 20.91667
Before 15 196 13.06667 39.35238
1 week 15 270 18 34.28571
2 weeks 15 278 18.53333 44.69524
3 weeks 15 361 24.06667 26.35238
Figure 61: ANOVA descriptive statistics

96
ANOVA
Source of Variation SS df MS F P-value F crit
Rows 1702.833 14 121.631 15.82722 1.95E-12 1.935009
Columns 910.9833 3 303.6611 39.51389 2.69E-12 2.827049
Error 322.7667 42 7.684921
Total 2936.583 59
Figure 62: Two-Factor ANOVA with repeated measures
For ANOVA with repeated measures, we aren’t interested in the analysis of the rows, only with
the columns, which correspond to variations by time. Since the test statistic F = 29.13 > 2.83 =
F-crit (or p-value = 2.69E-12 < .05 = α) in Figure 62, we reject the null hypothesis. We can
conclude that there are significant differences in mean anxiety levels by week.
Assumptions
The requirements for ANOVA with matched samples are as follows:
 Subjects are independent and randomly selected from the population
 Normality (or at least symmetry)–although overly restrictive, it is sufficient for this
condition to be met for each treatment (that is, for each column in the data range)
 Sphericity
Data satisfies the sphericity property when the pairwise differences in variance between the
samples are all equal. This assumption is stronger than the homogeneity assumption for the
other versions of ANOVA previously described.
We perform the test for sphericity in Figure 63:
Figure 63: Sphericity

97
Although the values in column U of Figure 63 are not equal, they aren’t too far apart. Also,
especially since the p-value obtained in Figure 62 is so low, we don’t expect sphericity to be
much of a problem for this example.
When sphericity is a problem, we can use a variety of correction factors to improve the ANOVA
results. Unfortunately, explaining this will take us too far afield for our purposes here. If you
would like more information, click here.
Follow-up Analysis
Since the mean score of 24.07 in Week 3 is much higher than the mean Before score of 13.07
(see Figure 61), it appears that there is a positive effect to be had from taking the program. We
can confirm this by using a paired sample t test comparing Before with three weeks after, as
shown in Figure 64:
Before 3 weeks
Mean 13.06667 24.06667
Variance 39.35238 26.35238
Observations 15 15
Df 14
t Stat -9.66536
P(T<=t) two-tail 1.42E-07
Figure 64: Paired t test: Before vs. three weeks after
We see that p-value = 1.42E-07, which is considerably less than .05, thus confirming that there
is a significant difference between Before and three weeks after.
In fact, we can see from Figure 65 that there is a significant improvement in the scores even
after one week:

98
Before 1 week
Mean 13.06667 18
Variance 39.35238 34.28571
Observations 15 15
Df 14
t Stat -5.09437
Figure 65: Paired t test: Before vs. one week
In fact, even if we apply Bonferroni’s correction factor using .05/2 = .025 (because we have now
performed two t tests), we can see that p-value = .000163 < .025) which confirms that there is a
significant improvement after only one week.

99
Chapter 10 Correlation and Covariance
Basic Definitions
The sample covariance between two sample random variables 𝑥 and 𝑦 is a measure of the
linear association between the two variables and is defined by the formula:
𝑐𝑜𝑣(𝑥, 𝑦) = ∑(𝑥𝑖 − 𝑥̅)(𝑦𝑖 − 𝑦̅)
𝑛
𝑖=1
(𝑛 − 1)⁄
The covariance is similar to the variance, except that the covariance is defined for two variables
( 𝑥 and 𝑦) whereas the variance is defined for only one variable. In fact, 𝑐𝑜𝑣(𝑥, 𝑥) = 𝑣𝑎𝑟(𝑥).
The correlation coefficient between two sample variables 𝑥 and 𝑦 is a scale-free measure of
linear association between the two variables and is defined by the formula:
𝑟 = 𝑐𝑜𝑣(𝑥, 𝑦) 𝑠 𝑥 𝑠 𝑦⁄
If necessary, we can write 𝑟 as 𝑟𝑥𝑦 to explicitly show the two variables. We also use the term
coefficient of determination for 𝑟2
.
If 𝑟 is close to 1, then 𝑥 and 𝑦 are positively correlated. A positive linear correlation means
that high values of 𝑦 are associated with high values of 𝑥, and low values of 𝑦 are associated
with low values of 𝑥.
If 𝑟 is close to -1, then 𝑥 and 𝑦 are negatively correlated. A negative linear correlation means
that high values of 𝑦 are associated with low values of 𝑥, and low values of 𝑦 are associated
with high values of 𝑥.
When 𝑟 is close to 0, there is little linear relationship between 𝑥 and 𝑦, that is, 𝑥 and 𝑦 are
independent.
Scatter Diagram
To better visualize the association between two data sets {𝑥1, … , 𝑥 𝑛} and {𝑦1, … , 𝑦𝑛}, we can
employ a chart called a scatter diagram (also called a scatter plot).
Example: A study is designed to check the relationship between smoking and longevity. A
sample of 15 men 50 years or older was taken. The average number of cigarettes smoked per
day and the age at death were recorded, as summarized on the left side of Figure 66. Create a
scatter diagram to show the level of association between the number of cigarettes smoked and
a person’s longevity.

100
Figure 66: Scatter diagram
This is done as follows:
1. Highlight the range A4:B18.and select Insert > Charts|Scatter
2. The Design, Layout, and Format ribbons now appear to help you refine the chart. At
any time, you can click the chart to get access to these ribbons.
The resulting chart is shown on the right side of Figure 66 although, initially, the chart does not
contain a chart title or axes titles. We add these as follows:
3. To add a chart title, click the chart, select Layout > Labels|Chart Title, then choose
Above Chart, and enter the title Longevity v. Cigarettes.
4. The title of the horizontal axis can be added in a similar manner by selecting Layout >
Labels|Axis Titles > Primary Horizontal Axis Title > Title Below Axis and entering
Number of Cigarettes Smoked (per day).
5. The title of the vertical axis is added by selecting Layout > Labels|Axis Titles >
Primary Vertical Axis Title > Rotated Title and entering the Life Expectancy (years)
6. The default legend is not particularly useful and so you can delete it by selecting Layout
> Labels|Legend > None
Since the smallest value of Longevity is 53, the chart will display better if we raise the horizontal
axis. This is done as follows:
7. Double-click the vertical axis (0 to 100). Format Axis will appear as shown in Figure 67.
Make sure the Axis Options tab is selected and change the Minimum field from Auto
to Manual and then enter the value 50. Click Close:

101
Figure 67: Format axis dialog box
Finally, we add the straight line that best fits the data points, as follows:
8. Select Layout > Analysis|Trendline > Linear Trendline
If desired, you can move the chart within the worksheet by left-clicking the chart and dragging it
to the desired location. You can also resize the chart, making it a little smaller or bigger, by
clicking one of the corners of the chart and dragging the corner to change the dimensions. To
ensure that the aspect ratio (that is, the ratio of the length to the width) doesn’t change, it is
important to hold Shift down while dragging the corner.
It appears from the scatter diagram in Figure 66 that there is a negative linear correlation
between longevity and cigarette smoking.

102
Excel Functions
Excel supplies the following worksheet functions for ranges R1 and R2 with the same number of
elements:
COVARIANCE.P(R1, R2) = the population covariance between the data in range R1 and
the data in range R2
COVARIANCE.S(R1, R2) = the sample covariance between the data in range R1 and the
data in range R2
CORREL(R1, R2) = the correlation coefficient between the data in range R1 and the data in
range R2
RSQ(R1, R2) = the coefficient of determination between the data in ranges R1 and R2; this
is equivalent to the formula =CORREL(R1, R2) ^ 2
Note that no distinction needs to be made between the population and sample correlation
coefficients since they have the same values.
COVARIANCE.P and COVARIANCE.S are only valid in versions of Excel starting with Excel
2010. For previous releases, Excel only supplies a function for the population covariance,
namely, COVAR(R1, R2) which is equivalent to COVARIANCE.P(R1, R2). The sample
covariance can be calculated using the formula:
=n * COVAR(R1, R2) / (n – 1)
where n = the number of elements in R1 (or R2)
For the data shown in Figure 66, the covariance is -97.44 and the correlation coefficient is
-0.71343. These values are calculated by the formulas:
=COVARIANCE.S(A4:A18,B4:B18)
=CORREL(A4:A18,B4:B18)
Excel Functions with Missing Data
Note that in computing CORREL(R1, R2) or any of the covariance functions described in this
chapter, any pairs for which one or both of the elements in the pair are non-numeric are simply
ignored.
Example: Calculate the correlation between smoking and longevity based on the data in Figure
68:

103
Figure 68: Correlation coefficient with missing data
The data in Figure 68 is identical to that in Figure 66 except that two sample elements have
been added (in rows 11 and 15) where one of the elements in the pair is blank. Since these
elements are ignored, =CORREL(A4:B20) yields the value -.71343, the same as we calculated
from the data in Figure 66.
Hypothesis Testing
It would be useful to know whether or not there is a significant correlation between smoking and
longevity. A correlation of 0 would indicate that smoking and longevity are independent (that is,
no association exists). Since the observed correlation coefficient of -.71343 is pretty far from 0,
it appears from the data that there is a negative correlation. But we would like to know whether
or not this is significant or due to random effects.
It turns out that, when the population correlation is 0, written ρ = 0 (that is, 𝑥 and 𝑦 are
independent), then for sufficiently large values of the sample size n, the sample correlation 𝑟 will
be approximately normally distributed. So we can use a t test with n – 2 degrees of freedom
where:
𝑡 =
𝑟
𝑠 𝑟
and 𝑠 𝑟 = √
1−𝑟2
𝑛−2

104
Here 𝑠 𝑟 is the standard error of 𝑟.
Example: Determine whether or not there is a significant correlation between smoking and
longevity.
The analysis is shown in Figure 69:
Figure 69: Hypothesis testing of the correlation coefficient
Based on the results in Figure 69 (p-value = 0.002822 < .05 = α), we reject the null hypothesis
that ρ = 0 and conclude that there is a significant correlation between smoking and longevity. It
is important to note that such a correlation doesn’t imply that we can claim that smoking causes
lower longevity. For example, perhaps there is some other factor that causes a person to smoke
and to have a shorter life expectancy.
If ρ ≠ 0, then the approach previously described won’t work since, even for large values of n, the
sample correlation 𝑟 is not approximately normally distributed. Fortunately, there is a
transformation of the correlation coefficient which will address this defect.
For any 𝑟, define the Fisher transformation of 𝑟 as follows:
𝑟′ =
1
2
ln
1 + 𝑟
1 − 𝑟
= (ln(1 + 𝑟) − ln(1 − 𝑟)) 2⁄

105
For n is sufficiently large, the Fisher transformation 𝑟′ of the correlation coefficient 𝑟 for samples
of size n has a normal distribution with standard deviation:
𝑠 𝑟′ =
1
√𝑛 − 3
Excel provides functions that calculate the Fisher transformation and its inverse.
FISHER (𝑟) = .5 * LN((1 + 𝑟) / (1 – 𝑟))
FISHERINV(𝑧) = (EXP(2 * 𝑧) – 1) / (EXP(2 * 𝑧) + 1)
Example: Test whether or not the correlation coefficient for the data in Figure 66 is significantly
different from -0.5:
Figure 70: Hypothesis testing using the Fisher transformation
From Figure 70, we see that there is no significant difference between the population correlation
coefficient and -.5 (p-value = 0.098485 > .05 = α). While this is not particularly useful
information, the calculation of the 95 percent confidence interval for the population correlation
coefficient is. This confidence interval is calculated as follows:
𝜌′
= 𝑟′
± 𝑧 𝑐𝑟𝑖𝑡 ∙ 𝑠. 𝑒.
𝜌𝑙𝑜𝑤𝑒𝑟 = FISHERINV(𝜌𝑙𝑜𝑤𝑒𝑟
′ ) 𝜌 𝑢𝑝𝑝𝑒𝑟 = FISHERINV(𝜌 𝑢𝑝𝑝𝑒𝑟
′
)

106
We see that the 95 percent confidence interval is (-.889,-.354). Note that 0 doesn’t lay in this
interval, once again demonstrating that the correlation coefficient is significantly different from 0
(as previously shown). Since -.5 is in this interval, we can’t rule out that the actual population
correlation coefficient is -.5 (with 95 percent confidence).
Data Analysis Tools
Excel provides the Correlation data analysis tool to calculate the pairwise correlation
coefficients for a number of variables.
Example: Calculate the pairwise correlation coefficients and covariance for the data in Figure
71:
We use Excel’s Correlation data analysis tool by selecting Data > Analysis | Data Analysis
and selecting Correlation from the menu that is displayed. The dialog box shown in Figure 72
appears:
Figure 72: Dialog box for the correlation data analysis tool

107
Using Excel’s Correlation data analysis tool, we can compute the pairwise correlation
coefficients for the various variables shown in the table in Figure 71. The results are shown in
Figure 73:
Figure 73: Pairwise correlation coefficients
Similarly, we use Excel’s Covariance data analysis tool to obtain the covariances as shown in
Figure 74:
Figure 74: Pairwise covariances

108
Chapter 11 Linear Regression
Linear Regression Analysis
The goal of linear regression is to create a model from observed data which captures the
relationship between an independent variable 𝑥 and a dependent 𝑦, and to use this model to
predict the values of the dependent variable based on values of the independent variable
(especially for values of the independent variable that were not originally observed). Even when
we can make such predictions, this doesn’t imply that we can claim any causal relationship
between the independent and dependent variables.
Regression Line
Essentially, we are looking for a straight line that best fits the observed data (𝑥1, 𝑦1), … , (𝑥 𝑛, 𝑦𝑛).
Recall that the equation for a straight line is 𝑦 = 𝑏𝑥 + 𝑎, where:
𝑏 = the slope of the line
𝑎 = y-intercept, that is, the value of 𝑥 where the line intersects with the y-axis
Using a technique called ordinary least squares, it turns out that:
𝑏 = 𝑐𝑜𝑣(𝑥, 𝑦) / 𝑣𝑎𝑟(𝑥)
𝑎 = 𝑦̅ − 𝑏𝑥̅
Excel provides the following functions where R1 is a range containing 𝑦 data values and R2 is a
range with 𝑥 data values:
SLOPE(R1, R2) = slope of the regression line as previously described
INTERCEPT(R1, R2) = y-intercept of the regression line as previously described
From our previous observation:
𝑏 = SLOPE(R1, R2) = COVARIANCE.S(R1, R2) / VAR.S(R2)
𝑎 = INTERCEPT(R1, R2) = AVERAGE(R1) – 𝑏 * AVERAGE(R2)

109
Example: Find the regression line for the data in Figure 66 (repeated on the left side of Figure
75). Based on this model, what is the life expectancy of someone who smokes 10, 20, or 30
cigarettes a day?
The results are shown on the right side of Figure 75:
Figure 75: Regression line and predictions
Therefore, the regression line is:
y = -0.6282x + 85.7204
For any value of x (= # of cigarettes smoked), the regression equation will generate a value for y
which is a prediction of the life expectancy of a person who smokes that number of cigarettes.
For example, if a person smokes 30 cigarettes a day, the model forecasts that the person will
live 66.87 years (cell E10 in Figure 75).
Excel provides the following functions to help carry out these forecasts, where R1 is a range
containing 𝑦 data values and R2 is a range with 𝑥 data values:
FORECAST(x, R1, R2) = calculates the predicted value 𝑦 for the given value x of 𝑥
TREND(R1, R2, R3) = array function which predicts the 𝑦 values corresponding to the 𝑥
values in R3
We show how to use FORECAST in cell E9 of Figure 75 and TREND in cell E8 of Figure 75.

110
Actually, TREND is an array function and so it can be used to carry out multiple predictions. In
fact, if we highlight the range E8:E10, enter the formula =TREND(B4:B18,A4:A18,D8;D10), and
then press Ctrl-Shift-Enter, we will get the same results as shown in Figure 75.
Residuals
If range R3 is omitted, then TREND will generate the forecasted values of 𝑦 for the various 𝑥
values in R1. We can use this function to see how well the regression model predicts the 𝑦
values which we have observed. This will help us determine how good the regression model
really is. The result is shown in Figure 76:
Figure 76: Residuals
Here we calculate the forecasted values in the range D4:D18 using the array formula:
=TREND(B4:B18,A4:A18)
The residuals, shown in the range E4:E18, are the differences between the actual 𝑦 values and
the forecasted 𝑦 values. These are the error terms. They can be calculated by inserting the
formula =B4-D4 in cell E4, highlighting the range E4:E18 and pressing Ctrl-D.

111
Although the forecasts are not 100 percent accurate, the mean of the forecasted values (cell
D19) is the same as the mean of the observed values (cell B19), and the mean of the error
terms (cell E19) is zero.
Model Fit
The sample variance of the observed 𝑦 values (Longevity) can be calculated by the formula
=VAR.S(B4:B18) and has a value 120.27. For our purposes, we will call this variance the total
mean square (abbreviated 𝑀𝑆 𝑇). As we have seen in Chapter 3, the sample variance can be
calculated by the formula:
𝑀𝑆 𝑇 =
∑ (𝑦𝑖 − 𝑦̅)2𝑛
𝑖=1
𝑛 − 1
=
DEVSQ(B4: B18)
COUNT(B4: B18) − 1
=
1683.73
14
= 120.27
The numerator is called the total sum of squares (abbreviated 𝑆𝑆 𝑇) and the denominator is the
total degrees of freedom (abbreviated 𝑑𝑓𝑇). The only thing new is the terminology.
It turns out that we can segregate the value of 𝑆𝑆 𝑇 (which represents the total variability of the
observed 𝑦 values) into the portion of variability explained by the model (the regression sum of
squares 𝑆𝑆 𝑅𝑒𝑔) and the portion not explained by the model (the residual sum of squares
𝑆𝑆 𝑅𝑒𝑠) and, similarly, for the degrees of freedom. Thus, we have:
𝑆𝑆 𝑇 = 𝑆𝑆 𝑅𝑒𝑔 + 𝑆𝑆 𝑅𝑒𝑠
𝑑𝑓𝑇 = 𝑑𝑓𝑅𝑒𝑔 + 𝑑𝑓𝑅𝑒𝑠
Now the values of each of these terms are:
𝑆𝑆 𝑇 = ∑ (𝑦𝑖 − 𝑦̅)2
𝑛
𝑖=1
𝑆𝑆 𝑅𝑒𝑔 = ∑ (𝑦̂𝑖 − 𝑦̅)2
𝑛
𝑖=1
𝑆𝑆 𝑅𝑒𝑠 = ∑ (𝑦𝑖 − 𝑦̂)2
𝑛
𝑖=1
𝑑𝑓𝑇 = 𝑛 − 1 𝑑𝑓𝑅𝑒𝑔 = 𝑘 𝑑𝑓𝑅𝑒𝑠 = 𝑛 − 𝑘 − 1
where the 𝑦𝑖 are the observed values for 𝑦, the 𝑦̂𝑖 are values of 𝑦 predicted by the model and k
= the number of independent variables in the model (k = 1 for now). We can also introduce
mean squares values:
𝑀𝑆 𝑇 =
𝑆𝑆 𝑇
𝑑𝑓𝑇
𝑀𝑆 𝑅𝑒𝑔 =
𝑆𝑆 𝑅𝑒𝑔
𝑑𝑓𝑅𝑒𝑔
𝑀𝑆 𝑅𝑒𝑠 =
𝑆𝑆 𝑅𝑒𝑠
𝑑𝑓𝑅𝑒𝑠

112
Clearly the model fits the observed data well if 𝑆𝑆 𝑅𝑒𝑠 is small in comparison to 𝑆𝑆 𝑅𝑒𝑔. Since
𝑆𝑆 𝑇 = 𝑆𝑆 𝑅𝑒𝑔 + 𝑆𝑆 𝑅𝑒𝑠, this is equivalent to saying that
𝑆𝑆 𝑇
is relatively close to 1. It turns out
that this fraction is equal to the coefficient of determination, that is, 𝑟2
. It is common to use a
capital 𝑅 instead of a small 𝑟 and so we have:
𝑅2
=
𝑆𝑆 𝑇
All these statistics can be calculated in Excel as shown in Figure 77:
Figure 77: Calculation of regression parameters in Excel
The formulas used are shown in Table 13 (based on the data Figure 76):
Cell Entity Formula
H4 𝑆𝑆 𝑇 =DEVSQ(B4:B18)
H5 𝑆𝑆 𝑅𝑒𝑔 =DEVSQ(D4:D18)
H6 𝑆𝑆 𝑅𝑒𝑠 =DEVSQ(E4:E18)
I4 𝑑𝑓𝑇 =COUNT(B4:B18)-1
I5 𝑑𝑓𝑅𝑒𝑔 =I5
I6 𝑑𝑓𝑅𝑒𝑠 =I4-I5

113
Cell Entity Formula
J4 𝑀𝑆 𝑇 =H4/I4
J5 𝑀𝑆 𝑅𝑒𝑔 =H5/I5
J6 𝑀𝑆 𝑅𝑒𝑠 =H6/I6
M4 𝑅2
=H5/H4
M5 𝑅 =SQRT(M4)
Table 13: Regression formulas
The residual mean squares term 𝑀𝑆 𝑅𝑒𝑠 is a sort of variance for the error not captured by the
model. The standard error of the estimate is simply the square root of 𝑀𝑆 𝑅𝑒𝑠. For the current
example, this is the square root of 63.596 (cell J6 of Figure 77), that is, 7.975. Excel provides
the following function to calculate the standard error of the estimate:
STEYX(R1, R2) = standard error of the estimate for the regression model where R1 is a
range containing the observed 𝑦 data values and R2 is a range containing the observed 𝑥
data values
Multiple Regression Analysis
We now extend the linear regression model previously introduced with one independent variable
to the case where there are one or more independent variables:
If 𝑦 is a dependent variable and 𝑥1, … , 𝑥 𝑘 are independent variables, then the general multiple
regression model provides a prediction of 𝑦 from the 𝑥𝑖 of the form:
𝑦 = 𝛽0 + 𝛽1 𝑥1 + 𝛽2 𝑥2 + ⋯ + 𝛽 𝑘 𝑥 𝑘 + 𝜀
where 𝛽0, 𝛽1, … , 𝛽 𝑘 are the regression coefficients of the model and 𝜀 is the random error.
We further assume that for any given values of the 𝑥𝑖 the random error 𝜀 is normally and
independently distributed with mean zero. Essentially, the last of these conditions means that
the model is a good fit for the data even if we eliminate the random error term.

114
Example: A jeweler prices diamonds on the basis of quality (with values from 0 to 8, with 8
being flawless and 0 containing numerous imperfections) and color (with values from 1 to 10,
with 10 being pure white and 1 being yellow). Based on the price per carat of the 11 diamonds,
each weighing between 1.0 and 1.5 carats as shown in Figure 78, build a regression model
which captures the relationship between quality, color, and price.
Based on this model, estimate the price for a diamond whose color is 4 and whose quality is 6:
We use Excel’s Regression analysis tool by pressing Data > Analysis | Data Analysis and
selecting Regression from the menu that is displayed. The dialog box shown in Figure 79
appears:

115
Figure 79: Dialog box for regression data analysis
After filling in the fields as shown in the figure and selecting OK, the output shown in Figure 80
will be displayed.
The calculations of the values in Figure 80 are extensions of the approaches described earlier in
this chapter for the case where there is only one independent variable (see especially Table 13).
Adjusted R Square (cell B6) is an attempt to create a better estimate of the true value of the
population coefficient of determination, and is calculated by the formula:
=1-(1-B5)*(B8-1)/(B8-B12-1)
The calculations of the regression coefficients (B17:B19) and the corresponding standard errors
(C17:C19) use the method of least squares but is beyond the scope of this book. The
calculation of the p-values and confidence intervals for these coefficients is as described in
Chapter 7.
One of the key results of the analysis in Figure 80 is that the regression model is given by the
formula:
Price = 1.7514 + 4.8953 ∙ Color + 3.7584 ∙ Quality

116
Figure 80: Regression data analysis
Since the p-value (cell F12) = 0.00497 < .05, we reject the null hypothesis. We conclude that the
regression model is a good fit for the data.
Note that the Color and Quality coefficients are significant (p-value < .05), that is, they are
significantly different from zero. Meanwhile, the intercept coefficient is not significantly different
from zero (p-value = .8077 > .05).
That R square = .8507 indicates that a good deal of the variability of price is captured by the
model, which supports the case that the regression model is a good fit for the data.
If a diamond has color 4 and quality 6, then we can use the model to estimate the price as
follows:
Price = 1.75 + 4.90 ∙ Color + 3.76 ∙ Quality = 1.75 + 4.90 ∙ 4 + 3.76 ∙ 6 = 43.88305
We get the same result using the TREND function as shown in Figure 81:
Figure 81: Prediction using regression model

Statistics using excel_succinctly

More Related Content

Similar to Statistics using excel_succinctly (20)

Recently uploaded (20)

Statistics using excel_succinctly