“Science is advanced by proposing and testing a hypothesis, not by declaring questions unsolvable” – Nick Matzke
Let’s start with a case study. I want you to think of your favorite restaurant right now. Let’s say you can predict a certain number of people arriving for lunch five days a week. At the end of the week, you observe that the expected footfall was different from the actual footfall.
Sounds like a prime statistics problem? That’s the idea!
So, how will you check the statistical significance between the observed and the expected footfall values? Remember this is a categorical variable – ‘Days of the week’ – with 5 categories [Monday, Tuesday, Wednesday, Thursday, Friday].
One of the best ways to deal with this is by using the Chi-Square test.
We can always opt for z-tests, t-tests or ANOVA when we’re dealing with continuous variables. But the situation becomes tricky when working with categorical features (as most data scientists will attest to!). I’ve found the chi-square test to be quite helpful in my own projects.
So let’s dive into the article to understand all about the chi-square test, what it is, how it works and how we can implement it in R.
If you are new statistics and data science, I would recommend the below resources to get a comprehensive overview of the two broad topics:
I’m sure you’ve encountered categorial variables before, even if you might not have intuitively recognized them. They can be tricky to deal with in the data science world so let’s first define them.
Categorical variables fall into a particular category of those variables that can be divided into finite categories. These categories are generally names or labels. These variables are also called qualitative variables as they depict the quality or characteristics of that particular variable.
For example, the category “Movie Genre” in a list of movies could contain the categorical variables – “Action”, “Fantasy”, “Comedy”, “Romance”, etc.
There are broadly two types of categorical variables:
When the data we want to analyze contains this type of variable, we turn to the chi-square test, denoted by χ², to test our hypothesis.
A Chi-Square test is a test of statistical significance for categorical variables.
Let’s learn the use of chi-square with an intuitive example.
A research scholar is interested in the relationship between the placement of students in the statistics department of a reputed University and their C.G.P.A (their final assessment score).
He obtains the placement records of the past five years from the placement cell database (at random). He records how many students who got placed fell into each of the following C.G.P.A. categories – 9-10, 8-9, 7-8, 6-7, and below 6.
Source: Anibrain School of Media Design
If there is no relationship between the placement rate and the C.G.P.A., then the placed students should be equally spread across the different C.G.P.A. categories (i.e. there should be similar numbers of placed students in each category).
However, if students having C.G.P.A more than 8 are more likely to get placed, then there would be a large number of placed students in the higher C.G.P.A. categories as compared to the lower C.G.P.A. categories. In this case, the data collected would make up the observed frequencies.
So the question is, are these frequencies being observed by chance or do they follow some pattern?
Here enters the chi-square test! The chi-square test helps us answer the above question by comparing the observed frequencies to the frequencies that we might expect to obtain purely by chance.
Chi-square test in hypothesis testing is used to test the hypothesis about the distribution of observations/frequencies in different categories.
Note: I strongly recommend going through the below article if you need to brush up your hypothesis testing concepts:
We are almost at the implementing aspect of chi-square tests but there’s one more thing we need to learn before we get there.
Just like any other statistical test, the chi-square test comes with a few assumptions of its own:
This is a non-parametric test. We typically use it to find how the observed value of a given event is significantly different from the expected value. In this case, we have categorical data for one independent variable, and we want to check whether the distribution of the data is similar or different from that of the expected distribution.
Let’s consider the above example where the research scholar was interested in the relationship between the placement of students in the statistics department of a reputed University and their C.G.P.A.
In this case, the independent variable is C.G.P.A with the categories 9-10, 8-9, 7-8, 6-7, and below 6.
The statistical question here is: whether or not the observed frequencies of placed students are equally distributed for different C.G.P.A categories (so that our theoretical frequency distribution contains the same number of students in each of the C.G.P.A categories).
We will arrange this data by using the contingency table which will consist of both the observed and expected values as below:
C.G.P.A | ||||||
10-9 | 9-8 | 8-7 | 7-6 | Below 6 | Total | |
Observed Frequency of Placed students |
30 | 35 | 20 | 10 | 5 | 100 |
Expected Frequency of Placed students |
20 | 20 | 20 | 20 | 20 | 100 |
After constructing the contingency table, the next task is to compute the value of the chi-square statistic. The formula for chi-square is given as:
where,
Let us look at the step-by-step approach to calculate the chi-square value:
In this case, the degrees of freedom are 5-1 = 4. So, the critical value at 5% level of significance is 9.49.
Our obtained value of 32.5 is much larger than the critical value of 9.49. Therefore, we can say that the observed frequencies are significantly different from the expected frequencies. In other words, C.G.P.A is related to the number of placements that occur in the department of statistics.
Let’s further solidify our understanding by performing the Chi-Square test in R.
Let’s implement the chi-square goodness of fit test in R. Time to fire up RStudio!
Problem Statement
Let’s understand the problem statement before we dive into R.
An organization claims that the experience of the employees of different departments is distributed in the following categories:
A random sample of 1470 employees is collected. Does this random sample provide evidence against the organization’s claim?
You can download the data here.
Setting up the hypothesis
Let’s begin!
Step 1: First, import the data
Step 2: Validate it for correctness in R:
Output:
#Count of Rows and columns [1] 1470 2 #View top 10 rows of the dataset age.intervals Experience.intervals 1 41 - 50 6 - 10 Years 2 41 - 50 6 - 10 Years 3 31 - 40 6 - 10 Years 4 31 - 40 6 - 10 Years 5 18 - 30 6 - 10 Years 6 31 - 40 6 - 10 Years 7 51 - 60 11 - 20 Years 8 18 - 30 Upto 5 Years 9 31 - 40 6 - 10 Years 10 31 - 40 11 - 20 Years
Step 3: Create a proportion table for expected frequencies:
Output:
11 - 20 Years 21 - 40 Years 6 - 10 Years Upto 5 Years 0.2312925 0.1408163 0.4129252 0.2149660
Step – 4: Calculate the chi-square value:
Output:
Chi-square test for given probabilities data: table(data$Experience.intervals) X-squared = 14.762, df = 3, p-value = 0.002032
The p-value here is less than 0.05. Therefore, we will reject our null hypothesis. Hence, the distribution of experience of the employees of different departments differs from what the organization states.
The second type of chi-square test is the Pearson’s chi-square test of association. This test is used when we have categorical data for two independent variables and we want to see if there is any relationship between the variables.
Let’s take another example to understand this. A teacher wants to know the answer to whether the outcome of a mathematics test is related to the gender of the person taking the test. Or in other words, she wants to know if males show a different pattern of pass/fail rates than females.
So, here are two categorical variables: Gender (Male and Female) and mathematics test outcome (Pass or Fail). Let us now look at the contingency table:
Boys | Girls | |
Pass | 17 | 20 |
Fail | 8 | 5 |
By looking at the above contingency table, we can see that the girls have a comparatively higher pass rate than boys. However, to test whether this observed difference is significant or not, we will carry out the chi-square test.
The steps to calculate the chi-square value are as follows:
Step 1: Calculate the row and column total of the above contingency table:
Boys | Girls | Total | |
Pass | 17 | 20 | 37 |
Fail | 8 | 5 | 13 |
Total | 25 | 25 | 50 |
Step 2: Calculate the expected frequency for each individual cell by multiplying row sum by column sum and dividing by total number:
Expected Frequency = (Row Total x Column Total)/Grand Total
For the first cell, the expected frequency would be (37*25)/50 = 18.5. Now, write them below the observed frequencies in brackets:
Boys | Girls | Total | |
Pass | 17 (18.5) |
20 (18.5) |
37 |
Fail | 8 (6.5) |
5 (6.5) |
13 |
Total | 25 | 25 | 50 |
Step 3: Calculate the value of chi-square using the formula:
Calculate the right-hand side part of each cell. For example, for the first cell, ((17-18.5)^2)/18.5 = 0.1216.
Step 4: Then, add all the values obtained for each cell. In this case, the values are:
0.1216+0.1216+0.3461+0.3461 = 0.9354
Step 5: Calculate the degrees of freedom, i.e. (Number of rows-1)*(Number of columns-1) = 1*1 = 1
The next task is to compare it with the critical chi-square value from the table we saw above.
The Chi-Square calculated value is 0.9354 which is less than the critical value of 3.84. So in this case, we fail to reject the null hypothesis. This means there is no significant association between the two variables, i.e, boys and girls have a statistically similar pattern of pass/fail rates on their mathematics tests.
Let’s further solidify our understanding by performing the chi-square test in R.
Problem statement
A Human Resources department of an organization wants to check whether age and experience of the employees are dependent on each other. For this purpose, a random sample of 1470 employees is collected with their age and experience. You can download the data here.
Setting up the hypothesis
Null hypothesis: Age and Experience are two independent variables
Alternative hypothesis: Age and Experience are two dependent variables
Let’s begin!
Step 1: First, import the data
Step 2: Validate it for correctness in R:
Output:
#Count of Rows and columns [1] 1470 2 > #View top 10 rows of the dataset age.intervals Experience.intervals 1 41 - 50 6 - 10 Years 2 41 - 50 6 - 10 Years 3 31 - 40 6 - 10 Years 4 31 - 40 6 - 10 Years 5 18 - 30 6 - 10 Years 6 31 - 40 6 - 10 Years 7 51 - 60 11 - 20 Years 8 18 - 30 Upto 5 Years 9 31 - 40 6 - 10 Years 10 31 - 40 11 - 20 Years
Step 3: Construct a contingency table and calculate the chi-square value:
Output:
ct<-table(data$age.intervals,data$Experience.intervals) > ct 11 - 20 Years 21 - 40 Years 6 - 10 Years Upto 5 Years 18 - 30 22 0 172 192 31 - 40 190 20 308 101 41 - 50 85 112 110 15 51 - 60 43 75 17 8 > chisq.test(ct) Pearson's Chi-squared test data: ct X-squared = 679.97, df = 9, p-value < 2.2e-16
The p-value here is less than 0.05. Therefore, we will reject our null hypothesis. We can conclude that age and experience are two dependent variables, aka as the experience increases, the age also increases (and vice versa).
In this article, we learned how to analyze the significant difference between data that contains categorical measures in it with the help of chi-square tests. We enhanced our knowledge on the use of chi-square, assumptions involved in carrying out the test, and how to conduct different types of chi-square tests both manually and in R.
If you are new to statistics, want to cover your basics, and also want to get a start in data science, I recommend taking the Introduction to Data Science course. It gives you a comprehensive overview of both descriptive and inferential statistics before diving into data science techniques.
Did you find this article useful? Can you think of any other applications of the chi-square test? Let me know in the comments section below and we can come up with more ideas!
Lorem ipsum dolor sit amet, consectetur adipiscing elit,
10 Nov 23 • 08:00pm
Very simple yet good explanation of Chi Square test. Thanks for the post