PDToolbox is a collection of methods helpful for doing probability distribution computations in Python. The aim of the PDToolbox Python module is to provide a set of features, based on simple probability distributions, that are not available from the scipy.stats module. This includes fast batch computations of (weighted) maximum likelihood estimates, computation of critical empirical distribution statistics, and more niche probability distributions or related code in the pdtoolbox.special module. The module contains code that is described in (ADD citations of the two articles).