Numba-accelerated implementations of common probability distributions
Project description
numba-stats
We provide numba-accelerated implementations of statistical functions for common probability distributions
- Uniform
- (Truncated) Normal
- Log-normal
- Poisson
- (Truncated) Exponential
- Student's t
- Voigtian
- Crystal Ball
- Generalised double-sided Crystal Ball
- Tsallis-Hagedorn, a model for the minimum bias pT distribution
- Q-Gaussian
- Bernstein density (not normalised to unity, use this in extended likelihood fits)
with more to come. The speed gains are huge, up to a factor of 100 compared to scipy. Benchmarks are included in the repository and are run by pytest.
Usage
Each distribution is implemented in a submodule. Import the submodule that you need.
from numba_stats import norm
import numpy as np
x = np.linspace(-10, 10)
mu = 2
sigma = 3
dp = norm.pdf(x, mu, sigma)
p = norm.cdf(x, mu, sigma)
The functions are vectorised on the variate x, but not on the shape parameters of the distribution. Ideally, the following functions are implemented for each distribution:
logpdfpdfcdfppf
cdf and ppf are missing for some distributions (e.g. voigt), if there is currently no fast implementation available. logpdf is only implemented if it is more efficient and accurate compared to computing log(dist.pdf(...)).
Documentation
To get documentation, please use help() in the Python interpreter.
Functions with equivalents in scipy.stats follow the Scipy calling conventions exactly. These conventions are sometimes a bit unusual, for example, in case of the exponential, the log-normal or the uniform distribution. See the SciPy docs for details.
Contributions
You can help with adding more distributions, patches are very welcome. Implementing a probability distribution is easy. You need to write it in simple Python that numba can understand. Special functions from scipy.special can be used after some wrapping, see submodule numba_stats._special.py how it is done.
numba-stats and numba-scipy
numba-scipy is the official package and repository for fast numba-accelerated scipy functions, are we reinventing the wheel?
Ideally, the functionality in this package should be in numba-scipy and we hope that eventually this will be case. In this package, we don't offer overloads for scipy functions and classes like numba-scipy does. This simplifies the implementation dramatically. numba-stats is intended as a temporary solution until fast statistical functions are included in numba-scipy. numba-stats currently does not depend on numba-scipy, only on numba and scipy.
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
File details
Details for the file numba-stats-1.1.1.tar.gz.
File metadata
- Download URL: numba-stats-1.1.1.tar.gz
- Upload date:
- Size: 18.7 kB
- Tags: Source
- Uploaded using Trusted Publishing? No
- Uploaded via: twine/4.0.2 CPython/3.9.16
File hashes
| Algorithm | Hash digest | |
|---|---|---|
| SHA256 | e1932a349c0d3f6a5988f903a9895ceac6537afb271afe5360aa1ea6005c7c25 |
|
| MD5 | ede01fd0433f34792efb6f4ee6cc9487 |
|
| BLAKE2b-256 | 58f368d9fac0374598b8a136739d1e61ed038548f83926236b37021dd60d89e5 |
File details
Details for the file numba_stats-1.1.1-py3-none-any.whl.
File metadata
- Download URL: numba_stats-1.1.1-py3-none-any.whl
- Upload date:
- Size: 17.6 kB
- Tags: Python 3
- Uploaded using Trusted Publishing? No
- Uploaded via: twine/4.0.2 CPython/3.9.16
File hashes
| Algorithm | Hash digest | |
|---|---|---|
| SHA256 | 3a7d6da1547ebb401a65682c550ce78e48f175af48b7423b9cef68d8d3ba425d |
|
| MD5 | d21a7f7bcc8f853ba12b1e69ae8ca01b |
|
| BLAKE2b-256 | eafabe441940c7932642e71f5b40fbf9ea56fc1d409c043d9679c398f82642ba |
