lakutata 2.0.81 → 2.0.82

package/vendor/Package.17.cjs

@@ -1,563 +1,563 @@
 "use strict";
 
-function r(r) {
-    var t;
+function linearRegression(r) {
     var e;
+    var t;
     var n = r.length;
     if (n === 1) {
-        t = 0;
-        e = r[0][1];
+        e = 0;
+        t = r[0][1];
     } else {
         var a = 0;
-        var o = 0;
         var i = 0;
+        var o = 0;
         var u = 0;
+        var s;
+        var l;
         var f;
-        var v;
-        var h;
-        for (var s = 0; s < n; s++) {
-            f = r[s];
-            v = f[0];
-            h = f[1];
-            a += v;
-            o += h;
-            i += v * v;
-            u += v * h;
-        }
-        t = (n * u - a * o) / (n * i - a * a);
-        e = o / n - t * a / n;
+        for (var v = 0; v < n; v++) {
+            s = r[v];
+            l = s[0];
+            f = s[1];
+            a += l;
+            i += f;
+            o += l * l;
+            u += l * f;
+        }
+        e = (n * u - a * i) / (n * o - a * a);
+        t = i / n - e * a / n;
     }
     return {
-        m: t,
-        b: e
+        m: e,
+        b: t
     };
 }
 
-function t(r) {
-    return function(t) {
-        return r.b + r.m * t;
+function linearRegressionLine(r) {
+    return function(e) {
+        return r.b + r.m * e;
     };
 }
 
-function e(r) {
+function sum(r) {
     if (r.length === 0) {
         return 0;
     }
-    var t = r[0];
-    var e = 0;
+    var e = r[0];
+    var t = 0;
     var n;
-    if (typeof t !== "number") {
+    if (typeof e !== "number") {
         return Number.NaN;
     }
     for (var a = 1; a < r.length; a++) {
         if (typeof r[a] !== "number") {
             return Number.NaN;
         }
-        n = t + r[a];
-        if (Math.abs(t) >= Math.abs(r[a])) {
-            e += t - n + r[a];
+        n = e + r[a];
+        if (Math.abs(e) >= Math.abs(r[a])) {
+            t += e - n + r[a];
         } else {
-            e += r[a] - n + t;
+            t += r[a] - n + e;
         }
-        t = n;
+        e = n;
     }
-    return t + e;
+    return e + t;
 }
 
-function n(r) {
+function mean(r) {
     if (r.length === 0) {
         throw new Error("mean requires at least one data point");
     }
-    return e(r) / r.length;
+    return sum(r) / r.length;
 }
 
-function a(r, t) {
-    var e = n(r);
-    var a = 0;
-    var o;
+function sumNthPowerDeviations(r, e) {
+    var t = mean(r);
+    var n = 0;
+    var a;
     var i;
-    if (t === 2) {
+    if (e === 2) {
         for (i = 0; i < r.length; i++) {
-            o = r[i] - e;
-            a += o * o;
+            a = r[i] - t;
+            n += a * a;
         }
     } else {
         for (i = 0; i < r.length; i++) {
-            a += Math.pow(r[i] - e, t);
+            n += Math.pow(r[i] - t, e);
         }
     }
-    return a;
+    return n;
 }
 
-function o(r) {
+function variance(r) {
     if (r.length === 0) {
         throw new Error("variance requires at least one data point");
     }
-    return a(r, 2) / r.length;
+    return sumNthPowerDeviations(r, 2) / r.length;
 }
 
-function i(r) {
+function standardDeviation(r) {
     if (r.length === 1) {
         return 0;
     }
-    var t = o(r);
-    return Math.sqrt(t);
+    var e = variance(r);
+    return Math.sqrt(e);
 }
 
-function u(r, t) {
+function rSquared(r, e) {
     if (r.length < 2) {
         return 1;
     }
-    var e = 0;
+    var t = 0;
     for (var n = 0; n < r.length; n++) {
-        e += r[n][1];
+        t += r[n][1];
     }
-    var a = e / r.length;
-    var o = 0;
-    for (var i = 0; i < r.length; i++) {
-        o += Math.pow(a - r[i][1], 2);
+    var a = t / r.length;
+    var i = 0;
+    for (var o = 0; o < r.length; o++) {
+        i += Math.pow(a - r[o][1], 2);
     }
     var u = 0;
-    for (var f = 0; f < r.length; f++) {
-        u += Math.pow(r[f][1] - t(r[f][0]), 2);
+    for (var s = 0; s < r.length; s++) {
+        u += Math.pow(r[s][1] - e(r[s][0]), 2);
     }
-    return 1 - u / o;
+    return 1 - u / i;
 }
 
-function f(r) {
+function modeSorted(r) {
     if (r.length === 0) {
         throw new Error("mode requires at least one data point");
     }
     if (r.length === 1) {
         return r[0];
     }
-    var t = r[0];
-    var e = Number.NaN;
+    var e = r[0];
+    var t = Number.NaN;
     var n = 0;
     var a = 1;
-    for (var o = 1; o < r.length + 1; o++) {
-        if (r[o] !== t) {
+    for (var i = 1; i < r.length + 1; i++) {
+        if (r[i] !== e) {
             if (a > n) {
                 n = a;
-                e = t;
+                t = e;
             }
             a = 1;
-            t = r[o];
+            e = r[i];
         } else {
             a++;
         }
     }
-    return e;
+    return t;
 }
 
-function v(r) {
-    return r.slice().sort((function(r, t) {
-        return r - t;
+function numericSort(r) {
+    return r.slice().sort((function(r, e) {
+        return r - e;
     }));
 }
 
-function h(r) {
-    return f(v(r));
+function mode(r) {
+    return modeSorted(numericSort(r));
 }
 
-function s(r) {
-    var t = new Map;
-    var e;
+function modeFast(r) {
+    var e = new Map;
+    var t;
     var n = 0;
     for (var a = 0; a < r.length; a++) {
-        var o = t.get(r[a]);
-        if (o === undefined) {
-            o = 1;
+        var i = e.get(r[a]);
+        if (i === undefined) {
+            i = 1;
         } else {
-            o++;
+            i++;
         }
-        if (o > n) {
-            e = r[a];
-            n = o;
+        if (i > n) {
+            t = r[a];
+            n = i;
         }
-        t.set(r[a], o);
+        e.set(r[a], i);
     }
     if (n === 0) {
         throw new Error("mode requires at last one data point");
     }
-    return e;
+    return t;
 }
 
-function l(r) {
+function min(r) {
     if (r.length === 0) {
         throw new Error("min requires at least one data point");
     }
-    var t = r[0];
-    for (var e = 1; e < r.length; e++) {
-        if (r[e] < t) {
-            t = r[e];
+    var e = r[0];
+    for (var t = 1; t < r.length; t++) {
+        if (r[t] < e) {
+            e = r[t];
         }
     }
-    return t;
+    return e;
 }
 
-function c(r) {
+function max(r) {
     if (r.length === 0) {
         throw new Error("max requires at least one data point");
     }
-    var t = r[0];
-    for (var e = 1; e < r.length; e++) {
-        if (r[e] > t) {
-            t = r[e];
+    var e = r[0];
+    for (var t = 1; t < r.length; t++) {
+        if (r[t] > e) {
+            e = r[t];
         }
     }
-    return t;
+    return e;
 }
 
-function p(r) {
+function extent(r) {
     if (r.length === 0) {
         throw new Error("extent requires at least one data point");
     }
-    var t = r[0];
     var e = r[0];
+    var t = r[0];
     for (var n = 1; n < r.length; n++) {
-        if (r[n] > e) {
-            e = r[n];
-        }
-        if (r[n] < t) {
+        if (r[n] > t) {
             t = r[n];
         }
+        if (r[n] < e) {
+            e = r[n];
+        }
     }
-    return [ t, e ];
+    return [ e, t ];
 }
 
-function g(r) {
+function minSorted(r) {
     return r[0];
 }
 
-function w(r) {
+function maxSorted(r) {
     return r[r.length - 1];
 }
 
-function M(r) {
+function extentSorted(r) {
     return [ r[0], r[r.length - 1] ];
 }
 
-function m(r) {
-    var t = 0;
-    for (var e = 0; e < r.length; e++) {
-        if (typeof r[e] !== "number") {
+function sumSimple(r) {
+    var e = 0;
+    for (var t = 0; t < r.length; t++) {
+        if (typeof r[t] !== "number") {
             return Number.NaN;
         }
-        t += r[e];
+        e += r[t];
     }
-    return t;
+    return e;
 }
 
-function d(r) {
-    var t = 1;
-    for (var e = 0; e < r.length; e++) {
-        t *= r[e];
+function product(r) {
+    var e = 1;
+    for (var t = 0; t < r.length; t++) {
+        e *= r[t];
     }
-    return t;
+    return e;
 }
 
-function b(r, t) {
-    var e = r.length * t;
+function quantileSorted(r, e) {
+    var t = r.length * e;
     if (r.length === 0) {
         throw new Error("quantile requires at least one data point.");
-    } else if (t < 0 || t > 1) {
+    } else if (e < 0 || e > 1) {
         throw new Error("quantiles must be between 0 and 1");
-    } else if (t === 1) {
+    } else if (e === 1) {
         return r[r.length - 1];
-    } else if (t === 0) {
+    } else if (e === 0) {
         return r[0];
-    } else if (e % 1 !== 0) {
-        return r[Math.ceil(e) - 1];
+    } else if (t % 1 !== 0) {
+        return r[Math.ceil(t) - 1];
     } else if (r.length % 2 === 0) {
-        return (r[e - 1] + r[e]) / 2;
+        return (r[t - 1] + r[t]) / 2;
     } else {
-        return r[e];
+        return r[t];
     }
 }
 
-function x(r, t, e, n) {
-    e = e || 0;
+function quickselect(r, e, t, n) {
+    t = t || 0;
     n = n || r.length - 1;
-    while (n > e) {
-        if (n - e > 600) {
-            var a = n - e + 1;
-            var o = t - e + 1;
-            var i = Math.log(a);
-            var u = .5 * Math.exp(2 * i / 3);
-            var f = .5 * Math.sqrt(i * u * (a - u) / a);
-            if (o - a / 2 < 0) {
-                f *= -1;
+    while (n > t) {
+        if (n - t > 600) {
+            var a = n - t + 1;
+            var i = e - t + 1;
+            var o = Math.log(a);
+            var u = .5 * Math.exp(2 * o / 3);
+            var s = .5 * Math.sqrt(o * u * (a - u) / a);
+            if (i - a / 2 < 0) {
+                s *= -1;
             }
-            var v = Math.max(e, Math.floor(t - o * u / a + f));
-            var h = Math.min(n, Math.floor(t + (a - o) * u / a + f));
-            x(r, t, v, h);
+            var l = Math.max(t, Math.floor(e - i * u / a + s));
+            var f = Math.min(n, Math.floor(e + (a - i) * u / a + s));
+            quickselect(r, e, l, f);
         }
-        var s = r[t];
-        var l = e;
+        var v = r[e];
+        var h = t;
         var c = n;
-        q(r, e, t);
-        if (r[n] > s) {
-            q(r, e, n);
+        swap(r, t, e);
+        if (r[n] > v) {
+            swap(r, t, n);
         }
-        while (l < c) {
-            q(r, l, c);
-            l++;
+        while (h < c) {
+            swap(r, h, c);
+            h++;
             c--;
-            while (r[l] < s) {
-                l++;
+            while (r[h] < v) {
+                h++;
             }
-            while (r[c] > s) {
+            while (r[c] > v) {
                 c--;
             }
         }
-        if (r[e] === s) {
-            q(r, e, c);
+        if (r[t] === v) {
+            swap(r, t, c);
         } else {
             c++;
-            q(r, c, n);
+            swap(r, c, n);
         }
-        if (c <= t) {
-            e = c + 1;
+        if (c <= e) {
+            t = c + 1;
         }
-        if (t <= c) {
+        if (e <= c) {
             n = c - 1;
         }
     }
 }
 
-function q(r, t, e) {
-    var n = r[t];
-    r[t] = r[e];
-    r[e] = n;
+function swap(r, e, t) {
+    var n = r[e];
+    r[e] = r[t];
+    r[t] = n;
 }
 
-function E(r, t) {
-    var e = r.slice();
-    if (Array.isArray(t)) {
-        N(e, t);
+function quantile(r, e) {
+    var t = r.slice();
+    if (Array.isArray(e)) {
+        multiQuantileSelect(t, e);
         var n = [];
-        for (var a = 0; a < t.length; a++) {
-            n[a] = b(e, t[a]);
+        for (var a = 0; a < e.length; a++) {
+            n[a] = quantileSorted(t, e[a]);
         }
         return n;
     } else {
-        var o = k(e.length, t);
-        y(e, o, 0, e.length - 1);
-        return b(e, t);
+        var i = quantileIndex(t.length, e);
+        quantileSelect(t, i, 0, t.length - 1);
+        return quantileSorted(t, e);
     }
 }
 
-function y(r, t, e, n) {
-    if (t % 1 === 0) {
-        x(r, t, e, n);
+function quantileSelect(r, e, t, n) {
+    if (e % 1 === 0) {
+        quickselect(r, e, t, n);
     } else {
-        t = Math.floor(t);
-        x(r, t, e, n);
-        x(r, t + 1, t + 1, n);
+        e = Math.floor(e);
+        quickselect(r, e, t, n);
+        quickselect(r, e + 1, e + 1, n);
     }
 }
 
-function N(r, t) {
-    var e = [ 0 ];
-    for (var n = 0; n < t.length; n++) {
-        e.push(k(r.length, t[n]));
+function multiQuantileSelect(r, e) {
+    var t = [ 0 ];
+    for (var n = 0; n < e.length; n++) {
+        t.push(quantileIndex(r.length, e[n]));
     }
-    e.push(r.length - 1);
-    e.sort(I);
-    var a = [ 0, e.length - 1 ];
+    t.push(r.length - 1);
+    t.sort(compare);
+    var a = [ 0, t.length - 1 ];
     while (a.length) {
-        var o = Math.ceil(a.pop());
-        var i = Math.floor(a.pop());
-        if (o - i <= 1) {
+        var i = Math.ceil(a.pop());
+        var o = Math.floor(a.pop());
+        if (i - o <= 1) {
             continue;
         }
-        var u = Math.floor((i + o) / 2);
-        y(r, e[u], Math.floor(e[i]), Math.ceil(e[o]));
-        a.push(i, u, u, o);
+        var u = Math.floor((o + i) / 2);
+        quantileSelect(r, t[u], Math.floor(t[o]), Math.ceil(t[i]));
+        a.push(o, u, u, i);
     }
 }
 
-function I(r, t) {
-    return r - t;
+function compare(r, e) {
+    return r - e;
 }
 
-function k(r, t) {
-    var e = r * t;
-    if (t === 1) {
+function quantileIndex(r, e) {
+    var t = r * e;
+    if (e === 1) {
         return r - 1;
-    } else if (t === 0) {
+    } else if (e === 0) {
         return 0;
-    } else if (e % 1 !== 0) {
-        return Math.ceil(e) - 1;
+    } else if (t % 1 !== 0) {
+        return Math.ceil(t) - 1;
     } else if (r % 2 === 0) {
-        return e - .5;
+        return t - .5;
     } else {
-        return e;
+        return t;
     }
 }
 
-function A(r, t) {
-    if (t < r[0]) {
+function quantileRankSorted(r, e) {
+    if (e < r[0]) {
         return 0;
     }
-    if (t > r[r.length - 1]) {
+    if (e > r[r.length - 1]) {
         return 1;
     }
-    var e = P(r, t);
-    if (r[e] !== t) {
-        return e / r.length;
+    var t = lowerBound(r, e);
+    if (r[t] !== e) {
+        return t / r.length;
     }
-    e++;
-    var n = S(r, t);
-    if (n === e) {
-        return e / r.length;
+    t++;
+    var n = upperBound(r, e);
+    if (n === t) {
+        return t / r.length;
     }
-    var a = n - e + 1;
-    var o = a * (n + e) / 2;
-    var i = o / a;
-    return i / r.length;
+    var a = n - t + 1;
+    var i = a * (n + t) / 2;
+    var o = i / a;
+    return o / r.length;
 }
 
-function P(r, t) {
-    var e = 0;
+function lowerBound(r, e) {
+    var t = 0;
     var n = 0;
     var a = r.length;
     while (n < a) {
-        e = n + a >>> 1;
-        if (t <= r[e]) {
-            a = e;
+        t = n + a >>> 1;
+        if (e <= r[t]) {
+            a = t;
         } else {
-            n = -~e;
+            n = -~t;
         }
     }
     return n;
 }
 
-function S(r, t) {
-    var e = 0;
+function upperBound(r, e) {
+    var t = 0;
     var n = 0;
     var a = r.length;
     while (n < a) {
-        e = n + a >>> 1;
-        if (t >= r[e]) {
-            n = -~e;
+        t = n + a >>> 1;
+        if (e >= r[t]) {
+            n = -~t;
         } else {
-            a = e;
+            a = t;
         }
     }
     return n;
 }
 
-function C(r, t) {
-    var e = v(r);
-    return A(e, t);
+function quantileRank(r, e) {
+    var t = numericSort(r);
+    return quantileRankSorted(t, e);
 }
 
-function _(r) {
-    var t = E(r, .75);
-    var e = E(r, .25);
-    if (typeof t === "number" && typeof e === "number") {
-        return t - e;
+function interquartileRange(r) {
+    var e = quantile(r, .75);
+    var t = quantile(r, .25);
+    if (typeof e === "number" && typeof t === "number") {
+        return e - t;
     }
 }
 
-function T(r) {
-    return +E(r, .5);
+function median(r) {
+    return +quantile(r, .5);
 }
 
-function V(r) {
-    var t = T(r);
-    var e = [];
+function medianAbsoluteDeviation(r) {
+    var e = median(r);
+    var t = [];
     for (var n = 0; n < r.length; n++) {
-        e.push(Math.abs(r[n] - t));
+        t.push(Math.abs(r[n] - e));
     }
-    return T(e);
+    return median(t);
 }
 
-function D(r, t) {
-    var e = [];
-    if (t < 1) {
+function chunk(r, e) {
+    var t = [];
+    if (e < 1) {
         throw new Error("chunk size must be a positive number");
     }
-    if (Math.floor(t) !== t) {
+    if (Math.floor(e) !== e) {
         throw new Error("chunk size must be an integer");
     }
-    for (var n = 0; n < r.length; n += t) {
-        e.push(r.slice(n, n + t));
+    for (var n = 0; n < r.length; n += e) {
+        t.push(r.slice(n, n + e));
     }
-    return e;
+    return t;
 }
 
-function L(r, t, e) {
+function sampleWithReplacement(r, e, t) {
     if (r.length === 0) {
         return [];
     }
-    e = e || Math.random;
+    t = t || Math.random;
     var n = r.length;
     var a = [];
-    for (var o = 0; o < t; o++) {
-        var i = Math.floor(e() * n);
-        a.push(r[i]);
+    for (var i = 0; i < e; i++) {
+        var o = Math.floor(t() * n);
+        a.push(r[o]);
     }
     return a;
 }
 
-function R(r, t) {
-    t = t || Math.random;
-    var e = r.length;
+function shuffleInPlace(r, e) {
+    e = e || Math.random;
+    var t = r.length;
     var n;
     var a;
-    while (e > 0) {
-        a = Math.floor(t() * e--);
-        n = r[e];
-        r[e] = r[a];
+    while (t > 0) {
+        a = Math.floor(e() * t--);
+        n = r[t];
+        r[t] = r[a];
         r[a] = n;
     }
     return r;
 }
 
-function U(r, t) {
-    var e = r.slice();
-    return R(e, t);
+function shuffle(r, e) {
+    var t = r.slice();
+    return shuffleInPlace(t, e);
 }
 
-function z(r, t, e) {
-    var n = U(r, e);
-    return n.slice(0, t);
+function sample(r, e, t) {
+    var n = shuffle(r, t);
+    return n.slice(0, e);
 }
 
-function F(r, t) {
-    var e = [];
+function makeMatrix(r, e) {
+    var t = [];
     for (var n = 0; n < r; n++) {
         var a = [];
-        for (var o = 0; o < t; o++) {
+        for (var i = 0; i < e; i++) {
             a.push(0);
         }
-        e.push(a);
+        t.push(a);
     }
-    return e;
+    return t;
 }
 
-function O(r) {
-    var t = 0;
-    var e;
+function uniqueCountSorted(r) {
+    var e = 0;
+    var t;
     for (var n = 0; n < r.length; n++) {
-        if (n === 0 || r[n] !== e) {
-            e = r[n];
-            t++;
+        if (n === 0 || r[n] !== t) {
+            t = r[n];
+            e++;
         }
     }
-    return t;
+    return e;
 }
 
-function X(r, t, e, n) {
+function ssq(r, e, t, n) {
     var a;
     if (r > 0) {
-        var o = (e[t] - e[r - 1]) / (t - r + 1);
-        a = n[t] - n[r - 1] - (t - r + 1) * o * o;
+        var i = (t[e] - t[r - 1]) / (e - r + 1);
+        a = n[e] - n[r - 1] - (e - r + 1) * i * i;
     } else {
-        a = n[t] - e[t] * e[t] / (t + 1);
+        a = n[e] - t[e] * t[e] / (e + 1);
     }
     if (a < 0) {
         return 0;
@@ -565,573 +565,573 @@ function X(r, t, e, n) {
     return a;
 }
 
-function Y(r, t, e, n, a, o, i) {
-    if (r > t) {
+function fillMatrixColumn(r, e, t, n, a, i, o) {
+    if (r > e) {
         return;
     }
-    var u = Math.floor((r + t) / 2);
-    n[e][u] = n[e - 1][u - 1];
-    a[e][u] = u;
-    var f = e;
-    if (r > e) {
-        f = Math.max(f, a[e][r - 1] || 0);
+    var u = Math.floor((r + e) / 2);
+    n[t][u] = n[t - 1][u - 1];
+    a[t][u] = u;
+    var s = t;
+    if (r > t) {
+        s = Math.max(s, a[t][r - 1] || 0);
     }
-    f = Math.max(f, a[e - 1][u] || 0);
-    var v = u - 1;
-    if (t < n[0].length - 1) {
-        v = Math.min(v, a[e][t + 1] || 0);
+    s = Math.max(s, a[t - 1][u] || 0);
+    var l = u - 1;
+    if (e < n[0].length - 1) {
+        l = Math.min(l, a[t][e + 1] || 0);
     }
+    var f;
+    var v;
     var h;
-    var s;
-    var l;
     var c;
-    for (var p = v; p >= f; --p) {
-        h = X(p, u, o, i);
-        if (h + n[e - 1][f - 1] >= n[e][u]) {
+    for (var p = l; p >= s; --p) {
+        f = ssq(p, u, i, o);
+        if (f + n[t - 1][s - 1] >= n[t][u]) {
             break;
         }
-        s = X(f, u, o, i);
-        l = s + n[e - 1][f - 1];
-        if (l < n[e][u]) {
-            n[e][u] = l;
-            a[e][u] = f;
+        v = ssq(s, u, i, o);
+        h = v + n[t - 1][s - 1];
+        if (h < n[t][u]) {
+            n[t][u] = h;
+            a[t][u] = s;
         }
-        f++;
-        c = h + n[e - 1][p - 1];
-        if (c < n[e][u]) {
-            n[e][u] = c;
-            a[e][u] = p;
+        s++;
+        c = f + n[t - 1][p - 1];
+        if (c < n[t][u]) {
+            n[t][u] = c;
+            a[t][u] = p;
         }
     }
-    Y(r, u - 1, e, n, a, o, i);
-    Y(u + 1, t, e, n, a, o, i);
+    fillMatrixColumn(r, u - 1, t, n, a, i, o);
+    fillMatrixColumn(u + 1, e, t, n, a, i, o);
 }
 
-function j(r, t, e) {
-    var n = t[0].length;
+function fillMatrices(r, e, t) {
+    var n = e[0].length;
     var a = r[Math.floor(n / 2)];
-    var o = [];
     var i = [];
-    for (var u = 0, f = void 0; u < n; ++u) {
-        f = r[u] - a;
+    var o = [];
+    for (var u = 0, s = void 0; u < n; ++u) {
+        s = r[u] - a;
         if (u === 0) {
-            o.push(f);
-            i.push(f * f);
+            i.push(s);
+            o.push(s * s);
         } else {
-            o.push(o[u - 1] + f);
-            i.push(i[u - 1] + f * f);
+            i.push(i[u - 1] + s);
+            o.push(o[u - 1] + s * s);
         }
-        t[0][u] = X(0, u, o, i);
-        e[0][u] = 0;
+        e[0][u] = ssq(0, u, i, o);
+        t[0][u] = 0;
     }
-    var v;
-    for (var h = 1; h < t.length; ++h) {
-        if (h < t.length - 1) {
-            v = h;
+    var l;
+    for (var f = 1; f < e.length; ++f) {
+        if (f < e.length - 1) {
+            l = f;
         } else {
-            v = n - 1;
+            l = n - 1;
         }
-        Y(v, n - 1, h, t, e, o, i);
+        fillMatrixColumn(l, n - 1, f, e, t, i, o);
     }
 }
 
-function K(r, t) {
-    if (t > r.length) {
+function ckmeans(r, e) {
+    if (e > r.length) {
         throw new Error("cannot generate more classes than there are data values");
     }
-    var e = v(r);
-    var n = O(e);
+    var t = numericSort(r);
+    var n = uniqueCountSorted(t);
     if (n === 1) {
-        return [ e ];
+        return [ t ];
     }
-    var a = F(t, e.length);
-    var o = F(t, e.length);
-    j(e, a, o);
-    var i = [];
-    var u = o[0].length - 1;
-    for (var f = o.length - 1; f >= 0; f--) {
-        var h = o[f][u];
-        i[f] = e.slice(h, u + 1);
-        if (f > 0) {
-            u = h - 1;
+    var a = makeMatrix(e, t.length);
+    var i = makeMatrix(e, t.length);
+    fillMatrices(t, a, i);
+    var o = [];
+    var u = i[0].length - 1;
+    for (var s = i.length - 1; s >= 0; s--) {
+        var l = i[s][u];
+        o[s] = t.slice(l, u + 1);
+        if (s > 0) {
+            u = l - 1;
         }
     }
-    return i;
+    return o;
 }
 
-function W(r, t, e) {
+function jenksBreaks(r, e, t) {
     var n = r.length;
     var a = [];
-    var o = e;
-    a[e] = r[r.length - 1];
-    while (o > 0) {
-        a[o - 1] = r[t[n][o] - 1];
-        n = t[n][o] - 1;
-        o--;
+    var i = t;
+    a[t] = r[r.length - 1];
+    while (i > 0) {
+        a[i - 1] = r[e[n][i] - 1];
+        n = e[n][i] - 1;
+        i--;
     }
     return a;
 }
 
-function B(r, t) {
-    var e = [];
+function jenksMatrices(r, e) {
+    var t = [];
     var n = [];
     var a;
-    var o;
-    var i = 0;
+    var i;
+    var o = 0;
     for (a = 0; a < r.length + 1; a++) {
         var u = [];
-        var f = [];
-        for (o = 0; o < t + 1; o++) {
+        var s = [];
+        for (i = 0; i < e + 1; i++) {
             u.push(0);
-            f.push(0);
+            s.push(0);
         }
-        e.push(u);
-        n.push(f);
+        t.push(u);
+        n.push(s);
     }
-    for (a = 1; a < t + 1; a++) {
-        e[1][a] = 1;
+    for (a = 1; a < e + 1; a++) {
+        t[1][a] = 1;
         n[1][a] = 0;
-        for (o = 2; o < r.length + 1; o++) {
-            n[o][a] = Number.POSITIVE_INFINITY;
+        for (i = 2; i < r.length + 1; i++) {
+            n[i][a] = Number.POSITIVE_INFINITY;
         }
     }
-    for (var v = 2; v < r.length + 1; v++) {
+    for (var l = 2; l < r.length + 1; l++) {
+        var f = 0;
+        var v = 0;
         var h = 0;
-        var s = 0;
-        var l = 0;
         var c = 0;
-        for (var p = 1; p < v + 1; p++) {
-            var g = v - p + 1;
-            var w = r[g - 1];
-            l++;
-            h += w;
-            s += w * w;
-            i = s - h * h / l;
-            c = g - 1;
+        for (var p = 1; p < l + 1; p++) {
+            var m = l - p + 1;
+            var g = r[m - 1];
+            h++;
+            f += g;
+            v += g * g;
+            o = v - f * f / h;
+            c = m - 1;
             if (c !== 0) {
-                for (o = 2; o < t + 1; o++) {
-                    if (n[v][o] >= i + n[c][o - 1]) {
-                        e[v][o] = g;
-                        n[v][o] = i + n[c][o - 1];
+                for (i = 2; i < e + 1; i++) {
+                    if (n[l][i] >= o + n[c][i - 1]) {
+                        t[l][i] = m;
+                        n[l][i] = o + n[c][i - 1];
                     }
                 }
             }
         }
-        e[v][1] = 1;
-        n[v][1] = i;
+        t[l][1] = 1;
+        n[l][1] = o;
     }
     return {
-        lowerClassLimits: e,
+        lowerClassLimits: t,
         varianceCombinations: n
     };
 }
 
-function G(r, t) {
-    if (t > r.length) {
+function jenks(r, e) {
+    if (e > r.length) {
         return null;
     }
-    r = r.slice().sort((function(r, t) {
-        return r - t;
+    r = r.slice().sort((function(r, e) {
+        return r - e;
     }));
-    var e = B(r, t);
-    var n = e.lowerClassLimits;
-    return W(r, n, t);
+    var t = jenksMatrices(r, e);
+    var n = t.lowerClassLimits;
+    return jenksBreaks(r, n, e);
 }
 
-function H(r, t) {
+function equalIntervalBreaks(r, e) {
     if (r.length < 2) {
         return r;
     }
-    var e = l(r);
-    var n = c(r);
-    var a = [ e ];
-    var o = (n - e) / t;
-    for (var i = 1; i < t; i++) {
-        a.push(a[0] + o * i);
+    var t = min(r);
+    var n = max(r);
+    var a = [ t ];
+    var i = (n - t) / e;
+    for (var o = 1; o < e; o++) {
+        a.push(a[0] + i * o);
     }
     a.push(n);
     return a;
 }
 
-function J(r, t) {
-    if (r.length !== t.length) {
+function sampleCovariance(r, e) {
+    if (r.length !== e.length) {
         throw new Error("sampleCovariance requires samples with equal lengths");
     }
     if (r.length < 2) {
         throw new Error("sampleCovariance requires at least two data points in each sample");
     }
-    var e = n(r);
-    var a = n(t);
-    var o = 0;
+    var t = mean(r);
+    var n = mean(e);
+    var a = 0;
     for (var i = 0; i < r.length; i++) {
-        o += (r[i] - e) * (t[i] - a);
+        a += (r[i] - t) * (e[i] - n);
     }
-    var u = r.length - 1;
-    return o / u;
+    var o = r.length - 1;
+    return a / o;
 }
 
-function Q(r) {
+function sampleVariance(r) {
     if (r.length < 2) {
         throw new Error("sampleVariance requires at least two data points");
     }
-    var t = a(r, 2);
-    var e = r.length - 1;
-    return t / e;
+    var e = sumNthPowerDeviations(r, 2);
+    var t = r.length - 1;
+    return e / t;
 }
 
-function Z(r) {
-    var t = Q(r);
-    return Math.sqrt(t);
+function sampleStandardDeviation(r) {
+    var e = sampleVariance(r);
+    return Math.sqrt(e);
 }
 
-function $(r, t) {
-    var e = J(r, t);
-    var n = Z(r);
-    var a = Z(t);
-    return e / n / a;
+function sampleCorrelation(r, e) {
+    var t = sampleCovariance(r, e);
+    var n = sampleStandardDeviation(r);
+    var a = sampleStandardDeviation(e);
+    return t / n / a;
 }
 
-function rr(r, t) {
-    var e = r.map((function(r, t) {
-        return [ r, t ];
-    })).sort((function(r, t) {
-        return r[0] - t[0];
+function sampleRankCorrelation(r, e) {
+    var t = r.map((function(r, e) {
+        return [ r, e ];
+    })).sort((function(r, e) {
+        return r[0] - e[0];
     })).map((function(r) {
         return r[1];
     }));
-    var n = t.map((function(r, t) {
-        return [ r, t ];
-    })).sort((function(r, t) {
-        return r[0] - t[0];
+    var n = e.map((function(r, e) {
+        return [ r, e ];
+    })).sort((function(r, e) {
+        return r[0] - e[0];
     })).map((function(r) {
         return r[1];
     }));
-    var a = Array(e.length);
-    var o = Array(e.length);
-    for (var i = 0; i < e.length; i++) {
-        a[e[i]] = i;
-        o[n[i]] = i;
+    var a = Array(t.length);
+    var i = Array(t.length);
+    for (var o = 0; o < t.length; o++) {
+        a[t[o]] = o;
+        i[n[o]] = o;
     }
-    return $(a, o);
+    return sampleCorrelation(a, i);
 }
 
-function tr(r) {
+function sampleSkewness(r) {
     if (r.length < 3) {
         throw new Error("sampleSkewness requires at least three data points");
     }
-    var t = n(r);
-    var e;
+    var e = mean(r);
+    var t;
+    var n = 0;
     var a = 0;
-    var o = 0;
     for (var i = 0; i < r.length; i++) {
-        e = r[i] - t;
-        a += e * e;
-        o += e * e * e;
+        t = r[i] - e;
+        n += t * t;
+        a += t * t * t;
     }
-    var u = r.length - 1;
-    var f = Math.sqrt(a / u);
-    var v = r.length;
-    var h = Math.pow(f, 3);
-    return v * o / ((v - 1) * (v - 2) * h);
+    var o = r.length - 1;
+    var u = Math.sqrt(n / o);
+    var s = r.length;
+    var l = Math.pow(u, 3);
+    return s * a / ((s - 1) * (s - 2) * l);
 }
 
-function er(r) {
-    var t = r.length;
-    if (t < 4) {
+function sampleKurtosis(r) {
+    var e = r.length;
+    if (e < 4) {
         throw new Error("sampleKurtosis requires at least four data points");
     }
-    var e = n(r);
-    var a;
-    var o = 0;
+    var t = mean(r);
+    var n;
+    var a = 0;
     var i = 0;
-    for (var u = 0; u < t; u++) {
-        a = r[u] - e;
-        o += a * a;
-        i += a * a * a * a;
+    for (var o = 0; o < e; o++) {
+        n = r[o] - t;
+        a += n * n;
+        i += n * n * n * n;
     }
-    return (t - 1) / ((t - 2) * (t - 3)) * (t * (t + 1) * i / (o * o) - 3 * (t - 1));
+    return (e - 1) / ((e - 2) * (e - 3)) * (e * (e + 1) * i / (a * a) - 3 * (e - 1));
 }
 
-function nr(r) {
-    var t = new Array(r.length);
-    var e = [ r.slice() ];
+function permutationsHeap(r) {
+    var e = new Array(r.length);
+    var t = [ r.slice() ];
     for (var n = 0; n < r.length; n++) {
-        t[n] = 0;
+        e[n] = 0;
     }
     for (var a = 0; a < r.length; ) {
-        if (t[a] < a) {
-            var o = 0;
+        if (e[a] < a) {
+            var i = 0;
             if (a % 2 !== 0) {
-                o = t[a];
+                i = e[a];
             }
-            var i = r[o];
-            r[o] = r[a];
-            r[a] = i;
-            e.push(r.slice());
-            t[a]++;
+            var o = r[i];
+            r[i] = r[a];
+            r[a] = o;
+            t.push(r.slice());
+            e[a]++;
             a = 0;
         } else {
-            t[a] = 0;
+            e[a] = 0;
             a++;
         }
     }
-    return e;
+    return t;
 }
 
-function ar(r, t) {
-    var e;
+function combinations(r, e) {
+    var t;
     var n;
     var a = [];
-    var o;
     var i;
-    for (e = 0; e < r.length; e++) {
-        if (t === 1) {
-            a.push([ r[e] ]);
+    var o;
+    for (t = 0; t < r.length; t++) {
+        if (e === 1) {
+            a.push([ r[t] ]);
         } else {
-            o = ar(r.slice(e + 1, r.length), t - 1);
-            for (n = 0; n < o.length; n++) {
-                i = o[n];
-                i.unshift(r[e]);
-                a.push(i);
+            i = combinations(r.slice(t + 1, r.length), e - 1);
+            for (n = 0; n < i.length; n++) {
+                o = i[n];
+                o.unshift(r[t]);
+                a.push(o);
             }
         }
     }
     return a;
 }
 
-function or(r, t) {
-    var e = [];
+function combinationsReplacement(r, e) {
+    var t = [];
     for (var n = 0; n < r.length; n++) {
-        if (t === 1) {
-            e.push([ r[n] ]);
+        if (e === 1) {
+            t.push([ r[n] ]);
         } else {
-            var a = or(r.slice(n, r.length), t - 1);
-            for (var o = 0; o < a.length; o++) {
-                e.push([ r[n] ].concat(a[o]));
+            var a = combinationsReplacement(r.slice(n, r.length), e - 1);
+            for (var i = 0; i < a.length; i++) {
+                t.push([ r[n] ].concat(a[i]));
             }
         }
     }
-    return e;
+    return t;
 }
 
-function ir(r, t, e) {
-    return r + (e - r) / (t + 1);
+function addToMean(r, e, t) {
+    return r + (t - r) / (e + 1);
 }
 
-function ur(r, t, e, n) {
-    return (r * t + e * n) / (t + n);
+function combineMeans(r, e, t, n) {
+    return (r * e + t * n) / (e + n);
 }
 
-function fr(r, t, e, n, a, o) {
-    var i = ur(t, e, a, o);
-    return (e * (r + Math.pow(t - i, 2)) + o * (n + Math.pow(a - i, 2))) / (e + o);
+function combineVariances(r, e, t, n, a, i) {
+    var o = combineMeans(e, t, a, i);
+    return (t * (r + Math.pow(e - o, 2)) + i * (n + Math.pow(a - o, 2))) / (t + i);
 }
 
-function vr(r) {
+function geometricMean(r) {
     if (r.length === 0) {
         throw new Error("geometricMean requires at least one data point");
     }
-    var t = 1;
-    for (var e = 0; e < r.length; e++) {
-        if (r[e] < 0) {
+    var e = 1;
+    for (var t = 0; t < r.length; t++) {
+        if (r[t] < 0) {
             throw new Error("geometricMean requires only non-negative numbers as input");
         }
-        t *= r[e];
+        e *= r[t];
     }
-    return Math.pow(t, 1 / r.length);
+    return Math.pow(e, 1 / r.length);
 }
 
-function hr(r) {
+function logAverage(r) {
     if (r.length === 0) {
         throw new Error("logAverage requires at least one data point");
     }
-    var t = 0;
-    for (var e = 0; e < r.length; e++) {
-        if (r[e] < 0) {
+    var e = 0;
+    for (var t = 0; t < r.length; t++) {
+        if (r[t] < 0) {
             throw new Error("logAverage requires only non-negative numbers as input");
         }
-        t += Math.log(r[e]);
+        e += Math.log(r[t]);
     }
-    return Math.exp(t / r.length);
+    return Math.exp(e / r.length);
 }
 
-function sr(r) {
+function harmonicMean(r) {
     if (r.length === 0) {
         throw new Error("harmonicMean requires at least one data point");
     }
-    var t = 0;
-    for (var e = 0; e < r.length; e++) {
-        if (r[e] <= 0) {
+    var e = 0;
+    for (var t = 0; t < r.length; t++) {
+        if (r[t] <= 0) {
             throw new Error("harmonicMean requires only positive numbers as input");
         }
-        t += 1 / r[e];
+        e += 1 / r[t];
     }
-    return r.length / t;
+    return r.length / e;
 }
 
-function lr(r) {
+function meanSimple(r) {
     if (r.length === 0) {
         throw new Error("meanSimple requires at least one data point");
     }
-    return m(r) / r.length;
+    return sumSimple(r) / r.length;
 }
 
-function cr(r) {
-    return b(r, .5);
+function medianSorted(r) {
+    return quantileSorted(r, .5);
 }
 
-function pr(r, t, e) {
-    return (r * t - e) / (t - 1);
+function subtractFromMean(r, e, t) {
+    return (r * e - t) / (e - 1);
 }
 
-function gr(r) {
+function rootMeanSquare(r) {
     if (r.length === 0) {
         throw new Error("rootMeanSquare requires at least one data point");
     }
-    var t = 0;
-    for (var e = 0; e < r.length; e++) {
-        t += Math.pow(r[e], 2);
+    var e = 0;
+    for (var t = 0; t < r.length; t++) {
+        e += Math.pow(r[t], 2);
     }
-    return Math.sqrt(t / r.length);
+    return Math.sqrt(e / r.length);
 }
 
-function wr(r) {
-    return Z(r) / n(r);
+function coefficientOfVariation(r) {
+    return sampleStandardDeviation(r) / mean(r);
 }
 
-function Mr(r, t) {
-    var e = n(r);
-    var a = i(r);
-    var o = Math.sqrt(r.length);
-    return (e - t) / (a / o);
+function tTest(r, e) {
+    var t = mean(r);
+    var n = standardDeviation(r);
+    var a = Math.sqrt(r.length);
+    return (t - e) / (n / a);
 }
 
-function mr(r, t, e) {
-    var a = r.length;
-    var o = t.length;
-    if (!a || !o) {
+function tTestTwoSample(r, e, t) {
+    var n = r.length;
+    var a = e.length;
+    if (!n || !a) {
         return null;
     }
-    if (!e) {
-        e = 0;
+    if (!t) {
+        t = 0;
     }
-    var i = n(r);
-    var u = n(t);
-    var f = Q(r);
-    var v = Q(t);
-    if (typeof i === "number" && typeof u === "number" && typeof f === "number" && typeof v === "number") {
-        var h = ((a - 1) * f + (o - 1) * v) / (a + o - 2);
-        return (i - u - e) / Math.sqrt(h * (1 / a + 1 / o));
+    var i = mean(r);
+    var o = mean(e);
+    var u = sampleVariance(r);
+    var s = sampleVariance(e);
+    if (typeof i === "number" && typeof o === "number" && typeof u === "number" && typeof s === "number") {
+        var l = ((n - 1) * u + (a - 1) * s) / (n + a - 2);
+        return (i - o - t) / Math.sqrt(l * (1 / n + 1 / a));
     }
 }
 
-function dr(r, t) {
-    if (!r.length || !t.length) {
+function wilcoxonRankSum(r, e) {
+    if (!r.length || !e.length) {
         throw new Error("Neither sample can be empty");
     }
-    var e = r.map((function(r) {
+    var t = r.map((function(r) {
         return {
             label: "x",
             value: r
         };
-    })).concat(t.map((function(r) {
+    })).concat(e.map((function(r) {
         return {
             label: "y",
             value: r
         };
-    }))).sort((function(r, t) {
-        return r.value - t.value;
+    }))).sort((function(r, e) {
+        return r.value - e.value;
     }));
-    for (var n = 0; n < e.length; n++) {
-        e[n].rank = n;
-    }
-    var a = [ e[0].rank ];
-    for (var o = 1; o < e.length; o++) {
-        if (e[o].value === e[o - 1].value) {
-            a.push(e[o].rank);
-            if (o === e.length - 1) {
-                i(e, a);
+    for (var n = 0; n < t.length; n++) {
+        t[n].rank = n;
+    }
+    var a = [ t[0].rank ];
+    for (var i = 1; i < t.length; i++) {
+        if (t[i].value === t[i - 1].value) {
+            a.push(t[i].rank);
+            if (i === t.length - 1) {
+                replaceRanksInPlace(t, a);
             }
         } else if (a.length > 1) {
-            i(e, a);
+            replaceRanksInPlace(t, a);
         } else {
-            a = [ e[o].rank ];
+            a = [ t[i].rank ];
         }
     }
-    function i(r, t) {
-        var e = (t[0] + t[t.length - 1]) / 2;
-        for (var n = 0; n < t.length; n++) {
-            r[t[n]].rank = e;
+    function replaceRanksInPlace(r, e) {
+        var t = (e[0] + e[e.length - 1]) / 2;
+        for (var n = 0; n < e.length; n++) {
+            r[e[n]].rank = t;
         }
     }
-    var u = 0;
-    for (var f = 0; f < e.length; f++) {
-        var v = e[f];
-        if (v.label === "x") {
-            u += v.rank + 1;
+    var o = 0;
+    for (var u = 0; u < t.length; u++) {
+        var s = t[u];
+        if (s.label === "x") {
+            o += s.rank + 1;
         }
     }
-    return u;
+    return o;
 }
 
-var br = function r() {
+var r = function BayesianClassifier() {
     this.totalCount = 0;
     this.data = {};
 };
 
-br.prototype.train = function r(t, e) {
+r.prototype.train = function train(r, e) {
     if (!this.data[e]) {
         this.data[e] = {};
     }
-    for (var n in t) {
-        var a = t[n];
-        if (this.data[e][n] === undefined) {
-            this.data[e][n] = {};
+    for (var t in r) {
+        var n = r[t];
+        if (this.data[e][t] === undefined) {
+            this.data[e][t] = {};
         }
-        if (this.data[e][n][a] === undefined) {
-            this.data[e][n][a] = 0;
+        if (this.data[e][t][n] === undefined) {
+            this.data[e][t][n] = 0;
         }
-        this.data[e][n][a]++;
+        this.data[e][t][n]++;
     }
     this.totalCount++;
 };
 
-br.prototype.score = function r(t) {
+r.prototype.score = function score(r) {
     var e = {};
-    var n;
-    for (var a in t) {
-        var o = t[a];
-        for (n in this.data) {
-            e[n] = {};
-            if (this.data[n][a]) {
-                e[n][a + "_" + o] = (this.data[n][a][o] || 0) / this.totalCount;
+    var t;
+    for (var n in r) {
+        var a = r[n];
+        for (t in this.data) {
+            e[t] = {};
+            if (this.data[t][n]) {
+                e[t][n + "_" + a] = (this.data[t][n][a] || 0) / this.totalCount;
             } else {
-                e[n][a + "_" + o] = 0;
+                e[t][n + "_" + a] = 0;
             }
         }
     }
     var i = {};
-    for (n in e) {
-        i[n] = 0;
-        for (var u in e[n]) {
-            i[n] += e[n][u];
+    for (t in e) {
+        i[t] = 0;
+        for (var o in e[t]) {
+            i[t] += e[t][o];
         }
     }
     return i;
 };
 
-var xr = function r() {
+var e = function PerceptronModel() {
     this.weights = [];
     this.bias = 0;
 };
 
-xr.prototype.predict = function r(t) {
-    if (t.length !== this.weights.length) {
+e.prototype.predict = function predict(r) {
+    if (r.length !== this.weights.length) {
         return null;
     }
     var e = 0;
-    for (var n = 0; n < this.weights.length; n++) {
-        e += this.weights[n] * t[n];
+    for (var t = 0; t < this.weights.length; t++) {
+        e += this.weights[t] * r[t];
     }
     e += this.bias;
     if (e > 0) {
@@ -1141,121 +1141,121 @@ xr.prototype.predict = function r(t) {
     }
 };
 
-xr.prototype.train = function r(t, e) {
+e.prototype.train = function train(r, e) {
     if (e !== 0 && e !== 1) {
         return null;
     }
-    if (t.length !== this.weights.length) {
-        this.weights = t;
+    if (r.length !== this.weights.length) {
+        this.weights = r;
         this.bias = 1;
     }
-    var n = this.predict(t);
-    if (typeof n === "number" && n !== e) {
-        var a = e - n;
-        for (var o = 0; o < this.weights.length; o++) {
-            this.weights[o] += a * t[o];
+    var t = this.predict(r);
+    if (typeof t === "number" && t !== e) {
+        var n = e - t;
+        for (var a = 0; a < this.weights.length; a++) {
+            this.weights[a] += n * r[a];
         }
-        this.bias += a;
+        this.bias += n;
     }
     return this;
 };
 
-var qr = 1e-4;
+var t = 1e-4;
 
-function Er(r) {
+function factorial(r) {
     if (r < 0) {
         throw new Error("factorial requires a non-negative value");
     }
     if (Math.floor(r) !== r) {
         throw new Error("factorial requires an integer input");
     }
-    var t = 1;
-    for (var e = 2; e <= r; e++) {
-        t *= e;
+    var e = 1;
+    for (var t = 2; t <= r; t++) {
+        e *= t;
     }
-    return t;
+    return e;
 }
 
-function yr(r) {
+function gamma(r) {
     if (Number.isInteger(r)) {
         if (r <= 0) {
             return Number.NaN;
         } else {
-            return Er(r - 1);
+            return factorial(r - 1);
         }
     }
     r--;
     if (r < 0) {
-        return Math.PI / (Math.sin(Math.PI * -r) * yr(-r));
+        return Math.PI / (Math.sin(Math.PI * -r) * gamma(-r));
     } else {
-        var t = Math.pow(r / Math.E, r) * Math.sqrt(2 * Math.PI * (r + 1 / 6));
-        var e = r + 1 / 4;
-        var n = 1 + 1 / 144 / Math.pow(e, 2) - 1 / 12960 / Math.pow(e, 3) - 257 / 207360 / Math.pow(e, 4) - 52 / 2612736 / Math.pow(e, 5) + 5741173 / 9405849600 / Math.pow(e, 6) + 37529 / 18811699200 / Math.pow(e, 7);
-        return t * n;
+        var e = Math.pow(r / Math.E, r) * Math.sqrt(2 * Math.PI * (r + 1 / 6));
+        var t = r + 1 / 4;
+        var n = 1 + 1 / 144 / Math.pow(t, 2) - 1 / 12960 / Math.pow(t, 3) - 257 / 207360 / Math.pow(t, 4) - 52 / 2612736 / Math.pow(t, 5) + 5741173 / 9405849600 / Math.pow(t, 6) + 37529 / 18811699200 / Math.pow(t, 7);
+        return e * n;
     }
 }
 
-var Nr = [ .9999999999999971, 57.15623566586292, -59.59796035547549, 14.136097974741746, -.4919138160976202, 3399464998481189e-20, 4652362892704858e-20, -9837447530487956e-20, .0001580887032249125, -.00021026444172410488, .00021743961811521265, -.0001643181065367639, 8441822398385275e-20, -26190838401581408e-21, 36899182659531625e-22 ];
+var n = [ .9999999999999971, 57.15623566586292, -59.59796035547549, 14.136097974741746, -.4919138160976202, 3399464998481189e-20, 4652362892704858e-20, -9837447530487956e-20, .0001580887032249125, -.00021026444172410488, .00021743961811521265, -.0001643181065367639, 8441822398385275e-20, -26190838401581408e-21, 36899182659531625e-22 ];
 
-var Ir = 607 / 128;
+var a = 607 / 128;
 
-var kr = Math.log(Math.sqrt(2 * Math.PI));
+var i = Math.log(Math.sqrt(2 * Math.PI));
 
-function Ar(r) {
+function gammaln(r) {
     if (r <= 0) {
         return Number.POSITIVE_INFINITY;
     }
     r--;
-    var t = Nr[0];
-    for (var e = 1; e < 15; e++) {
-        t += Nr[e] / (r + e);
+    var e = n[0];
+    for (var t = 1; t < 15; t++) {
+        e += n[t] / (r + t);
     }
-    var n = Ir + .5 + r;
-    return kr + Math.log(t) - n + (r + .5) * Math.log(n);
+    var o = a + .5 + r;
+    return i + Math.log(e) - o + (r + .5) * Math.log(o);
 }
 
-function Pr(r) {
+function bernoulliDistribution(r) {
     if (r < 0 || r > 1) {
         throw new Error("bernoulliDistribution requires probability to be between 0 and 1 inclusive");
     }
     return [ 1 - r, r ];
 }
 
-function Sr(r, t) {
-    if (t < 0 || t > 1 || r <= 0 || r % 1 !== 0) {
+function binomialDistribution(r, e) {
+    if (e < 0 || e > 1 || r <= 0 || r % 1 !== 0) {
         return undefined;
     }
-    var e = 0;
     var n = 0;
-    var a = [];
+    var a = 0;
+    var i = [];
     var o = 1;
     do {
-        a[e] = o * Math.pow(t, e) * Math.pow(1 - t, r - e);
-        n += a[e];
-        e++;
-        o = o * (r - e + 1) / e;
-    } while (n < 1 - qr);
-    return a;
+        i[n] = o * Math.pow(e, n) * Math.pow(1 - e, r - n);
+        a += i[n];
+        n++;
+        o = o * (r - n + 1) / n;
+    } while (a < 1 - t);
+    return i;
 }
 
-function Cr(r) {
+function poissonDistribution(r) {
     if (r <= 0) {
         return undefined;
     }
-    var t = 0;
     var e = 0;
-    var n = [];
-    var a = 1;
+    var n = 0;
+    var a = [];
+    var i = 1;
     do {
-        n[t] = Math.exp(-r) * Math.pow(r, t) / a;
-        e += n[t];
-        t++;
-        a *= t;
-    } while (e < 1 - qr);
-    return n;
+        a[e] = Math.exp(-r) * Math.pow(r, e) / i;
+        n += a[e];
+        e++;
+        i *= e;
+    } while (n < 1 - t);
+    return a;
 }
 
-var _r = {
+var o = {
     1: {
         .995: 0,
         .99: 0,
@@ -1739,215 +1739,215 @@ var _r = {
     }
 };
 
-function Tr(r, t, e) {
-    var a = n(r);
-    var o = 0;
+function chiSquaredGoodnessOfFit(r, e, t) {
+    var n = mean(r);
+    var a = 0;
     var i = 1;
-    var u = t(a);
-    var f = [];
-    var v = [];
-    for (var h = 0; h < r.length; h++) {
-        if (f[r[h]] === undefined) {
-            f[r[h]] = 0;
+    var u = e(n);
+    var s = [];
+    var l = [];
+    for (var f = 0; f < r.length; f++) {
+        if (s[r[f]] === undefined) {
+            s[r[f]] = 0;
         }
-        f[r[h]]++;
+        s[r[f]]++;
     }
-    for (var s = 0; s < f.length; s++) {
-        if (f[s] === undefined) {
-            f[s] = 0;
+    for (var v = 0; v < s.length; v++) {
+        if (s[v] === undefined) {
+            s[v] = 0;
         }
     }
-    for (var l in u) {
-        if (l in f) {
-            v[+l] = u[l] * r.length;
+    for (var h in u) {
+        if (h in s) {
+            l[+h] = u[h] * r.length;
         }
     }
-    for (var c = v.length - 1; c >= 0; c--) {
-        if (v[c] < 3) {
-            v[c - 1] += v[c];
-            v.pop();
-            f[c - 1] += f[c];
-            f.pop();
+    for (var c = l.length - 1; c >= 0; c--) {
+        if (l[c] < 3) {
+            l[c - 1] += l[c];
+            l.pop();
+            s[c - 1] += s[c];
+            s.pop();
         }
     }
-    for (var p = 0; p < f.length; p++) {
-        o += Math.pow(f[p] - v[p], 2) / v[p];
+    for (var p = 0; p < s.length; p++) {
+        a += Math.pow(s[p] - l[p], 2) / l[p];
     }
-    var g = f.length - i - 1;
-    return _r[g][e] < o;
+    var m = s.length - i - 1;
+    return o[m][t] < a;
 }
 
-var Vr = Math.sqrt(2 * Math.PI);
+var u = Math.sqrt(2 * Math.PI);
 
-var Dr = {
+var s = {
     gaussian: function(r) {
-        return Math.exp(-.5 * r * r) / Vr;
+        return Math.exp(-.5 * r * r) / u;
     }
 };
 
-var Lr = {
+var l = {
     nrd: function(r) {
-        var t = Z(r);
-        var e = _(r);
-        if (typeof e === "number") {
-            t = Math.min(t, e / 1.34);
+        var e = sampleStandardDeviation(r);
+        var t = interquartileRange(r);
+        if (typeof t === "number") {
+            e = Math.min(e, t / 1.34);
         }
-        return 1.06 * t * Math.pow(r.length, -.2);
+        return 1.06 * e * Math.pow(r.length, -.2);
     }
 };
 
-function Rr(r, t, e) {
+function kernelDensityEstimation(r, e, t) {
     var n;
-    if (t === undefined) {
-        n = Dr.gaussian;
-    } else if (typeof t === "string") {
-        if (!Dr[t]) {
-            throw new Error('Unknown kernel "' + t + '"');
+    if (e === undefined) {
+        n = s.gaussian;
+    } else if (typeof e === "string") {
+        if (!s[e]) {
+            throw new Error('Unknown kernel "' + e + '"');
         }
-        n = Dr[t];
+        n = s[e];
     } else {
-        n = t;
+        n = e;
     }
     var a;
-    if (typeof e === "undefined") {
-        a = Lr.nrd(r);
-    } else if (typeof e === "string") {
-        if (!Lr[e]) {
-            throw new Error('Unknown bandwidth method "' + e + '"');
+    if (typeof t === "undefined") {
+        a = l.nrd(r);
+    } else if (typeof t === "string") {
+        if (!l[t]) {
+            throw new Error('Unknown bandwidth method "' + t + '"');
         }
-        a = Lr[e](r);
+        a = l[t](r);
     } else {
-        a = e;
+        a = t;
     }
-    return function(t) {
-        var e = 0;
-        var o = 0;
-        for (e = 0; e < r.length; e++) {
-            o += n((t - r[e]) / a);
+    return function(e) {
+        var t = 0;
+        var i = 0;
+        for (t = 0; t < r.length; t++) {
+            i += n((e - r[t]) / a);
         }
-        return o / a / r.length;
+        return i / a / r.length;
     };
 }
 
-function Ur(r, t, e) {
-    return (r - t) / e;
+function zScore(r, e, t) {
+    return (r - e) / t;
 }
 
-var zr = Math.sqrt(2 * Math.PI);
+var f = Math.sqrt(2 * Math.PI);
 
-function Fr(r) {
-    var t = r;
+function cumulativeDistribution(r) {
     var e = r;
+    var t = r;
     for (var n = 1; n < 15; n++) {
-        e *= r * r / (2 * n + 1);
-        t += e;
+        t *= r * r / (2 * n + 1);
+        e += t;
     }
-    return Math.round((.5 + t / zr * Math.exp(-r * r / 2)) * 1e4) / 1e4;
+    return Math.round((.5 + e / f * Math.exp(-r * r / 2)) * 1e4) / 1e4;
 }
 
-var Or = [];
+var v = [];
 
-for (var Xr = 0; Xr <= 3.09; Xr += .01) {
-    Or.push(Fr(Xr));
+for (var h = 0; h <= 3.09; h += .01) {
+    v.push(cumulativeDistribution(h));
 }
 
-function Yr(r) {
-    var t = Math.abs(r);
-    var e = Math.min(Math.round(t * 100), Or.length - 1);
+function cumulativeStdNormalProbability(r) {
+    var e = Math.abs(r);
+    var t = Math.min(Math.round(e * 100), v.length - 1);
     if (r >= 0) {
-        return Or[e];
+        return v[t];
     } else {
-        return Math.round((1 - Or[e]) * 1e4) / 1e4;
+        return Math.round((1 - v[t]) * 1e4) / 1e4;
     }
 }
 
-function jr(r) {
+function cumulativeStdLogisticProbability(r) {
     return 1 / (Math.exp(-r) + 1);
 }
 
-function Kr(r) {
-    var t = 1 / (1 + .5 * Math.abs(r));
-    var e = t * Math.exp(-r * r + ((((((((.17087277 * t - .82215223) * t + 1.48851587) * t - 1.13520398) * t + .27886807) * t - .18628806) * t + .09678418) * t + .37409196) * t + 1.00002368) * t - 1.26551223);
+function errorFunction(r) {
+    var e = 1 / (1 + .5 * Math.abs(r));
+    var t = e * Math.exp(-r * r + ((((((((.17087277 * e - .82215223) * e + 1.48851587) * e - 1.13520398) * e + .27886807) * e - .18628806) * e + .09678418) * e + .37409196) * e + 1.00002368) * e - 1.26551223);
     if (r >= 0) {
-        return 1 - e;
+        return 1 - t;
     } else {
-        return e - 1;
+        return t - 1;
     }
 }
 
-function Wr(r) {
-    var t = 8 * (Math.PI - 3) / (3 * Math.PI * (4 - Math.PI));
-    var e = Math.sqrt(Math.sqrt(Math.pow(2 / (Math.PI * t) + Math.log(1 - r * r) / 2, 2) - Math.log(1 - r * r) / t) - (2 / (Math.PI * t) + Math.log(1 - r * r) / 2));
+function inverseErrorFunction(r) {
+    var e = 8 * (Math.PI - 3) / (3 * Math.PI * (4 - Math.PI));
+    var t = Math.sqrt(Math.sqrt(Math.pow(2 / (Math.PI * e) + Math.log(1 - r * r) / 2, 2) - Math.log(1 - r * r) / e) - (2 / (Math.PI * e) + Math.log(1 - r * r) / 2));
     if (r >= 0) {
-        return e;
+        return t;
     } else {
-        return -e;
+        return -t;
     }
 }
 
-function Br(r) {
+function probit(r) {
     if (r === 0) {
-        r = qr;
+        r = t;
     } else if (r >= 1) {
-        r = 1 - qr;
+        r = 1 - t;
     }
-    return Math.sqrt(2) * Wr(2 * r - 1);
+    return Math.sqrt(2) * inverseErrorFunction(2 * r - 1);
 }
 
-function Gr(r) {
+function logit(r) {
     if (r <= 0 || r >= 1) {
         throw new Error("p must be strictly between zero and one");
     }
     return Math.log(r / (1 - r));
 }
 
-function Hr(r, t, e, a, o) {
-    if (a === undefined) {
-        a = 1e4;
+function permutationTest(r, e, t, n, a) {
+    if (n === undefined) {
+        n = 1e4;
     }
-    if (e === undefined) {
-        e = "two_side";
+    if (t === undefined) {
+        t = "two_side";
     }
-    if (e !== "two_side" && e !== "greater" && e !== "less") {
+    if (t !== "two_side" && t !== "greater" && t !== "less") {
         throw new Error("`alternative` must be either 'two_side', 'greater', or 'less'.");
     }
-    var i = n(r);
-    var u = n(t);
-    var f = i - u;
-    var v = new Array(a);
-    var h = r.concat(t);
-    var s = Math.floor(h.length / 2);
-    for (var l = 0; l < a; l++) {
-        R(h, o);
-        var c = h.slice(0, s);
-        var p = h.slice(s, h.length);
-        var g = n(c) - n(p);
-        v[l] = g;
-    }
-    var w = 0;
-    if (e === "two_side") {
-        for (var M = 0; M <= a; M++) {
-            if (Math.abs(v[M]) >= Math.abs(f)) {
-                w += 1;
+    var i = mean(r);
+    var o = mean(e);
+    var u = i - o;
+    var s = new Array(n);
+    var l = r.concat(e);
+    var f = Math.floor(l.length / 2);
+    for (var v = 0; v < n; v++) {
+        shuffleInPlace(l, a);
+        var h = l.slice(0, f);
+        var c = l.slice(f, l.length);
+        var p = mean(h) - mean(c);
+        s[v] = p;
+    }
+    var m = 0;
+    if (t === "two_side") {
+        for (var g = 0; g <= n; g++) {
+            if (Math.abs(s[g]) >= Math.abs(u)) {
+                m += 1;
             }
         }
-    } else if (e === "greater") {
-        for (var m = 0; m <= a; m++) {
-            if (v[m] >= f) {
-                w += 1;
+    } else if (t === "greater") {
+        for (var d = 0; d <= n; d++) {
+            if (s[d] >= u) {
+                m += 1;
             }
         }
     } else {
-        for (var d = 0; d <= a; d++) {
-            if (v[d] <= f) {
-                w += 1;
+        for (var w = 0; w <= n; w++) {
+            if (s[w] <= u) {
+                m += 1;
             }
         }
     }
-    return w / a;
+    return m / n;
 }
 
-function Jr(r) {
+function sign(r) {
     if (typeof r === "number") {
         if (r < 0) {
             return -1;
@@ -1961,59 +1961,59 @@ function Jr(r) {
     }
 }
 
-function Qr(r, t, e, n, a) {
+function bisect(r, e, t, n, a) {
     if (typeof r !== "function") {
         throw new TypeError("func must be a function");
     }
-    for (var o = 0; o < n; o++) {
-        var i = (t + e) / 2;
-        if (r(i) === 0 || Math.abs((e - t) / 2) < a) {
-            return i;
+    for (var i = 0; i < n; i++) {
+        var o = (e + t) / 2;
+        if (r(o) === 0 || Math.abs((t - e) / 2) < a) {
+            return o;
         }
-        if (Jr(r(i)) === Jr(r(t))) {
-            t = i;
+        if (sign(r(o)) === sign(r(e))) {
+            e = o;
         } else {
-            e = i;
+            t = o;
         }
     }
     throw new Error("maximum number of iterations exceeded");
 }
 
-function Zr(r, t) {
-    var e = 0;
+function euclideanDistance(r, e) {
+    var t = 0;
     for (var n = 0; n < r.length; n++) {
-        var a = r[n] - t[n];
-        e += a * a;
+        var a = r[n] - e[n];
+        t += a * a;
     }
-    return Math.sqrt(e);
+    return Math.sqrt(t);
 }
 
-function $r(r, t, e) {
-    if (e === void 0) e = Math.random;
+function kMeansCluster(r, e, t) {
+    if (t === void 0) t = Math.random;
     var n = null;
-    var a = z(r, t, e);
-    var o = null;
-    var i = Number.MAX_VALUE;
-    while (i !== 0) {
-        o = rt(r, a);
+    var a = sample(r, e, t);
+    var i = null;
+    var o = Number.MAX_VALUE;
+    while (o !== 0) {
+        i = labelPoints(r, a);
         n = a;
-        a = tt(r, o, t);
-        i = et(a, n);
+        a = calculateCentroids(r, i, e);
+        o = calculateChange(a, n);
     }
     return {
-        labels: o,
+        labels: i,
         centroids: a
     };
 }
 
-function rt(r, t) {
+function labelPoints(r, e) {
     return r.map((function(r) {
-        var e = Number.MAX_VALUE;
+        var t = Number.MAX_VALUE;
         var n = -1;
-        for (var a = 0; a < t.length; a++) {
-            var o = Zr(r, t[a]);
-            if (o < e) {
-                e = o;
+        for (var a = 0; a < e.length; a++) {
+            var i = euclideanDistance(r, e[a]);
+            if (i < t) {
+                t = i;
                 n = a;
             }
         }
@@ -2021,229 +2021,229 @@ function rt(r, t) {
     }));
 }
 
-function tt(r, t, e) {
+function calculateCentroids(r, e, t) {
     var n = r[0].length;
-    var a = F(e, n);
-    var o = Array(e).fill(0);
-    var i = r.length;
-    for (var u = 0; u < i; u++) {
-        var f = r[u];
-        var v = t[u];
-        var h = a[v];
-        for (var s = 0; s < n; s++) {
-            h[s] += f[s];
-        }
-        o[v] += 1;
-    }
-    for (var l = 0; l < e; l++) {
-        if (o[l] === 0) {
-            throw new Error("Centroid " + l + " has no friends");
-        }
-        var c = a[l];
+    var a = makeMatrix(t, n);
+    var i = Array(t).fill(0);
+    var o = r.length;
+    for (var u = 0; u < o; u++) {
+        var s = r[u];
+        var l = e[u];
+        var f = a[l];
+        for (var v = 0; v < n; v++) {
+            f[v] += s[v];
+        }
+        i[l] += 1;
+    }
+    for (var h = 0; h < t; h++) {
+        if (i[h] === 0) {
+            throw new Error("Centroid " + h + " has no friends");
+        }
+        var c = a[h];
         for (var p = 0; p < n; p++) {
-            c[p] /= o[l];
+            c[p] /= i[h];
         }
     }
     return a;
 }
 
-function et(r, t) {
-    var e = 0;
+function calculateChange(r, e) {
+    var t = 0;
     for (var n = 0; n < r.length; n++) {
-        e += Zr(r[n], t[n]);
+        t += euclideanDistance(r[n], e[n]);
     }
-    return e;
+    return t;
 }
 
-function nt(r, t) {
-    if (r.length !== t.length) {
+function silhouette(r, e) {
+    if (r.length !== e.length) {
         throw new Error("must have exactly as many labels as points");
     }
-    var e = at(t);
-    var n = ot(r);
+    var t = createGroups(e);
+    var n = calculateAllDistances(r);
     var a = [];
-    for (var o = 0; o < r.length; o++) {
-        var i = 0;
-        if (e[t[o]].length > 1) {
-            var u = ut(o, e[t[o]], n);
-            var f = it(o, t, e, n);
-            i = (f - u) / Math.max(u, f);
+    for (var i = 0; i < r.length; i++) {
+        var o = 0;
+        if (t[e[i]].length > 1) {
+            var u = meanDistanceFromPointToGroup(i, t[e[i]], n);
+            var s = meanDistanceToNearestGroup(i, e, t, n);
+            o = (s - u) / Math.max(u, s);
         }
-        a.push(i);
+        a.push(o);
     }
     return a;
 }
 
-function at(r) {
-    var t = 1 + c(r);
-    var e = Array(t);
+function createGroups(r) {
+    var e = 1 + max(r);
+    var t = Array(e);
     for (var n = 0; n < r.length; n++) {
         var a = r[n];
-        if (e[a] === undefined) {
-            e[a] = [];
+        if (t[a] === undefined) {
+            t[a] = [];
         }
-        e[a].push(n);
+        t[a].push(n);
     }
-    return e;
+    return t;
 }
 
-function ot(r) {
-    var t = r.length;
-    var e = F(t, t);
-    for (var n = 0; n < t; n++) {
+function calculateAllDistances(r) {
+    var e = r.length;
+    var t = makeMatrix(e, e);
+    for (var n = 0; n < e; n++) {
         for (var a = 0; a < n; a++) {
-            e[n][a] = Zr(r[n], r[a]);
-            e[a][n] = e[n][a];
+            t[n][a] = euclideanDistance(r[n], r[a]);
+            t[a][n] = t[n][a];
         }
     }
-    return e;
+    return t;
 }
 
-function it(r, t, e, n) {
-    var a = t[r];
-    var o = Number.MAX_VALUE;
-    for (var i = 0; i < e.length; i++) {
-        if (i !== a) {
-            var u = ut(r, e[i], n);
-            if (u < o) {
-                o = u;
+function meanDistanceToNearestGroup(r, e, t, n) {
+    var a = e[r];
+    var i = Number.MAX_VALUE;
+    for (var o = 0; o < t.length; o++) {
+        if (o !== a) {
+            var u = meanDistanceFromPointToGroup(r, t[o], n);
+            if (u < i) {
+                i = u;
             }
         }
     }
-    return o;
+    return i;
 }
 
-function ut(r, t, e) {
+function meanDistanceFromPointToGroup(r, e, t) {
     var n = 0;
-    for (var a = 0; a < t.length; a++) {
-        n += e[r][t[a]];
+    for (var a = 0; a < e.length; a++) {
+        n += t[r][e[a]];
     }
-    return n / t.length;
+    return n / e.length;
 }
 
-function ft(r, t) {
-    var e = nt(r, t);
-    return c(e);
+function silhouetteMetric(r, e) {
+    var t = silhouette(r, e);
+    return max(t);
 }
 
-function vt(r, t) {
-    if (r === 0 && t === 0) {
+function relativeError(r, e) {
+    if (r === 0 && e === 0) {
         return 0;
     }
-    return Math.abs((r - t) / t);
+    return Math.abs((r - e) / e);
 }
 
-function ht(r, t, e) {
-    if (e === void 0) e = qr;
-    return vt(r, t) <= e;
+function approxEqual(r, e, n) {
+    if (n === void 0) n = t;
+    return relativeError(r, e) <= n;
 }
 
-exports.addToMean = ir;
+exports.addToMean = addToMean;
 
-exports.approxEqual = ht;
+exports.approxEqual = approxEqual;
 
-exports.bernoulliDistribution = Pr;
+exports.bernoulliDistribution = bernoulliDistribution;
 
-exports.binomialDistribution = Sr;
+exports.binomialDistribution = binomialDistribution;
 
-exports.bisect = Qr;
+exports.bisect = bisect;
 
-exports.chunk = D;
+exports.chunk = chunk;
 
-exports.coefficientOfVariation = wr;
+exports.coefficientOfVariation = coefficientOfVariation;
 
-exports.combinations = ar;
+exports.combinations = combinations;
 
-exports.combinationsReplacement = or;
+exports.combinationsReplacement = combinationsReplacement;
 
-exports.combineMeans = ur;
+exports.combineMeans = combineMeans;
 
-exports.combineVariances = fr;
+exports.combineVariances = combineVariances;
 
-exports.cumulativeStdLogisticProbability = jr;
+exports.cumulativeStdLogisticProbability = cumulativeStdLogisticProbability;
 
-exports.epsilon = qr;
+exports.epsilon = t;
 
-exports.extent = p;
+exports.extent = extent;
 
-exports.factorial = Er;
+exports.factorial = factorial;
 
-exports.gamma = yr;
+exports.gamma = gamma;
 
-exports.gammaln = Ar;
+exports.gammaln = gammaln;
 
-exports.geometricMean = vr;
+exports.geometricMean = geometricMean;
 
-exports.harmonicMean = sr;
+exports.harmonicMean = harmonicMean;
 
-exports.interquartileRange = _;
+exports.interquartileRange = interquartileRange;
 
-exports.jenks = G;
+exports.jenks = jenks;
 
-exports.linearRegression = r;
+exports.linearRegression = linearRegression;
 
-exports.linearRegressionLine = t;
+exports.linearRegressionLine = linearRegressionLine;
 
-exports.logAverage = hr;
+exports.logAverage = logAverage;
 
-exports.logit = Gr;
+exports.logit = logit;
 
-exports.max = c;
+exports.max = max;
 
-exports.mean = n;
+exports.mean = mean;
 
-exports.median = T;
+exports.median = median;
 
-exports.medianAbsoluteDeviation = V;
+exports.medianAbsoluteDeviation = medianAbsoluteDeviation;
 
-exports.min = l;
+exports.min = min;
 
-exports.mode = h;
+exports.mode = mode;
 
-exports.poissonDistribution = Cr;
+exports.poissonDistribution = poissonDistribution;
 
-exports.probit = Br;
+exports.probit = probit;
 
-exports.product = d;
+exports.product = product;
 
-exports.quantile = E;
+exports.quantile = quantile;
 
-exports.quantileRank = C;
+exports.quantileRank = quantileRank;
 
-exports.quickselect = x;
+exports.quickselect = quickselect;
 
-exports.rSquared = u;
+exports.rSquared = rSquared;
 
-exports.rootMeanSquare = gr;
+exports.rootMeanSquare = rootMeanSquare;
 
-exports.sample = z;
+exports.sample = sample;
 
-exports.sampleCorrelation = $;
+exports.sampleCorrelation = sampleCorrelation;
 
-exports.sampleCovariance = J;
+exports.sampleCovariance = sampleCovariance;
 
-exports.sampleSkewness = tr;
+exports.sampleSkewness = sampleSkewness;
 
-exports.sampleStandardDeviation = Z;
+exports.sampleStandardDeviation = sampleStandardDeviation;
 
-exports.sampleVariance = Q;
+exports.sampleVariance = sampleVariance;
 
-exports.sampleWithReplacement = L;
+exports.sampleWithReplacement = sampleWithReplacement;
 
-exports.shuffle = U;
+exports.shuffle = shuffle;
 
-exports.standardDeviation = i;
+exports.standardDeviation = standardDeviation;
 
-exports.subtractFromMean = pr;
+exports.subtractFromMean = subtractFromMean;
 
-exports.sum = e;
+exports.sum = sum;
 
-exports.sumNthPowerDeviations = a;
+exports.sumNthPowerDeviations = sumNthPowerDeviations;
 
-exports.tTest = Mr;
+exports.tTest = tTest;
 
-exports.tTestTwoSample = mr;
+exports.tTestTwoSample = tTestTwoSample;
 
-exports.variance = o;
+exports.variance = variance;
 
-exports.zScore = Ur;
+exports.zScore = zScore;