4.5 C
London
Friday, December 20, 2024
HomePythonChi-Square Tests in PythonHow to Perform a Chi-Square Goodness of Fit Test in Python

How to Perform a Chi-Square Goodness of Fit Test in Python

Related stories

Learn About Opening an Automobile Repair Shop in India

Starting a car repair shop is quite a good...

Unlocking the Power: Embracing the Benefits of Tax-Free Investing

  Unlocking the Power: Embracing the Benefits of Tax-Free Investing For...

Income Splitting in Canada for 2023

  Income Splitting in Canada for 2023 The federal government’s expanded...

Can I Deduct Home Office Expenses on my Tax Return 2023?

Can I Deduct Home Office Expenses on my Tax...

Canadian Tax – Personal Tax Deadline 2022

  Canadian Tax – Personal Tax Deadline 2022 Resources and Tools...

A Chi-Square Goodness of Fit Test is used to determine whether or not a categorical variable follows a hypothesized distribution.

This tutorial explains how to perform a Chi-Square Goodness of Fit Test in Python.

Example: Chi-Square Goodness of Fit Test in Python

A shop owner claims that an equal number of customers come into his shop each weekday. To test this hypothesis, a researcher records the number of customers that come into the shop in a given week and finds the following:

  • Monday: 50 customers
  • Tuesday: 60 customers
  • Wednesday: 40 customers
  • Thursday: 47 customers
  • Friday: 53 customers

Use the following steps to perform a Chi-Square goodness of fit test in Python to determine if the data is consistent with the shop owner’s claim.

Step 1: Create the data.

First, we will create two arrays to hold our observed and expected number of customers for each day:

expected = [50, 50, 50, 50, 50]
observed = [50, 60, 40, 47, 53]

Step 2: Perform the Chi-Square Goodness of Fit Test.

Next, we can perform the Chi-Square Goodness of Fit Test using the chisquare function from the SciPy library, which uses the following syntax:

chisquare(f_obs, f_exp) 

where:

  • f_obs: An array of observed counts.
  • f_exp: An array of expected counts. By default, each category is assumed to be equally likely.

The following code shows how to use this function in our specific example:

import scipy.stats as stats

#perform Chi-Square Goodness of Fit Test
stats.chisquare(f_obs=observed, f_exp=expected)

(statistic=4.36, pvalue=0.35947)

The Chi-Square test statistic is found to be 4.36 and the corresponding p-value is 0.35947.

Note that the p-value corresponds to a Chi-Square value with n-1 degrees of freedom (dof), where n is the number of different categories. In this case, dof = 5-1 = 4. You can use the Chi-Square to P Value Calculator to confirm that the p-value that corresponds to X2 = 4.36 with dof = 4 is 0.35947.

Recall that a Chi-Square Goodness of Fit Test uses the following null and alternative hypotheses:

  • H0: (null hypothesis) A variable follows a hypothesized distribution.
  • H1: (alternative hypothesis) A variable does not follow a hypothesized distribution.

Since the p-value (.35947) is not less than 0.05, we fail to reject the null hypothesis. This means we do not have sufficient evidence to say that the true distribution of customers is different from the distribution that the shop owner claimed.

Subscribe

- Never miss a story with notifications

- Gain full access to our premium content

- Browse free from up to 5 devices at once

Latest stories