@milaboratories/miplots4 1.0.157 → 1.0.158

package/dist/scatterplot/linearRegression.js

@@ -1,40 +1,40 @@
-import nn from "../node_modules/@stdlib/stats-base-dists-f-cdf/lib/index.js";
-import tn from "../node_modules/d3-array/src/extent.js";
-import en from "../node_modules/d3-shape/src/area.js";
-import g from "../node_modules/d3-array/src/sum.js";
-import w from "../node_modules/d3-array/src/mean.js";
-const D = 0.05;
+import nn from "@stdlib/stats-base-dists-f-cdf";
+import tn from "../node_modules/.pnpm/d3-array@3.2.4/node_modules/d3-array/src/extent.js";
+import en from "../node_modules/.pnpm/d3-shape@3.2.0/node_modules/d3-shape/src/area.js";
+import g from "../node_modules/.pnpm/d3-array@3.2.4/node_modules/d3-array/src/sum.js";
+import w from "../node_modules/.pnpm/d3-array@3.2.4/node_modules/d3-array/src/mean.js";
+const I = 0.05;
 function on(n) {
-  const e = 3.3871328727963665, o = 133.14166789178438, a = 1971.5909503065513, s = 13731.69376550946, h = 45921.95393154987, f = 67265.7709270087, u = 33430.57558358813, i = 2509.0809287301227, r = 42.31333070160091, c = 687.1870074920579, l = 5394.196021424751, x = 21213.794301586597, m = 39307.89580009271, N = 28729.085735721943, v = 5226.495278852854, X = 1.4234371107496835, p = 4.630337846156546, b = 5.769497221460691, q = 3.6478483247632045, y = 1.2704582524523684, V = 0.2417807251774506, Y = 0.022723844989269184, F = 7745450142783414e-19, R = 2.053191626637759, S = 1.6763848301838038, _ = 0.6897673349851, L = 0.14810397642748008, P = 0.015198666563616457, z = 5475938084995345e-19, A = 10507500716444169e-25, B = 6.657904643501103, O = 5.463784911164114, T = 1.7848265399172913, j = 0.29656057182850487, k = 0.026532189526576124, G = 0.0012426609473880784, H = 27115555687434876e-21, J = 20103343992922881e-23, Q = 0.599832206555888, U = 0.1369298809227358, W = 0.014875361290850615, Z = 7868691311456133e-19, $ = 18463183175100548e-21, K = 1421511758316446e-22, d = 20442631033899397e-31, C = n - 0.5;
+  const e = 3.3871328727963665, o = 133.14166789178438, a = 1971.5909503065513, s = 13731.69376550946, m = 45921.95393154987, f = 67265.7709270087, u = 33430.57558358813, i = 2509.0809287301227, r = 42.31333070160091, c = 687.1870074920579, l = 5394.196021424751, x = 21213.794301586597, h = 39307.89580009271, N = 28729.085735721943, v = 5226.495278852854, X = 1.4234371107496835, p = 4.630337846156546, b = 5.769497221460691, q = 3.6478483247632045, y = 1.2704582524523684, Y = 0.2417807251774506, D = 0.022723844989269184, R = 7745450142783414e-19, _ = 2.053191626637759, F = 1.6763848301838038, S = 0.6897673349851, L = 0.14810397642748008, P = 0.015198666563616457, $ = 5475938084995345e-19, z = 10507500716444169e-25, A = 6.657904643501103, B = 5.463784911164114, O = 1.7848265399172913, T = 0.29656057182850487, j = 0.026532189526576124, k = 0.0012426609473880784, G = 27115555687434876e-21, H = 20103343992922881e-23, J = 0.599832206555888, Q = 0.1369298809227358, U = 0.014875361290850615, W = 7868691311456133e-19, Z = 18463183175100548e-21, K = 1421511758316446e-22, d = 20442631033899397e-31, E = n - 0.5;
   let t, M;
-  return Math.abs(C) <= 0.425 ? (t = 0.180625 - C * C, M = C * (((((((i * t + u) * t + f) * t + h) * t + s) * t + a) * t + o) * t + e) / (((((((v * t + N) * t + m) * t + x) * t + l) * t + c) * t + r) * t + 1), M) : (t = n, t = Math.sqrt(-Math.log(t)), t <= 5 ? (t += -1.6, M = (((((((F * t + Y) * t + V) * t + y) * t + q) * t + b) * t + p) * t + X) / (((((((A * t + z) * t + P) * t + L) * t + _) * t + S) * t + R) * t + 1)) : (t += -5, M = (((((((J * t + H) * t + G) * t + k) * t + j) * t + T) * t + O) * t + B) / (((((((d * t + K) * t + $) * t + Z) * t + W) * t + U) * t + Q) * t + 1)), M = -M, M);
+  return Math.abs(E) <= 0.425 ? (t = 0.180625 - E * E, M = E * (((((((i * t + u) * t + f) * t + m) * t + s) * t + a) * t + o) * t + e) / (((((((v * t + N) * t + h) * t + x) * t + l) * t + c) * t + r) * t + 1), M) : (t = n, t = Math.sqrt(-Math.log(t)), t <= 5 ? (t += -1.6, M = (((((((R * t + D) * t + Y) * t + y) * t + q) * t + b) * t + p) * t + X) / (((((((z * t + $) * t + P) * t + L) * t + S) * t + F) * t + _) * t + 1)) : (t += -5, M = (((((((H * t + G) * t + k) * t + j) * t + T) * t + O) * t + B) * t + A) / (((((((d * t + K) * t + Z) * t + W) * t + U) * t + Q) * t + J) * t + 1)), M = -M, M);
 }
-function I(n, e) {
-  const { sin: o, cos: a, sqrt: s, pow: h, exp: f, PI: u } = Math;
+function V(n, e) {
+  const { sin: o, cos: a, sqrt: s, pow: m, exp: f, PI: u } = Math;
   let i, r;
   if (e == 1) return a(n * u / 2) / o(n * u / 2);
   if (e == 2) return s(2 / (n * (2 - n)) - 2);
   const c = 1 / (e - 0.5), l = 48 / (c * c);
   let x = ((20700 * c / l - 98) * c - 16) * c + 96.36;
-  const m = ((94.5 / (l + x) - 3) / l + 1) * s(c * u * 0.5) * e;
-  return i = m * n, r = h(i, 2 / e), r > 0.05 + c ? (i = on(n / 2), r = i * i, e < 5 && (x = x + 0.3 * (e - 4.5) * (i + 0.6)), x = (((0.05 * m * i - 5) * i - 7) * i - 2) * i + l + x, r = (((((0.4 * r + 6.3) * r + 36) * r + 94.5) / x - r - 3) / l + 1) * i, r = c * r * r, r > 2e-3 ? r = f(r) - 1 : r = 0.5 * r * r + r) : r = ((1 / (((e + 6) / (e * r) - 0.089 * m - 0.822) * (e + 2) * 3) + 0.5 / (e + 4)) * r - 1) * (e + 1) / (e + 2) + 1 / r, s(e * r);
+  const h = ((94.5 / (l + x) - 3) / l + 1) * s(c * u * 0.5) * e;
+  return i = h * n, r = m(i, 2 / e), r > 0.05 + c ? (i = on(n / 2), r = i * i, e < 5 && (x = x + 0.3 * (e - 4.5) * (i + 0.6)), x = (((0.05 * h * i - 5) * i - 7) * i - 2) * i + l + x, r = (((((0.4 * r + 6.3) * r + 36) * r + 94.5) / x - r - 3) / l + 1) * i, r = c * r * r, r > 2e-3 ? r = f(r) - 1 : r = 0.5 * r * r + r) : r = ((1 / (((e + 6) / (e * r) - 0.089 * h - 0.822) * (e + 2) * 3) + 0.5 / (e + 4)) * r - 1) * (e + 1) / (e + 2) + 1 / r, s(e * r);
 }
 function rn(n) {
-  let e = n.length, o = 0, a = 0, s = 0, h = 0;
+  let e = n.length, o = 0, a = 0, s = 0, m = 0;
   for (let i = 0; i < n.length; ++i) {
     const r = n[i].x, c = n[i].y;
     if (!Number.isFinite(r) || !Number.isFinite(c)) {
       --e;
       continue;
     }
-    o += r, a += c, s += r * c, h += r * r;
+    o += r, a += c, s += r * c, m += r * r;
   }
-  const f = (e * s - o * a) / (e * h - o * o), u = (a - f * o) / e;
+  const f = (e * s - o * a) / (e * m - o * o), u = (a - f * o) / e;
   return [f, u];
 }
 function cn(n, e, o) {
-  const a = (l) => l * e + o, s = n.length, h = g(n, (l) => l.x) / n.length, f = I(D, s - 2), u = Math.sqrt(
-    g(n, (l) => (l.y - a(l.x)) ** 2) / (s - 2) / g(n, (l) => (l.x - h) ** 2)
+  const a = (l) => l * e + o, s = n.length, m = g(n, (l) => l.x) / n.length, f = V(I, s - 2), u = Math.sqrt(
+    g(n, (l) => (l.y - a(l.x)) ** 2) / (s - 2) / g(n, (l) => (l.x - m) ** 2)
   ), i = u * f, c = u * Math.sqrt(g(n, (l) => l.x ** 2) / s) * f;
   return { slope: e, intercept: o, slopeError: i, interceptError: c };
 }
@@ -43,9 +43,9 @@ function sn(n, e) {
   let a = 0, s = 0;
   for (let u = 0; u < n.length; ++u)
     a += Math.pow(n[u].x - o, 2), s += Math.pow(n[u].y - e(n[u].x), 2);
-  const h = Math.sqrt(s / (n.length - 2)), f = I(D, n.length - 2);
+  const m = Math.sqrt(s / (n.length - 2)), f = V(I, n.length - 2);
   return function(u) {
-    const i = e(u), r = h * Math.sqrt(1 / n.length + Math.pow(u - o, 2) / a), c = i - f * r, l = i + f * r;
+    const i = e(u), r = m * Math.sqrt(1 / n.length + Math.pow(u - o, 2) / a), c = i - f * r, l = i + f * r;
     return { x: u, y: i, left: isNaN(c) ? i : c, right: isNaN(l) ? i : l };
   };
 }
@@ -70,48 +70,48 @@ function ln(n, e) {
   const o = w(n, e);
   return w(n.map((a) => (e(a) - o) ** 2));
 }
-function E(n, e) {
+function C(n, e) {
   return Math.sqrt(ln(n, e));
 }
 function an(n, e, o) {
-  const a = w(n, e), s = w(n, o), h = n.length;
-  return g(n, (f) => (f.x - a) * (f.y - s)) / h;
+  const a = w(n, e), s = w(n, o), m = n.length;
+  return g(n, (f) => (f.x - a) * (f.y - s)) / m;
 }
 function fn(n) {
   const e = an(
     n,
     (s) => s.x,
     (s) => s.y
-  ), o = E(n, (s) => s.x), a = E(n, (s) => s.y);
+  ), o = C(n, (s) => s.x), a = C(n, (s) => s.y);
   return e / (o * a);
 }
-function hn(n, e) {
-  const o = n.length, a = 1, s = 1, h = o - a - 1, f = w(n, (m) => m.y), u = g(n, (m) => (e(m.x) - f) ** 2), i = g(n, (m) => (e(m.x) - m.y) ** 2), r = u / s, c = i / h, l = r / c;
-  return 1 - nn(l, s, h);
+function mn(n, e) {
+  const o = n.length, a = 1, s = 1, m = o - a - 1, f = w(n, (h) => h.y), u = g(n, (h) => (e(h.x) - f) ** 2), i = g(n, (h) => (e(h.x) - h.y) ** 2), r = u / s, c = i / m, l = r / c;
+  return 1 - nn(l, s, m);
 }
-function Mn(n, e, o, a, s, h) {
-  return h ? o.reduce((f, u) => {
+function Mn(n, e, o, a, s, m) {
+  return m ? o.reduce((f, u) => {
     const i = e[u].dots;
     return f[u] = a.map((r) => {
       const c = i.filter(
         (p) => s.every((b, q) => String(n.getColumnValue(b.value, p.idx)) === r[q])
       );
       if (c.length === 0) return null;
-      const [l = 0, x = 0] = tn(c, (p) => p.x), [m, N] = rn(c), v = (p) => m * p + N, X = sn(c, v);
+      const [l = 0, x = 0] = tn(c, (p) => p.x), [h, N] = rn(c), v = (p) => h * p + N, X = sn(c, v);
       return {
         idx: c[0].idx,
         predict: v,
         getInterval: X,
         getArea: (p, b, q) => en().x((y) => p(y.x)).y0((y) => b(y.left)).y1((y) => b(y.right))(q.map(X)) ?? "",
-        color: h.color,
-        bounded: h.bounded,
+        color: m.color,
+        bounded: m.bounded,
         showStats: !0,
         xBounds: [l, x],
         stats: {
           r2: un(c),
           r: fn(c),
-          pValue: hn(c, v),
-          coefficients: cn(c, m, N)
+          pValue: mn(c, v),
+          coefficients: cn(c, h, N)
         }
       };
     }).filter((r) => r !== null), f;

package/dist/scatterplot/linearRegression.js.map

@@ -1 +1 @@
-{"version":3,"file":"linearRegression.js","sources":["../../src/scatterplot/linearRegression.ts"],"sourcesContent":["import fCDF from '@stdlib/stats-base-dists-f-cdf';\nimport { extent, mean, sum } from 'd3-array';\nimport type { ScaleLinear } from 'd3-scale';\nimport { area } from 'd3-shape';\nimport type { DataFrame } from '../DataFrame';\nimport type { ColumnName, ContinuousAesFromColumn, InheritAesScatterplot } from '../types';\nimport type { Dot, GroupedDots } from './dots';\nimport type { ScatterplotSettingsImpl } from './ScatterplotSettingsImpl';\n\ntype NumericDot = {x: number, y:number};\nconst CONFIDENCE_LEVEL = 0.05;\n\n/*\n Returns the normal deviate Z corresponding to a given lower tail area of P; Z is accurate to about 1 part in 10**16.\n\n Wichura, M. J. (1988) Algorithm AS 241: The percentage points of the normal distribution. Applied Statistics, 37, 477–484.\n\n extension of:\n Beasley, J. D./ Springer, S. G. (1977), The percentage points of the NormalDistribution, Applied Statistics. 26, 118–121.\n\n Excel: STANDNORMINV(); R: qnorm()\n */\n\nfunction normdev(p: number) {\n    if (p < 0 || p > 1) throw Error('P p < 0 || p > 1 in linear regression calculation');\n    if (p == 0) return -Infinity;\n    if (p == 1) return Infinity;\n\n    /* eslint-disable @typescript-eslint/no-loss-of-precision */\n    const a0 = 3.387132872796366608,\n        a1 = 1.3314166789178437745e2,\n        a2 = 1.9715909503065514427e3,\n        a3 = 1.3731693765509461125e4,\n        a4 = 4.5921953931549871457e4,\n        a5 = 6.7265770927008700853e4,\n        a6 = 3.3430575583588128105e4,\n        a7 = 2.5090809287301226727e3,\n        b1 = 4.2313330701600911252e1,\n        b2 = 6.871870074920579083e2,\n        b3 = 5.3941960214247511077e3,\n        b4 = 2.1213794301586595867e4,\n        b5 = 3.930789580009271061e4,\n        b6 = 2.8729085735721942674e4,\n        b7 = 5.226495278852854561e3,\n        c0 = 1.42343711074968357734,\n        c1 = 4.6303378461565452959,\n        c2 = 5.7694972214606914055,\n        c3 = 3.64784832476320460504,\n        c4 = 1.27045825245236838258,\n        c5 = 2.4178072517745061177e-1,\n        c6 = 2.27238449892691845833e-2,\n        c7 = 7.7454501427834140764e-4,\n        d1 = 2.05319162663775882187,\n        d2 = 1.6763848301838038494,\n        d3 = 6.8976733498510000455e-1,\n        d4 = 1.4810397642748007459e-1,\n        d5 = 1.51986665636164571966e-2,\n        d6 = 5.475938084995344946e-4,\n        d7 = 1.05075007164441684324e-9,\n        e0 = 6.6579046435011037772,\n        e1 = 5.4637849111641143699,\n        e2 = 1.7848265399172913358,\n        e3 = 2.9656057182850489123e-1,\n        e4 = 2.6532189526576123093e-2,\n        e5 = 1.2426609473880784386e-3,\n        e6 = 2.71155556874348757815e-5,\n        e7 = 2.01033439929228813265e-7,\n        f1 = 5.9983220655588793769e-1,\n        f2 = 1.3692988092273580531e-1,\n        f3 = 1.48753612908506148525e-2,\n        f4 = 7.868691311456132591e-4,\n        f5 = 1.8463183175100546818e-5,\n        f6 = 1.4215117583164458887e-7,\n        f7 = 2.04426310338993978564e-15;\n    /* eslint-enable @typescript-eslint/no-loss-of-precision */\n\n    const q = p - 0.5;\n    let r, z;\n\n    // p close to 0.5\n    if (Math.abs(q) <= 0.425) {\n        r = 0.180625 - q * q;\n        z =\n            (q * (((((((a7 * r + a6) * r + a5) * r + a4) * r + a3) * r + a2) * r + a1) * r + a0)) /\n            (((((((b7 * r + b6) * r + b5) * r + b4) * r + b3) * r + b2) * r + b1) * r + 1);\n        return z;\n    }\n\n    if (q > 0) r = 1 - p;\n    else r = p;\n    r = Math.sqrt(-Math.log(r));\n\n    // p neither close to 0.5 nor 0 or 1\n    if (r <= 5) {\n        r += -1.6;\n        z =\n            (((((((c7 * r + c6) * r + c5) * r + c4) * r + c3) * r + c2) * r + c1) * r + c0) /\n            (((((((d7 * r + d6) * r + d5) * r + d4) * r + d3) * r + d2) * r + d1) * r + 1);\n    }\n    // p near 0 or 1\n    else {\n        r += -5;\n        z =\n            (((((((e7 * r + e6) * r + e5) * r + e4) * r + e3) * r + e2) * r + e1) * r + e0) /\n            (((((((f7 * r + f6) * r + f5) * r + f4) * r + f3) * r + f2) * r + f1) * r + 1);\n    }\n\n    if (q < 0.0) z = -z;\n    return z;\n}\n\n/*\nHill's approximated inverse t-distribution calculates t given df and two-tail probability:\nHill, G. W. (1970), Algorithm 396: Student's t-quantiles. Communications of the ACM, 13(10), 619–620.\n\nResult should be \"correct to at least 6 significant digits even for the analytic continuation through noninteger values of n > 5\". For higher precision (used in R, ...) see:\nHill, G. W. (1981) Remark on Algorithm 396, ACM Transactions on Mathematical Software, 7, 250–1.\n\nExcel: TINV(); R: qt()\n*/\nfunction inverseT(p: number, df: number) {\n    const {sin, cos, sqrt, pow, exp, PI} = Math;\n    let x: number, y: number;\n\n    if (df == 1) return cos((p * PI) / 2) / sin((p * PI) / 2);\n    if (df == 2) return sqrt(2 / (p * (2 - p)) - 2);\n\n    const a = 1 / (df - 0.5);\n    const b = 48 / (a * a);\n    let c = (((20700 * a) / b - 98) * a - 16) * a + 96.36;\n    const d = ((94.5 / (b + c) - 3) / b + 1) * sqrt(a * PI * 0.5) * df;\n    x = d * p;\n    y = pow(x, 2 / df);\n\n    if (y > 0.05 + a) {\n        // The procedure normdev(p) is assumed to return a negative normal\n        // deviate at the lower tail probability level p, e.g. -2.32 for p = 0.01.\n        x = normdev(p / 2);\n        y = x * x;\n        if (df < 5) c = c + 0.3 * (df - 4.5) * (x + 0.6);\n        c = (((0.05 * d * x - 5) * x - 7) * x - 2) * x + b + c;\n        y = (((((0.4 * y + 6.3) * y + 36) * y + 94.5) / c - y - 3) / b + 1) * x;\n        y = a * y * y;\n        if (y > 0.002) y = exp(y) - 1;\n        else y = 0.5 * y * y + y;\n    } else {\n        y =\n            (((1 / (((df + 6) / (df * y) - 0.089 * d - 0.822) * (df + 2) * 3) + 0.5 / (df + 4)) * y - 1) * (df + 1)) /\n                (df + 2) +\n            1 / y;\n    }\n\n    return sqrt(df * y);\n}\n\nexport function linearRegression(data: NumericDot[]): [number, number] {\n    let n = data.length,\n        sumX = 0,\n        sumY = 0,\n        sumXY = 0,\n        sumX2 = 0;\n    for (let i = 0; i < data.length; ++i) {\n        const xi = data[i].x;\n        const yi = data[i].y;\n        if (!Number.isFinite(xi) || !Number.isFinite(yi)) {\n            --n;\n            continue;\n        }\n        sumX += xi;\n        sumY += yi;\n        sumXY += xi * yi;\n        sumX2 += xi * xi;\n    }\n    const slope = (n * sumXY - sumX * sumY) / (n * sumX2 - sumX * sumX);\n    const intercept = (sumY - slope * sumX) / n;\n\n    return [slope, intercept];\n}\n\nfunction getCoefficientsInterval(dots: NumericDot[], slope: number, intercept: number) {\n    const predict = (x: number) => x * slope + intercept;\n    const n = dots.length;\n    const meanX = sum(dots, d => d.x) / dots.length;\n    const t = inverseT(CONFIDENCE_LEVEL, n - 2);\n\n    const se_b = Math.sqrt(\n        sum(dots, dot => (dot.y - predict(dot.x)) ** 2) / (n - 2) / sum(dots, dot => (dot.x - meanX) ** 2)\n    );\n    const slopeError = se_b * t;\n\n    const se_a = se_b * Math.sqrt(sum(dots, dot => dot.x ** 2) / n);\n    const interceptError = se_a * t;\n\n    return {slope, intercept, slopeError, interceptError};\n}\n\nexport type PointIntervalInfo = {\n    x: number; y: number, left: number; right: number\n}\n// https://stats.stackexchange.com/questions/101318/understanding-shape-and-calculation-of-confidence-bands-in-linear-regression\nexport function confidenceInterval(\n    data: NumericDot[],\n    predict: (x: number) => number\n): (x: number) => PointIntervalInfo {\n    const mean = sum(data, d => d.x) / data.length;\n    let a = 0,\n        b = 0;\n    for (let i = 0; i < data.length; ++i) {\n        a += Math.pow(data[i].x - mean, 2);\n        b += Math.pow(data[i].y - predict(data[i].x), 2);\n    }\n    const sy = Math.sqrt(b / (data.length - 2));\n    const t = inverseT(CONFIDENCE_LEVEL, data.length - 2);\n    return function (x: number) {\n        const Y = predict(x);\n        const se = sy * Math.sqrt(1 / data.length + Math.pow(x - mean, 2) / a);\n        const left = Y - t * se;\n        const right = Y + t * se;\n        return {x, y: Y, left: isNaN(left) ? Y : left, right: isNaN(right) ? Y : right};\n    };\n}\n\nfunction getR2(dots: NumericDot[]) {\n    const n = dots.length;\n    return Math.pow(\n        (n * sum(dots, (d: NumericDot) => d.x * d.y) - sum(dots, (d: NumericDot) => d.x) * sum(dots, (d: NumericDot) => d.y)) /\n            (Math.sqrt(\n                n * sum(dots, d => d.x * d.x) -\n                    Math.pow(\n                        sum(dots, d => d.x),\n                        2\n                    )\n            ) *\n                Math.sqrt(\n                    n * sum(dots, d => d.y * d.y) -\n                        Math.pow(\n                            sum(dots, d => d.y),\n                            2\n                        )\n                )),\n        2\n    );\n}\n\nfunction getVariance(arr: NumericDot[], getter: (d: NumericDot) => number) {\n    const vMean = mean(arr, getter) as number;\n\n    return mean(arr.map(v => (getter(v) - vMean) ** 2)) as number;\n}\n\nfunction getStandardDeviation(arr: NumericDot[], getter: (d: NumericDot) => number) {\n    return Math.sqrt(getVariance(arr, getter));\n}\n\nfunction getCovariance(arr: NumericDot[], getterX: (d: NumericDot) => number, getterY: (d: NumericDot) => number) {\n    const meanX = mean(arr, getterX) as number;\n    const meanY = mean(arr, getterY) as number;\n    const n = arr.length;\n\n    return sum(arr, d => (d.x - meanX) * (d.y - meanY)) / n;\n}\n\nfunction getR(dots: NumericDot[]) {\n    const cov = getCovariance(\n        dots,\n        d => d.x,\n        d => d.y\n    );\n    const sigmaX = getStandardDeviation(dots, d => d.x);\n    const sigmaY = getStandardDeviation(dots, d => d.y);\n\n    return cov / (sigmaX * sigmaY);\n}\n\nfunction getPValue(dots: NumericDot[], predict: (x: number) => number) {\n    const n = dots.length;\n    const p = 1; // number of coefficients\n    const dfm = 1; // degrees of freedom model\n    const dfe = n - p - 1; // degrees of freedom error\n    const meanY = mean(dots, d => d.y) as number;\n    const ssr = sum(dots, d => (predict(d.x) - meanY) ** 2); // sum of squares regression\n    const sse = sum(dots, d => (predict(d.x) - d.y) ** 2); // sum of squares errors\n\n    const msr = ssr / dfm;\n    const mse = sse / dfe;\n\n    const f = msr / mse;\n\n    const pval = fCDF(f, dfm, dfe);\n\n    return 1.0 - pval; // \"greater\" p-value\n}\n\ntype TrendInfo = {\n    idx: number;\n    predict: (x: number) => number;\n    getInterval: (x: number) => PointIntervalInfo;\n    getArea: (x: ScaleLinear<number, number>, y: ScaleLinear<number, number>, range: number[]) => string;\n    color: string | InheritAesScatterplot | ContinuousAesFromColumn<string>;\n    bounded: boolean;\n    showStats: boolean;\n    xBounds: [number, number];\n    stats: {\n        r2?: number;\n        r?: number;\n        pValue?: number;\n        coefficients?: {slope: number, intercept: number, slopeError:number, interceptError: number};\n    };\n};\nexport type TrendsData = Record<string, TrendInfo[]>;\n\nexport function getRegressionData(\n    dataFrame: DataFrame,\n    dotsByFacets: GroupedDots,\n    facetKeys: string[],\n    groupingKeysValues: string[][],\n    grouping: ColumnName[],\n    trend: ScatterplotSettingsImpl['trend']\n): TrendsData | null {\n    if (!trend) {\n        return null;\n    }\n    return facetKeys.reduce((res: TrendsData, facetKey) => {\n        const facetDots = dotsByFacets[facetKey].dots;\n        \n        res[facetKey] = groupingKeysValues\n            .map(values => {\n                const dots = facetDots.filter(dot =>\n                    grouping.every((column, idx) => String(dataFrame.getColumnValue(column.value, dot.idx)) === values[idx])\n                ) as NumericDot[];\n\n                if (dots.length === 0) return null;\n\n                const [minX = 0, maxX = 0] = extent(dots, (dot: NumericDot) => dot.x);\n                const [slope, intercept] = linearRegression(dots);\n                const predict = (x: number) => slope * x + intercept;\n                const getInterval = confidenceInterval(dots, predict);\n\n                return {\n                    idx: (dots[0] as Dot).idx,\n                    predict,\n                    getInterval,\n                    getArea: (scaleX: ScaleLinear<number, number>, scaleY: ScaleLinear<number, number>, range: number[]) =>\n                        area<PointIntervalInfo>()\n                            .x(d => scaleX(d.x))\n                            .y0(d => scaleY(d.left))\n                            .y1(d => scaleY(d.right))(range.map(getInterval)) ?? '',\n                    color: trend.color,\n                    bounded: trend.bounded,\n                    showStats: true,\n                    xBounds: [minX, maxX],\n                    stats: {\n                        r2: getR2(dots),\n                        r: getR(dots),\n                        pValue: getPValue(dots, predict),\n                        coefficients: getCoefficientsInterval(dots, slope, intercept),\n                    },\n                };\n            })\n            .filter((item) => item !== null) as TrendInfo[];\n        \n        return res;\n    }, {});\n}\n"],"names":["CONFIDENCE_LEVEL","normdev","p","a0","a1","a2","a3","a4","a5","a6","a7","b1","b2","b3","b4","b5","b6","b7","c0","c1","c2","c3","c4","c5","c6","c7","d1","d2","d3","d4","d5","d6","d7","e0","e1","e2","e3","e4","e5","e6","e7","f1","f2","f3","f4","f5","f6","f7","q","r","z","inverseT","df","sin","cos","sqrt","pow","exp","PI","x","y","a","b","c","d","linearRegression","data","n","sumX","sumY","sumXY","sumX2","xi","yi","slope","intercept","getCoefficientsInterval","dots","predict","meanX","sum","t","se_b","dot","slopeError","interceptError","confidenceInterval","mean","i","sy","Y","se","left","right","getR2","getVariance","arr","getter","vMean","v","getStandardDeviation","getCovariance","getterX","getterY","meanY","getR","cov","sigmaX","sigmaY","getPValue","dfm","dfe","ssr","sse","msr","mse","f","fCDF","getRegressionData","dataFrame","dotsByFacets","facetKeys","groupingKeysValues","grouping","trend","res","facetKey","facetDots","values","column","idx","minX","maxX","extent","getInterval","scaleX","scaleY","range","area","item"],"mappings":";;;;;AAUA,MAAMA,IAAmB;AAazB,SAASC,GAAQC,GAAW;AAMxB,QAAMC,IAAK,oBACPC,IAAK,oBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,kBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,oBACLC,IAAK,oBACLC,IAAK,sBACLC,IAAK,sBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,iBACLC,IAAK,qBACLC,IAAK,sBACLC,IAAK,sBACLC,IAAK,uBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,qBACLC,IAAK,sBACLC,IAAK,uBACLC,IAAK,uBACLC,IAAK,uBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,sBACLC,IAAK,sBACLC,IAAK,uBACLC,IAAK,sBACLC,IAAK,uBAGHC,IAAI9C,IAAI;AACd,MAAI+C,GAAGC;AAGP,SAAI,KAAK,IAAIF,CAAC,KAAK,SACfC,IAAI,WAAWD,IAAIA,GACnBE,IACKF,WAAWtC,IAAKuC,IAAIxC,KAAMwC,IAAIzC,KAAMyC,IAAI1C,KAAM0C,IAAI3C,KAAM2C,IAAI5C,KAAM4C,IAAI7C,KAAM6C,IAAI9C,YAC1Ec,IAAKgC,IAAIjC,KAAMiC,IAAIlC,KAAMkC,IAAInC,KAAMmC,IAAIpC,KAAMoC,IAAIrC,KAAMqC,IAAItC,KAAMsC,IAAI,IACzEC,MAIND,IAAI/C,GACT+C,IAAI,KAAK,KAAK,CAAC,KAAK,IAAIA,CAAC,CAAC,GAGtBA,KAAK,KACLA,KAAK,MACLC,WACWzB,IAAKwB,IAAIzB,KAAMyB,IAAI1B,KAAM0B,IAAI3B,KAAM2B,IAAI5B,KAAM4B,IAAI7B,KAAM6B,IAAI9B,KAAM8B,IAAI/B,YACrEc,IAAKiB,IAAIlB,KAAMkB,IAAInB,KAAMmB,IAAIpB,KAAMoB,IAAIrB,KAAMqB,IAAItB,KAAMsB,IAAIvB,KAAMuB,IAAI,OAIhFA,KAAK,IACLC,WACWV,IAAKS,IAAIV,KAAMU,IAAIX,KAAMW,IAAIZ,KAAMY,IAAIb,KAAMa,IAAId,KAAMc,IAAIf,KAAMe,IAAIhB,YACrEc,IAAKE,IAAIH,KAAMG,IAAIJ,KAAMI,IAAIL,KAAMK,IAAIN,KAAMM,IAAIP,KAAMO,IAAIR,KAAMQ,IAAI,KAGvEC,IAAI,CAACA,GACXA;AACX;AAWA,SAASC,EAASjD,GAAWkD,GAAY;AACrC,QAAM,EAAC,KAAAC,GAAK,KAAAC,GAAK,MAAAC,GAAM,KAAAC,GAAK,KAAAC,GAAK,IAAAC,MAAM;AACvC,MAAIC,GAAWC;AAEf,MAAIR,KAAM,EAAG,QAAOE,EAAKpD,IAAIwD,IAAM,CAAC,IAAIL,EAAKnD,IAAIwD,IAAM,CAAC;AACxD,MAAIN,KAAM,EAAG,QAAOG,EAAK,KAAKrD,KAAK,IAAIA,MAAM,CAAC;AAE9C,QAAM2D,IAAI,KAAKT,IAAK,MACdU,IAAI,MAAMD,IAAIA;AACpB,MAAIE,MAAO,QAAQF,IAAKC,IAAI,MAAMD,IAAI,MAAMA,IAAI;AAChD,QAAMG,MAAM,QAAQF,IAAIC,KAAK,KAAKD,IAAI,KAAKP,EAAKM,IAAIH,IAAK,GAAG,IAAIN;AAChE,SAAAO,IAAIK,IAAI9D,GACR0D,IAAIJ,EAAIG,GAAG,IAAIP,CAAE,GAEbQ,IAAI,OAAOC,KAGXF,IAAI1D,GAAQC,IAAI,CAAC,GACjB0D,IAAID,IAAIA,GACJP,IAAK,MAAGW,IAAIA,IAAI,OAAOX,IAAK,QAAQO,IAAI,OAC5CI,OAAO,OAAOC,IAAIL,IAAI,KAAKA,IAAI,KAAKA,IAAI,KAAKA,IAAIG,IAAIC,GACrDH,SAAS,MAAMA,IAAI,OAAOA,IAAI,MAAMA,IAAI,QAAQG,IAAIH,IAAI,KAAKE,IAAI,KAAKH,GACtEC,IAAIC,IAAID,IAAIA,GACRA,IAAI,OAAOA,IAAIH,EAAIG,CAAC,IAAI,IACvBA,IAAI,MAAMA,IAAIA,IAAIA,KAEvBA,MACO,OAAOR,IAAK,MAAMA,IAAKQ,KAAK,QAAQI,IAAI,UAAUZ,IAAK,KAAK,KAAK,OAAOA,IAAK,MAAMQ,IAAI,MAAMR,IAAK,MAChGA,IAAK,KACV,IAAIQ,GAGLL,EAAKH,IAAKQ,CAAC;AACtB;AAEO,SAASK,GAAiBC,GAAsC;AACnE,MAAIC,IAAID,EAAK,QACTE,IAAO,GACPC,IAAO,GACPC,IAAQ,GACRC,IAAQ;AACZ,WAAS,IAAI,GAAG,IAAIL,EAAK,QAAQ,EAAE,GAAG;AAClC,UAAMM,IAAKN,EAAK,CAAC,EAAE,GACbO,IAAKP,EAAK,CAAC,EAAE;AACnB,QAAI,CAAC,OAAO,SAASM,CAAE,KAAK,CAAC,OAAO,SAASC,CAAE,GAAG;AAC9C,QAAEN;AACF;AAAA,IACJ;AACA,IAAAC,KAAQI,GACRH,KAAQI,GACRH,KAASE,IAAKC,GACdF,KAASC,IAAKA;AAAA,EAClB;AACA,QAAME,KAASP,IAAIG,IAAQF,IAAOC,MAASF,IAAII,IAAQH,IAAOA,IACxDO,KAAaN,IAAOK,IAAQN,KAAQD;AAE1C,SAAO,CAACO,GAAOC,CAAS;AAC5B;AAEA,SAASC,GAAwBC,GAAoBH,GAAeC,GAAmB;AACnF,QAAMG,IAAU,CAACnB,MAAcA,IAAIe,IAAQC,GACrCR,IAAIU,EAAK,QACTE,IAAQC,EAAIH,GAAM,OAAKb,EAAE,CAAC,IAAIa,EAAK,QACnCI,IAAI9B,EAASnD,GAAkBmE,IAAI,CAAC,GAEpCe,IAAO,KAAK;AAAA,IACdF,EAAIH,GAAM,CAAAM,OAAQA,EAAI,IAAIL,EAAQK,EAAI,CAAC,MAAM,CAAC,KAAKhB,IAAI,KAAKa,EAAIH,GAAM,QAAQM,EAAI,IAAIJ,MAAU,CAAC;AAAA,EAAA,GAE/FK,IAAaF,IAAOD,GAGpBI,IADOH,IAAO,KAAK,KAAKF,EAAIH,GAAM,CAAAM,MAAOA,EAAI,KAAK,CAAC,IAAIhB,CAAC,IAChCc;AAE9B,SAAO,EAAC,OAAAP,GAAO,WAAAC,GAAW,YAAAS,GAAY,gBAAAC,EAAA;AAC1C;AAMO,SAASC,GACZpB,GACAY,GACgC;AAChC,QAAMS,IAAOP,EAAId,GAAM,OAAKF,EAAE,CAAC,IAAIE,EAAK;AACxC,MAAI,IAAI,GACJJ,IAAI;AACR,WAAS0B,IAAI,GAAGA,IAAItB,EAAK,QAAQ,EAAEsB;AAC/B,SAAK,KAAK,IAAItB,EAAKsB,CAAC,EAAE,IAAID,GAAM,CAAC,GACjCzB,KAAK,KAAK,IAAII,EAAKsB,CAAC,EAAE,IAAIV,EAAQZ,EAAKsB,CAAC,EAAE,CAAC,GAAG,CAAC;AAEnD,QAAMC,IAAK,KAAK,KAAK3B,KAAKI,EAAK,SAAS,EAAE,GACpCe,IAAI9B,EAASnD,GAAkBkE,EAAK,SAAS,CAAC;AACpD,SAAO,SAAUP,GAAW;AACxB,UAAM+B,IAAIZ,EAAQnB,CAAC,GACbgC,IAAKF,IAAK,KAAK,KAAK,IAAIvB,EAAK,SAAS,KAAK,IAAIP,IAAI4B,GAAM,CAAC,IAAI,CAAC,GAC/DK,IAAOF,IAAIT,IAAIU,GACfE,IAAQH,IAAIT,IAAIU;AACtB,WAAO,EAAC,GAAAhC,GAAG,GAAG+B,GAAG,MAAM,MAAME,CAAI,IAAIF,IAAIE,GAAM,OAAO,MAAMC,CAAK,IAAIH,IAAIG,EAAA;AAAA,EAC7E;AACJ;AAEA,SAASC,GAAMjB,GAAoB;AAC/B,QAAMV,IAAIU,EAAK;AACf,SAAO,KAAK;AAAA,KACPV,IAAIa,EAAIH,GAAM,CAACb,MAAkBA,EAAE,IAAIA,EAAE,CAAC,IAAIgB,EAAIH,GAAM,CAACb,MAAkBA,EAAE,CAAC,IAAIgB,EAAIH,GAAM,CAACb,MAAkBA,EAAE,CAAC,MAC9G,KAAK;AAAA,MACFG,IAAIa,EAAIH,GAAM,CAAAb,MAAKA,EAAE,IAAIA,EAAE,CAAC,IACxB,KAAK;AAAA,QACDgB,EAAIH,GAAM,CAAAb,MAAKA,EAAE,CAAC;AAAA,QAClB;AAAA,MAAA;AAAA,IACJ,IAEJ,KAAK;AAAA,MACDG,IAAIa,EAAIH,GAAM,CAAAb,MAAKA,EAAE,IAAIA,EAAE,CAAC,IACxB,KAAK;AAAA,QACDgB,EAAIH,GAAM,CAAAb,MAAKA,EAAE,CAAC;AAAA,QAClB;AAAA,MAAA;AAAA,IACJ;AAAA,IAEhB;AAAA,EAAA;AAER;AAEA,SAAS+B,GAAYC,GAAmBC,GAAmC;AACvE,QAAMC,IAAQX,EAAKS,GAAKC,CAAM;AAE9B,SAAOV,EAAKS,EAAI,IAAI,CAAAG,OAAMF,EAAOE,CAAC,IAAID,MAAU,CAAC,CAAC;AACtD;AAEA,SAASE,EAAqBJ,GAAmBC,GAAmC;AAChF,SAAO,KAAK,KAAKF,GAAYC,GAAKC,CAAM,CAAC;AAC7C;AAEA,SAASI,GAAcL,GAAmBM,GAAoCC,GAAoC;AAC9G,QAAMxB,IAAQQ,EAAKS,GAAKM,CAAO,GACzBE,IAAQjB,EAAKS,GAAKO,CAAO,GACzBpC,IAAI6B,EAAI;AAEd,SAAOhB,EAAIgB,GAAK,CAAAhC,OAAMA,EAAE,IAAIe,MAAUf,EAAE,IAAIwC,EAAM,IAAIrC;AAC1D;AAEA,SAASsC,GAAK5B,GAAoB;AAC9B,QAAM6B,IAAML;AAAA,IACRxB;AAAA,IACA,OAAKb,EAAE;AAAA,IACP,OAAKA,EAAE;AAAA,EAAA,GAEL2C,IAASP,EAAqBvB,GAAM,CAAAb,MAAKA,EAAE,CAAC,GAC5C4C,IAASR,EAAqBvB,GAAM,CAAAb,MAAKA,EAAE,CAAC;AAElD,SAAO0C,KAAOC,IAASC;AAC3B;AAEA,SAASC,GAAUhC,GAAoBC,GAAgC;AACnE,QAAMX,IAAIU,EAAK,QACT3E,IAAI,GACJ4G,IAAM,GACNC,IAAM5C,IAAIjE,IAAI,GACdsG,IAAQjB,EAAKV,GAAM,CAAAb,MAAKA,EAAE,CAAC,GAC3BgD,IAAMhC,EAAIH,GAAM,CAAAb,OAAMc,EAAQd,EAAE,CAAC,IAAIwC,MAAU,CAAC,GAChDS,IAAMjC,EAAIH,GAAM,CAAAb,OAAMc,EAAQd,EAAE,CAAC,IAAIA,EAAE,MAAM,CAAC,GAE9CkD,IAAMF,IAAMF,GACZK,IAAMF,IAAMF,GAEZK,IAAIF,IAAMC;AAIhB,SAAO,IAFME,GAAKD,GAAGN,GAAKC,CAAG;AAGjC;AAoBO,SAASO,GACZC,GACAC,GACAC,GACAC,GACAC,GACAC,GACiB;AACjB,SAAKA,IAGEH,EAAU,OAAO,CAACI,GAAiBC,MAAa;AACnD,UAAMC,IAAYP,EAAaM,CAAQ,EAAE;AAEzC,WAAAD,EAAIC,CAAQ,IAAIJ,EACX,IAAI,CAAAM,MAAU;AACX,YAAMnD,IAAOkD,EAAU;AAAA,QAAO,OAC1BJ,EAAS,MAAM,CAACM,GAAQC,MAAQ,OAAOX,EAAU,eAAeU,EAAO,OAAO9C,EAAI,GAAG,CAAC,MAAM6C,EAAOE,CAAG,CAAC;AAAA,MAAA;AAG3G,UAAIrD,EAAK,WAAW,EAAG,QAAO;AAE9B,YAAM,CAACsD,IAAO,GAAGC,IAAO,CAAC,IAAIC,GAAOxD,GAAM,CAACM,MAAoBA,EAAI,CAAC,GAC9D,CAACT,GAAOC,CAAS,IAAIV,GAAiBY,CAAI,GAC1CC,IAAU,CAACnB,MAAce,IAAQf,IAAIgB,GACrC2D,IAAchD,GAAmBT,GAAMC,CAAO;AAEpD,aAAO;AAAA,QACH,KAAMD,EAAK,CAAC,EAAU;AAAA,QACtB,SAAAC;AAAA,QACA,aAAAwD;AAAA,QACA,SAAS,CAACC,GAAqCC,GAAqCC,MAChFC,GAAA,EACK,EAAE,CAAA1E,MAAKuE,EAAOvE,EAAE,CAAC,CAAC,EAClB,GAAG,CAAAA,MAAKwE,EAAOxE,EAAE,IAAI,CAAC,EACtB,GAAG,OAAKwE,EAAOxE,EAAE,KAAK,CAAC,EAAEyE,EAAM,IAAIH,CAAW,CAAC,KAAK;AAAA,QAC7D,OAAOV,EAAM;AAAA,QACb,SAASA,EAAM;AAAA,QACf,WAAW;AAAA,QACX,SAAS,CAACO,GAAMC,CAAI;AAAA,QACpB,OAAO;AAAA,UACH,IAAItC,GAAMjB,CAAI;AAAA,UACd,GAAG4B,GAAK5B,CAAI;AAAA,UACZ,QAAQgC,GAAUhC,GAAMC,CAAO;AAAA,UAC/B,cAAcF,GAAwBC,GAAMH,GAAOC,CAAS;AAAA,QAAA;AAAA,MAChE;AAAA,IAER,CAAC,EACA,OAAO,CAACgE,MAASA,MAAS,IAAI,GAE5Bd;AAAA,EACX,GAAG,CAAA,CAAE,IA1CM;AA2Cf;"}
+{"version":3,"file":"linearRegression.js","sources":["../../src/scatterplot/linearRegression.ts"],"sourcesContent":["/* eslint-disable no-loss-of-precision */\n// @ts-ignore\nimport fCDF from '@stdlib/stats-base-dists-f-cdf';\nimport { extent, mean, sum } from 'd3-array';\nimport type { ScaleLinear } from 'd3-scale';\nimport { area } from 'd3-shape';\nimport type { DataFrame } from '../DataFrame';\nimport type { ColumnName, ContinuousAesFromColumn, InheritAesScatterplot } from '../types';\nimport type { Dot, GroupedDots } from './dots';\nimport type { ScatterplotSettingsImpl } from './ScatterplotSettingsImpl';\n\ntype NumericDot = {x: number, y:number};\nconst CONFIDENCE_LEVEL = 0.05;\n\n/*\n Returns the normal deviate Z corresponding to a given lower tail area of P; Z is accurate to about 1 part in 10**16.\n\n Wichura, M. J. (1988) Algorithm AS 241: The percentage points of the normal distribution. Applied Statistics, 37, 477–484.\n\n extension of:\n Beasley, J. D./ Springer, S. G. (1977), The percentage points of the NormalDistribution, Applied Statistics. 26, 118–121.\n\n Excel: STANDNORMINV(); R: qnorm()\n */\n\nfunction normdev(p: number) {\n    if (p < 0 || p > 1) throw Error('P p < 0 || p > 1 in linear regression calculation');\n    if (p == 0) return -Infinity;\n    if (p == 1) return Infinity;\n\n    const a0 = 3.387132872796366608,\n        a1 = 1.3314166789178437745e2,\n        a2 = 1.9715909503065514427e3,\n        a3 = 1.3731693765509461125e4,\n        a4 = 4.5921953931549871457e4,\n        a5 = 6.7265770927008700853e4,\n        a6 = 3.3430575583588128105e4,\n        a7 = 2.5090809287301226727e3,\n        b1 = 4.2313330701600911252e1,\n        b2 = 6.871870074920579083e2,\n        b3 = 5.3941960214247511077e3,\n        b4 = 2.1213794301586595867e4,\n        b5 = 3.930789580009271061e4,\n        b6 = 2.8729085735721942674e4,\n        b7 = 5.226495278852854561e3,\n        c0 = 1.42343711074968357734,\n        c1 = 4.6303378461565452959,\n        c2 = 5.7694972214606914055,\n        c3 = 3.64784832476320460504,\n        c4 = 1.27045825245236838258,\n        c5 = 2.4178072517745061177e-1,\n        c6 = 2.27238449892691845833e-2,\n        c7 = 7.7454501427834140764e-4,\n        d1 = 2.05319162663775882187,\n        d2 = 1.6763848301838038494,\n        d3 = 6.8976733498510000455e-1,\n        d4 = 1.4810397642748007459e-1,\n        d5 = 1.51986665636164571966e-2,\n        d6 = 5.475938084995344946e-4,\n        d7 = 1.05075007164441684324e-9,\n        e0 = 6.6579046435011037772,\n        e1 = 5.4637849111641143699,\n        e2 = 1.7848265399172913358,\n        e3 = 2.9656057182850489123e-1,\n        e4 = 2.6532189526576123093e-2,\n        e5 = 1.2426609473880784386e-3,\n        e6 = 2.71155556874348757815e-5,\n        e7 = 2.01033439929228813265e-7,\n        f1 = 5.9983220655588793769e-1,\n        f2 = 1.3692988092273580531e-1,\n        f3 = 1.48753612908506148525e-2,\n        f4 = 7.868691311456132591e-4,\n        f5 = 1.8463183175100546818e-5,\n        f6 = 1.4215117583164458887e-7,\n        f7 = 2.04426310338993978564e-15;\n\n    const q = p - 0.5;\n    let r, z;\n\n    // p close to 0.5\n    if (Math.abs(q) <= 0.425) {\n        r = 0.180625 - q * q;\n        z =\n            (q * (((((((a7 * r + a6) * r + a5) * r + a4) * r + a3) * r + a2) * r + a1) * r + a0)) /\n            (((((((b7 * r + b6) * r + b5) * r + b4) * r + b3) * r + b2) * r + b1) * r + 1);\n        return z;\n    }\n\n    if (q > 0) r = 1 - p;\n    else r = p;\n    r = Math.sqrt(-Math.log(r));\n\n    // p neither close to 0.5 nor 0 or 1\n    if (r <= 5) {\n        r += -1.6;\n        z =\n            (((((((c7 * r + c6) * r + c5) * r + c4) * r + c3) * r + c2) * r + c1) * r + c0) /\n            (((((((d7 * r + d6) * r + d5) * r + d4) * r + d3) * r + d2) * r + d1) * r + 1);\n    }\n    // p near 0 or 1\n    else {\n        r += -5;\n        z =\n            (((((((e7 * r + e6) * r + e5) * r + e4) * r + e3) * r + e2) * r + e1) * r + e0) /\n            (((((((f7 * r + f6) * r + f5) * r + f4) * r + f3) * r + f2) * r + f1) * r + 1);\n    }\n\n    if (q < 0.0) z = -z;\n    return z;\n}\n\n/*\nHill's approximated inverse t-distribution calculates t given df and two-tail probability:\nHill, G. W. (1970), Algorithm 396: Student's t-quantiles. Communications of the ACM, 13(10), 619–620.\n\nResult should be \"correct to at least 6 significant digits even for the analytic continuation through noninteger values of n > 5\". For higher precision (used in R, ...) see:\nHill, G. W. (1981) Remark on Algorithm 396, ACM Transactions on Mathematical Software, 7, 250–1.\n\nExcel: TINV(); R: qt()\n*/\nfunction inverseT(p: number, df: number) {\n    const {sin, cos, sqrt, pow, exp, PI} = Math;\n    let x: number, y: number;\n\n    if (df == 1) return cos((p * PI) / 2) / sin((p * PI) / 2);\n    if (df == 2) return sqrt(2 / (p * (2 - p)) - 2);\n\n    const a = 1 / (df - 0.5);\n    const b = 48 / (a * a);\n    let c = (((20700 * a) / b - 98) * a - 16) * a + 96.36;\n    const d = ((94.5 / (b + c) - 3) / b + 1) * sqrt(a * PI * 0.5) * df;\n    x = d * p;\n    y = pow(x, 2 / df);\n\n    if (y > 0.05 + a) {\n        // The procedure normdev(p) is assumed to return a negative normal\n        // deviate at the lower tail probability level p, e.g. -2.32 for p = 0.01.\n        x = normdev(p / 2);\n        y = x * x;\n        if (df < 5) c = c + 0.3 * (df - 4.5) * (x + 0.6);\n        c = (((0.05 * d * x - 5) * x - 7) * x - 2) * x + b + c;\n        y = (((((0.4 * y + 6.3) * y + 36) * y + 94.5) / c - y - 3) / b + 1) * x;\n        y = a * y * y;\n        if (y > 0.002) y = exp(y) - 1;\n        else y = 0.5 * y * y + y;\n    } else {\n        y =\n            (((1 / (((df + 6) / (df * y) - 0.089 * d - 0.822) * (df + 2) * 3) + 0.5 / (df + 4)) * y - 1) * (df + 1)) /\n                (df + 2) +\n            1 / y;\n    }\n\n    return sqrt(df * y);\n}\n\nexport function linearRegression(data: NumericDot[]): [number, number] {\n    let n = data.length,\n        sumX = 0,\n        sumY = 0,\n        sumXY = 0,\n        sumX2 = 0;\n    for (let i = 0; i < data.length; ++i) {\n        const xi = data[i].x;\n        const yi = data[i].y;\n        if (!Number.isFinite(xi) || !Number.isFinite(yi)) {\n            --n;\n            continue;\n        }\n        sumX += xi;\n        sumY += yi;\n        sumXY += xi * yi;\n        sumX2 += xi * xi;\n    }\n    const slope = (n * sumXY - sumX * sumY) / (n * sumX2 - sumX * sumX);\n    const intercept = (sumY - slope * sumX) / n;\n\n    return [slope, intercept];\n}\n\nfunction getCoefficientsInterval(dots: NumericDot[], slope: number, intercept: number) {\n    const predict = (x: number) => x * slope + intercept;\n    const n = dots.length;\n    const meanX = sum(dots, d => d.x) / dots.length;\n    const t = inverseT(CONFIDENCE_LEVEL, n - 2);\n\n    const se_b = Math.sqrt(\n        sum(dots, dot => (dot.y - predict(dot.x)) ** 2) / (n - 2) / sum(dots, dot => (dot.x - meanX) ** 2)\n    );\n    const slopeError = se_b * t;\n\n    const se_a = se_b * Math.sqrt(sum(dots, dot => dot.x ** 2) / n);\n    const interceptError = se_a * t;\n\n    return {slope, intercept, slopeError, interceptError};\n}\n\nexport type PointIntervalInfo = {\n    x: number; y: number, left: number; right: number\n}\n// https://stats.stackexchange.com/questions/101318/understanding-shape-and-calculation-of-confidence-bands-in-linear-regression\nexport function confidenceInterval(\n    data: NumericDot[],\n    predict: (x: number) => number\n): (x: number) => PointIntervalInfo {\n    const mean = sum(data, d => d.x) / data.length;\n    let a = 0,\n        b = 0;\n    for (let i = 0; i < data.length; ++i) {\n        a += Math.pow(data[i].x - mean, 2);\n        b += Math.pow(data[i].y - predict(data[i].x), 2);\n    }\n    const sy = Math.sqrt(b / (data.length - 2));\n    const t = inverseT(CONFIDENCE_LEVEL, data.length - 2);\n    return function (x: number) {\n        const Y = predict(x);\n        const se = sy * Math.sqrt(1 / data.length + Math.pow(x - mean, 2) / a);\n        const left = Y - t * se;\n        const right = Y + t * se;\n        return {x, y: Y, left: isNaN(left) ? Y : left, right: isNaN(right) ? Y : right};\n    };\n}\n\nfunction getR2(dots: NumericDot[]) {\n    const n = dots.length;\n    return Math.pow(\n        (n * sum(dots, (d: NumericDot) => d.x * d.y) - sum(dots, (d: NumericDot) => d.x) * sum(dots, (d: NumericDot) => d.y)) /\n            (Math.sqrt(\n                n * sum(dots, d => d.x * d.x) -\n                    Math.pow(\n                        sum(dots, d => d.x),\n                        2\n                    )\n            ) *\n                Math.sqrt(\n                    n * sum(dots, d => d.y * d.y) -\n                        Math.pow(\n                            sum(dots, d => d.y),\n                            2\n                        )\n                )),\n        2\n    );\n}\n\nfunction getVariance(arr: NumericDot[], getter: (d: NumericDot) => number) {\n    const vMean = mean(arr, getter) as number;\n\n    return mean(arr.map(v => (getter(v) - vMean) ** 2)) as number;\n}\n\nfunction getStandardDeviation(arr: NumericDot[], getter: (d: NumericDot) => number) {\n    return Math.sqrt(getVariance(arr, getter));\n}\n\nfunction getCovariance(arr: NumericDot[], getterX: (d: NumericDot) => number, getterY: (d: NumericDot) => number) {\n    const meanX = mean(arr, getterX) as number;\n    const meanY = mean(arr, getterY) as number;\n    const n = arr.length;\n\n    return sum(arr, d => (d.x - meanX) * (d.y - meanY)) / n;\n}\n\nfunction getR(dots: NumericDot[]) {\n    const cov = getCovariance(\n        dots,\n        d => d.x,\n        d => d.y\n    );\n    const sigmaX = getStandardDeviation(dots, d => d.x);\n    const sigmaY = getStandardDeviation(dots, d => d.y);\n\n    return cov / (sigmaX * sigmaY);\n}\n\nfunction getPValue(dots: NumericDot[], predict: (x: number) => number) {\n    const n = dots.length;\n    const p = 1; // number of coefficients\n    const dfm = 1; // degrees of freedom model\n    const dfe = n - p - 1; // degrees of freedom error\n    const meanY = mean(dots, d => d.y) as number;\n    const ssr = sum(dots, d => (predict(d.x) - meanY) ** 2); // sum of squares regression\n    const sse = sum(dots, d => (predict(d.x) - d.y) ** 2); // sum of squares errors\n\n    const msr = ssr / dfm;\n    const mse = sse / dfe;\n\n    const f = msr / mse;\n\n    const pval = fCDF(f, dfm, dfe);\n\n    return 1.0 - pval; // \"greater\" p-value\n}\n\ntype TrendInfo = {\n    idx: number;\n    predict: (x: number) => number;\n    getInterval: (x: number) => PointIntervalInfo;\n    getArea: (x: ScaleLinear<number, number>, y: ScaleLinear<number, number>, range: number[]) => string;\n    color: string | InheritAesScatterplot | ContinuousAesFromColumn<string>;\n    bounded: boolean;\n    showStats: boolean;\n    xBounds: [number, number];\n    stats: {\n        r2?: number;\n        r?: number;\n        pValue?: number;\n        coefficients?: {slope: number, intercept: number, slopeError:number, interceptError: number};\n    };\n};\nexport type TrendsData = Record<string, TrendInfo[]>;\n\nexport function getRegressionData(\n    dataFrame: DataFrame,\n    dotsByFacets: GroupedDots,\n    facetKeys: string[],\n    groupingKeysValues: string[][],\n    grouping: ColumnName[],\n    trend: ScatterplotSettingsImpl['trend']\n): TrendsData | null {\n    if (!trend) {\n        return null;\n    }\n    return facetKeys.reduce((res: TrendsData, facetKey) => {\n        const facetDots = dotsByFacets[facetKey].dots;\n        \n        res[facetKey] = groupingKeysValues\n            .map(values => {\n                const dots = facetDots.filter(dot =>\n                    grouping.every((column, idx) => String(dataFrame.getColumnValue(column.value, dot.idx)) === values[idx])\n                ) as NumericDot[];\n\n                if (dots.length === 0) return null;\n\n                const [minX = 0, maxX = 0] = extent(dots, (dot: NumericDot) => dot.x);\n                const [slope, intercept] = linearRegression(dots);\n                const predict = (x: number) => slope * x + intercept;\n                const getInterval = confidenceInterval(dots, predict);\n\n                return {\n                    idx: (dots[0] as Dot).idx,\n                    predict,\n                    getInterval,\n                    getArea: (scaleX: ScaleLinear<number, number>, scaleY: ScaleLinear<number, number>, range: number[]) =>\n                        area<PointIntervalInfo>()\n                            .x(d => scaleX(d.x))\n                            .y0(d => scaleY(d.left))\n                            .y1(d => scaleY(d.right))(range.map(getInterval)) ?? '',\n                    color: trend.color,\n                    bounded: trend.bounded,\n                    showStats: true,\n                    xBounds: [minX, maxX],\n                    stats: {\n                        r2: getR2(dots),\n                        r: getR(dots),\n                        pValue: getPValue(dots, predict),\n                        coefficients: getCoefficientsInterval(dots, slope, intercept),\n                    },\n                };\n            })\n            .filter((item) => item !== null) as TrendInfo[];\n        \n        return res;\n    }, {});\n}\n"],"names":["CONFIDENCE_LEVEL","normdev","p","a0","a1","a2","a3","a4","a5","a6","a7","b1","b2","b3","b4","b5","b6","b7","c0","c1","c2","c3","c4","c5","c6","c7","d1","d2","d3","d4","d5","d6","d7","e0","e1","e2","e3","e4","e5","e6","e7","f1","f2","f3","f4","f5","f6","f7","q","r","z","inverseT","df","sin","cos","sqrt","pow","exp","PI","x","y","a","b","c","d","linearRegression","data","n","sumX","sumY","sumXY","sumX2","xi","yi","slope","intercept","getCoefficientsInterval","dots","predict","meanX","sum","t","se_b","dot","slopeError","interceptError","confidenceInterval","mean","i","sy","Y","se","left","right","getR2","getVariance","arr","getter","vMean","v","getStandardDeviation","getCovariance","getterX","getterY","meanY","getR","cov","sigmaX","sigmaY","getPValue","dfm","dfe","ssr","sse","msr","mse","f","fCDF","getRegressionData","dataFrame","dotsByFacets","facetKeys","groupingKeysValues","grouping","trend","res","facetKey","facetDots","values","column","idx","minX","maxX","extent","getInterval","scaleX","scaleY","range","area","item"],"mappings":";;;;;AAYA,MAAMA,IAAmB;AAazB,SAASC,GAAQC,GAAW;AAKxB,QAAMC,IAAK,oBACPC,IAAK,oBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,kBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,oBACLC,IAAK,oBACLC,IAAK,sBACLC,IAAK,sBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,iBACLC,IAAK,qBACLC,IAAK,sBACLC,IAAK,sBACLC,IAAK,uBACLC,IAAK,mBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,qBACLC,IAAK,sBACLC,IAAK,uBACLC,IAAK,uBACLC,IAAK,uBACLC,IAAK,mBACLC,IAAK,oBACLC,IAAK,sBACLC,IAAK,sBACLC,IAAK,uBACLC,IAAK,sBACLC,IAAK,uBAEHC,IAAI9C,IAAI;AACd,MAAI+C,GAAGC;AAGP,SAAI,KAAK,IAAIF,CAAC,KAAK,SACfC,IAAI,WAAWD,IAAIA,GACnBE,IACKF,WAAWtC,IAAKuC,IAAIxC,KAAMwC,IAAIzC,KAAMyC,IAAI1C,KAAM0C,IAAI3C,KAAM2C,IAAI5C,KAAM4C,IAAI7C,KAAM6C,IAAI9C,YAC1Ec,IAAKgC,IAAIjC,KAAMiC,IAAIlC,KAAMkC,IAAInC,KAAMmC,IAAIpC,KAAMoC,IAAIrC,KAAMqC,IAAItC,KAAMsC,IAAI,IACzEC,MAIND,IAAI/C,GACT+C,IAAI,KAAK,KAAK,CAAC,KAAK,IAAIA,CAAC,CAAC,GAGtBA,KAAK,KACLA,KAAK,MACLC,WACWzB,IAAKwB,IAAIzB,KAAMyB,IAAI1B,KAAM0B,IAAI3B,KAAM2B,IAAI5B,KAAM4B,IAAI7B,KAAM6B,IAAI9B,KAAM8B,IAAI/B,YACrEc,IAAKiB,IAAIlB,KAAMkB,IAAInB,KAAMmB,IAAIpB,KAAMoB,IAAIrB,KAAMqB,IAAItB,KAAMsB,IAAIvB,KAAMuB,IAAI,OAIhFA,KAAK,IACLC,WACWV,IAAKS,IAAIV,KAAMU,IAAIX,KAAMW,IAAIZ,KAAMY,IAAIb,KAAMa,IAAId,KAAMc,IAAIf,KAAMe,IAAIhB,YACrEc,IAAKE,IAAIH,KAAMG,IAAIJ,KAAMI,IAAIL,KAAMK,IAAIN,KAAMM,IAAIP,KAAMO,IAAIR,KAAMQ,IAAI,KAGvEC,IAAI,CAACA,GACXA;AACX;AAWA,SAASC,EAASjD,GAAWkD,GAAY;AACrC,QAAM,EAAC,KAAAC,GAAK,KAAAC,GAAK,MAAAC,GAAM,KAAAC,GAAK,KAAAC,GAAK,IAAAC,MAAM;AACvC,MAAIC,GAAWC;AAEf,MAAIR,KAAM,EAAG,QAAOE,EAAKpD,IAAIwD,IAAM,CAAC,IAAIL,EAAKnD,IAAIwD,IAAM,CAAC;AACxD,MAAIN,KAAM,EAAG,QAAOG,EAAK,KAAKrD,KAAK,IAAIA,MAAM,CAAC;AAE9C,QAAM2D,IAAI,KAAKT,IAAK,MACdU,IAAI,MAAMD,IAAIA;AACpB,MAAIE,MAAO,QAAQF,IAAKC,IAAI,MAAMD,IAAI,MAAMA,IAAI;AAChD,QAAMG,MAAM,QAAQF,IAAIC,KAAK,KAAKD,IAAI,KAAKP,EAAKM,IAAIH,IAAK,GAAG,IAAIN;AAChE,SAAAO,IAAIK,IAAI9D,GACR0D,IAAIJ,EAAIG,GAAG,IAAIP,CAAE,GAEbQ,IAAI,OAAOC,KAGXF,IAAI1D,GAAQC,IAAI,CAAC,GACjB0D,IAAID,IAAIA,GACJP,IAAK,MAAGW,IAAIA,IAAI,OAAOX,IAAK,QAAQO,IAAI,OAC5CI,OAAO,OAAOC,IAAIL,IAAI,KAAKA,IAAI,KAAKA,IAAI,KAAKA,IAAIG,IAAIC,GACrDH,SAAS,MAAMA,IAAI,OAAOA,IAAI,MAAMA,IAAI,QAAQG,IAAIH,IAAI,KAAKE,IAAI,KAAKH,GACtEC,IAAIC,IAAID,IAAIA,GACRA,IAAI,OAAOA,IAAIH,EAAIG,CAAC,IAAI,IACvBA,IAAI,MAAMA,IAAIA,IAAIA,KAEvBA,MACO,OAAOR,IAAK,MAAMA,IAAKQ,KAAK,QAAQI,IAAI,UAAUZ,IAAK,KAAK,KAAK,OAAOA,IAAK,MAAMQ,IAAI,MAAMR,IAAK,MAChGA,IAAK,KACV,IAAIQ,GAGLL,EAAKH,IAAKQ,CAAC;AACtB;AAEO,SAASK,GAAiBC,GAAsC;AACnE,MAAIC,IAAID,EAAK,QACTE,IAAO,GACPC,IAAO,GACPC,IAAQ,GACRC,IAAQ;AACZ,WAAS,IAAI,GAAG,IAAIL,EAAK,QAAQ,EAAE,GAAG;AAClC,UAAMM,IAAKN,EAAK,CAAC,EAAE,GACbO,IAAKP,EAAK,CAAC,EAAE;AACnB,QAAI,CAAC,OAAO,SAASM,CAAE,KAAK,CAAC,OAAO,SAASC,CAAE,GAAG;AAC9C,QAAEN;AACF;AAAA,IACJ;AACA,IAAAC,KAAQI,GACRH,KAAQI,GACRH,KAASE,IAAKC,GACdF,KAASC,IAAKA;AAAA,EAClB;AACA,QAAME,KAASP,IAAIG,IAAQF,IAAOC,MAASF,IAAII,IAAQH,IAAOA,IACxDO,KAAaN,IAAOK,IAAQN,KAAQD;AAE1C,SAAO,CAACO,GAAOC,CAAS;AAC5B;AAEA,SAASC,GAAwBC,GAAoBH,GAAeC,GAAmB;AACnF,QAAMG,IAAU,CAACnB,MAAcA,IAAIe,IAAQC,GACrCR,IAAIU,EAAK,QACTE,IAAQC,EAAIH,GAAM,OAAKb,EAAE,CAAC,IAAIa,EAAK,QACnCI,IAAI9B,EAASnD,GAAkBmE,IAAI,CAAC,GAEpCe,IAAO,KAAK;AAAA,IACdF,EAAIH,GAAM,CAAAM,OAAQA,EAAI,IAAIL,EAAQK,EAAI,CAAC,MAAM,CAAC,KAAKhB,IAAI,KAAKa,EAAIH,GAAM,QAAQM,EAAI,IAAIJ,MAAU,CAAC;AAAA,EAAA,GAE/FK,IAAaF,IAAOD,GAGpBI,IADOH,IAAO,KAAK,KAAKF,EAAIH,GAAM,CAAAM,MAAOA,EAAI,KAAK,CAAC,IAAIhB,CAAC,IAChCc;AAE9B,SAAO,EAAC,OAAAP,GAAO,WAAAC,GAAW,YAAAS,GAAY,gBAAAC,EAAA;AAC1C;AAMO,SAASC,GACZpB,GACAY,GACgC;AAChC,QAAMS,IAAOP,EAAId,GAAM,OAAKF,EAAE,CAAC,IAAIE,EAAK;AACxC,MAAI,IAAI,GACJJ,IAAI;AACR,WAAS0B,IAAI,GAAGA,IAAItB,EAAK,QAAQ,EAAEsB;AAC/B,SAAK,KAAK,IAAItB,EAAKsB,CAAC,EAAE,IAAID,GAAM,CAAC,GACjCzB,KAAK,KAAK,IAAII,EAAKsB,CAAC,EAAE,IAAIV,EAAQZ,EAAKsB,CAAC,EAAE,CAAC,GAAG,CAAC;AAEnD,QAAMC,IAAK,KAAK,KAAK3B,KAAKI,EAAK,SAAS,EAAE,GACpCe,IAAI9B,EAASnD,GAAkBkE,EAAK,SAAS,CAAC;AACpD,SAAO,SAAUP,GAAW;AACxB,UAAM+B,IAAIZ,EAAQnB,CAAC,GACbgC,IAAKF,IAAK,KAAK,KAAK,IAAIvB,EAAK,SAAS,KAAK,IAAIP,IAAI4B,GAAM,CAAC,IAAI,CAAC,GAC/DK,IAAOF,IAAIT,IAAIU,GACfE,IAAQH,IAAIT,IAAIU;AACtB,WAAO,EAAC,GAAAhC,GAAG,GAAG+B,GAAG,MAAM,MAAME,CAAI,IAAIF,IAAIE,GAAM,OAAO,MAAMC,CAAK,IAAIH,IAAIG,EAAA;AAAA,EAC7E;AACJ;AAEA,SAASC,GAAMjB,GAAoB;AAC/B,QAAMV,IAAIU,EAAK;AACf,SAAO,KAAK;AAAA,KACPV,IAAIa,EAAIH,GAAM,CAACb,MAAkBA,EAAE,IAAIA,EAAE,CAAC,IAAIgB,EAAIH,GAAM,CAACb,MAAkBA,EAAE,CAAC,IAAIgB,EAAIH,GAAM,CAACb,MAAkBA,EAAE,CAAC,MAC9G,KAAK;AAAA,MACFG,IAAIa,EAAIH,GAAM,CAAAb,MAAKA,EAAE,IAAIA,EAAE,CAAC,IACxB,KAAK;AAAA,QACDgB,EAAIH,GAAM,CAAAb,MAAKA,EAAE,CAAC;AAAA,QAClB;AAAA,MAAA;AAAA,IACJ,IAEJ,KAAK;AAAA,MACDG,IAAIa,EAAIH,GAAM,CAAAb,MAAKA,EAAE,IAAIA,EAAE,CAAC,IACxB,KAAK;AAAA,QACDgB,EAAIH,GAAM,CAAAb,MAAKA,EAAE,CAAC;AAAA,QAClB;AAAA,MAAA;AAAA,IACJ;AAAA,IAEhB;AAAA,EAAA;AAER;AAEA,SAAS+B,GAAYC,GAAmBC,GAAmC;AACvE,QAAMC,IAAQX,EAAKS,GAAKC,CAAM;AAE9B,SAAOV,EAAKS,EAAI,IAAI,CAAAG,OAAMF,EAAOE,CAAC,IAAID,MAAU,CAAC,CAAC;AACtD;AAEA,SAASE,EAAqBJ,GAAmBC,GAAmC;AAChF,SAAO,KAAK,KAAKF,GAAYC,GAAKC,CAAM,CAAC;AAC7C;AAEA,SAASI,GAAcL,GAAmBM,GAAoCC,GAAoC;AAC9G,QAAMxB,IAAQQ,EAAKS,GAAKM,CAAO,GACzBE,IAAQjB,EAAKS,GAAKO,CAAO,GACzBpC,IAAI6B,EAAI;AAEd,SAAOhB,EAAIgB,GAAK,CAAAhC,OAAMA,EAAE,IAAIe,MAAUf,EAAE,IAAIwC,EAAM,IAAIrC;AAC1D;AAEA,SAASsC,GAAK5B,GAAoB;AAC9B,QAAM6B,IAAML;AAAA,IACRxB;AAAA,IACA,OAAKb,EAAE;AAAA,IACP,OAAKA,EAAE;AAAA,EAAA,GAEL2C,IAASP,EAAqBvB,GAAM,CAAAb,MAAKA,EAAE,CAAC,GAC5C4C,IAASR,EAAqBvB,GAAM,CAAAb,MAAKA,EAAE,CAAC;AAElD,SAAO0C,KAAOC,IAASC;AAC3B;AAEA,SAASC,GAAUhC,GAAoBC,GAAgC;AACnE,QAAMX,IAAIU,EAAK,QACT3E,IAAI,GACJ4G,IAAM,GACNC,IAAM5C,IAAIjE,IAAI,GACdsG,IAAQjB,EAAKV,GAAM,CAAAb,MAAKA,EAAE,CAAC,GAC3BgD,IAAMhC,EAAIH,GAAM,CAAAb,OAAMc,EAAQd,EAAE,CAAC,IAAIwC,MAAU,CAAC,GAChDS,IAAMjC,EAAIH,GAAM,CAAAb,OAAMc,EAAQd,EAAE,CAAC,IAAIA,EAAE,MAAM,CAAC,GAE9CkD,IAAMF,IAAMF,GACZK,IAAMF,IAAMF,GAEZK,IAAIF,IAAMC;AAIhB,SAAO,IAFME,GAAKD,GAAGN,GAAKC,CAAG;AAGjC;AAoBO,SAASO,GACZC,GACAC,GACAC,GACAC,GACAC,GACAC,GACiB;AACjB,SAAKA,IAGEH,EAAU,OAAO,CAACI,GAAiBC,MAAa;AACnD,UAAMC,IAAYP,EAAaM,CAAQ,EAAE;AAEzC,WAAAD,EAAIC,CAAQ,IAAIJ,EACX,IAAI,CAAAM,MAAU;AACX,YAAMnD,IAAOkD,EAAU;AAAA,QAAO,OAC1BJ,EAAS,MAAM,CAACM,GAAQC,MAAQ,OAAOX,EAAU,eAAeU,EAAO,OAAO9C,EAAI,GAAG,CAAC,MAAM6C,EAAOE,CAAG,CAAC;AAAA,MAAA;AAG3G,UAAIrD,EAAK,WAAW,EAAG,QAAO;AAE9B,YAAM,CAACsD,IAAO,GAAGC,IAAO,CAAC,IAAIC,GAAOxD,GAAM,CAACM,MAAoBA,EAAI,CAAC,GAC9D,CAACT,GAAOC,CAAS,IAAIV,GAAiBY,CAAI,GAC1CC,IAAU,CAACnB,MAAce,IAAQf,IAAIgB,GACrC2D,IAAchD,GAAmBT,GAAMC,CAAO;AAEpD,aAAO;AAAA,QACH,KAAMD,EAAK,CAAC,EAAU;AAAA,QACtB,SAAAC;AAAA,QACA,aAAAwD;AAAA,QACA,SAAS,CAACC,GAAqCC,GAAqCC,MAChFC,GAAA,EACK,EAAE,CAAA1E,MAAKuE,EAAOvE,EAAE,CAAC,CAAC,EAClB,GAAG,CAAAA,MAAKwE,EAAOxE,EAAE,IAAI,CAAC,EACtB,GAAG,OAAKwE,EAAOxE,EAAE,KAAK,CAAC,EAAEyE,EAAM,IAAIH,CAAW,CAAC,KAAK;AAAA,QAC7D,OAAOV,EAAM;AAAA,QACb,SAASA,EAAM;AAAA,QACf,WAAW;AAAA,QACX,SAAS,CAACO,GAAMC,CAAI;AAAA,QACpB,OAAO;AAAA,UACH,IAAItC,GAAMjB,CAAI;AAAA,UACd,GAAG4B,GAAK5B,CAAI;AAAA,UACZ,QAAQgC,GAAUhC,GAAMC,CAAO;AAAA,UAC/B,cAAcF,GAAwBC,GAAMH,GAAOC,CAAS;AAAA,QAAA;AAAA,MAChE;AAAA,IAER,CAAC,EACA,OAAO,CAACgE,MAASA,MAAS,IAAI,GAE5Bd;AAAA,EACX,GAAG,CAAA,CAAE,IA1CM;AA2Cf;"}

package/dist/scatterplot/utils/createAesGetter.d.ts

@@ -1,4 +1,3 @@
 import { DataFrame } from '../../DataFrame';
 import { AesItem, ContinuousAesFromColumn, InheritAesScatterplot, LineShape, PointShape, ScatterplotLegendInfo } from '../../types';
 export declare function createAesGetter<ValueType extends string | number | PointShape | LineShape>(dataFrame: DataFrame, legendInfo: ScatterplotLegendInfo, aesItem: ValueType | ContinuousAesFromColumn<ValueType> | InheritAesScatterplot, field: keyof AesItem): (rowIndex: number) => ValueType;
-//# sourceMappingURL=createAesGetter.d.ts.map

package/dist/scatterplot/utils/createAesGetter.js

@@ -7,7 +7,7 @@ import "../../types/histogram.js";
 import "../../types/bubble.js";
 import { getContinuousColorScale as a } from "../../utils/getContinuousColorScale.js";
 import { DEFAULT_DOT_AES as g, DEFAULT_COMMON_AES as s } from "../constants.js";
-import { sqrt as f } from "../../node_modules/d3-scale/src/pow.js";
+import { sqrt as f } from "../../node_modules/.pnpm/d3-scale@4.0.2/node_modules/d3-scale/src/pow.js";
 function C(o) {
   return typeof o == "object" && o !== null && "type" in o && "value" in o && o.type === "grouping";
 }

package/dist/scatterplot/utils/createLegendInfo.d.ts

@@ -3,4 +3,3 @@ import { AesRecord, AesItem, Category, ColumnName, ScatterplotLayer, Scatterplot
 import { ScatterplotSettingsImpl } from '../ScatterplotSettingsImpl';
 export declare function addPalettesToAesMapping(columnName: ColumnName, usedAesFromPalettes: Partial<Record<keyof AesItem, string[]>>, inheritedAes?: AesRecord): (category: Category, field: keyof AesItem) => string | number | undefined;
 export declare function createLegendInfo(data: DataFrame, grouping: ScatterplotSettingsImpl['grouping'], layers: ScatterplotLayer[]): ScatterplotLegendInfo;
-//# sourceMappingURL=createLegendInfo.d.ts.map

package/dist/scatterplot/utils/formatColumnValue.d.ts

@@ -2,4 +2,3 @@ import { DataFrame } from '../../DataFrame';
 import { ColumnName } from '../../types';
 import { numberFormat } from '../../utils/numberFormat';
 export declare function formatColumnValue(data: DataFrame, rowIdx: number, column: ColumnName | null | undefined, numberFormatter?: typeof numberFormat): string | null;
-//# sourceMappingURL=formatColumnValue.d.ts.map

package/dist/scatterplot/utils/getTicksAndFormat.d.ts

@@ -5,4 +5,3 @@ export declare function getTicksAndFormat(scale: Scale, onlyPositive: boolean, h
     format: (v: number) => string;
 };
 export {};
-//# sourceMappingURL=getTicksAndFormat.d.ts.map

package/dist/scatterplot/utils/getVisibleLabels.d.ts

@@ -22,4 +22,3 @@ export declare function getLabelMinX(label: Label): number;
 export declare function getLabelMinY(label: Label): number;
 export declare function createLabelPositioner(chartW: number, chartH: number, yOffset?: number, xOffset?: number): (x: number, y: number, w: number, h: number, padding: number) => null | [XPosition, YPosition];
 export {};
-//# sourceMappingURL=getVisibleLabels.d.ts.map

package/dist/scatterplot/utils/getVisibleLabels.js

@@ -1,4 +1,4 @@
-import L from "../../node_modules/rbush/index.js";
+import L from "../../node_modules/.pnpm/rbush@4.0.1/node_modules/rbush/index.js";
 function d(n, r, t, o, i) {
   return o !== "middle" && i !== "middle" && (t *= 0.7), o === "right" ? n + t : o === "middle" ? n - r / 2 : n - r - t;
 }

package/dist/scatterplot/utils/getVisibleLabels.js.map

@@ -1 +1 @@
-{"version":3,"file":"getVisibleLabels.js","sources":["../../../src/scatterplot/utils/getVisibleLabels.ts"],"sourcesContent":["import RBush from 'rbush';\nimport type { Dot } from '../dots';\n\ntype XPosition = 'right' | 'middle' | 'left';\ntype YPosition = 'top' | 'middle' | 'bottom';\nexport type Label = {\n    x: number;\n    y: number;\n    dot: Dot;\n    name: string;\n    width: number;\n    height: number;\n    padding: number;\n    xPosition: XPosition;\n    yPosition: YPosition;\n};\nexport type LabelBBox = {\n    minX: number;\n    maxX: number;\n    minY: number;\n    maxY: number;\n}\n\nfunction getX(x: number, width: number, padding: number, xPos: XPosition, yPos: YPosition) {\n    const isCorner = xPos !== 'middle' && yPos !== 'middle';\n    isCorner && (padding *= 0.7);\n    return xPos === 'right'\n        ? x + padding\n        : xPos === 'middle'\n            ? x - width / 2\n            : x - width - padding;\n}\nfunction getY(y: number, height: number, padding: number, xPos: XPosition, yPos: YPosition) {\n    const isCorner = xPos !== 'middle' && yPos !== 'middle';\n    isCorner && (padding *= 0.7);\n    return yPos === 'middle'\n        ? y - height / 2\n        : yPos === 'top'\n            ? y - padding - height\n            : y + padding;\n}\n\nfunction getXOuter(x: number, width: number, padding: number, xPos: XPosition, yPos: YPosition, offsetH:number) {\n    return getX(x, width, padding, xPos, yPos) - offsetH;\n}\nfunction getYOuter(y: number, height: number, padding: number, xPos: XPosition, yPos: YPosition, offsetV:number) {\n    return getY(y, height, padding, xPos, yPos) - offsetV;\n}\nfunction getLabelOuterBBox(\n    x: number,\n    y: number,\n    width: number,\n    height: number,\n    padding: number,\n    xPos: XPosition,\n    yPos: YPosition,\n    yOffset: number,\n    xOffset: number,\n    out = {minX: 0, maxX: 0, minY: 0, maxY: 0}\n):LabelBBox {\n    out.minX = getXOuter(x, width, padding, xPos, yPos, xOffset); // consider padding for dot\n    out.maxX = getX(x, width, padding, xPos, yPos) + width + xOffset;\n    out.minY = getYOuter(y, height, padding, xPos, yPos, yOffset);\n    out.maxY = out.minY + height + 2 * yOffset;\n\n    return out;\n}\nfunction inBounds (bbox: LabelBBox, width: number, height: number) {\n    return bbox.minX >= 0 && bbox.maxX <= width && bbox.minY >= 0 && bbox.maxY <= height;\n}\n\nexport function getLabelMinX(label: Label) {\n    return getX(label.x, label.width, label.padding, label.xPosition, label.yPosition);\n}\nexport function getLabelMinY(label: Label) {\n    return getY(label.y, label.height, label.padding, label.xPosition, label.yPosition);\n}\n\nconst positionOptions: [YPosition, XPosition][] = [\n    ['middle', 'right'],\n    ['top', 'right'],\n    ['bottom', 'right'],\n    ['middle', 'left'],\n    ['top', 'left'],\n    ['bottom', 'left'],\n    ['top', 'middle'],\n    ['bottom', 'middle'],\n];\nconst LABEL_OFFSET = 2;\nconst MUT_BBOX: LabelBBox = {minX: 0, maxX: 0, minY: 0, maxY: 0};\nexport function createLabelPositioner(chartW: number, chartH: number, yOffset = LABEL_OFFSET, xOffset = LABEL_OFFSET) {\n    const outerBBox = MUT_BBOX;\n    const bboxTree = new RBush<LabelBBox>();\n    \n    return (x: number, y: number, w: number, h: number, padding: number): null | [XPosition, YPosition] => {\n        for (const [yPosition, xPosition] of positionOptions) {\n            getLabelOuterBBox(x, y, w, h, padding, xPosition, yPosition, yOffset, xOffset, outerBBox);\n\n            if (inBounds(outerBBox, chartW, chartH) && !bboxTree.collides(outerBBox)) {\n                bboxTree.insert({...outerBBox});\n                return [xPosition, yPosition];\n            }\n        }\n\n        return null;\n    };\n}"],"names":["getX","x","width","padding","xPos","yPos","getY","y","height","getXOuter","offsetH","getYOuter","offsetV","getLabelOuterBBox","yOffset","xOffset","out","inBounds","bbox","getLabelMinX","label","getLabelMinY","positionOptions","LABEL_OFFSET","MUT_BBOX","createLabelPositioner","chartW","chartH","outerBBox","bboxTree","RBush","w","h","yPosition","xPosition"],"mappings":";AAuBA,SAASA,EAAKC,GAAWC,GAAeC,GAAiBC,GAAiBC,GAAiB;AAEvF,SADiBD,MAAS,YAAYC,MAAS,aAClCF,KAAW,MACjBC,MAAS,UACVH,IAAIE,IACJC,MAAS,WACLH,IAAIC,IAAQ,IACZD,IAAIC,IAAQC;AAC1B;AACA,SAASG,EAAKC,GAAWC,GAAgBL,GAAiBC,GAAiBC,GAAiB;AAExF,SADiBD,MAAS,YAAYC,MAAS,aAClCF,KAAW,MACjBE,MAAS,WACVE,IAAIC,IAAS,IACbH,MAAS,QACLE,IAAIJ,IAAUK,IACdD,IAAIJ;AAClB;AAEA,SAASM,EAAUR,GAAWC,GAAeC,GAAiBC,GAAiBC,GAAiBK,GAAgB;AAC5G,SAAOV,EAAKC,GAAGC,GAAOC,GAASC,GAAMC,CAAI,IAAIK;AACjD;AACA,SAASC,EAAUJ,GAAWC,GAAgBL,GAAiBC,GAAiBC,GAAiBO,GAAgB;AAC7G,SAAON,EAAKC,GAAGC,GAAQL,GAASC,GAAMC,CAAI,IAAIO;AAClD;AACA,SAASC,EACLZ,GACAM,GACAL,GACAM,GACAL,GACAC,GACAC,GACAS,GACAC,GACAC,IAAM,EAAC,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,KAChC;AACR,SAAAA,EAAI,OAAOP,EAAUR,GAAGC,GAAOC,GAASC,GAAMC,GAAMU,CAAO,GAC3DC,EAAI,OAAOhB,EAAKC,GAAGC,GAAOC,GAASC,GAAMC,CAAI,IAAIH,IAAQa,GACzDC,EAAI,OAAOL,EAAUJ,GAAGC,GAAQL,GAASC,GAAMC,GAAMS,CAAO,GAC5DE,EAAI,OAAOA,EAAI,OAAOR,IAAS,IAAIM,GAE5BE;AACX;AACA,SAASC,EAAUC,GAAiBhB,GAAeM,GAAgB;AAC/D,SAAOU,EAAK,QAAQ,KAAKA,EAAK,QAAQhB,KAASgB,EAAK,QAAQ,KAAKA,EAAK,QAAQV;AAClF;AAEO,SAASW,EAAaC,GAAc;AACvC,SAAOpB,EAAKoB,EAAM,GAAGA,EAAM,OAAOA,EAAM,SAASA,EAAM,WAAWA,EAAM,SAAS;AACrF;AACO,SAASC,EAAaD,GAAc;AACvC,SAAOd,EAAKc,EAAM,GAAGA,EAAM,QAAQA,EAAM,SAASA,EAAM,WAAWA,EAAM,SAAS;AACtF;AAEA,MAAME,IAA4C;AAAA,EAC9C,CAAC,UAAU,OAAO;AAAA,EAClB,CAAC,OAAO,OAAO;AAAA,EACf,CAAC,UAAU,OAAO;AAAA,EAClB,CAAC,UAAU,MAAM;AAAA,EACjB,CAAC,OAAO,MAAM;AAAA,EACd,CAAC,UAAU,MAAM;AAAA,EACjB,CAAC,OAAO,QAAQ;AAAA,EAChB,CAAC,UAAU,QAAQ;AACvB,GACMC,IAAe,GACfC,IAAsB,EAAC,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,EAAA;AACvD,SAASC,EAAsBC,GAAgBC,GAAgBb,IAAUS,GAAcR,IAAUQ,GAAc;AAClH,QAAMK,IAAYJ,GACZK,IAAW,IAAIC,EAAA;AAErB,SAAO,CAAC7B,GAAWM,GAAWwB,GAAWC,GAAW7B,MAAmD;AACnG,eAAW,CAAC8B,GAAWC,CAAS,KAAKZ;AAGjC,UAFAT,EAAkBZ,GAAGM,GAAGwB,GAAGC,GAAG7B,GAAS+B,GAAWD,GAAWnB,GAASC,GAASa,CAAS,GAEpFX,EAASW,GAAWF,GAAQC,CAAM,KAAK,CAACE,EAAS,SAASD,CAAS;AACnE,eAAAC,EAAS,OAAO,EAAC,GAAGD,GAAU,GACvB,CAACM,GAAWD,CAAS;AAIpC,WAAO;AAAA,EACX;AACJ;"}
+{"version":3,"file":"getVisibleLabels.js","sources":["../../../src/scatterplot/utils/getVisibleLabels.ts"],"sourcesContent":["import RBush from 'rbush';\nimport type { Dot } from '../dots';\n\ntype XPosition = 'right' | 'middle' | 'left';\ntype YPosition = 'top' | 'middle' | 'bottom';\nexport type Label = {\n    x: number;\n    y: number;\n    dot: Dot;\n    name: string;\n    width: number;\n    height: number;\n    padding: number;\n    xPosition: XPosition;\n    yPosition: YPosition;\n};\nexport type LabelBBox = {\n    minX: number;\n    maxX: number;\n    minY: number;\n    maxY: number;\n}\n\nfunction getX(x: number, width: number, padding: number, xPos: XPosition, yPos: YPosition) {\n    const isCorner = xPos !== 'middle' && yPos !== 'middle';\n    // eslint-disable-next-line @typescript-eslint/no-unused-expressions\n    isCorner && (padding *= 0.7);\n    return xPos === 'right'\n        ? x + padding\n        : xPos === 'middle'\n            ? x - width / 2\n            : x - width - padding;\n}\nfunction getY(y: number, height: number, padding: number, xPos: XPosition, yPos: YPosition) {\n    const isCorner = xPos !== 'middle' && yPos !== 'middle';\n    // eslint-disable-next-line @typescript-eslint/no-unused-expressions\n    isCorner && (padding *= 0.7);\n    return yPos === 'middle'\n        ? y - height / 2\n        : yPos === 'top'\n            ? y - padding - height\n            : y + padding;\n}\n\nfunction getXOuter(x: number, width: number, padding: number, xPos: XPosition, yPos: YPosition, offsetH:number) {\n    return getX(x, width, padding, xPos, yPos) - offsetH;\n}\nfunction getYOuter(y: number, height: number, padding: number, xPos: XPosition, yPos: YPosition, offsetV:number) {\n    return getY(y, height, padding, xPos, yPos) - offsetV;\n}\nfunction getLabelOuterBBox(\n    x: number,\n    y: number,\n    width: number,\n    height: number,\n    padding: number,\n    xPos: XPosition,\n    yPos: YPosition,\n    yOffset: number,\n    xOffset: number,\n    out = {minX: 0, maxX: 0, minY: 0, maxY: 0}\n):LabelBBox {\n    out.minX = getXOuter(x, width, padding, xPos, yPos, xOffset); // consider padding for dot\n    out.maxX = getX(x, width, padding, xPos, yPos) + width + xOffset;\n    out.minY = getYOuter(y, height, padding, xPos, yPos, yOffset);\n    out.maxY = out.minY + height + 2 * yOffset;\n\n    return out;\n}\nfunction inBounds (bbox: LabelBBox, width: number, height: number) {\n    return bbox.minX >= 0 && bbox.maxX <= width && bbox.minY >= 0 && bbox.maxY <= height;\n}\n\nexport function getLabelMinX(label: Label) {\n    return getX(label.x, label.width, label.padding, label.xPosition, label.yPosition);\n}\nexport function getLabelMinY(label: Label) {\n    return getY(label.y, label.height, label.padding, label.xPosition, label.yPosition);\n}\n\nconst positionOptions: [YPosition, XPosition][] = [\n    ['middle', 'right'],\n    ['top', 'right'],\n    ['bottom', 'right'],\n    ['middle', 'left'],\n    ['top', 'left'],\n    ['bottom', 'left'],\n    ['top', 'middle'],\n    ['bottom', 'middle'],\n];\nconst LABEL_OFFSET = 2;\nconst MUT_BBOX: LabelBBox = {minX: 0, maxX: 0, minY: 0, maxY: 0};\nexport function createLabelPositioner(chartW: number, chartH: number, yOffset = LABEL_OFFSET, xOffset = LABEL_OFFSET) {\n    const outerBBox = MUT_BBOX;\n    const bboxTree = new RBush<LabelBBox>();\n    \n    return (x: number, y: number, w: number, h: number, padding: number): null | [XPosition, YPosition] => {\n        for (const [yPosition, xPosition] of positionOptions) {\n            getLabelOuterBBox(x, y, w, h, padding, xPosition, yPosition, yOffset, xOffset, outerBBox);\n\n            if (inBounds(outerBBox, chartW, chartH) && !bboxTree.collides(outerBBox)) {\n                bboxTree.insert({...outerBBox});\n                return [xPosition, yPosition];\n            }\n        }\n\n        return null;\n    };\n}"],"names":["getX","x","width","padding","xPos","yPos","getY","y","height","getXOuter","offsetH","getYOuter","offsetV","getLabelOuterBBox","yOffset","xOffset","out","inBounds","bbox","getLabelMinX","label","getLabelMinY","positionOptions","LABEL_OFFSET","MUT_BBOX","createLabelPositioner","chartW","chartH","outerBBox","bboxTree","RBush","w","h","yPosition","xPosition"],"mappings":";AAuBA,SAASA,EAAKC,GAAWC,GAAeC,GAAiBC,GAAiBC,GAAiB;AAGvF,SAFiBD,MAAS,YAAYC,MAAS,aAElCF,KAAW,MACjBC,MAAS,UACVH,IAAIE,IACJC,MAAS,WACLH,IAAIC,IAAQ,IACZD,IAAIC,IAAQC;AAC1B;AACA,SAASG,EAAKC,GAAWC,GAAgBL,GAAiBC,GAAiBC,GAAiB;AAGxF,SAFiBD,MAAS,YAAYC,MAAS,aAElCF,KAAW,MACjBE,MAAS,WACVE,IAAIC,IAAS,IACbH,MAAS,QACLE,IAAIJ,IAAUK,IACdD,IAAIJ;AAClB;AAEA,SAASM,EAAUR,GAAWC,GAAeC,GAAiBC,GAAiBC,GAAiBK,GAAgB;AAC5G,SAAOV,EAAKC,GAAGC,GAAOC,GAASC,GAAMC,CAAI,IAAIK;AACjD;AACA,SAASC,EAAUJ,GAAWC,GAAgBL,GAAiBC,GAAiBC,GAAiBO,GAAgB;AAC7G,SAAON,EAAKC,GAAGC,GAAQL,GAASC,GAAMC,CAAI,IAAIO;AAClD;AACA,SAASC,EACLZ,GACAM,GACAL,GACAM,GACAL,GACAC,GACAC,GACAS,GACAC,GACAC,IAAM,EAAC,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,KAChC;AACR,SAAAA,EAAI,OAAOP,EAAUR,GAAGC,GAAOC,GAASC,GAAMC,GAAMU,CAAO,GAC3DC,EAAI,OAAOhB,EAAKC,GAAGC,GAAOC,GAASC,GAAMC,CAAI,IAAIH,IAAQa,GACzDC,EAAI,OAAOL,EAAUJ,GAAGC,GAAQL,GAASC,GAAMC,GAAMS,CAAO,GAC5DE,EAAI,OAAOA,EAAI,OAAOR,IAAS,IAAIM,GAE5BE;AACX;AACA,SAASC,EAAUC,GAAiBhB,GAAeM,GAAgB;AAC/D,SAAOU,EAAK,QAAQ,KAAKA,EAAK,QAAQhB,KAASgB,EAAK,QAAQ,KAAKA,EAAK,QAAQV;AAClF;AAEO,SAASW,EAAaC,GAAc;AACvC,SAAOpB,EAAKoB,EAAM,GAAGA,EAAM,OAAOA,EAAM,SAASA,EAAM,WAAWA,EAAM,SAAS;AACrF;AACO,SAASC,EAAaD,GAAc;AACvC,SAAOd,EAAKc,EAAM,GAAGA,EAAM,QAAQA,EAAM,SAASA,EAAM,WAAWA,EAAM,SAAS;AACtF;AAEA,MAAME,IAA4C;AAAA,EAC9C,CAAC,UAAU,OAAO;AAAA,EAClB,CAAC,OAAO,OAAO;AAAA,EACf,CAAC,UAAU,OAAO;AAAA,EAClB,CAAC,UAAU,MAAM;AAAA,EACjB,CAAC,OAAO,MAAM;AAAA,EACd,CAAC,UAAU,MAAM;AAAA,EACjB,CAAC,OAAO,QAAQ;AAAA,EAChB,CAAC,UAAU,QAAQ;AACvB,GACMC,IAAe,GACfC,IAAsB,EAAC,MAAM,GAAG,MAAM,GAAG,MAAM,GAAG,MAAM,EAAA;AACvD,SAASC,EAAsBC,GAAgBC,GAAgBb,IAAUS,GAAcR,IAAUQ,GAAc;AAClH,QAAMK,IAAYJ,GACZK,IAAW,IAAIC,EAAA;AAErB,SAAO,CAAC7B,GAAWM,GAAWwB,GAAWC,GAAW7B,MAAmD;AACnG,eAAW,CAAC8B,GAAWC,CAAS,KAAKZ;AAGjC,UAFAT,EAAkBZ,GAAGM,GAAGwB,GAAGC,GAAG7B,GAAS+B,GAAWD,GAAWnB,GAASC,GAASa,CAAS,GAEpFX,EAASW,GAAWF,GAAQC,CAAM,KAAK,CAACE,EAAS,SAASD,CAAS;AACnE,eAAAC,EAAS,OAAO,EAAC,GAAGD,GAAU,GACvB,CAACM,GAAWD,CAAS;AAIpC,WAAO;AAAA,EACX;AACJ;"}

package/dist/scatterplot/utils/isNumericScale.d.ts

@@ -3,4 +3,3 @@ import { ChartScales } from '../components/types';
 import { AxisSettings, AxisSettingsContinuous } from '../../types';
 export declare function isNumericScale(scale: ChartScales['x' | 'y']): scale is Exclude<ChartScales['x' | 'y'], ScalePoint<string | number>>;
 export declare function isContinuousAxis(axis: AxisSettings): axis is AxisSettingsContinuous;
-//# sourceMappingURL=isNumericScale.d.ts.map

package/dist/scatterplot/utils/sortDotsByGrouping.d.ts

@@ -5,4 +5,3 @@ export declare function sortDotsByGrouping(dataFrame: DataFrame, dots: Dot[], gr
     columnName: ColumnName;
     order?: DataValue[];
 }[], highlight: ColumnName | null): void;
-//# sourceMappingURL=sortDotsByGrouping.d.ts.map

package/dist/scatterplot-umap/ChartRenderer.d.ts

@@ -19,14 +19,14 @@ declare class ChartRenderer {
     svgLayerComponent: ReactElement;
     aesGetters: AesGetters;
     tools: {
-        webglChart: CartesianChart<any, any>;
+        webglChart: CartesianChart<unknown, unknown>;
         qt?: KDBush;
         fillColor: () => {
             value: (setter: (d: Dot) => [number, number, number, number]) => void;
             data: (setter: () => Dot[]) => void;
         };
-        pointSeries: () => any;
-        zoom: ZoomBehavior<any, any>;
+        pointSeries: () => unknown;
+        zoom: ZoomBehavior<Element, unknown>;
     } | null;
     margins: Margins;
     chartSizes: ChartSizes;
@@ -69,4 +69,3 @@ declare class ChartRenderer {
     private updateLabelsPosition;
 }
 export default ChartRenderer;
-//# sourceMappingURL=ChartRenderer.d.ts.map

package/dist/scatterplot-umap/ChartRenderer.js

@@ -1,23 +1,19 @@
 var B = Object.defineProperty;
 var k = (S, t, e) => t in S ? B(S, t, { enumerable: !0, configurable: !0, writable: !0, value: e }) : S[t] = e;
 var h = (S, t, e) => k(S, typeof t != "symbol" ? t + "" : t, e);
-import { j as T } from "../node_modules/react/jsx-runtime.js";
-import K from "../node_modules/@d3fc/d3fc-chart/src/cartesian.js";
-import "../node_modules/@d3fc/d3fc-webgl/src/scale/scaleMapper.js";
-import X from "../node_modules/@d3fc/d3fc-webgl/src/style/fillColor.js";
-import Y from "../node_modules/@d3fc/d3fc-series/src/webgl/point.js";
-import j from "../node_modules/@d3fc/d3fc-series/src/svg/multi.js";
-import q from "../node_modules/@d3fc/d3fc-series/src/webgl/multi.js";
-import $ from "../node_modules/d3-zoom/src/zoom.js";
-import "../node_modules/d3-zoom/src/transform.js";
-import { createRoot as Q } from "../node_modules/react-dom/client.js";
-import { Error as J } from "../common/Error.js";
-import { BLACK as G, DEFAULT_HEIGHT as w, DEFAULT_WIDTH as D, TITLE_LINE_HEIGHT as Z, MAX_SHOWN_UNIQUE_VALUES_IN_LABELS as tt, TITLE_MARGIN as et } from "../constants.js";
-import { splitTextByWidth as it } from "../utils/splitTextByWidth.js";
+import { j as T } from "../_virtual/jsx-runtime.js";
+import K from "../node_modules/.pnpm/@d3fc_d3fc-chart@5.1.9_d3-array@3.2.4_d3-path@3.1.0_d3-scale-chromatic@3.1.0_d3-scale@4.0.2_d_kyhwsugamw24x7v3obvyiw5jfe/node_modules/@d3fc/d3fc-chart/src/cartesian.js";
+import { seriesWebglPoint as X, seriesSvgMulti as Y, seriesWebglMulti as j } from "@d3fc/d3fc-series";
+import q from "../node_modules/.pnpm/d3-zoom@3.0.0/node_modules/d3-zoom/src/zoom.js";
+import "../node_modules/.pnpm/d3-zoom@3.0.0/node_modules/d3-zoom/src/transform.js";
+import { c as $ } from "../_virtual/client.js";
+import { Error as Q } from "../common/Error.js";
+import { BLACK as G, DEFAULT_HEIGHT as w, DEFAULT_WIDTH as D, TITLE_LINE_HEIGHT as J, MAX_SHOWN_UNIQUE_VALUES_IN_LABELS as Z, TITLE_MARGIN as tt } from "../constants.js";
+import { splitTextByWidth as et } from "../utils/splitTextByWidth.js";
 import { MIN_MARGIN as L, DEFAULT_DOT_AES as F } from "../scatterplot/constants.js";
 import { createAesGetter as H } from "../scatterplot/utils/createAesGetter.js";
-import { getTicksAndFormat as st } from "../scatterplot/utils/getTicksAndFormat.js";
-import { createLabelPositioner as ot } from "../scatterplot/utils/getVisibleLabels.js";
+import { getTicksAndFormat as it } from "../scatterplot/utils/getTicksAndFormat.js";
+import { createLabelPositioner as st } from "../scatterplot/utils/getVisibleLabels.js";
 import { isContinuousAes as A } from "../types/common.js";
 import "../types/discrete.js";
 import "../types/scatterplot.js";
@@ -25,28 +21,30 @@ import "../types/heatmap.js";
 import "../types/dendro.js";
 import "../types/histogram.js";
 import "../types/bubble.js";
-import { getContinuousColorScale as at } from "../utils/getContinuousColorScale.js";
-import { getContinuousLegendTicks as nt } from "../utils/getContinuousLegendTicks.js";
+import { getContinuousColorScale as ot } from "../utils/getContinuousColorScale.js";
+import { getContinuousLegendTicks as at } from "../utils/getContinuousLegendTicks.js";
 import { TextMeasurer as R } from "../utils/TextMeasurer/TextMeasurer.js";
-import { SVGLayer as rt } from "./components/SVGLayer.js";
-import { PADDINGS as W, TITLE_LINE as I, LEGEND_OFFSET as lt } from "./constants.js";
-import ht from "../node_modules/@d3fc/d3fc-pointer/src/pointer.js";
-import ct from "../node_modules/kdbush/index.js";
+import { SVGLayer as nt } from "./components/SVGLayer.js";
+import { PADDINGS as E, TITLE_LINE as I, LEGEND_OFFSET as rt } from "./constants.js";
+import "../node_modules/.pnpm/@d3fc_d3fc-webgl@3.2.1_d3-scale@4.0.2_d3-shape@3.2.0/node_modules/@d3fc/d3fc-webgl/src/scale/scaleMapper.js";
+import lt from "../node_modules/.pnpm/@d3fc_d3fc-webgl@3.2.1_d3-scale@4.0.2_d3-shape@3.2.0/node_modules/@d3fc/d3fc-webgl/src/style/fillColor.js";
+import ht from "../node_modules/.pnpm/@d3fc_d3fc-pointer@3.0.3_d3-dispatch@3.0.1_d3-selection@3.0.0/node_modules/@d3fc/d3fc-pointer/src/pointer.js";
+import ct from "../node_modules/.pnpm/kdbush@4.0.2/node_modules/kdbush/index.js";
 import { DataFrameProvider as dt } from "../common/useDataFrame.js";
 import { DataFrame as mt } from "../DataFrame.js";
 import { arrangeLegendParts as pt } from "../utils/arrangeLegendParts.js";
 import { stringToNumberRgba as M, getColorWithFakeOpacity as P, GET_BLACK as O } from "./colors.js";
-import C from "../node_modules/d3-scale/src/linear.js";
-import U from "../node_modules/d3-scale/src/log.js";
-import N from "../node_modules/d3-selection/src/select.js";
-import gt from "../node_modules/d3-scale/src/ordinal.js";
-import ut from "../node_modules/d3-scale/src/symlog.js";
-import { sqrt as ft } from "../node_modules/d3-scale/src/pow.js";
-import St from "../node_modules/d3-polygon/src/contains.js";
+import C from "../node_modules/.pnpm/d3-scale@4.0.2/node_modules/d3-scale/src/linear.js";
+import U from "../node_modules/.pnpm/d3-scale@4.0.2/node_modules/d3-scale/src/log.js";
+import N from "../node_modules/.pnpm/d3-selection@3.0.0/node_modules/d3-selection/src/select.js";
+import gt from "../node_modules/.pnpm/d3-scale@4.0.2/node_modules/d3-scale/src/ordinal.js";
+import ut from "../node_modules/.pnpm/d3-scale@4.0.2/node_modules/d3-scale/src/symlog.js";
+import { sqrt as ft } from "../node_modules/.pnpm/d3-scale@4.0.2/node_modules/d3-scale/src/pow.js";
+import St from "../node_modules/.pnpm/d3-polygon@3.0.1/node_modules/d3-polygon/src/contains.js";
 function _(S) {
   return S / 2.35 * (S / 2.35) * 64;
 }
-function E(S, t, e) {
+function W(S, t, e) {
   return e.length === 0 ? !0 : e.some((i) => i.closed && i.points.length > 2 && St(i.points, [S, t]));
 }
 function V(S, t, e) {
@@ -57,7 +55,7 @@ function V(S, t, e) {
     lineType: H(S, t, (l == null ? void 0 : l.aes.lineShape) ?? "solid", "lineShape")
   };
 }
-class re {
+class ae {
   constructor() {
     h(this, "reactRoot", null);
     h(this, "parentNode", null);
@@ -116,7 +114,7 @@ class re {
     h(this, "dataFrame", new mt("", {}));
     h(this, "onPolygonUpdate", (t) => {
       this.polygons = t, this.updateByLasso();
-      const e = this.allDots.filter((i) => E(i.x, i.y, this.polygons)).map((i) => i.idx);
+      const e = this.allDots.filter((i) => W(i.x, i.y, this.polygons)).map((i) => i.idx);
       this.onPolygonUpdateOutside(e, t);
     });
   }
@@ -128,7 +126,7 @@ class re {
     });
   }
   init(t) {
-    this.parentNode === null && (this.parentNode = t, this.rootNode = document.createElement("div"), this.canvasNode = document.createElement("div"), this.parentNode.appendChild(this.rootNode), this.parentNode.appendChild(this.canvasNode), this.reactRoot = Q(this.rootNode));
+    this.parentNode === null && (this.parentNode = t, this.rootNode = document.createElement("div"), this.canvasNode = document.createElement("div"), this.parentNode.appendChild(this.rootNode), this.parentNode.appendChild(this.canvasNode), this.reactRoot = $.createRoot(this.rootNode));
   }
   updateChartSizes(t) {
     (t.width !== this.chartSizes.chartWidth || t.height !== this.chartSizes.chartHeight) && (this.chartSizes.chartWidth = t.width, this.chartSizes.chartHeight = t.height, this.scales.x.range([0, this.chartSizes.chartWidth]), this.scales.y.range([this.chartSizes.chartHeight, 0]));
@@ -150,9 +148,9 @@ class re {
         c = r - f, u = r + f;
       }
     }
-    const y = [0, this.chartSizes.chartWidth], x = l.copy().domain([o, m]).range([W.LEFT, this.chartSizes.chartWidth - W.RIGHT]);
+    const y = [0, this.chartSizes.chartWidth], x = l.copy().domain([o, m]).range([E.LEFT, this.chartSizes.chartWidth - E.RIGHT]);
     l.domain([x.invert(0), x.invert(this.chartSizes.chartWidth)]).range(y).nice(), this.scales.x = l, this.scales.xOriginal = l.copy();
-    const d = e.scale === "log" ? U() : C(), n = [this.chartSizes.chartHeight, 0], g = d.copy().domain([c, u]).range([this.chartSizes.chartHeight - W.BOTTOM, W.TOP]);
+    const d = e.scale === "log" ? U() : C(), n = [this.chartSizes.chartHeight, 0], g = d.copy().domain([c, u]).range([this.chartSizes.chartHeight - E.BOTTOM, E.TOP]);
     d.domain([g.invert(this.chartSizes.chartHeight), g.invert(0)]).range(n).nice(), this.scales.y = d, this.scales.yOriginal = d.copy();
   }
   updateCaptionsSize() {
@@ -160,21 +158,21 @@ class re {
     function e(c) {
       return Math.max(...c.map((m) => t.getTextWidth(m)));
     }
-    const { ticks: i, format: l } = st(this.scales.y, !1), o = e(i.map(l));
+    const { ticks: i, format: l } = it(this.scales.y, !1), o = e(i.map(l));
     this.captionsSizes = {
       xAxisCaptionsWidth: 20,
       yAxisCaptionsWidth: o
     };
   }
   createMainTitle(t) {
-    this.mainTitle = it(
+    this.mainTitle = et(
       t.name,
       this.chartSizes.totalWidth - this.margins.left - this.margins.right,
       20
     );
   }
   updateMargins() {
-    const t = Z * this.mainTitle.length, e = t > 0 ? t + et * 2 : 0;
+    const t = J * this.mainTitle.length, e = t > 0 ? t + tt * 2 : 0;
     this.margins = {
       top: Math.max(e, L),
       bottom: L + this.captionsSizes.xAxisCaptionsWidth + I,
@@ -192,7 +190,7 @@ class re {
       const n = e[d.value];
       if (n.usedAes.length === 0 || !n.aesMap)
         return;
-      if (n.values.length > tt) {
+      if (n.values.length > Z) {
         console.warn(`Too many values for discrete legend (${n.values.length})`);
         return;
       }
@@ -206,7 +204,7 @@ class re {
       o.push({ ...c, id: d.value, type: "dots", title: r, scale: f, values: n.values, labels: n.labels });
     }), l.forEach((d) => {
       if (d.type === "dots" && A(d.aes.dotFill)) {
-        const { domain: n, range: g, columnName: r, type: f = "linear" } = d.aes.dotFill, p = r.label ?? r.value, s = at(g, n, "linear"), a = (f === "log" ? ut() : C()).domain(n).range([this.chartSizes.chartHeight, 0]), z = nt(a, n);
+        const { domain: n, range: g, columnName: r, type: f = "linear" } = d.aes.dotFill, p = r.label ?? r.value, s = ot(g, n, "linear"), a = (f === "log" ? ut() : C()).domain(n).range([this.chartSizes.chartHeight, 0]), z = at(a, n);
         o.push({ ...c, id: "dotFill", type: "continuous", title: p, scale: s, tickPositionScale: a, values: z });
       }
       if (d.type === "dots" && A(d.aes.dotSize)) {
@@ -217,7 +215,7 @@ class re {
       this.legend = { width: 0, height: 0, items: [] };
       return;
     }
-    const m = pt(o, this.chartSizes.chartHeight), u = m[m.length - 1], y = u.left + u.width + lt, x = this.chartSizes.chartHeight;
+    const m = pt(o, this.chartSizes.chartHeight), u = m[m.length - 1], y = u.left + u.width + rt, x = this.chartSizes.chartHeight;
     this.legend = {
       width: y,
       height: x,
@@ -227,10 +225,10 @@ class re {
   initSettings(t, e, i, l, o, c) {
     const { xAxis: m, yAxis: u, size: y, title: x, legend: d } = e;
     this.updateChartSizes(y), this.updateViewport(m, u, i), this.updateCaptionsSize(), this.updateLegendSize(d, l, c, o), this.createMainTitle(x), this.updateMargins(), this.aesGetters = V(this.dataFrame, l, o);
-    const n = X().data(t).value((s) => {
+    const n = lt().data(t).value((s) => {
       const a = M(this.aesGetters.dotFill(s.idx));
-      return a[3] = !s.dimmed && E(s.x, s.y, this.polygons) ? 1 : 0.4, P(a, a);
-    }), g = Y().equals((s, a) => s === a).size((s) => _(this.aesGetters.dotSize(s.idx))).crossValue((s) => s.x).mainValue((s) => s.y).decorate((s) => n(s)), r = $().scaleExtent([0.1, 1e3]).filter((s) => s.type === "wheel" || s.metaKey || s.ctrlKey).on("start", () => {
+      return a[3] = !s.dimmed && W(s.x, s.y, this.polygons) ? 1 : 0.4, P(a, a);
+    }), g = X().equals((s, a) => s === a).size((s) => _(this.aesGetters.dotSize(s.idx))).crossValue((s) => s.x).mainValue((s) => s.y).decorate((s) => n(s)), r = q().scaleExtent([0.1, 1e3]).filter((s) => s.type === "wheel" || s.metaKey || s.ctrlKey).on("start", () => {
       this.selectedDot = null;
     }).on("zoom", (s) => {
       const a = s.transform, z = a.rescaleX(this.scales.xOriginal), v = a.rescaleY(this.scales.yOriginal);
@@ -246,9 +244,9 @@ class re {
     }), p = K({
       xScale: this.scales.x,
       yScale: this.scales.y
-    }).xAxisHeight(() => "0").yAxisWidth(() => "0").svgPlotArea(j()).webglPlotArea(
+    }).xAxisHeight(() => "0").yAxisWidth(() => "0").svgPlotArea(Y()).webglPlotArea(
       // only render the point series on the WebGL layer
-      q().series([g]).mapping((s) => s)
+      j().series([g]).mapping((s) => s)
     ).decorate((s) => {
       const a = s.enter();
       a && !this.svgLayerElement && (a.style("grid-template-columns", "0 auto 1fr auto 0"), a.style("grid-template-rows", "0 auto 1fr auto 0"), this.svgLayerElement = a.select(".svg-plot-area.plot-area svg").node()), a.select(".svg-plot-area.plot-area").attr("class", "svg-plot-area plot-area miplots-scatterplot-plot-area").call(r).call(f);
@@ -263,7 +261,7 @@ class re {
   updateAes(t, e) {
     this.tools !== null && (this.aesGetters = V(this.dataFrame, t, e), this.tools.fillColor.value((i) => {
       const l = this.aesGetters.dotFill ?? O, o = M(l(i.idx));
-      return o[3] = !i.dimmed && E(i.x, i.y, this.polygons) ? 1 : 0.4, P(o, o);
+      return o[3] = !i.dimmed && W(i.x, i.y, this.polygons) ? 1 : 0.4, P(o, o);
     }));
   }
   updatePointSize() {
@@ -273,7 +271,7 @@ class re {
     var t;
     this.tools !== null && ((t = this.tools) == null || t.fillColor.value((e) => {
       const i = this.aesGetters.dotFill ?? O, l = M(i(e.idx));
-      return l[3] = !e.dimmed && E(e.x, e.y, this.polygons) ? 1 : 0.4, P(l, l);
+      return l[3] = !e.dimmed && W(e.x, e.y, this.polygons) ? 1 : 0.4, P(l, l);
     }), this.renderWebglLayer());
   }
   createQuadtree(t) {
@@ -350,14 +348,14 @@ class re {
       onTooltipHintSwitch: this.onTooltipHintSwitch,
       onLassoControlsStateUpdate: this.onLassoControlsStateUpdate
     };
-    this.svgLayerComponent = /* @__PURE__ */ T.jsx(dt, { dataFrame: this.dataFrame, children: /* @__PURE__ */ T.jsx(rt, { ...t }) }), (e = this.reactRoot) == null || e.render(this.svgLayerComponent);
+    this.svgLayerComponent = /* @__PURE__ */ T.jsx(dt, { dataFrame: this.dataFrame, children: /* @__PURE__ */ T.jsx(nt, { ...t }) }), (e = this.reactRoot) == null || e.render(this.svgLayerComponent);
   }
   renderError(t) {
     var e;
-    (e = this.reactRoot) == null || e.render(/* @__PURE__ */ T.jsx(J, { message: t }));
+    (e = this.reactRoot) == null || e.render(/* @__PURE__ */ T.jsx(Q, { message: t }));
   }
   computeLabels(t) {
-    const e = new R("16px Manrope"), i = ot(this.chartSizes.chartWidth, this.chartSizes.chartHeight), l = [];
+    const e = new R("16px Manrope"), i = st(this.chartSizes.chartWidth, this.chartSizes.chartHeight), l = [];
     for (let o = 0; o < t.length; o++) {
       const c = t[o];
       if (c.label == null)
@@ -385,6 +383,6 @@ class re {
   }
 }
 export {
-  re as default
+  ae as default
 };
 //# sourceMappingURL=ChartRenderer.js.map