Skip to main content

L-Moments for robust statistics.

Project description

jorenham/lmo

Lmo - Trimmed L-moments and L-comoments

GitHub Workflow Status (with branch) PyPI versions license


Is your tail too heavy? 
Can't find a moment? 
Are the swans black? 
The distribution pathological?

... then look no further: Lmo's got you covered!

Uniform or multi-dimensional, Lmo can summarize it all with one quick glance!

Unlike the legacy moments, L-moments uniquely describe a probability distribution, and are more robust and efficient. The "L" stands for Linear; it is a linear combination of order statistics. So Lmo is as fast as sorting your samples (in terms of time-complexity).

Key Features

  • Calculates trimmed L-moments and L-comoments, from samples or any scipy.stats distribution.
  • Full support for trimmed L-moment (TL-moments), e.g. lmo.l_moment(..., trim=(1/137, 3.1416)).
  • Method of L-moments for robust distribution fitting.
  • Fast estimation of L-comoment matrices from your multidimensional data or multivariate distribution.
  • Goodness-of-fit test, using L-moment or L-moment ratio's.
  • Non-parametric estimation of continuous distributions with lmo.l_rv_nonparametric
  • Exact (co)variance structure of the sample- and population L-moments.
  • Theoretical & empirical influence functions of L-moments & L-ratio's.
  • Complete docs, including detailed API reference with usage examples and with mathematical $\TeX$ definitions.
  • Clean Pythonic syntax for ease of use.
  • Vectorized functions for very fast fitting.
  • Fully typed, tested, and tickled.

Quick example

Even if your data is pathological like Cauchy, and the L-moments are not defined, the trimmed L-moments (TL-moments) can be used instead. Let's calculate the TL-location and TL-scale of a small amount of samples:

>>> import numpy as np
>>> import lmo
>>> rng = np.random.default_rng(1980)
>>> x = rng.standard_cauchy(96)  # pickle me, Lmo
>>> lmo.l_moment(x, [1, 2], trim=(1, 1)).
array([-0.17937038,  0.68287665])

Now compare with the theoretical standard Cauchy TL-moments:

>>> from scipy.stats import cauchy
>>> cauchy.l_moment([1, 2], trim=(1, 1))
array([0.        , 0.69782723])

See the documentation for more examples and the API reference.

Roadmap

  • Automatic trim-length selection.
  • Plotting utilities (deps optional), e.g. for L-moment ratio diagrams.

Installation

Lmo is on PyPI, so you can do something like:

pip install lmo

Dependencies

  • python >= 3.10
  • numpy >= 1.22
  • scipy >= 1.9

Foundational Literature

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

lmo-0.10.2.tar.gz (57.9 kB view details)

Uploaded Source

Built Distribution

lmo-0.10.2-py3-none-any.whl (62.4 kB view details)

Uploaded Python 3

File details

Details for the file lmo-0.10.2.tar.gz.

File metadata

  • Download URL: lmo-0.10.2.tar.gz
  • Upload date:
  • Size: 57.9 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: poetry/1.6.1 CPython/3.11.5 Linux/6.4.6-76060406-generic

File hashes

Hashes for lmo-0.10.2.tar.gz
Algorithm Hash digest
SHA256 118adc6a59630a3196757c5eb6799c77bcc9b7cdc6fcaa10b8090a2222a4805f
MD5 ccaf58230130181c48e81b896018b4b2
BLAKE2b-256 c38d538a1895e62e9bd86eccaa156611d68856b1ee2063b716e956fa9ca3fb3a

See more details on using hashes here.

File details

Details for the file lmo-0.10.2-py3-none-any.whl.

File metadata

  • Download URL: lmo-0.10.2-py3-none-any.whl
  • Upload date:
  • Size: 62.4 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: poetry/1.6.1 CPython/3.11.5 Linux/6.4.6-76060406-generic

File hashes

Hashes for lmo-0.10.2-py3-none-any.whl
Algorithm Hash digest
SHA256 35424552820dfa495b6f3640f9c9ccc60f74c94ed9f433e142e27389e2127a5a
MD5 144c3ed45ba0b048a8f1c5dbfe0026af
BLAKE2b-256 48f377e590541195da57059cf90d02d5dd63cb2524f358fe7e0155f55d89f90f

See more details on using hashes here.

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page