pava

Greatest convex majorants and least concave minorants https://github.com/dschrempf/pava#readme

 Stackage Nightly 2020-09-29: 0.1.0.0 Latest on Hackage: 0.1.0.0

See all snapshots pava appears in

Module documentation for 0.1.0.0

This version can be pinned in stack with:pava-0.1.0.0@sha256:a069338665384de1f85541c3b881a187d755a7e4f8b5d75833d0491c352ac474,1805

Compute greatest convex majorants and least concave minorants using the pool adjacent violators algorithm.

The pool adjacent violators algorithm (PAVA) is an iterative algorithm for solving monotonic regression problems. In particular, (antitonic) regression is the computation of a non-decreasing (non-increasing) sequence of values such that a given problem is optimized. PAVA can also be used to compute the greatest convex minorant and the least concave majorant of a given set of observables.

At the moment, greatest convex majorants and least concave minorants can be computed efficiently. More general isotonic regression is not yet supported, but may be in future releases.

Run times

Run times with random exponential vectors of different length:

stack bench 2>&1

pava> benchmarks
Running 1 benchmarks...
Benchmark pava-bench: RUNNING...
benchmarking Greatest convex minorant/Vector of length 1e3
time                 397.8 μs   (384.9 μs .. 421.6 μs)
0.984 R²   (0.959 R² .. 0.999 R²)
mean                 392.7 μs   (387.4 μs .. 406.6 μs)
std dev              30.48 μs   (10.82 μs .. 55.79 μs)
variance introduced by outliers: 67% (severely inflated)

benchmarking Greatest convex minorant/Vector of length 1e4
time                 3.806 ms   (3.689 ms .. 3.892 ms)
0.959 R²   (0.869 R² .. 0.999 R²)
mean                 4.001 ms   (3.859 ms .. 4.631 ms)
std dev              801.4 μs   (123.5 μs .. 1.807 ms)
variance introduced by outliers: 88% (severely inflated)

benchmarking Greatest convex minorant/Vector of length 1e5
time                 39.91 ms   (38.41 ms .. 41.42 ms)
0.996 R²   (0.993 R² .. 0.998 R²)
mean                 40.01 ms   (39.19 ms .. 41.10 ms)
std dev              1.906 ms   (1.265 ms .. 2.777 ms)
variance introduced by outliers: 13% (moderately inflated)

Benchmark pava-bench: FINISH
Depends on 2 packages(full list with versions):