This project is python code, with examples, for computing rate spectra from noisy time series data. While numerical inverse Laplace transforms are typically ill-posed, this code provides a straightforward method of regularization, by sampling from a posterior distribution of possible models.