@voila.dev/chart 4.0.0

package/dist/core/marks/stats/hexbin.js

@@ -0,0 +1,113 @@
+import { channelLabel, isPlaceable, readChannel } from "../../channel.js";
+import { markId, valueOn } from "../shared.js";
+import { polygonPath } from "../../paths.js";
+const DEFAULT_RADIUS = 10;
+const SQRT3 = Math.sqrt(3);
+function hexCenter(x, y, radius) {
+  const dx = radius * SQRT3;
+  const dy = radius * 1.5;
+  const row = Math.round(y / dy);
+  const neighbourRow = y / dy > row ? row + 1 : row - 1;
+  const centres = [row, neighbourRow].map((candidate) => {
+    const shift = (candidate & 1) / 2;
+    return { x: (Math.round(x / dx - shift) + shift) * dx, y: candidate * dy };
+  });
+  return Math.hypot(centres[0].x - x, centres[0].y - y) <= Math.hypot(centres[1].x - x, centres[1].y - y) ? centres[0] : centres[1];
+}
+function hexagon(cx, cy, radius) {
+  return Array.from({ length: 6 }, (_unused, index) => {
+    const angle = Math.PI / 3 * index + Math.PI / 6;
+    return {
+      x: cx + radius * Math.cos(angle),
+      y: cy + radius * Math.sin(angle)
+    };
+  });
+}
+function hexbin(data, options) {
+  const xs = readChannel(data, options.x);
+  const ys = readChannel(data, options.y);
+  const reach = options.radius ?? DEFAULT_RADIUS;
+  return {
+    id: options.id,
+    coordinate: "cartesian",
+    focusOrder: "point",
+    channels: {
+      x: {
+        values: xs.filter(isPlaceable),
+        label: channelLabel(options.x),
+        inset: reach
+      },
+      y: {
+        values: ys.filter(isPlaceable),
+        label: channelLabel(options.y),
+        inset: reach
+      }
+    },
+    render(context) {
+      const id = markId(options.id, "hexbin", context);
+      const { x: xScale, y: yScale } = context.scales;
+      if (xScale === void 0 || yScale === void 0) {
+        return { nodes: [] };
+      }
+      const radius = options.radius ?? DEFAULT_RADIUS;
+      const bins = /* @__PURE__ */ new Map();
+      for (const [index, raw] of xs.entries()) {
+        const x = valueOn(xScale, raw);
+        const y = valueOn(yScale, ys[index]);
+        if (x === void 0 || y === void 0) continue;
+        const center = hexCenter(xScale.center(x), yScale.center(y), radius);
+        const key = `${Math.round(center.x)}:${Math.round(center.y)}`;
+        const bin = bins.get(key) ?? { ...center, members: [] };
+        bin.members.push(index);
+        bins.set(key, bin);
+      }
+      const fullest = Math.max(
+        1,
+        ...[...bins.values()].map((bin) => bin.members.length)
+      );
+      const color = options.fill ?? context.paletteColor(context.markIndex);
+      const nodes = [];
+      const points = [];
+      for (const [key, bin] of bins) {
+        const share = bin.members.length / fullest;
+        const fill = `color-mix(in oklab, ${color} ${Math.round(15 + 85 * share)}%, ${context.theme.background})`;
+        nodes.push({
+          kind: "path",
+          key: `${id}:${key}`,
+          role: "mark",
+          d: polygonPath(hexagon(bin.x, bin.y, radius - 0.5)),
+          paint: { fill }
+        });
+        points.push({
+          key: `${id}:${key}`,
+          markId: id,
+          index: bin.members[0],
+          datum: bin.members.map((index) => data[index]),
+          x: bin.x,
+          y: bin.y,
+          xValue: xScale.invert(bin.x),
+          yValue: yScale.invert(bin.y),
+          seriesLabel: options.label,
+          color: fill,
+          title: `${context.formatX(xScale.invert(bin.x))} \xB7 ${context.formatY(yScale.invert(bin.y))}`,
+          value: String(bin.members.length),
+          hit: {
+            kind: "rect",
+            rect: {
+              x: bin.x - radius,
+              y: bin.y - radius,
+              width: radius * 2,
+              height: radius * 2
+            }
+          }
+        });
+      }
+      return { nodes, points };
+    }
+  };
+}
+export {
+  hexCenter,
+  hexbin
+};
+//# sourceMappingURL=hexbin.js.map

package/dist/core/marks/stats/hexbin.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../../../../src/core/marks/stats/hexbin.ts"],
+  "sourcesContent": ["import { channelLabel, isPlaceable, readChannel } from \"#/core/channel.ts\";\nimport { markId, valueOn } from \"#/core/marks/shared.ts\";\nimport { polygonPath } from \"#/core/paths.ts\";\nimport type {\n\tChartAccessor,\n\tChartMark,\n\tChartPoint,\n\tChartValue,\n\tSceneNode,\n} from \"#/core/types.ts\";\n\nexport interface HexbinOptions<TDatum> {\n\treadonly id?: string;\n\treadonly x: ChartAccessor<TDatum, ChartValue>;\n\treadonly y: ChartAccessor<TDatum, ChartValue>;\n\t/** Pixel radius of a hexagon. */\n\treadonly radius?: number;\n\t/** Colour of the fullest bin; emptier bins fade towards the background. */\n\treadonly fill?: string;\n\treadonly label?: string;\n}\n\nconst DEFAULT_RADIUS = 10;\nconst SQRT3 = Math.sqrt(3);\n\n/**\n * The centre of the hexagon a pixel falls in, on a pointy-top grid whose odd\n * rows are offset by half a hexagon. Near a slanted edge the nearest row is\n * not always the right one, so the closer of the two candidate centres wins.\n */\nexport function hexCenter(\n\tx: number,\n\ty: number,\n\tradius: number,\n): { x: number; y: number } {\n\tconst dx = radius * SQRT3;\n\tconst dy = radius * 1.5;\n\tconst row = Math.round(y / dy);\n\tconst neighbourRow = y / dy > row ? row + 1 : row - 1;\n\tconst centres = [row, neighbourRow].map((candidate) => {\n\t\tconst shift = (candidate & 1) / 2;\n\t\treturn { x: (Math.round(x / dx - shift) + shift) * dx, y: candidate * dy };\n\t});\n\treturn Math.hypot(centres[0].x - x, centres[0].y - y) <=\n\t\tMath.hypot(centres[1].x - x, centres[1].y - y)\n\t\t? centres[0]\n\t\t: centres[1];\n}\n\nfunction hexagon(cx: number, cy: number, radius: number) {\n\treturn Array.from({ length: 6 }, (_unused, index) => {\n\t\tconst angle = (Math.PI / 3) * index + Math.PI / 6;\n\t\treturn {\n\t\t\tx: cx + radius * Math.cos(angle),\n\t\t\ty: cy + radius * Math.sin(angle),\n\t\t};\n\t});\n}\n\n/**\n * Points counted into hexagons: a scatter too dense to read as dots. The\n * bins are found in pixels, so they stay hexagonal whatever the scales.\n */\nexport function hexbin<TDatum>(\n\tdata: ReadonlyArray<TDatum>,\n\toptions: HexbinOptions<TDatum>,\n): ChartMark {\n\tconst xs = readChannel(data, options.x);\n\tconst ys = readChannel(data, options.y);\n\tconst reach = options.radius ?? DEFAULT_RADIUS;\n\treturn {\n\t\tid: options.id,\n\t\tcoordinate: \"cartesian\",\n\t\tfocusOrder: \"point\",\n\t\tchannels: {\n\t\t\tx: {\n\t\t\t\tvalues: xs.filter(isPlaceable),\n\t\t\t\tlabel: channelLabel(options.x),\n\t\t\t\tinset: reach,\n\t\t\t},\n\t\t\ty: {\n\t\t\t\tvalues: ys.filter(isPlaceable),\n\t\t\t\tlabel: channelLabel(options.y),\n\t\t\t\tinset: reach,\n\t\t\t},\n\t\t},\n\t\trender(context) {\n\t\t\tconst id = markId(options.id, \"hexbin\", context);\n\t\t\tconst { x: xScale, y: yScale } = context.scales;\n\t\t\tif (xScale === undefined || yScale === undefined) {\n\t\t\t\treturn { nodes: [] };\n\t\t\t}\n\t\t\tconst radius = options.radius ?? DEFAULT_RADIUS;\n\t\t\tconst bins = new Map<\n\t\t\t\tstring,\n\t\t\t\t{ x: number; y: number; members: number[] }\n\t\t\t>();\n\t\t\tfor (const [index, raw] of xs.entries()) {\n\t\t\t\tconst x = valueOn(xScale, raw);\n\t\t\t\tconst y = valueOn(yScale, ys[index]);\n\t\t\t\tif (x === undefined || y === undefined) continue;\n\t\t\t\tconst center = hexCenter(xScale.center(x), yScale.center(y), radius);\n\t\t\t\tconst key = `${Math.round(center.x)}:${Math.round(center.y)}`;\n\t\t\t\tconst bin = bins.get(key) ?? { ...center, members: [] };\n\t\t\t\tbin.members.push(index);\n\t\t\t\tbins.set(key, bin);\n\t\t\t}\n\t\t\tconst fullest = Math.max(\n\t\t\t\t1,\n\t\t\t\t...[...bins.values()].map((bin) => bin.members.length),\n\t\t\t);\n\t\t\tconst color = options.fill ?? context.paletteColor(context.markIndex);\n\t\t\tconst nodes: SceneNode[] = [];\n\t\t\tconst points: ChartPoint[] = [];\n\t\t\tfor (const [key, bin] of bins) {\n\t\t\t\tconst share = bin.members.length / fullest;\n\t\t\t\tconst fill = `color-mix(in oklab, ${color} ${Math.round(15 + 85 * share)}%, ${context.theme.background})`;\n\t\t\t\tnodes.push({\n\t\t\t\t\tkind: \"path\",\n\t\t\t\t\tkey: `${id}:${key}`,\n\t\t\t\t\trole: \"mark\",\n\t\t\t\t\td: polygonPath(hexagon(bin.x, bin.y, radius - 0.5)),\n\t\t\t\t\tpaint: { fill },\n\t\t\t\t});\n\t\t\t\tpoints.push({\n\t\t\t\t\tkey: `${id}:${key}`,\n\t\t\t\t\tmarkId: id,\n\t\t\t\t\tindex: bin.members[0],\n\t\t\t\t\tdatum: bin.members.map((index) => data[index]),\n\t\t\t\t\tx: bin.x,\n\t\t\t\t\ty: bin.y,\n\t\t\t\t\txValue: xScale.invert(bin.x),\n\t\t\t\t\tyValue: yScale.invert(bin.y),\n\t\t\t\t\tseriesLabel: options.label,\n\t\t\t\t\tcolor: fill,\n\t\t\t\t\ttitle: `${context.formatX(xScale.invert(bin.x))} \u00B7 ${context.formatY(yScale.invert(bin.y))}`,\n\t\t\t\t\tvalue: String(bin.members.length),\n\t\t\t\t\thit: {\n\t\t\t\t\t\tkind: \"rect\",\n\t\t\t\t\t\trect: {\n\t\t\t\t\t\t\tx: bin.x - radius,\n\t\t\t\t\t\t\ty: bin.y - radius,\n\t\t\t\t\t\t\twidth: radius * 2,\n\t\t\t\t\t\t\theight: radius * 2,\n\t\t\t\t\t\t},\n\t\t\t\t\t},\n\t\t\t\t});\n\t\t\t}\n\t\t\treturn { nodes, points };\n\t\t},\n\t};\n}\n"],
+  "mappings": "AAAA,SAAS,cAAc,aAAa,mBAAmB;AACvD,SAAS,QAAQ,eAAe;AAChC,SAAS,mBAAmB;AAoB5B,MAAM,iBAAiB;AACvB,MAAM,QAAQ,KAAK,KAAK,CAAC;AAOlB,SAAS,UACf,GACA,GACA,QAC2B;AAC3B,QAAM,KAAK,SAAS;AACpB,QAAM,KAAK,SAAS;AACpB,QAAM,MAAM,KAAK,MAAM,IAAI,EAAE;AAC7B,QAAM,eAAe,IAAI,KAAK,MAAM,MAAM,IAAI,MAAM;AACpD,QAAM,UAAU,CAAC,KAAK,YAAY,EAAE,IAAI,CAAC,cAAc;AACtD,UAAM,SAAS,YAAY,KAAK;AAChC,WAAO,EAAE,IAAI,KAAK,MAAM,IAAI,KAAK,KAAK,IAAI,SAAS,IAAI,GAAG,YAAY,GAAG;AAAA,EAC1E,CAAC;AACD,SAAO,KAAK,MAAM,QAAQ,CAAC,EAAE,IAAI,GAAG,QAAQ,CAAC,EAAE,IAAI,CAAC,KACnD,KAAK,MAAM,QAAQ,CAAC,EAAE,IAAI,GAAG,QAAQ,CAAC,EAAE,IAAI,CAAC,IAC3C,QAAQ,CAAC,IACT,QAAQ,CAAC;AACb;AAEA,SAAS,QAAQ,IAAY,IAAY,QAAgB;AACxD,SAAO,MAAM,KAAK,EAAE,QAAQ,EAAE,GAAG,CAAC,SAAS,UAAU;AACpD,UAAM,QAAS,KAAK,KAAK,IAAK,QAAQ,KAAK,KAAK;AAChD,WAAO;AAAA,MACN,GAAG,KAAK,SAAS,KAAK,IAAI,KAAK;AAAA,MAC/B,GAAG,KAAK,SAAS,KAAK,IAAI,KAAK;AAAA,IAChC;AAAA,EACD,CAAC;AACF;AAMO,SAAS,OACf,MACA,SACY;AACZ,QAAM,KAAK,YAAY,MAAM,QAAQ,CAAC;AACtC,QAAM,KAAK,YAAY,MAAM,QAAQ,CAAC;AACtC,QAAM,QAAQ,QAAQ,UAAU;AAChC,SAAO;AAAA,IACN,IAAI,QAAQ;AAAA,IACZ,YAAY;AAAA,IACZ,YAAY;AAAA,IACZ,UAAU;AAAA,MACT,GAAG;AAAA,QACF,QAAQ,GAAG,OAAO,WAAW;AAAA,QAC7B,OAAO,aAAa,QAAQ,CAAC;AAAA,QAC7B,OAAO;AAAA,MACR;AAAA,MACA,GAAG;AAAA,QACF,QAAQ,GAAG,OAAO,WAAW;AAAA,QAC7B,OAAO,aAAa,QAAQ,CAAC;AAAA,QAC7B,OAAO;AAAA,MACR;AAAA,IACD;AAAA,IACA,OAAO,SAAS;AACf,YAAM,KAAK,OAAO,QAAQ,IAAI,UAAU,OAAO;AAC/C,YAAM,EAAE,GAAG,QAAQ,GAAG,OAAO,IAAI,QAAQ;AACzC,UAAI,WAAW,UAAa,WAAW,QAAW;AACjD,eAAO,EAAE,OAAO,CAAC,EAAE;AAAA,MACpB;AACA,YAAM,SAAS,QAAQ,UAAU;AACjC,YAAM,OAAO,oBAAI,IAGf;AACF,iBAAW,CAAC,OAAO,GAAG,KAAK,GAAG,QAAQ,GAAG;AACxC,cAAM,IAAI,QAAQ,QAAQ,GAAG;AAC7B,cAAM,IAAI,QAAQ,QAAQ,GAAG,KAAK,CAAC;AACnC,YAAI,MAAM,UAAa,MAAM,OAAW;AACxC,cAAM,SAAS,UAAU,OAAO,OAAO,CAAC,GAAG,OAAO,OAAO,CAAC,GAAG,MAAM;AACnE,cAAM,MAAM,GAAG,KAAK,MAAM,OAAO,CAAC,CAAC,IAAI,KAAK,MAAM,OAAO,CAAC,CAAC;AAC3D,cAAM,MAAM,KAAK,IAAI,GAAG,KAAK,EAAE,GAAG,QAAQ,SAAS,CAAC,EAAE;AACtD,YAAI,QAAQ,KAAK,KAAK;AACtB,aAAK,IAAI,KAAK,GAAG;AAAA,MAClB;AACA,YAAM,UAAU,KAAK;AAAA,QACpB;AAAA,QACA,GAAG,CAAC,GAAG,KAAK,OAAO,CAAC,EAAE,IAAI,CAAC,QAAQ,IAAI,QAAQ,MAAM;AAAA,MACtD;AACA,YAAM,QAAQ,QAAQ,QAAQ,QAAQ,aAAa,QAAQ,SAAS;AACpE,YAAM,QAAqB,CAAC;AAC5B,YAAM,SAAuB,CAAC;AAC9B,iBAAW,CAAC,KAAK,GAAG,KAAK,MAAM;AAC9B,cAAM,QAAQ,IAAI,QAAQ,SAAS;AACnC,cAAM,OAAO,uBAAuB,KAAK,IAAI,KAAK,MAAM,KAAK,KAAK,KAAK,CAAC,MAAM,QAAQ,MAAM,UAAU;AACtG,cAAM,KAAK;AAAA,UACV,MAAM;AAAA,UACN,KAAK,GAAG,EAAE,IAAI,GAAG;AAAA,UACjB,MAAM;AAAA,UACN,GAAG,YAAY,QAAQ,IAAI,GAAG,IAAI,GAAG,SAAS,GAAG,CAAC;AAAA,UAClD,OAAO,EAAE,KAAK;AAAA,QACf,CAAC;AACD,eAAO,KAAK;AAAA,UACX,KAAK,GAAG,EAAE,IAAI,GAAG;AAAA,UACjB,QAAQ;AAAA,UACR,OAAO,IAAI,QAAQ,CAAC;AAAA,UACpB,OAAO,IAAI,QAAQ,IAAI,CAAC,UAAU,KAAK,KAAK,CAAC;AAAA,UAC7C,GAAG,IAAI;AAAA,UACP,GAAG,IAAI;AAAA,UACP,QAAQ,OAAO,OAAO,IAAI,CAAC;AAAA,UAC3B,QAAQ,OAAO,OAAO,IAAI,CAAC;AAAA,UAC3B,aAAa,QAAQ;AAAA,UACrB,OAAO;AAAA,UACP,OAAO,GAAG,QAAQ,QAAQ,OAAO,OAAO,IAAI,CAAC,CAAC,CAAC,SAAM,QAAQ,QAAQ,OAAO,OAAO,IAAI,CAAC,CAAC,CAAC;AAAA,UAC1F,OAAO,OAAO,IAAI,QAAQ,MAAM;AAAA,UAChC,KAAK;AAAA,YACJ,MAAM;AAAA,YACN,MAAM;AAAA,cACL,GAAG,IAAI,IAAI;AAAA,cACX,GAAG,IAAI,IAAI;AAAA,cACX,OAAO,SAAS;AAAA,cAChB,QAAQ,SAAS;AAAA,YAClB;AAAA,UACD;AAAA,QACD,CAAC;AAAA,MACF;AACA,aAAO,EAAE,OAAO,OAAO;AAAA,IACxB;AAAA,EACD;AACD;",
+  "names": []
+}

package/dist/core/marks/stats/histogram.d.ts

@@ -0,0 +1,27 @@
+import type { ChartMarkOptions } from "../shared.js";
+import type { ChartAccessor, ChartMark, ChartValue } from "../../types.js";
+export interface HistogramOptions<TDatum> extends ChartMarkOptions<TDatum> {
+    readonly x: ChartAccessor<TDatum, ChartValue>;
+    /** About how many bins: the edges land on round numbers, so the count is approximate. */
+    readonly bins?: number;
+    readonly fill?: string;
+    /** Each bin's share of the total instead of its count. */
+    readonly normalize?: boolean;
+}
+export interface Bin {
+    readonly x1: number;
+    readonly x2: number;
+    readonly count: number;
+    readonly series?: ChartValue;
+}
+/** Round-edged bins over the values' extent, every value in exactly one. */
+export declare function binValues(values: ReadonlyArray<number>, target: number): ReadonlyArray<{
+    x1: number;
+    x2: number;
+}>;
+/**
+ * Counts per bin, as touching rectangles over a continuous x. With a color
+ * channel the series stack inside each bin.
+ */
+export declare function histogram<TDatum>(data: ReadonlyArray<TDatum>, options: HistogramOptions<TDatum>): ChartMark;
+//# sourceMappingURL=histogram.d.ts.map

package/dist/core/marks/stats/histogram.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"histogram.d.ts","sourceRoot":"","sources":["../../../../src/core/marks/stats/histogram.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,gBAAgB,EAAE,MAAM,wBAAwB,CAAC;AAE/D,OAAO,KAAK,EAAE,aAAa,EAAE,SAAS,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAE5E,MAAM,WAAW,gBAAgB,CAAC,MAAM,CAAE,SAAQ,gBAAgB,CAAC,MAAM,CAAC;IACzE,QAAQ,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAC9C,yFAAyF;IACzF,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB,0DAA0D;IAC1D,QAAQ,CAAC,SAAS,CAAC,EAAE,OAAO,CAAC;CAC7B;AAED,MAAM,WAAW,GAAG;IACnB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,MAAM,CAAC,EAAE,UAAU,CAAC;CAC7B;AAED,4EAA4E;AAC5E,wBAAgB,SAAS,CACxB,MAAM,EAAE,aAAa,CAAC,MAAM,CAAC,EAC7B,MAAM,EAAE,MAAM,GACZ,aAAa,CAAC;IAAE,EAAE,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,MAAM,CAAA;CAAE,CAAC,CAiB3C;AAID;;;GAGG;AACH,wBAAgB,SAAS,CAAC,MAAM,EAC/B,IAAI,EAAE,aAAa,CAAC,MAAM,CAAC,EAC3B,OAAO,EAAE,gBAAgB,CAAC,MAAM,CAAC,GAC/B,SAAS,CAqDX"}

package/dist/core/marks/stats/histogram.js

@@ -0,0 +1,70 @@
+import { isPlaceable, numericValue, readChannel } from "../../channel.js";
+import { rectY } from "../rect.js";
+import { tickStep } from "../../ticks.js";
+function binValues(values, target) {
+  const finite = values.filter(Number.isFinite);
+  if (finite.length === 0) {
+    return [];
+  }
+  const low = Math.min(...finite);
+  const high = Math.max(...finite);
+  const step = tickStep(low, high, target) || 1;
+  const first = Math.floor(low / step) * step;
+  const count = Math.max(
+    1,
+    Math.ceil((high - first) / step + (high === first ? 1 : 0))
+  );
+  return Array.from({ length: count }, (_unused, index) => ({
+    x1: first + index * step,
+    x2: first + (index + 1) * step
+  }));
+}
+const DEFAULT_BINS = 20;
+function histogram(data, options) {
+  const xs = readChannel(data, options.x).map(numericValue);
+  const colors = readChannel(data, options.color);
+  const edges = binValues(
+    xs.filter((value) => value !== void 0),
+    options.bins ?? DEFAULT_BINS
+  );
+  const seriesValues = options.color === void 0 ? [void 0] : [
+    ...new Map(
+      colors.filter(isPlaceable).map((value) => [String(value), value])
+    ).values()
+  ];
+  const total = xs.filter((value) => value !== void 0).length || 1;
+  const bins = [];
+  for (const series of seriesValues) {
+    for (const [binIndex, edge] of edges.entries()) {
+      const last = binIndex === edges.length - 1;
+      let count = 0;
+      for (const [index, x] of xs.entries()) {
+        const inSeries = series === void 0 || String(colors[index]) === String(series);
+        if (inSeries && x !== void 0 && x >= edge.x1 && (x < edge.x2 || last && x <= edge.x2)) {
+          count += 1;
+        }
+      }
+      bins.push({
+        ...edge,
+        count: options.normalize ? count / total : count,
+        series
+      });
+    }
+  }
+  return rectY(bins, {
+    id: options.id,
+    label: options.label,
+    tip: options.tip,
+    opacity: options.opacity,
+    fill: options.fill,
+    x1: (bin) => bin.x1,
+    x2: (bin) => bin.x2,
+    y: (bin) => bin.count,
+    ...options.color === void 0 ? {} : { color: (bin) => bin.series }
+  });
+}
+export {
+  binValues,
+  histogram
+};
+//# sourceMappingURL=histogram.js.map

package/dist/core/marks/stats/histogram.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../../../../src/core/marks/stats/histogram.ts"],
+  "sourcesContent": ["import { isPlaceable, numericValue, readChannel } from \"#/core/channel.ts\";\nimport { rectY } from \"#/core/marks/rect.ts\";\nimport type { ChartMarkOptions } from \"#/core/marks/shared.ts\";\nimport { tickStep } from \"#/core/ticks.ts\";\nimport type { ChartAccessor, ChartMark, ChartValue } from \"#/core/types.ts\";\n\nexport interface HistogramOptions<TDatum> extends ChartMarkOptions<TDatum> {\n\treadonly x: ChartAccessor<TDatum, ChartValue>;\n\t/** About how many bins: the edges land on round numbers, so the count is approximate. */\n\treadonly bins?: number;\n\treadonly fill?: string;\n\t/** Each bin's share of the total instead of its count. */\n\treadonly normalize?: boolean;\n}\n\nexport interface Bin {\n\treadonly x1: number;\n\treadonly x2: number;\n\treadonly count: number;\n\treadonly series?: ChartValue;\n}\n\n/** Round-edged bins over the values' extent, every value in exactly one. */\nexport function binValues(\n\tvalues: ReadonlyArray<number>,\n\ttarget: number,\n): ReadonlyArray<{ x1: number; x2: number }> {\n\tconst finite = values.filter(Number.isFinite);\n\tif (finite.length === 0) {\n\t\treturn [];\n\t}\n\tconst low = Math.min(...finite);\n\tconst high = Math.max(...finite);\n\tconst step = tickStep(low, high, target) || 1;\n\tconst first = Math.floor(low / step) * step;\n\tconst count = Math.max(\n\t\t1,\n\t\tMath.ceil((high - first) / step + (high === first ? 1 : 0)),\n\t);\n\treturn Array.from({ length: count }, (_unused, index) => ({\n\t\tx1: first + index * step,\n\t\tx2: first + (index + 1) * step,\n\t}));\n}\n\nconst DEFAULT_BINS = 20;\n\n/**\n * Counts per bin, as touching rectangles over a continuous x. With a color\n * channel the series stack inside each bin.\n */\nexport function histogram<TDatum>(\n\tdata: ReadonlyArray<TDatum>,\n\toptions: HistogramOptions<TDatum>,\n): ChartMark {\n\tconst xs = readChannel(data, options.x).map(numericValue);\n\tconst colors = readChannel(data, options.color);\n\tconst edges = binValues(\n\t\txs.filter((value): value is number => value !== undefined),\n\t\toptions.bins ?? DEFAULT_BINS,\n\t);\n\tconst seriesValues =\n\t\toptions.color === undefined\n\t\t\t? [undefined]\n\t\t\t: [\n\t\t\t\t\t...new Map(\n\t\t\t\t\t\tcolors.filter(isPlaceable).map((value) => [String(value), value]),\n\t\t\t\t\t).values(),\n\t\t\t\t];\n\tconst total = xs.filter((value) => value !== undefined).length || 1;\n\tconst bins: Bin[] = [];\n\tfor (const series of seriesValues) {\n\t\tfor (const [binIndex, edge] of edges.entries()) {\n\t\t\tconst last = binIndex === edges.length - 1;\n\t\t\tlet count = 0;\n\t\t\tfor (const [index, x] of xs.entries()) {\n\t\t\t\tconst inSeries =\n\t\t\t\t\tseries === undefined || String(colors[index]) === String(series);\n\t\t\t\tif (\n\t\t\t\t\tinSeries &&\n\t\t\t\t\tx !== undefined &&\n\t\t\t\t\tx >= edge.x1 &&\n\t\t\t\t\t(x < edge.x2 || (last && x <= edge.x2))\n\t\t\t\t) {\n\t\t\t\t\tcount += 1;\n\t\t\t\t}\n\t\t\t}\n\t\t\tbins.push({\n\t\t\t\t...edge,\n\t\t\t\tcount: options.normalize ? count / total : count,\n\t\t\t\tseries,\n\t\t\t});\n\t\t}\n\t}\n\treturn rectY(bins, {\n\t\tid: options.id,\n\t\tlabel: options.label,\n\t\ttip: options.tip,\n\t\topacity: options.opacity,\n\t\tfill: options.fill,\n\t\tx1: (bin) => bin.x1,\n\t\tx2: (bin) => bin.x2,\n\t\ty: (bin) => bin.count,\n\t\t...(options.color === undefined\n\t\t\t? {}\n\t\t\t: { color: (bin: Bin) => bin.series as ChartValue }),\n\t});\n}\n"],
+  "mappings": "AAAA,SAAS,aAAa,cAAc,mBAAmB;AACvD,SAAS,aAAa;AAEtB,SAAS,gBAAgB;AAoBlB,SAAS,UACf,QACA,QAC4C;AAC5C,QAAM,SAAS,OAAO,OAAO,OAAO,QAAQ;AAC5C,MAAI,OAAO,WAAW,GAAG;AACxB,WAAO,CAAC;AAAA,EACT;AACA,QAAM,MAAM,KAAK,IAAI,GAAG,MAAM;AAC9B,QAAM,OAAO,KAAK,IAAI,GAAG,MAAM;AAC/B,QAAM,OAAO,SAAS,KAAK,MAAM,MAAM,KAAK;AAC5C,QAAM,QAAQ,KAAK,MAAM,MAAM,IAAI,IAAI;AACvC,QAAM,QAAQ,KAAK;AAAA,IAClB;AAAA,IACA,KAAK,MAAM,OAAO,SAAS,QAAQ,SAAS,QAAQ,IAAI,EAAE;AAAA,EAC3D;AACA,SAAO,MAAM,KAAK,EAAE,QAAQ,MAAM,GAAG,CAAC,SAAS,WAAW;AAAA,IACzD,IAAI,QAAQ,QAAQ;AAAA,IACpB,IAAI,SAAS,QAAQ,KAAK;AAAA,EAC3B,EAAE;AACH;AAEA,MAAM,eAAe;AAMd,SAAS,UACf,MACA,SACY;AACZ,QAAM,KAAK,YAAY,MAAM,QAAQ,CAAC,EAAE,IAAI,YAAY;AACxD,QAAM,SAAS,YAAY,MAAM,QAAQ,KAAK;AAC9C,QAAM,QAAQ;AAAA,IACb,GAAG,OAAO,CAAC,UAA2B,UAAU,MAAS;AAAA,IACzD,QAAQ,QAAQ;AAAA,EACjB;AACA,QAAM,eACL,QAAQ,UAAU,SACf,CAAC,MAAS,IACV;AAAA,IACA,GAAG,IAAI;AAAA,MACN,OAAO,OAAO,WAAW,EAAE,IAAI,CAAC,UAAU,CAAC,OAAO,KAAK,GAAG,KAAK,CAAC;AAAA,IACjE,EAAE,OAAO;AAAA,EACV;AACH,QAAM,QAAQ,GAAG,OAAO,CAAC,UAAU,UAAU,MAAS,EAAE,UAAU;AAClE,QAAM,OAAc,CAAC;AACrB,aAAW,UAAU,cAAc;AAClC,eAAW,CAAC,UAAU,IAAI,KAAK,MAAM,QAAQ,GAAG;AAC/C,YAAM,OAAO,aAAa,MAAM,SAAS;AACzC,UAAI,QAAQ;AACZ,iBAAW,CAAC,OAAO,CAAC,KAAK,GAAG,QAAQ,GAAG;AACtC,cAAM,WACL,WAAW,UAAa,OAAO,OAAO,KAAK,CAAC,MAAM,OAAO,MAAM;AAChE,YACC,YACA,MAAM,UACN,KAAK,KAAK,OACT,IAAI,KAAK,MAAO,QAAQ,KAAK,KAAK,KAClC;AACD,mBAAS;AAAA,QACV;AAAA,MACD;AACA,WAAK,KAAK;AAAA,QACT,GAAG;AAAA,QACH,OAAO,QAAQ,YAAY,QAAQ,QAAQ;AAAA,QAC3C;AAAA,MACD,CAAC;AAAA,IACF;AAAA,EACD;AACA,SAAO,MAAM,MAAM;AAAA,IAClB,IAAI,QAAQ;AAAA,IACZ,OAAO,QAAQ;AAAA,IACf,KAAK,QAAQ;AAAA,IACb,SAAS,QAAQ;AAAA,IACjB,MAAM,QAAQ;AAAA,IACd,IAAI,CAAC,QAAQ,IAAI;AAAA,IACjB,IAAI,CAAC,QAAQ,IAAI;AAAA,IACjB,GAAG,CAAC,QAAQ,IAAI;AAAA,IAChB,GAAI,QAAQ,UAAU,SACnB,CAAC,IACD,EAAE,OAAO,CAAC,QAAa,IAAI,OAAqB;AAAA,EACpD,CAAC;AACF;",
+  "names": []
+}

package/dist/core/marks/stats/regression.d.ts

@@ -0,0 +1,13 @@
+import type { ChartAccessor, ChartMark, ChartValue } from "../../types.js";
+export interface RegressionOptions<TDatum> {
+    readonly id?: string;
+    readonly x: ChartAccessor<TDatum, ChartValue>;
+    readonly y: ChartAccessor<TDatum, ChartValue>;
+    readonly stroke?: string;
+    /** Shade the 95 % confidence band of the fitted mean. */
+    readonly band?: boolean;
+    readonly label?: string;
+}
+/** The least-squares line through x and y, and how sure it is. Drawn over a scatter. */
+export declare function regressionY<TDatum>(data: ReadonlyArray<TDatum>, options: RegressionOptions<TDatum>): ChartMark;
+//# sourceMappingURL=regression.d.ts.map

package/dist/core/marks/stats/regression.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"regression.d.ts","sourceRoot":"","sources":["../../../../src/core/marks/stats/regression.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EACX,aAAa,EACb,SAAS,EACT,UAAU,EAEV,MAAM,iBAAiB,CAAC;AAEzB,MAAM,WAAW,iBAAiB,CAAC,MAAM;IACxC,QAAQ,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAC9C,QAAQ,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAC9C,QAAQ,CAAC,MAAM,CAAC,EAAE,MAAM,CAAC;IACzB,yDAAyD;IACzD,QAAQ,CAAC,IAAI,CAAC,EAAE,OAAO,CAAC;IACxB,QAAQ,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC;CACxB;AAID,wFAAwF;AACxF,wBAAgB,WAAW,CAAC,MAAM,EACjC,IAAI,EAAE,aAAa,CAAC,MAAM,CAAC,EAC3B,OAAO,EAAE,iBAAiB,CAAC,MAAM,CAAC,GAChC,SAAS,CA+EX"}

package/dist/core/marks/stats/regression.js

@@ -0,0 +1,80 @@
+import { isPlaceable, numericValue, readChannel } from "../../channel.js";
+import { markId, paint } from "../shared.js";
+import { linearFit } from "./statistics.js";
+import { areaPath, linePath } from "../../paths.js";
+import { toNumber } from "../../scales/continuous.js";
+const BAND_SAMPLES = 24;
+function regressionY(data, options) {
+  const pairs = readChannel(data, options.x).map((x, index) => ({
+    x,
+    y: numericValue(readChannel(data, options.y)[index])
+  })).filter(
+    (pair) => isPlaceable(pair.x) && pair.y !== void 0
+  );
+  const xs = pairs.map((pair) => toNumber(pair.x));
+  const fit = linearFit(
+    xs,
+    pairs.map((pair) => pair.y)
+  );
+  const dates = pairs.some((pair) => pair.x instanceof Date);
+  return {
+    id: options.id,
+    coordinate: "cartesian",
+    channels: {
+      x: { values: pairs.map((pair) => pair.x) },
+      y: { values: pairs.map((pair) => pair.y) }
+    },
+    render(context) {
+      const id = markId(options.id, "regressionY", context);
+      const { x: xScale, y: yScale } = context.scales;
+      if (fit === void 0 || xScale === void 0 || yScale === void 0 || xs.length === 0) {
+        return { nodes: [] };
+      }
+      const low = Math.min(...xs);
+      const high = Math.max(...xs);
+      const at = Array.from(
+        { length: BAND_SAMPLES },
+        (_unused, index) => low + (high - low) * index / (BAND_SAMPLES - 1)
+      );
+      const place = (x) => xScale.map(dates ? new Date(x) : x);
+      const stroke = options.stroke ?? context.theme.foreground;
+      const nodes = [];
+      if (options.band !== false) {
+        nodes.push({
+          kind: "path",
+          key: `${id}:band`,
+          role: "mark",
+          d: areaPath(
+            at.map((x) => ({
+              x: place(x),
+              y: yScale.map(fit.intercept + fit.slope * x + fit.margin(x))
+            })),
+            at.map((x) => ({
+              x: place(x),
+              y: yScale.map(fit.intercept + fit.slope * x - fit.margin(x))
+            }))
+          ),
+          paint: { fill: stroke, fillOpacity: 0.12 }
+        });
+      }
+      nodes.push({
+        kind: "path",
+        key: `${id}:line`,
+        role: "mark",
+        d: linePath([
+          { x: place(low), y: yScale.map(fit.intercept + fit.slope * low) },
+          { x: place(high), y: yScale.map(fit.intercept + fit.slope * high) }
+        ]),
+        paint: paint({ fill: "none", stroke, strokeWidth: 2 })
+      });
+      return {
+        nodes,
+        legend: options.label ? [{ key: id, label: options.label, color: stroke, shape: "line" }] : []
+      };
+    }
+  };
+}
+export {
+  regressionY
+};
+//# sourceMappingURL=regression.js.map

package/dist/core/marks/stats/regression.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../../../../src/core/marks/stats/regression.ts"],
+  "sourcesContent": ["import { isPlaceable, numericValue, readChannel } from \"#/core/channel.ts\";\nimport { markId, paint } from \"#/core/marks/shared.ts\";\nimport { linearFit } from \"#/core/marks/stats/statistics.ts\";\nimport { areaPath, linePath } from \"#/core/paths.ts\";\nimport { toNumber } from \"#/core/scales/continuous.ts\";\nimport type {\n\tChartAccessor,\n\tChartMark,\n\tChartValue,\n\tSceneNode,\n} from \"#/core/types.ts\";\n\nexport interface RegressionOptions<TDatum> {\n\treadonly id?: string;\n\treadonly x: ChartAccessor<TDatum, ChartValue>;\n\treadonly y: ChartAccessor<TDatum, ChartValue>;\n\treadonly stroke?: string;\n\t/** Shade the 95 % confidence band of the fitted mean. */\n\treadonly band?: boolean;\n\treadonly label?: string;\n}\n\nconst BAND_SAMPLES = 24;\n\n/** The least-squares line through x and y, and how sure it is. Drawn over a scatter. */\nexport function regressionY<TDatum>(\n\tdata: ReadonlyArray<TDatum>,\n\toptions: RegressionOptions<TDatum>,\n): ChartMark {\n\tconst pairs = readChannel(data, options.x)\n\t\t.map((x, index) => ({\n\t\t\tx,\n\t\t\ty: numericValue(readChannel(data, options.y)[index]),\n\t\t}))\n\t\t.filter(\n\t\t\t(pair): pair is { x: ChartValue; y: number } =>\n\t\t\t\tisPlaceable(pair.x) && pair.y !== undefined,\n\t\t);\n\tconst xs = pairs.map((pair) => toNumber(pair.x));\n\tconst fit = linearFit(\n\t\txs,\n\t\tpairs.map((pair) => pair.y),\n\t);\n\tconst dates = pairs.some((pair) => pair.x instanceof Date);\n\treturn {\n\t\tid: options.id,\n\t\tcoordinate: \"cartesian\",\n\t\tchannels: {\n\t\t\tx: { values: pairs.map((pair) => pair.x) },\n\t\t\ty: { values: pairs.map((pair) => pair.y) },\n\t\t},\n\t\trender(context) {\n\t\t\tconst id = markId(options.id, \"regressionY\", context);\n\t\t\tconst { x: xScale, y: yScale } = context.scales;\n\t\t\tif (\n\t\t\t\tfit === undefined ||\n\t\t\t\txScale === undefined ||\n\t\t\t\tyScale === undefined ||\n\t\t\t\txs.length === 0\n\t\t\t) {\n\t\t\t\treturn { nodes: [] };\n\t\t\t}\n\t\t\tconst low = Math.min(...xs);\n\t\t\tconst high = Math.max(...xs);\n\t\t\tconst at = Array.from(\n\t\t\t\t{ length: BAND_SAMPLES },\n\t\t\t\t(_unused, index) => low + ((high - low) * index) / (BAND_SAMPLES - 1),\n\t\t\t);\n\t\t\tconst place = (x: number) => xScale.map(dates ? new Date(x) : x);\n\t\t\tconst stroke = options.stroke ?? context.theme.foreground;\n\t\t\tconst nodes: SceneNode[] = [];\n\t\t\tif (options.band !== false) {\n\t\t\t\tnodes.push({\n\t\t\t\t\tkind: \"path\",\n\t\t\t\t\tkey: `${id}:band`,\n\t\t\t\t\trole: \"mark\",\n\t\t\t\t\td: areaPath(\n\t\t\t\t\t\tat.map((x) => ({\n\t\t\t\t\t\t\tx: place(x),\n\t\t\t\t\t\t\ty: yScale.map(fit.intercept + fit.slope * x + fit.margin(x)),\n\t\t\t\t\t\t})),\n\t\t\t\t\t\tat.map((x) => ({\n\t\t\t\t\t\t\tx: place(x),\n\t\t\t\t\t\t\ty: yScale.map(fit.intercept + fit.slope * x - fit.margin(x)),\n\t\t\t\t\t\t})),\n\t\t\t\t\t),\n\t\t\t\t\tpaint: { fill: stroke, fillOpacity: 0.12 },\n\t\t\t\t});\n\t\t\t}\n\t\t\tnodes.push({\n\t\t\t\tkind: \"path\",\n\t\t\t\tkey: `${id}:line`,\n\t\t\t\trole: \"mark\",\n\t\t\t\td: linePath([\n\t\t\t\t\t{ x: place(low), y: yScale.map(fit.intercept + fit.slope * low) },\n\t\t\t\t\t{ x: place(high), y: yScale.map(fit.intercept + fit.slope * high) },\n\t\t\t\t]),\n\t\t\t\tpaint: paint({ fill: \"none\", stroke, strokeWidth: 2 }),\n\t\t\t});\n\t\t\treturn {\n\t\t\t\tnodes,\n\t\t\t\tlegend: options.label\n\t\t\t\t\t? [{ key: id, label: options.label, color: stroke, shape: \"line\" }]\n\t\t\t\t\t: [],\n\t\t\t};\n\t\t},\n\t};\n}\n"],
+  "mappings": "AAAA,SAAS,aAAa,cAAc,mBAAmB;AACvD,SAAS,QAAQ,aAAa;AAC9B,SAAS,iBAAiB;AAC1B,SAAS,UAAU,gBAAgB;AACnC,SAAS,gBAAgB;AAkBzB,MAAM,eAAe;AAGd,SAAS,YACf,MACA,SACY;AACZ,QAAM,QAAQ,YAAY,MAAM,QAAQ,CAAC,EACvC,IAAI,CAAC,GAAG,WAAW;AAAA,IACnB;AAAA,IACA,GAAG,aAAa,YAAY,MAAM,QAAQ,CAAC,EAAE,KAAK,CAAC;AAAA,EACpD,EAAE,EACD;AAAA,IACA,CAAC,SACA,YAAY,KAAK,CAAC,KAAK,KAAK,MAAM;AAAA,EACpC;AACD,QAAM,KAAK,MAAM,IAAI,CAAC,SAAS,SAAS,KAAK,CAAC,CAAC;AAC/C,QAAM,MAAM;AAAA,IACX;AAAA,IACA,MAAM,IAAI,CAAC,SAAS,KAAK,CAAC;AAAA,EAC3B;AACA,QAAM,QAAQ,MAAM,KAAK,CAAC,SAAS,KAAK,aAAa,IAAI;AACzD,SAAO;AAAA,IACN,IAAI,QAAQ;AAAA,IACZ,YAAY;AAAA,IACZ,UAAU;AAAA,MACT,GAAG,EAAE,QAAQ,MAAM,IAAI,CAAC,SAAS,KAAK,CAAC,EAAE;AAAA,MACzC,GAAG,EAAE,QAAQ,MAAM,IAAI,CAAC,SAAS,KAAK,CAAC,EAAE;AAAA,IAC1C;AAAA,IACA,OAAO,SAAS;AACf,YAAM,KAAK,OAAO,QAAQ,IAAI,eAAe,OAAO;AACpD,YAAM,EAAE,GAAG,QAAQ,GAAG,OAAO,IAAI,QAAQ;AACzC,UACC,QAAQ,UACR,WAAW,UACX,WAAW,UACX,GAAG,WAAW,GACb;AACD,eAAO,EAAE,OAAO,CAAC,EAAE;AAAA,MACpB;AACA,YAAM,MAAM,KAAK,IAAI,GAAG,EAAE;AAC1B,YAAM,OAAO,KAAK,IAAI,GAAG,EAAE;AAC3B,YAAM,KAAK,MAAM;AAAA,QAChB,EAAE,QAAQ,aAAa;AAAA,QACvB,CAAC,SAAS,UAAU,OAAQ,OAAO,OAAO,SAAU,eAAe;AAAA,MACpE;AACA,YAAM,QAAQ,CAAC,MAAc,OAAO,IAAI,QAAQ,IAAI,KAAK,CAAC,IAAI,CAAC;AAC/D,YAAM,SAAS,QAAQ,UAAU,QAAQ,MAAM;AAC/C,YAAM,QAAqB,CAAC;AAC5B,UAAI,QAAQ,SAAS,OAAO;AAC3B,cAAM,KAAK;AAAA,UACV,MAAM;AAAA,UACN,KAAK,GAAG,EAAE;AAAA,UACV,MAAM;AAAA,UACN,GAAG;AAAA,YACF,GAAG,IAAI,CAAC,OAAO;AAAA,cACd,GAAG,MAAM,CAAC;AAAA,cACV,GAAG,OAAO,IAAI,IAAI,YAAY,IAAI,QAAQ,IAAI,IAAI,OAAO,CAAC,CAAC;AAAA,YAC5D,EAAE;AAAA,YACF,GAAG,IAAI,CAAC,OAAO;AAAA,cACd,GAAG,MAAM,CAAC;AAAA,cACV,GAAG,OAAO,IAAI,IAAI,YAAY,IAAI,QAAQ,IAAI,IAAI,OAAO,CAAC,CAAC;AAAA,YAC5D,EAAE;AAAA,UACH;AAAA,UACA,OAAO,EAAE,MAAM,QAAQ,aAAa,KAAK;AAAA,QAC1C,CAAC;AAAA,MACF;AACA,YAAM,KAAK;AAAA,QACV,MAAM;AAAA,QACN,KAAK,GAAG,EAAE;AAAA,QACV,MAAM;AAAA,QACN,GAAG,SAAS;AAAA,UACX,EAAE,GAAG,MAAM,GAAG,GAAG,GAAG,OAAO,IAAI,IAAI,YAAY,IAAI,QAAQ,GAAG,EAAE;AAAA,UAChE,EAAE,GAAG,MAAM,IAAI,GAAG,GAAG,OAAO,IAAI,IAAI,YAAY,IAAI,QAAQ,IAAI,EAAE;AAAA,QACnE,CAAC;AAAA,QACD,OAAO,MAAM,EAAE,MAAM,QAAQ,QAAQ,aAAa,EAAE,CAAC;AAAA,MACtD,CAAC;AACD,aAAO;AAAA,QACN;AAAA,QACA,QAAQ,QAAQ,QACb,CAAC,EAAE,KAAK,IAAI,OAAO,QAAQ,OAAO,OAAO,QAAQ,OAAO,OAAO,CAAC,IAChE,CAAC;AAAA,MACL;AAAA,IACD;AAAA,EACD;AACD;",
+  "names": []
+}

package/dist/core/marks/stats/statistics.d.ts

@@ -0,0 +1,34 @@
+/**
+ * The few statistics the statistical marks need, written out: quantiles,
+ * kernel density, least squares. Each is a handful of lines, which is less to
+ * carry than a dependency.
+ */
+/** Linear-interpolated quantile of sorted values (the R-7 / Excel definition). */
+export declare function quantile(sorted: ReadonlyArray<number>, p: number): number;
+export interface BoxSummary {
+    readonly min: number;
+    readonly q1: number;
+    readonly median: number;
+    readonly q3: number;
+    readonly max: number;
+    /** Whisker ends: the furthest values within 1.5 IQR of the box. */
+    readonly low: number;
+    readonly high: number;
+    readonly outliers: ReadonlyArray<number>;
+    readonly count: number;
+}
+export declare function boxSummary(values: ReadonlyArray<number>): BoxSummary | undefined;
+/** Silverman's rule of thumb: a bandwidth that suits most unimodal data. */
+export declare function silvermanBandwidth(values: ReadonlyArray<number>): number;
+/** Gaussian kernel density of `values`, sampled at `at`. */
+export declare function kernelDensity(values: ReadonlyArray<number>, at: ReadonlyArray<number>, bandwidth?: number): ReadonlyArray<number>;
+/** `count` evenly spaced samples over `[low, high]`. */
+export declare function samples(low: number, high: number, count: number): number[];
+export interface LinearFit {
+    readonly slope: number;
+    readonly intercept: number;
+    /** Half-width of the 95 % confidence band of the mean at `x`. */
+    readonly margin: (x: number) => number;
+}
+export declare function linearFit(xs: ReadonlyArray<number>, ys: ReadonlyArray<number>): LinearFit | undefined;
+//# sourceMappingURL=statistics.d.ts.map

package/dist/core/marks/stats/statistics.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"statistics.d.ts","sourceRoot":"","sources":["../../../../src/core/marks/stats/statistics.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,kFAAkF;AAClF,wBAAgB,QAAQ,CAAC,MAAM,EAAE,aAAa,CAAC,MAAM,CAAC,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAQzE;AAED,MAAM,WAAW,UAAU;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,mEAAmE;IACnE,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,QAAQ,EAAE,aAAa,CAAC,MAAM,CAAC,CAAC;IACzC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;CACvB;AAID,wBAAgB,UAAU,CACzB,MAAM,EAAE,aAAa,CAAC,MAAM,CAAC,GAC3B,UAAU,GAAG,SAAS,CA0BxB;AAED,4EAA4E;AAC5E,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,MAAM,CAexE;AAID,4DAA4D;AAC5D,wBAAgB,aAAa,CAC5B,MAAM,EAAE,aAAa,CAAC,MAAM,CAAC,EAC7B,EAAE,EAAE,aAAa,CAAC,MAAM,CAAC,EACzB,SAAS,SAA6B,GACpC,aAAa,CAAC,MAAM,CAAC,CAUvB;AAED,wDAAwD;AACxD,wBAAgB,OAAO,CAAC,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,EAAE,CAQ1E;AAED,MAAM,WAAW,SAAS;IACzB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,iEAAiE;IACjE,QAAQ,CAAC,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,CAAC;CACvC;AAKD,wBAAgB,SAAS,CACxB,EAAE,EAAE,aAAa,CAAC,MAAM,CAAC,EACzB,EAAE,EAAE,aAAa,CAAC,MAAM,CAAC,GACvB,SAAS,GAAG,SAAS,CA6BvB"}

package/dist/core/marks/stats/statistics.js

@@ -0,0 +1,111 @@
+function quantile(sorted, p) {
+  if (sorted.length === 0) {
+    return Number.NaN;
+  }
+  const position = (sorted.length - 1) * p;
+  const below = Math.floor(position);
+  const above = Math.ceil(position);
+  return sorted[below] + (sorted[above] - sorted[below]) * (position - below);
+}
+const TUKEY_FENCE = 1.5;
+function boxSummary(values) {
+  const sorted = values.filter(Number.isFinite).sort((left, right) => left - right);
+  if (sorted.length === 0) {
+    return void 0;
+  }
+  const q1 = quantile(sorted, 0.25);
+  const q3 = quantile(sorted, 0.75);
+  const fence = (q3 - q1) * TUKEY_FENCE;
+  const inside = sorted.filter(
+    (value) => value >= q1 - fence && value <= q3 + fence
+  );
+  return {
+    min: sorted[0],
+    q1,
+    median: quantile(sorted, 0.5),
+    q3,
+    max: sorted[sorted.length - 1],
+    low: inside[0] ?? q1,
+    high: inside[inside.length - 1] ?? q3,
+    outliers: sorted.filter(
+      (value) => value < q1 - fence || value > q3 + fence
+    ),
+    count: sorted.length
+  };
+}
+function silvermanBandwidth(values) {
+  const count = values.length;
+  if (count < 2) {
+    return 1;
+  }
+  const mean = values.reduce((sum, value) => sum + value, 0) / count;
+  const deviation = Math.sqrt(
+    values.reduce((sum, value) => sum + (value - mean) ** 2, 0) / (count - 1)
+  );
+  const sorted = [...values].sort((left, right) => left - right);
+  const spread = Math.min(
+    deviation,
+    (quantile(sorted, 0.75) - quantile(sorted, 0.25)) / 1.34
+  );
+  return 0.9 * (spread || deviation || 1) * count ** -0.2;
+}
+const GAUSSIAN = 1 / Math.sqrt(2 * Math.PI);
+function kernelDensity(values, at, bandwidth = silvermanBandwidth(values)) {
+  const count = values.length;
+  return at.map((x) => {
+    let sum = 0;
+    for (const value of values) {
+      const u = (x - value) / bandwidth;
+      sum += GAUSSIAN * Math.exp(-0.5 * u * u);
+    }
+    return count === 0 ? 0 : sum / (count * bandwidth);
+  });
+}
+function samples(low, high, count) {
+  if (count < 2 || low === high) {
+    return [low];
+  }
+  return Array.from(
+    { length: count },
+    (_unused, index) => low + (high - low) * index / (count - 1)
+  );
+}
+const Z_95 = 1.96;
+function linearFit(xs, ys) {
+  const count = xs.length;
+  if (count < 2) {
+    return void 0;
+  }
+  const meanX = xs.reduce((sum, x) => sum + x, 0) / count;
+  const meanY = ys.reduce((sum, y) => sum + y, 0) / count;
+  let sxx = 0;
+  let sxy = 0;
+  for (let index = 0; index < count; index += 1) {
+    sxx += (xs[index] - meanX) ** 2;
+    sxy += (xs[index] - meanX) * (ys[index] - meanY);
+  }
+  if (sxx === 0) {
+    return void 0;
+  }
+  const slope = sxy / sxx;
+  const intercept = meanY - slope * meanX;
+  let residuals = 0;
+  for (let index = 0; index < count; index += 1) {
+    residuals += (ys[index] - (intercept + slope * xs[index])) ** 2;
+  }
+  const standardError = count > 2 ? Math.sqrt(residuals / (count - 2)) : 0;
+  return {
+    slope,
+    intercept,
+    margin: (x) => Z_95 * standardError * Math.sqrt(1 / count + (x - meanX) ** 2 / sxx)
+  };
+}
+export {
+  boxSummary,
+  kernelDensity,
+  linearFit,
+  quantile,
+  samples,
+  silvermanBandwidth
+};
+//# sourceMappingURL=statistics.js.map

package/dist/core/marks/stats/statistics.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../../../../src/core/marks/stats/statistics.ts"],
+  "sourcesContent": ["/**\n * The few statistics the statistical marks need, written out: quantiles,\n * kernel density, least squares. Each is a handful of lines, which is less to\n * carry than a dependency.\n */\n\n/** Linear-interpolated quantile of sorted values (the R-7 / Excel definition). */\nexport function quantile(sorted: ReadonlyArray<number>, p: number): number {\n\tif (sorted.length === 0) {\n\t\treturn Number.NaN;\n\t}\n\tconst position = (sorted.length - 1) * p;\n\tconst below = Math.floor(position);\n\tconst above = Math.ceil(position);\n\treturn sorted[below] + (sorted[above] - sorted[below]) * (position - below);\n}\n\nexport interface BoxSummary {\n\treadonly min: number;\n\treadonly q1: number;\n\treadonly median: number;\n\treadonly q3: number;\n\treadonly max: number;\n\t/** Whisker ends: the furthest values within 1.5 IQR of the box. */\n\treadonly low: number;\n\treadonly high: number;\n\treadonly outliers: ReadonlyArray<number>;\n\treadonly count: number;\n}\n\nconst TUKEY_FENCE = 1.5;\n\nexport function boxSummary(\n\tvalues: ReadonlyArray<number>,\n): BoxSummary | undefined {\n\tconst sorted = values\n\t\t.filter(Number.isFinite)\n\t\t.sort((left, right) => left - right);\n\tif (sorted.length === 0) {\n\t\treturn undefined;\n\t}\n\tconst q1 = quantile(sorted, 0.25);\n\tconst q3 = quantile(sorted, 0.75);\n\tconst fence = (q3 - q1) * TUKEY_FENCE;\n\tconst inside = sorted.filter(\n\t\t(value) => value >= q1 - fence && value <= q3 + fence,\n\t);\n\treturn {\n\t\tmin: sorted[0],\n\t\tq1,\n\t\tmedian: quantile(sorted, 0.5),\n\t\tq3,\n\t\tmax: sorted[sorted.length - 1],\n\t\tlow: inside[0] ?? q1,\n\t\thigh: inside[inside.length - 1] ?? q3,\n\t\toutliers: sorted.filter(\n\t\t\t(value) => value < q1 - fence || value > q3 + fence,\n\t\t),\n\t\tcount: sorted.length,\n\t};\n}\n\n/** Silverman's rule of thumb: a bandwidth that suits most unimodal data. */\nexport function silvermanBandwidth(values: ReadonlyArray<number>): number {\n\tconst count = values.length;\n\tif (count < 2) {\n\t\treturn 1;\n\t}\n\tconst mean = values.reduce((sum, value) => sum + value, 0) / count;\n\tconst deviation = Math.sqrt(\n\t\tvalues.reduce((sum, value) => sum + (value - mean) ** 2, 0) / (count - 1),\n\t);\n\tconst sorted = [...values].sort((left, right) => left - right);\n\tconst spread = Math.min(\n\t\tdeviation,\n\t\t(quantile(sorted, 0.75) - quantile(sorted, 0.25)) / 1.34,\n\t);\n\treturn 0.9 * (spread || deviation || 1) * count ** -0.2;\n}\n\nconst GAUSSIAN = 1 / Math.sqrt(2 * Math.PI);\n\n/** Gaussian kernel density of `values`, sampled at `at`. */\nexport function kernelDensity(\n\tvalues: ReadonlyArray<number>,\n\tat: ReadonlyArray<number>,\n\tbandwidth = silvermanBandwidth(values),\n): ReadonlyArray<number> {\n\tconst count = values.length;\n\treturn at.map((x) => {\n\t\tlet sum = 0;\n\t\tfor (const value of values) {\n\t\t\tconst u = (x - value) / bandwidth;\n\t\t\tsum += GAUSSIAN * Math.exp(-0.5 * u * u);\n\t\t}\n\t\treturn count === 0 ? 0 : sum / (count * bandwidth);\n\t});\n}\n\n/** `count` evenly spaced samples over `[low, high]`. */\nexport function samples(low: number, high: number, count: number): number[] {\n\tif (count < 2 || low === high) {\n\t\treturn [low];\n\t}\n\treturn Array.from(\n\t\t{ length: count },\n\t\t(_unused, index) => low + ((high - low) * index) / (count - 1),\n\t);\n}\n\nexport interface LinearFit {\n\treadonly slope: number;\n\treadonly intercept: number;\n\t/** Half-width of the 95 % confidence band of the mean at `x`. */\n\treadonly margin: (x: number) => number;\n}\n\n/** 1.96: a normal approximation of Student's t, fine past a dozen points. */\nconst Z_95 = 1.96;\n\nexport function linearFit(\n\txs: ReadonlyArray<number>,\n\tys: ReadonlyArray<number>,\n): LinearFit | undefined {\n\tconst count = xs.length;\n\tif (count < 2) {\n\t\treturn undefined;\n\t}\n\tconst meanX = xs.reduce((sum, x) => sum + x, 0) / count;\n\tconst meanY = ys.reduce((sum, y) => sum + y, 0) / count;\n\tlet sxx = 0;\n\tlet sxy = 0;\n\tfor (let index = 0; index < count; index += 1) {\n\t\tsxx += (xs[index] - meanX) ** 2;\n\t\tsxy += (xs[index] - meanX) * (ys[index] - meanY);\n\t}\n\tif (sxx === 0) {\n\t\treturn undefined;\n\t}\n\tconst slope = sxy / sxx;\n\tconst intercept = meanY - slope * meanX;\n\tlet residuals = 0;\n\tfor (let index = 0; index < count; index += 1) {\n\t\tresiduals += (ys[index] - (intercept + slope * xs[index])) ** 2;\n\t}\n\tconst standardError = count > 2 ? Math.sqrt(residuals / (count - 2)) : 0;\n\treturn {\n\t\tslope,\n\t\tintercept,\n\t\tmargin: (x) =>\n\t\t\tZ_95 * standardError * Math.sqrt(1 / count + (x - meanX) ** 2 / sxx),\n\t};\n}\n"],
+  "mappings": "AAOO,SAAS,SAAS,QAA+B,GAAmB;AAC1E,MAAI,OAAO,WAAW,GAAG;AACxB,WAAO,OAAO;AAAA,EACf;AACA,QAAM,YAAY,OAAO,SAAS,KAAK;AACvC,QAAM,QAAQ,KAAK,MAAM,QAAQ;AACjC,QAAM,QAAQ,KAAK,KAAK,QAAQ;AAChC,SAAO,OAAO,KAAK,KAAK,OAAO,KAAK,IAAI,OAAO,KAAK,MAAM,WAAW;AACtE;AAeA,MAAM,cAAc;AAEb,SAAS,WACf,QACyB;AACzB,QAAM,SAAS,OACb,OAAO,OAAO,QAAQ,EACtB,KAAK,CAAC,MAAM,UAAU,OAAO,KAAK;AACpC,MAAI,OAAO,WAAW,GAAG;AACxB,WAAO;AAAA,EACR;AACA,QAAM,KAAK,SAAS,QAAQ,IAAI;AAChC,QAAM,KAAK,SAAS,QAAQ,IAAI;AAChC,QAAM,SAAS,KAAK,MAAM;AAC1B,QAAM,SAAS,OAAO;AAAA,IACrB,CAAC,UAAU,SAAS,KAAK,SAAS,SAAS,KAAK;AAAA,EACjD;AACA,SAAO;AAAA,IACN,KAAK,OAAO,CAAC;AAAA,IACb;AAAA,IACA,QAAQ,SAAS,QAAQ,GAAG;AAAA,IAC5B;AAAA,IACA,KAAK,OAAO,OAAO,SAAS,CAAC;AAAA,IAC7B,KAAK,OAAO,CAAC,KAAK;AAAA,IAClB,MAAM,OAAO,OAAO,SAAS,CAAC,KAAK;AAAA,IACnC,UAAU,OAAO;AAAA,MAChB,CAAC,UAAU,QAAQ,KAAK,SAAS,QAAQ,KAAK;AAAA,IAC/C;AAAA,IACA,OAAO,OAAO;AAAA,EACf;AACD;AAGO,SAAS,mBAAmB,QAAuC;AACzE,QAAM,QAAQ,OAAO;AACrB,MAAI,QAAQ,GAAG;AACd,WAAO;AAAA,EACR;AACA,QAAM,OAAO,OAAO,OAAO,CAAC,KAAK,UAAU,MAAM,OAAO,CAAC,IAAI;AAC7D,QAAM,YAAY,KAAK;AAAA,IACtB,OAAO,OAAO,CAAC,KAAK,UAAU,OAAO,QAAQ,SAAS,GAAG,CAAC,KAAK,QAAQ;AAAA,EACxE;AACA,QAAM,SAAS,CAAC,GAAG,MAAM,EAAE,KAAK,CAAC,MAAM,UAAU,OAAO,KAAK;AAC7D,QAAM,SAAS,KAAK;AAAA,IACnB;AAAA,KACC,SAAS,QAAQ,IAAI,IAAI,SAAS,QAAQ,IAAI,KAAK;AAAA,EACrD;AACA,SAAO,OAAO,UAAU,aAAa,KAAK,SAAS;AACpD;AAEA,MAAM,WAAW,IAAI,KAAK,KAAK,IAAI,KAAK,EAAE;AAGnC,SAAS,cACf,QACA,IACA,YAAY,mBAAmB,MAAM,GACb;AACxB,QAAM,QAAQ,OAAO;AACrB,SAAO,GAAG,IAAI,CAAC,MAAM;AACpB,QAAI,MAAM;AACV,eAAW,SAAS,QAAQ;AAC3B,YAAM,KAAK,IAAI,SAAS;AACxB,aAAO,WAAW,KAAK,IAAI,OAAO,IAAI,CAAC;AAAA,IACxC;AACA,WAAO,UAAU,IAAI,IAAI,OAAO,QAAQ;AAAA,EACzC,CAAC;AACF;AAGO,SAAS,QAAQ,KAAa,MAAc,OAAyB;AAC3E,MAAI,QAAQ,KAAK,QAAQ,MAAM;AAC9B,WAAO,CAAC,GAAG;AAAA,EACZ;AACA,SAAO,MAAM;AAAA,IACZ,EAAE,QAAQ,MAAM;AAAA,IAChB,CAAC,SAAS,UAAU,OAAQ,OAAO,OAAO,SAAU,QAAQ;AAAA,EAC7D;AACD;AAUA,MAAM,OAAO;AAEN,SAAS,UACf,IACA,IACwB;AACxB,QAAM,QAAQ,GAAG;AACjB,MAAI,QAAQ,GAAG;AACd,WAAO;AAAA,EACR;AACA,QAAM,QAAQ,GAAG,OAAO,CAAC,KAAK,MAAM,MAAM,GAAG,CAAC,IAAI;AAClD,QAAM,QAAQ,GAAG,OAAO,CAAC,KAAK,MAAM,MAAM,GAAG,CAAC,IAAI;AAClD,MAAI,MAAM;AACV,MAAI,MAAM;AACV,WAAS,QAAQ,GAAG,QAAQ,OAAO,SAAS,GAAG;AAC9C,YAAQ,GAAG,KAAK,IAAI,UAAU;AAC9B,YAAQ,GAAG,KAAK,IAAI,UAAU,GAAG,KAAK,IAAI;AAAA,EAC3C;AACA,MAAI,QAAQ,GAAG;AACd,WAAO;AAAA,EACR;AACA,QAAM,QAAQ,MAAM;AACpB,QAAM,YAAY,QAAQ,QAAQ;AAClC,MAAI,YAAY;AAChB,WAAS,QAAQ,GAAG,QAAQ,OAAO,SAAS,GAAG;AAC9C,kBAAc,GAAG,KAAK,KAAK,YAAY,QAAQ,GAAG,KAAK,OAAO;AAAA,EAC/D;AACA,QAAM,gBAAgB,QAAQ,IAAI,KAAK,KAAK,aAAa,QAAQ,EAAE,IAAI;AACvE,SAAO;AAAA,IACN;AAAA,IACA;AAAA,IACA,QAAQ,CAAC,MACR,OAAO,gBAAgB,KAAK,KAAK,IAAI,SAAS,IAAI,UAAU,IAAI,GAAG;AAAA,EACrE;AACD;",
+  "names": []
+}

package/dist/core/marks/stats/waffle.d.ts

@@ -0,0 +1,18 @@
+import { type ChartMarkOptions } from "../shared.js";
+import type { ChartAccessor, ChartMark, ChartValue } from "../../types.js";
+export interface WaffleOptions<TDatum> extends ChartMarkOptions<TDatum> {
+    readonly x: ChartAccessor<TDatum, ChartValue>;
+    readonly y: ChartAccessor<TDatum, ChartValue>;
+    /** What one cell stands for. */
+    readonly unit?: number;
+    /** Cells across a band. */
+    readonly columns?: number;
+    readonly fill?: string;
+}
+/**
+ * Bars made of countable cells: "12 missions" reads as twelve squares. A
+ * row of `columns` cells is `columns × unit` on the y axis, so the waffle
+ * still reads against it; series continue one another inside a column.
+ */
+export declare function waffleY<TDatum>(data: ReadonlyArray<TDatum>, options: WaffleOptions<TDatum>): ChartMark;
+//# sourceMappingURL=waffle.d.ts.map

package/dist/core/marks/stats/waffle.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"waffle.d.ts","sourceRoot":"","sources":["../../../../src/core/marks/stats/waffle.ts"],"names":[],"mappings":"AAMA,OAAO,EACN,KAAK,gBAAgB,EAKrB,MAAM,wBAAwB,CAAC;AAEhC,OAAO,KAAK,EACX,aAAa,EACb,SAAS,EAET,UAAU,EAEV,MAAM,iBAAiB,CAAC;AAEzB,MAAM,WAAW,aAAa,CAAC,MAAM,CAAE,SAAQ,gBAAgB,CAAC,MAAM,CAAC;IACtE,QAAQ,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAC9C,QAAQ,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAC9C,gCAAgC;IAChC,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB,2BAA2B;IAC3B,QAAQ,CAAC,OAAO,CAAC,EAAE,MAAM,CAAC;IAC1B,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;CACvB;AAKD;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,MAAM,EAC7B,IAAI,EAAE,aAAa,CAAC,MAAM,CAAC,EAC3B,OAAO,EAAE,aAAa,CAAC,MAAM,CAAC,GAC5B,SAAS,CAiGX"}