How to Calculate Standard Deviation of a Matrix in Python

Understanding data variability is fundamental to statistical analysis and machine learning. Calculating the standard deviation of a matrix helps you measure how spread out values are from the mean, enabling better feature engineering, outlier detection, and data normalization. In this guide, you'll learn how to efficiently compute standard deviation across entire matrices, rows, and columns using NumPy, the Python's most powerful numerical computing library.

Global Standard Deviation

By default, the np.std() function treats the entire matrix as a single list of numbers and returns a single "global" standard deviation value.

import numpy as np

matrix = np.array([
    [10, 20, 30],
    [40, 50, 60],
    [70, 80, 90]
])

global_sd = np.std(matrix)

print(f"Global Standard Deviation: {global_sd:.2f}")

Output:

Global Standard Deviation: 25.82

Row-wise and Column-wise Calculations

For advanced feature engineering, you often need to analyze the variance across specific dimensions. You can control this using the axis parameter:

• axis=0: Calculates the standard deviation of each column
• axis=1: Calculates the standard deviation of each row

import numpy as np

matrix = np.array([
    [10, 20, 30],
    [40, 50, 60],
    [70, 80, 90]
])

col_sd = np.std(matrix, axis=0)

row_sd = np.std(matrix, axis=1)

print(f"Column Standard Deviations: {col_sd}")
print(f"Row Standard Deviations:    {row_sd}")

Output:

Column Standard Deviations: [24.49489743 24.49489743 24.49489743]
Row Standard Deviations:    [8.16496581 8.16496581 8.16496581]

Understanding the Axis Parameter

Thinking about axes can be confusing at first. Use this mental model:

• Axis 0: Operations move down the rows, collapsing them into a single row of results
• Axis 1: Operations move across the columns, collapsing them into a single column of results

By mastering axis-aware calculations in NumPy, you can extract precise variability metrics from any multi-dimensional dataset with minimal code.