Standard deviation is a key concept in statistics, measuring the amount of variation or dispersion in a dataset’s values. In the built-in Python statistics module exists the statistics.stdev() function, which calculates the standard deviation of a sample dataset. This article aims to guide you through the use of this function, including its syntax, various examples, and handling common errors.

Using the statistics.stdev() Function

Before using the statistics.stdev() function, you need to import the statistics module.

import statistics

The statistics.stdev() function is designed to calculate the sample standard deviation of a dataset.

Recall the distinction between sample standard deviation and population standard deviation: the sample standard deviation is calculated using a subset (sample) of a population, whereas the population standard deviation uses all members of a population, the computational difference being that the formula for sample standard deviation divides by (n – 1) — where n is the number of observations — to account for the fact that we are estimating the population standard deviation from a sample. Population standard deviation uses n as the denominator.

Moving on, the syntax for the statistics.stdev() function is as follows:

statistics.stdev(data, xbar=None)

• data: This is the primary argument and must be an iterable (such a list or tuple) containing numeric values
• xbar (optional): This parameter allows you to specify the mean of the data; if not provided, the function will compute the mean internally

Let’s have a look at some examples of using the function. We will start with a simple example using a list of numbers.

import statistics

data = [10, 12, 23, 23, 16, 23, 21, 16]
stdev_value = statistics.stdev(data)
print(stdev_value)

# Output: 5.237229365663818

When you run this code, it will output the standard deviation of the given dataset. The output will be a numeric value representing the sample standard deviation.

As you know, outliers can significantly impact the standard deviation. Here’s an example where the dataset includes an outlier:

import statistics

data_with_outliers = [10, 12, 23, 23, 16, 100, 21, 16]
stdev_with_outliers = statistics.stdev(data_with_outliers)
print(stdev_with_outliers)

# Output: 29.64040437356703

In this case, the presence of the outlier (100) will increase the standard deviation compared to the previous example, demonstrating how outliers can distort the measure of variability in a dataset.

You may want to use the xbar parameter when you already know the mean of your dataset, which can save computation time. This can be especially useful in very large datasets, and when performing multiple statistical calculations in the same script which may rely on computing the same metrics multiple times. Here’s how we do it.

import statistics

data = [10, 12, 23, 23, 16, 23, 21, 16]
mean_value = sum(data) / len(data)
stdev_with_xbar = statistics.stdev(data, xbar=mean_value)
print(stdev_with_xbar)

# Output: 5.237229365663818

The statistics.stdev() can handle different iterable types, including lists and tuples. Here’s an example using a tuple.

import statistics

data_tuple = (10, 12, 23, 23, 16)
stdev_tuple = statistics.stdev(data_tuple)
print(stdev_tuple)

# Output: 6.058052492344384

This code will yield the same result as if you used a list, demonstrating the flexibility of the function.

While using statistics.stdev(), trying to calculate the standard deviation of an empty dataset, Python will raise a StatisticsError. Here’s how you can handle such errors.

import statistics
from statistics import StatisticsError

try:
    invalid_data = []
    statistics.stdev(invalid_data)
except StatisticsError as e:
    print(f"Error: {e}")
	
# Output: Error: stdev requires at least two data points

Wrapping Up

In this article, we have covered the essentials of using the statistics.stdev() function in Python. We discussed the concept of standard deviation, how to import the necessary module, the function’s syntax, and provided examples that demonstrate its application. For more information, consult the statistics module documentation.