prographs 0.12.0 → 0.13.0

package/dist/{chunk-7FG4JSVO.cjs → chunk-2GL2VNJU.cjs}

@@ -1,4 +1,4 @@
-"use strict";Object.defineProperty(exports, "__esModule", {value: true}); function _nullishCoalesce(lhs, rhsFn) { if (lhs != null) { return lhs; } else { return rhsFn(); } } function _optionalChain(ops) { let lastAccessLHS = undefined; let value = ops[0]; let i = 1; while (i < ops.length) { const op = ops[i]; const fn = ops[i + 1]; i += 2; if ((op === 'optionalAccess' || op === 'optionalCall') && value == null) { return undefined; } if (op === 'access' || op === 'optionalAccess') { lastAccessLHS = value; value = fn(value); } else if (op === 'call' || op === 'optionalCall') { value = fn((...args) => value.call(lastAccessLHS, ...args)); lastAccessLHS = undefined; } } return value; } var _class; var _class2; var _class3; var _class4; var _class5; var _class6; var _class7; var _class8; var _class9; var _class10;
+"use strict";Object.defineProperty(exports, "__esModule", {value: true}); function _nullishCoalesce(lhs, rhsFn) { if (lhs != null) { return lhs; } else { return rhsFn(); } } function _optionalChain(ops) { let lastAccessLHS = undefined; let value = ops[0]; let i = 1; while (i < ops.length) { const op = ops[i]; const fn = ops[i + 1]; i += 2; if ((op === 'optionalAccess' || op === 'optionalCall') && value == null) { return undefined; } if (op === 'access' || op === 'optionalAccess') { lastAccessLHS = value; value = fn(value); } else if (op === 'call' || op === 'optionalCall') { value = fn((...args) => value.call(lastAccessLHS, ...args)); lastAccessLHS = undefined; } } return value; } var _class; var _class2; var _class3; var _class4; var _class5; var _class6; var _class7; var _class8; var _class9; var _class10; var _class11;
 
 var _chunkXX2YTCORcjs = require('./chunk-XX2YTCOR.cjs');
 
@@ -6,7 +6,7 @@ var _chunkXX2YTCORcjs = require('./chunk-XX2YTCOR.cjs');
 var _chunkEEF5WRRZcjs = require('./chunk-EEF5WRRZ.cjs');
 
 
-var _chunkUMEKIIHJcjs = require('./chunk-UMEKIIHJ.cjs');
+var _chunkLNWTC72Ycjs = require('./chunk-LNWTC72Y.cjs');
 
 // src/charts/fieldcharts/fieldEnums.ts
 var FieldStyle = /* @__PURE__ */ ((FieldStyle2) => {
@@ -51,6 +51,12 @@ var RarefactionBasis = /* @__PURE__ */ ((RarefactionBasis2) => {
   RarefactionBasis2["Samples"] = "samples";
   return RarefactionBasis2;
 })(RarefactionBasis || {});
+var AccumulationOrder = /* @__PURE__ */ ((AccumulationOrder2) => {
+  AccumulationOrder2["Collector"] = "collector";
+  AccumulationOrder2["Randomised"] = "randomised";
+  AccumulationOrder2["Both"] = "collector-and-randomised";
+  return AccumulationOrder2;
+})(AccumulationOrder || {});
 var KiteScaling = /* @__PURE__ */ ((KiteScaling2) => {
   KiteScaling2["Shared"] = "shared";
   KiteScaling2["Own"] = "own";
@@ -151,19 +157,19 @@ var FieldKit = class _FieldKit {
   }
   /** Area-weighted centroid, falling back to the vertex mean for degenerate rings. */
   static centroid(polygon) {
-    let area3 = 0, cx = 0, cy = 0;
+    let area4 = 0, cx = 0, cy = 0;
     polygon.forEach(([x, y], i) => {
       const [nx, ny] = polygon[(i + 1) % polygon.length];
       const cross = x * ny - nx * y;
-      area3 += cross;
+      area4 += cross;
       cx += (x + nx) * cross;
       cy += (y + ny) * cross;
     });
-    if (Math.abs(area3) < 1e-9) {
+    if (Math.abs(area4) < 1e-9) {
       const n = polygon.length || 1;
       return [polygon.reduce((sum, p) => sum + p[0], 0) / n, polygon.reduce((sum, p) => sum + p[1], 0) / n];
     }
-    return [cx / (3 * area3), cy / (3 * area3)];
+    return [cx / (3 * area4), cy / (3 * area4)];
   }
   /** Width of a polygon along the horizontal line through `y`. */
   static span(polygon, y) {
@@ -230,7 +236,7 @@ var ConcordiaDiagram = (_class = class _ConcordiaDiagram extends _chunkXX2YTCORc
     description: "Plots U\u2013Pb isotope ratios with error ellipses against the concordia curve, in Wetherill or Tera-Wasserburg form, with an optional discordia line.",
     tags: ["geochronology", "u-pb", "zircon", "wetherill", "tera-wasserburg", "discordia", "isotopes"],
     sample: (seed) => {
-      const { u, n } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+      const { u, n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
       const upper = _ConcordiaDiagram.wetherill((1780 + u() * 160) * MA), lower = _ConcordiaDiagram.wetherill((380 + u() * 120) * MA);
       return Array.from({ length: 16 }, (_, i) => {
         const f = i < 4 ? 1 : 0.25 + u() * 0.72;
@@ -378,8 +384,8 @@ var ConcordiaDiagram = (_class = class _ConcordiaDiagram extends _chunkXX2YTCORc
     const px = (p) => [x(p[0]), y(p[1])];
     const span = this.#span(x, y);
     const curve = span ? Array.from({ length: 160 }, (_, i) => px(this.#curve(span[0] + i / 159 * (span[1] - span[0])))) : [];
-    const line4 = ctx.marks("concordia", "path", curve.length ? [curve] : [], () => "curve").attr("fill", "none").attr("stroke", ctx.theme.inkSecondary).attr("stroke-width", 2).attr("pointer-events", "none");
-    ctx.set(line4, { d: (d) => FieldKit.path(d, false) });
+    const line5 = ctx.marks("concordia", "path", curve.length ? [curve] : [], () => "curve").attr("fill", "none").attr("stroke", ctx.theme.inkSecondary).attr("stroke-width", 2).attr("pointer-events", "none");
+    ctx.set(line5, { d: (d) => FieldKit.path(d, false) });
     const tickAges = span ? _d3.ticks.call(void 0, span[0] / MA, span[1] / MA, Math.max(2, Math.min(7, Math.floor(width / 90)))).filter((age) => age > 0) : [];
     const ticks = tickAges.map((age) => ({ age, at: px(this.#curve(age * MA)) }));
     const knots = ctx.marks("age-ticks", "circle", ticks, (d) => d.age).attr("fill", ctx.theme.surface).attr("stroke", ctx.theme.inkSecondary).attr("stroke-width", 1.5).attr("pointer-events", "none");
@@ -463,7 +469,7 @@ var KiteDiagram = (_class2 = class extends _chunkXX2YTCORcjs.CartesianChart {
     description: "Shows species abundance along a transect as bands drawn symmetrically about each species\u2019 own axis, to reveal zonation.",
     tags: ["ecology", "transect", "zonation", "abundance", "field sampling", "rocky shore"],
     sample: (seed) => {
-      const { n } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+      const { n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
       return SHORE.flatMap(([species, centre, spread, peak]) => Array.from({ length: 22 }, (_, i) => {
         const at = i * 2, value = peak * Math.exp(-(((at - centre) / spread) ** 2) / 2) * (1 + n() * 0.12);
         return { position: at, species, abundance: Math.max(0, Math.round(value < 3 ? 0 : value)) };
@@ -569,7 +575,7 @@ var KobePlot = (_class3 = class _KobePlot extends _chunkXX2YTCORcjs.CartesianCha
     description: "Tracks a fish stock\u2019s B/BMSY against F/FMSY through the four Kobe status quadrants, as a yearly trajectory or the latest status.",
     tags: ["fisheries", "stock assessment", "msy", "kobe", "status", "ecology"],
     sample: (seed) => {
-      const { n } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+      const { n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
       return STOCKS.flatMap(([stock, start, peakF, cut]) => {
         let b = start;
         return Array.from({ length: 40 }, (_, i) => {
@@ -713,7 +719,7 @@ var PollenDiagram = (_class4 = class _PollenDiagram extends _chunkEEF5WRRZcjs.Ch
     description: "Plots each pollen taxon\u2019s percentage of the pollen sum against sediment depth in side-by-side panels, with \xD710 exaggeration and CONISS zones.",
     tags: ["palynology", "palaeoecology", "stratigraphy", "sediment core", "coniss", "holocene"],
     sample: (seed) => {
-      const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+      const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
       return Array.from({ length: 30 }, (_, i) => i * 10).flatMap((depth) => {
         const age = depth / 290, raw = TAXA.map(([, , profile]) => Math.max(0, profile(age) * (0.8 + u() * 0.4))), total = raw.reduce((a, b) => a + b, 0) || 1;
         return TAXA.map(([taxon, type], k) => ({ depth, taxon, type, count: Math.round(raw[k] / total * 400) }));
@@ -914,7 +920,7 @@ var QapfDiagram = (_class5 = class _QapfDiagram extends _chunkEEF5WRRZcjs.Chart2
     description: "Classifies igneous rocks by modal quartz, alkali feldspar, plagioclase and feldspathoid on the IUGS plutonic or volcanic double triangle.",
     tags: ["petrology", "igneous", "mineralogy", "iugs", "streckeisen", "ternary", "classification"],
     sample: (seed) => {
-      const { u } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+      const { u } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
       return SUITES.flatMap(([suite, q, f, ratio]) => Array.from({ length: 14 }, (_, i) => {
         const quartz = q[0] + u() * (q[1] - q[0]), foid = f[0] + u() * (f[1] - f[0]), r = (ratio[0] + u() * (ratio[1] - ratio[0])) / 100, rest = 100 - quartz - foid;
         return { sample: `${suite.split(" ")[0]} ${i + 1}`, suite, Q: +quartz.toFixed(1), A: +(rest * (1 - r)).toFixed(1), P: +(rest * r).toFixed(1), F: +foid.toFixed(1) };
@@ -1009,7 +1015,7 @@ var RarefactionCurve = (_class6 = class _RarefactionCurve extends _chunkXX2YTCOR
     description: "Plots expected species richness against individuals or samples drawn per site, with intervals and Chao richness estimates.",
     tags: ["ecology", "biodiversity", "species richness", "rarefaction", "chao", "accumulation"],
     sample: (seed) => {
-      const { u, n } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+      const { u, n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
       return [["Old-growth", 42], ["Secondary", 30], ["Plantation", 16]].flatMap(([site, pool]) => {
         const weights = Array.from({ length: +pool }, () => Math.exp(n() * 1.3) * 1.1);
         return Array.from({ length: 10 }, (_, j) => weights.flatMap((weight, s) => {
@@ -1116,6 +1122,217 @@ var RarefactionCurve = (_class6 = class _RarefactionCurve extends _chunkXX2YTCOR
   #title = "";
 }, _class6.__initStatic6(), _class6);
 
+// src/charts/fieldcharts/SpeciesAccumulationCurve.ts
+
+var STEPS2 = 48;
+var AHEAD = 24;
+var Z = 1.96;
+var diamond = _d3.symbol.call(void 0, _d3.symbolDiamond, 64)();
+var HABITATS = [["Riparian forest", 58, 12], ["Upland forest", 40, 10], ["Pasture", 20, 8]];
+var SpeciesAccumulationCurve = (_class7 = class _SpeciesAccumulationCurve extends _chunkXX2YTCORcjs.CartesianChart {
+  static __initStatic7() {this.meta = {
+    id: "species-accumulation-curve",
+    name: "Species accumulation curve",
+    family: "statistical" /* Statistical */,
+    dimension: 2 /* Two */,
+    shape: "table" /* Table */,
+    description: "Plots species found against sampling effort per site, as sampled and as the randomised mean with its band, extrapolated towards a Chao richness estimate.",
+    tags: ["ecology", "biodiversity", "species richness", "collector curve", "sampling effort", "extrapolation", "chao", "asymptote"],
+    sample: (seed) => {
+      const { u, n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
+      const poisson = (mean) => {
+        let k = 0;
+        for (let p = u(); p > Math.exp(-mean); p *= u()) k++;
+        return k;
+      };
+      return HABITATS.flatMap(([site, pool, traps]) => {
+        const rates = Array.from({ length: pool }, () => Math.exp(n() * 1.5 - 1.7));
+        return Array.from({ length: traps }, (_, j) => rates.flatMap((rate, s) => {
+          const count = poisson(rate);
+          return count > 0 ? [{ site, sample: `Trap ${j + 1}`, species: `sp${s + 1}`, count }] : [];
+        })).flat();
+      });
+    },
+    variants: {
+      order: ["collector-and-randomised" /* Both */, "collector" /* Collector */, "randomised" /* Randomised */],
+      effort: ["samples" /* Samples */, "individuals" /* Individuals */],
+      extrapolate: [true, false]
+    }
+  }}
+  defaults() {
+    return { site: "site", sample: "sample", species: "species", count: "count", order: "collector-and-randomised" /* Both */, effort: "samples" /* Samples */, extrapolate: true, reach: 2, yLabel: "Species" };
+  }
+  /** Species found so far after each sample, in the order given. Each sample maps a species to the individuals counted. */
+  static collector(samples) {
+    const seen = /* @__PURE__ */ new Set();
+    let individuals = 0;
+    return samples.map((found) => {
+      const before = seen.size;
+      for (const [name, count] of found) {
+        if (count <= 0) continue;
+        seen.add(name);
+        individuals += count;
+      }
+      return { species: seen.size, added: seen.size - before, individuals };
+    });
+  }
+  /** Chao's bias-corrected count of species present but not found, from those found once and twice. `scale` is (T - 1) / T for T samples and 1 for individuals. */
+  static unseen(ones, twos, scale = 1) {
+    return ones > 1 ? scale * ones * (ones - 1) / (2 * (twos + 1)) : 0;
+  }
+  /**
+   * Expected species `extra` units of effort past the `effort` observed, after Colwell et al. (2012): the curve closes
+   * on the Chao estimate at a rate set by the species found once. The variance treats the frequency counts as
+   * multinomial, so it meets the interpolated band at `extra = 0` and grows to that of the estimate itself.
+   */
+  static extrapolated(observed, ones, twos, effort, extra, scale = 1) {
+    const gain = (q1, q2) => {
+      const missing2 = _SpeciesAccumulationCurve.unseen(q1, q2, scale);
+      return missing2 > 0 ? missing2 * (1 - (1 - q1 / (effort * missing2 + q1)) ** extra) : 0;
+    };
+    const missing = _SpeciesAccumulationCurve.unseen(ones, twos, scale), total = observed + missing;
+    if (!(total > 0)) return { s: 0, variance: 0 };
+    const h = 1e-4;
+    const d1 = (gain(ones + h, twos) - gain(ones - h, twos)) / (2 * h), d2 = (gain(ones, twos + h) - gain(ones, twos - h)) / (2 * h);
+    const variance = missing * (1 - missing / total) + 2 * missing * (d1 * ones + d2 * twos) / total + d1 * d1 * ones * (1 - ones / total) + d2 * d2 * twos * (1 - twos / total) - 2 * d1 * d2 * ones * twos / total;
+    return { s: observed + gain(ones, twos), variance: Math.max(0, variance) };
+  }
+  get #bySamples() {
+    return this.options.effort !== "individuals" /* Individuals */;
+  }
+  get #chao() {
+    return this.#bySamples ? "Chao2" : "Chao1";
+  }
+  #curves() {
+    const { data, site, sample, species, count, extrapolate, reach } = this.options;
+    const sites = /* @__PURE__ */ new Map();
+    for (const row of Array.isArray(data) ? data : []) {
+      const value = FieldKit.number(_optionalChain([row, 'optionalAccess', _20 => _20[count]]));
+      if (value === null || value < 0) continue;
+      const key = FieldKit.text(row[site]) || "All", unit = FieldKit.text(row[sample]) || "All", name = FieldKit.text(row[species]);
+      const units = _nullishCoalesce(sites.get(key), () => ( sites.set(key, /* @__PURE__ */ new Map()).get(key)));
+      const found = _nullishCoalesce(units.get(unit), () => ( units.set(unit, /* @__PURE__ */ new Map()).get(unit)));
+      if (name && value > 0) found.set(name, (_nullishCoalesce(found.get(name), () => ( 0))) + value);
+    }
+    const bySamples = this.#bySamples, far = extrapolate ? Math.max(1, Number(reach) || 2) : 1;
+    const point = (n, { s, variance }) => {
+      const half = Z * Math.sqrt(variance);
+      return { n, s, low: Math.max(0, s - half), high: s + half };
+    };
+    return [...sites].map(([key, units]) => {
+      const names = [...units.keys()], samples = [...units.values()];
+      const tally = /* @__PURE__ */ new Map();
+      for (const found of samples) for (const [name, value] of found) {
+        const entry = _nullishCoalesce(tally.get(name), () => ( tally.set(name, { total: 0, units: 0 }).get(name)));
+        entry.total += value;
+        entry.units++;
+      }
+      const counts = [...tally.values()].map((entry) => bySamples ? entry.units : Math.round(entry.total)).filter((k) => k > 0);
+      const effort = bySamples ? samples.length : counts.reduce((sum, k) => sum + k, 0), observed = counts.length;
+      const ones = counts.filter((k) => k === 1).length, twos = counts.filter((k) => k === 2).length;
+      const scale = bySamples && effort ? (effort - 1) / effort : 1;
+      const estimate = observed + _SpeciesAccumulationCurve.unseen(ones, twos, scale);
+      const ahead = (extra) => _SpeciesAccumulationCurve.extrapolated(observed, ones, twos, effort, extra, scale);
+      const drawn = [...new Set(Array.from({ length: STEPS2 + 1 }, (_, i) => Math.max(1, Math.round(i / STEPS2 * effort))))];
+      const half = Z * Math.sqrt(ahead(Infinity).variance);
+      return {
+        site: key,
+        observed,
+        estimate,
+        effort,
+        ceiling: [Math.max(observed, estimate - half), estimate + half],
+        steps: _SpeciesAccumulationCurve.collector(samples).map((step, i) => ({ ...step, key: `${key}\0${i}`, site: key, sample: names[i], effort: bySamples ? i + 1 : step.individuals })),
+        mean: effort > 0 ? drawn.map((n) => point(n, RarefactionCurve.samples(counts, effort, n, estimate))) : [],
+        beyond: far > 1 ? Array.from({ length: AHEAD + 1 }, (_, i) => point(effort * (1 + i / AHEAD * (far - 1)), ahead(i / AHEAD * (far - 1) * effort))) : []
+      };
+    }).filter((curve) => curve.effort > 0);
+  }
+  get #collected() {
+    return this.options.order !== "randomised" /* Randomised */;
+  }
+  get #averaged() {
+    return this.options.order !== "collector" /* Collector */;
+  }
+  #label(curve) {
+    return `${this.#chao} ${curve.estimate.toFixed(1)}`;
+  }
+  /** Room on the right for the richness estimates the extrapolation is marked with. */
+  plotMargin() {
+    const margin = super.plotMargin(), marked = this.#curves().filter((curve) => curve.beyond.length);
+    return marked.length ? { ...margin, right: 16 + Math.max(...marked.map((curve) => FieldKit.width(this.#label(curve), 10))) } : margin;
+  }
+  scales(ctx) {
+    const curves = this.#curves();
+    const far = Math.max(1, ...curves.map((curve) => _nullishCoalesce(_optionalChain([curve, 'access', _21 => _21.beyond, 'access', _22 => _22.at, 'call', _23 => _23(-1), 'optionalAccess', _24 => _24.n]), () => ( curve.effort))));
+    const top = Math.max(1, ...curves.flatMap((curve) => [
+      curve.observed,
+      ...this.#averaged ? curve.mean.map((p) => p.high) : [],
+      ...curve.beyond.map((p) => p.high),
+      ...curve.beyond.length ? [curve.estimate] : []
+    ]));
+    return { x: _d3.scaleLinear.call(void 0, [0, far], [0, ctx.plot.width]).nice(), y: _d3.scaleLinear.call(void 0, [0, top * 1.05], [ctx.plot.height, 0]).nice() };
+  }
+  plotMarks(ctx, { x, y }) {
+    const curves = this.#curves(), { width, height } = ctx.plot;
+    const sx = (n) => _nullishCoalesce(x(n), () => ( 0)), sy = (s) => _nullishCoalesce(y(s), () => ( 0));
+    const tone = (d) => ctx.palette.color(d.site);
+    const unit = this.#bySamples ? "Samples" : "Individuals", chao = this.#chao;
+    const averaged = this.#averaged ? curves : [], ahead = curves.filter((curve) => curve.beyond.length);
+    const summary = (d) => [["Observed species", d.observed], [unit, d.effort], [`${chao} estimate`, d.estimate.toFixed(1)]];
+    const band = _d3.area.call(void 0, ).x((p) => sx(p.n)).y0((p) => sy(p.low)).y1((p) => sy(p.high));
+    const bands = ctx.marks("intervals", "path", averaged, (d) => d.site).attr("fill", tone).attr("fill-opacity", 0.16).attr("stroke", "none").attr("pointer-events", "none");
+    ctx.set(bands, { d: (d) => _nullishCoalesce(band(d.mean), () => ( "")) });
+    const cones = ctx.marks("extrapolated-intervals", "path", ahead, (d) => d.site).attr("fill", tone).attr("fill-opacity", 0.1).attr("stroke", "none").attr("pointer-events", "none");
+    ctx.set(cones, { d: (d) => _nullishCoalesce(band(d.beyond), () => ( "")) });
+    const trace = _d3.line.call(void 0, ).x((p) => sx(p.n)).y((p) => sy(p.s));
+    const dashed = ctx.marks("extrapolated", "path", ahead, (d) => d.site).attr("fill", "none").attr("stroke", tone).attr("stroke-width", 2).attr("stroke-dasharray", "6 4");
+    ctx.set(dashed, { d: (d) => _nullishCoalesce(trace(d.beyond), () => ( "")) });
+    ctx.hover(dashed, (d) => {
+      const end = d.beyond.at(-1);
+      return { title: `${d.site}, extrapolated`, color: tone(d), rows: [[`Expected at ${Math.round(end.n)} ${unit.toLowerCase()}`, end.s.toFixed(1)], ["95% interval", `${end.low.toFixed(1)} to ${end.high.toFixed(1)}`], ...summary(d)] };
+    });
+    const means = ctx.marks("curves", "path", averaged, (d) => d.site).attr("fill", "none").attr("stroke", tone).attr("stroke-width", 2).attr("stroke-linejoin", "round");
+    ctx.set(means, { d: (d) => _nullishCoalesce(trace(d.mean), () => ( "")) });
+    ctx.drawOn(means);
+    ctx.hover(means, (d) => ({ title: `${d.site}, randomised mean`, color: tone(d), rows: summary(d) }));
+    const collected = this.#collected ? curves : [];
+    const walk = _d3.line.call(void 0, ).x((step) => sx(step.effort)).y((step) => sy(step.species));
+    const walks = ctx.marks("collector", "path", collected, (d) => d.site).attr("fill", "none").attr("stroke", tone).attr("stroke-width", this.#averaged ? 1.5 : 2).attr("stroke-linejoin", "round");
+    ctx.set(walks, { d: (d) => _nullishCoalesce(walk(d.steps), () => ( "")) });
+    ctx.hover(walks, (d) => ({ title: `${d.site}, as sampled`, color: tone(d), rows: summary(d) }));
+    const roomy = collected.every((curve) => (sx(curve.steps.at(-1).effort) - sx(curve.steps[0].effort)) / curve.steps.length >= 12);
+    const dots = ctx.marks("samples", "circle", roomy ? collected.flatMap((curve) => curve.steps.slice(0, -1)) : [], (step) => step.key).attr("fill", tone).attr("stroke", ctx.theme.surface).attr("stroke-width", 1.5);
+    ctx.set(dots, { cx: (step) => sx(step.effort), cy: (step) => sy(step.species), r: 4 }, { r: 0 });
+    ctx.hover(dots, (step) => ({ title: `${step.site}, ${step.sample}`, color: tone(step), rows: [["Species so far", step.species], ["New in this sample", step.added], [unit, this.#bySamples ? step.effort : Math.round(step.effort)]] }));
+    const ends = ctx.marks("observed", "circle", curves, (d) => d.site).attr("fill", tone).attr("stroke", ctx.theme.surface).attr("stroke-width", 2);
+    ctx.set(ends, { cx: (d) => sx(d.effort), cy: (d) => sy(d.observed), r: 5 }, { r: 0 });
+    ctx.hover(ends, (d) => ({ title: d.site, color: tone(d), rows: summary(d) }));
+    const place = (d) => `translate(${width},${sy(d.estimate)})`;
+    const marks = ctx.marks("estimates", "path", ahead, (d) => d.site).attr("d", diamond).attr("fill", tone).attr("stroke", ctx.theme.surface).attr("stroke-width", 2).attr("paint-order", "stroke");
+    ctx.set(marks, { transform: place }, { transform: place });
+    ctx.hover(marks, (d) => ({
+      title: d.site,
+      color: tone(d),
+      rows: [[`${chao} estimate`, d.estimate.toFixed(1)], ["95% interval", `${d.ceiling[0].toFixed(1)} to ${d.ceiling[1].toFixed(1)}`], ["Observed species", d.observed], ["Share found", `${Math.round(100 * d.observed / Math.max(d.estimate, 1e-9))}%`]]
+    }));
+    const stacked = [...ahead].sort((a, b) => a.estimate - b.estimate).map((curve) => ({ curve, at: sy(curve.estimate) + 3.5 }));
+    stacked.forEach((label, i) => {
+      if (i) label.at = Math.min(label.at, stacked[i - 1].at - 11);
+    });
+    const lift = Math.max(0, 9 - (_nullishCoalesce(_optionalChain([stacked, 'access', _25 => _25.at, 'call', _26 => _26(-1), 'optionalAccess', _27 => _27.at]), () => ( height))));
+    const names = ctx.text(ctx.marks("estimate-labels", "text", stacked, (d) => d.curve.site), 10).attr("pointer-events", "none").text((d) => this.#label(d.curve));
+    ctx.set(names, { x: width + 9, y: (d) => d.at + lift }, { x: width + 9, y: (d) => d.at + lift });
+  }
+  /** The x title names the effort unless the caller set their own. */
+  paint() {
+    const title = this.#bySamples ? "Samples taken" : "Individuals counted";
+    if (this.options.xLabel === void 0 || this.options.xLabel === this.#title) this.options.xLabel = title;
+    this.#title = title;
+    super.paint();
+  }
+  #title = "";
+}, _class7.__initStatic7(), _class7);
+
 // src/charts/fieldcharts/SpiderDiagram.ts
 
 var CHONDRITE = { La: 0.237, Ce: 0.612, Pr: 0.095, Nd: 0.467, Sm: 0.153, Eu: 0.058, Gd: 0.2055, Tb: 0.0374, Dy: 0.254, Ho: 0.0566, Er: 0.1655, Tm: 0.0255, Yb: 0.17, Lu: 0.0254 };
@@ -1249,8 +1466,8 @@ var REFERENCES = {
   ["n-morb" /* NMorb */]: ["N-MORB", N_MORB]
 };
 var SUITES2 = [["OIB-like", OIB], ["E-MORB-like", E_MORB], ["N-MORB-like", N_MORB]];
-var SpiderDiagram = (_class7 = class extends _chunkXX2YTCORcjs.CartesianChart {
-  static __initStatic7() {this.meta = {
+var SpiderDiagram = (_class8 = class extends _chunkXX2YTCORcjs.CartesianChart {
+  static __initStatic8() {this.meta = {
     id: "trace-element-spider",
     name: "Spider diagram (normalised)",
     family: "scientific" /* Scientific */,
@@ -1259,7 +1476,7 @@ var SpiderDiagram = (_class7 = class extends _chunkXX2YTCORcjs.CartesianChart {
     description: "Plots element concentrations normalised to chondrite, primitive mantle or N-MORB on an ordered element axis, to compare rock suites.",
     tags: ["geochemistry", "trace elements", "ree", "normalised", "multi-element", "sun mcdonough"],
     sample: (seed) => {
-      const { n } = _chunkUMEKIIHJcjs.SampleData.rng(seed);
+      const { n } = _chunkLNWTC72Ycjs.SampleData.rng(seed);
       return SUITES2.flatMap(([group, reference]) => Array.from({ length: 4 }, (_, i) => {
         const shift = n() * 0.12;
         return { sample: `${group.split("-")[0]}-${i + 1}`, suite: group, ...Object.fromEntries(Object.entries(reference).map(([element, value]) => [element, +(value * Math.exp(shift + n() * 0.1)).toPrecision(3)])) };
@@ -1341,19 +1558,19 @@ var SpiderDiagram = (_class7 = class extends _chunkXX2YTCORcjs.CartesianChart {
     super.paint();
   }
   #title = "";
-}, _class7.__initStatic7(), _class7);
+}, _class8.__initStatic8(), _class8);
 
 // src/charts/fieldcharts/FieldDiagram.ts
 
-var FieldDiagram = (_class8 = class extends _chunkXX2YTCORcjs.CartesianChart {constructor(...args) { super(...args); _class8.prototype.__init.call(this); }
+var FieldDiagram = (_class9 = class extends _chunkXX2YTCORcjs.CartesianChart {constructor(...args) { super(...args); _class9.prototype.__init.call(this); }
   __init() {this.outside = "Outside the fields"}
   #samples() {
     const { data, x, y, group, label } = this.options;
     const shapes = this.fieldShapes();
     return (Array.isArray(data) ? data : []).flatMap((row, i) => {
-      const px = FieldKit.number(_optionalChain([row, 'optionalAccess', _20 => _20[x]])), py = FieldKit.number(_optionalChain([row, 'optionalAccess', _21 => _21[y]]));
+      const px = FieldKit.number(_optionalChain([row, 'optionalAccess', _28 => _28[x]])), py = FieldKit.number(_optionalChain([row, 'optionalAccess', _29 => _29[y]]));
       if (px === null || py === null) return [];
-      const field2 = _nullishCoalesce(_optionalChain([shapes, 'access', _22 => _22.find, 'call', _23 => _23((shape) => FieldKit.inside([px, py], shape.polygon)), 'optionalAccess', _24 => _24.name]), () => ( this.outside));
+      const field2 = _nullishCoalesce(_optionalChain([shapes, 'access', _30 => _30.find, 'call', _31 => _31((shape) => FieldKit.inside([px, py], shape.polygon)), 'optionalAccess', _32 => _32.name]), () => ( this.outside));
       return [{ key: `${i}`, row, x: px, y: py, group: FieldKit.text(row[group]) || "Samples", label: FieldKit.text(row[label]), field: field2 }];
     });
   }
@@ -1369,7 +1586,7 @@ var FieldDiagram = (_class8 = class extends _chunkXX2YTCORcjs.CartesianChart {co
     const toPixel = ([a, b]) => [_nullishCoalesce(x(a), () => ( 0)), _nullishCoalesce(y(b), () => ( 0))];
     const shapes = this.fieldShapes().map((shape) => ({ ...shape, pixels: FieldKit.clip(shape.polygon, extent.x, extent.y).map(toPixel) })).filter((shape) => shape.pixels.length > 2);
     const filled = fields === "filled" /* Filled */;
-    const columns = filled || this.options.grid === false ? [] : _nullishCoalesce(_optionalChain([x, 'access', _25 => _25.ticks, 'optionalCall', _26 => _26(Math.max(2, Math.floor(ctx.plot.width / 80)))]), () => ( []));
+    const columns = filled || this.options.grid === false ? [] : _nullishCoalesce(_optionalChain([x, 'access', _33 => _33.ticks, 'optionalCall', _34 => _34(Math.max(2, Math.floor(ctx.plot.width / 80)))]), () => ( []));
     const verticals = ctx.marks("axis-grid-x", "line", columns, (d) => d).attr("stroke", ctx.theme.grid).attr("stroke-width", 1);
     ctx.set(verticals, { x1: (d) => _nullishCoalesce(x(d), () => ( 0)), x2: (d) => _nullishCoalesce(x(d), () => ( 0)), y1: 0, y2: ctx.plot.height });
     const areas = ctx.marks("fields", "path", shapes, (d) => d.key).attr("stroke-linejoin", "round").attr("fill", (_, i) => filled ? ctx.palette.sequential(0.06 + i % 4 * 0.08) : "transparent").attr("stroke", filled ? ctx.theme.surface : ctx.theme.inkMuted).attr("stroke-width", filled ? 2 : 1.25);
@@ -1396,7 +1613,7 @@ var FieldDiagram = (_class8 = class extends _chunkXX2YTCORcjs.CartesianChart {co
     ctx.set(dots, { cx: (d) => d.px, cy: (d) => d.py, r: 4.5 }, { r: 0 }, Math.min(3, 300 / (shown.length || 1)));
     ctx.hover(dots, (d) => ({ title: d.label || d.field, color: ctx.palette.color(d.group), rows: [["Field", d.field], [xName, +d.x.toFixed(2)], [yName, +d.y.toFixed(2)], ["Group", d.group]] }));
   }
-}, _class8);
+}, _class9);
 
 // src/charts/fieldcharts/TasDiagram.ts
 var TAS = [
@@ -1419,8 +1636,8 @@ var TAS = [
 ];
 var SHAPES = TAS.map(([code, name, polygon]) => ({ key: code, code, name, polygon }));
 var UNITS = ["Lower unit", "Middle unit", "Upper unit"];
-var TasDiagram = (_class9 = class extends FieldDiagram {
-  static __initStatic8() {this.meta = {
+var TasDiagram = (_class10 = class extends FieldDiagram {
+  static __initStatic9() {this.meta = {
     id: "tas-diagram",
     name: "TAS diagram",
     family: "scientific" /* Scientific */,
@@ -1428,7 +1645,7 @@ var TasDiagram = (_class9 = class extends FieldDiagram {
     shape: "table" /* Table */,
     description: "Classifies volcanic rocks by total alkali (Na2O + K2O) against silica on the Le Bas et al. (1986) IUGS fields.",
     tags: ["geochemistry", "petrology", "igneous", "classification", "iugs", "volcanic"],
-    sample: (seed) => _chunkUMEKIIHJcjs.SampleData.points(seed, 90, 3).map((row, i) => {
+    sample: (seed) => _chunkLNWTC72Ycjs.SampleData.points(seed, 90, 3).map((row, i) => {
       const g = i % 3, silica = 41 + (+row.x - 5) * 0.48;
       const alkali = 1.2 + (silica - 41) * 0.2 + (+row.y - +row.x * 0.8 - 10) / 9 * 0.9 + (g === 1 ? 1.6 : 0);
       return { sample: `S-${i + 1}`, unit: UNITS[g], SiO2: +silica.toFixed(2), alkali: +Math.max(0.2, alkali).toFixed(2) };
@@ -1462,7 +1679,7 @@ var TasDiagram = (_class9 = class extends FieldDiagram {
   measures() {
     return ["SiO\u2082 wt%", "Na\u2082O + K\u2082O wt%"];
   }
-}, _class9.__initStatic8(), _class9);
+}, _class10.__initStatic9(), _class10);
 
 // src/charts/fieldcharts/WhittakerBiomeDiagram.ts
 var BIOMES = [
@@ -1478,8 +1695,8 @@ var BIOMES = [
 ];
 var SHAPES2 = BIOMES.map(([code, name, polygon]) => ({ key: code, code, name, polygon }));
 var NETWORKS = { Alpha: "Network A", Bravo: "Network B", Charlie: "Network C" };
-var WhittakerBiomeDiagram = (_class10 = class extends FieldDiagram {constructor(...args2) { super(...args2); _class10.prototype.__init2.call(this); }
-  static __initStatic9() {this.meta = {
+var WhittakerBiomeDiagram = (_class11 = class extends FieldDiagram {constructor(...args2) { super(...args2); _class11.prototype.__init2.call(this); }
+  static __initStatic10() {this.meta = {
     id: "whittaker-biome-diagram",
     name: "Whittaker biome diagram",
     family: "scientific" /* Scientific */,
@@ -1487,7 +1704,7 @@ var WhittakerBiomeDiagram = (_class10 = class extends FieldDiagram {constructor(
     shape: "table" /* Table */,
     description: "Places sites on Whittaker biome regions by mean annual temperature and precipitation, naming the biome each falls in.",
     tags: ["ecology", "biome", "climate", "vegetation", "whittaker", "ricklefs"],
-    sample: (seed) => _chunkUMEKIIHJcjs.SampleData.points(seed, 140, 3).map((row, i) => {
+    sample: (seed) => _chunkLNWTC72Ycjs.SampleData.points(seed, 140, 3).map((row, i) => {
       const x = +row.x, wet = +row.y - x * 0.8 - 10;
       return { site: `Site ${i + 1}`, network: _nullishCoalesce(NETWORKS[String(row.group)], () => ( "Network A")), temperature: +(-15 + (x - 5) * 0.58).toFixed(1), precipitation: Math.round(15 + (x - 5) * 3 + wet * 7) };
     }).filter((row) => SHAPES2.some((shape) => FieldKit.inside([row.temperature, row.precipitation], shape.polygon))).slice(0, 70),
@@ -1521,7 +1738,7 @@ var WhittakerBiomeDiagram = (_class10 = class extends FieldDiagram {constructor(
   measures() {
     return ["Temperature \xB0C", "Precipitation cm"];
   }
-}, _class10.__initStatic9(), _class10);
+}, _class11.__initStatic10(), _class11);
 
 // src/charts/fieldcharts/index.ts
 var fieldchartsCatalog = [
@@ -1533,7 +1750,8 @@ var fieldchartsCatalog = [
   KiteDiagram,
   PollenDiagram,
   WhittakerBiomeDiagram,
-  KobePlot
+  KobePlot,
+  SpeciesAccumulationCurve
 ];
 
 
@@ -1559,4 +1777,6 @@ var fieldchartsCatalog = [
 
 
 
-exports.FieldStyle = FieldStyle; exports.FieldSamples = FieldSamples; exports.ConcordiaProjection = ConcordiaProjection; exports.ConcordiaMarks = ConcordiaMarks; exports.SpiderNormaliser = SpiderNormaliser; exports.SpiderDisplay = SpiderDisplay; exports.QapfRocks = QapfRocks; exports.RarefactionBasis = RarefactionBasis; exports.KiteScaling = KiteScaling; exports.PollenStyle = PollenStyle; exports.PollenWidth = PollenWidth; exports.KobePath = KobePath; exports.ConcordiaDiagram = ConcordiaDiagram; exports.KiteDiagram = KiteDiagram; exports.KobePlot = KobePlot; exports.PollenDiagram = PollenDiagram; exports.QapfDiagram = QapfDiagram; exports.RarefactionCurve = RarefactionCurve; exports.SpiderDiagram = SpiderDiagram; exports.TasDiagram = TasDiagram; exports.WhittakerBiomeDiagram = WhittakerBiomeDiagram; exports.fieldchartsCatalog = fieldchartsCatalog;
+
+
+exports.FieldStyle = FieldStyle; exports.FieldSamples = FieldSamples; exports.ConcordiaProjection = ConcordiaProjection; exports.ConcordiaMarks = ConcordiaMarks; exports.SpiderNormaliser = SpiderNormaliser; exports.SpiderDisplay = SpiderDisplay; exports.QapfRocks = QapfRocks; exports.RarefactionBasis = RarefactionBasis; exports.AccumulationOrder = AccumulationOrder; exports.KiteScaling = KiteScaling; exports.PollenStyle = PollenStyle; exports.PollenWidth = PollenWidth; exports.KobePath = KobePath; exports.ConcordiaDiagram = ConcordiaDiagram; exports.KiteDiagram = KiteDiagram; exports.KobePlot = KobePlot; exports.PollenDiagram = PollenDiagram; exports.QapfDiagram = QapfDiagram; exports.RarefactionCurve = RarefactionCurve; exports.SpeciesAccumulationCurve = SpeciesAccumulationCurve; exports.SpiderDiagram = SpiderDiagram; exports.TasDiagram = TasDiagram; exports.WhittakerBiomeDiagram = WhittakerBiomeDiagram; exports.fieldchartsCatalog = fieldchartsCatalog;