Bonsai

Statistical benchmarking in J.

Installation

Using pacman

install 'github:jitwit/bonsai@a'

Examples

Monad

Inverting a random 50 by 50 matrix.

   load 'stats/bonsai'
   bonsai '%. ? 50 50 $ 0'
┌────────┬─────────┬──────────┬──────────┐
│N = 2000│lower    │estimate  │upper     │
├────────┼─────────┼──────────┼──────────┤
│μ       │9.3306e_5│9.3794e_5 │9.42814e_5│
├────────┼─────────┼──────────┼──────────┤
│σ       │9.0145e_6│1.09376e_5│1.26696e_5│
└────────┴─────────┴──────────┴──────────┘

Dyad

Comparing two ways to get the 1000 by 1000 identity matrix.

   'e. i. 1000' bonsai '= i. 1000'
┌──────────┬───────────┬───────────┬───────────┐
│comparison│lower      │estimate   │upper      │
├──────────┼───────────┼───────────┼───────────┤
│- & μ     │0.00147287 │0.00148659 │0.00150041 │
├──────────┼───────────┼───────────┼───────────┤
│μ(x)      │0.00165997 │0.00167314 │0.00168693 │
├──────────┼───────────┼───────────┼───────────┤
│μ(y)      │0.000185358│0.000186555│0.000187817│
├──────────┼───────────┼───────────┼───────────┤
│x (-~%[) y│_0.888336  │_0.8885    │_0.888663  │
└──────────┴───────────┴───────────┴───────────┘