@sjcrh/proteinpaint-client 2.207.1 → 2.209.0

package/dist/chunk-V3SOBDIT.js

@@ -0,0 +1,255 @@
+import {
+  first_genetrack_tolist
+} from "./chunk-C3HEDQPT.js";
+import {
+  HYPER_COLOR,
+  HYPO_COLOR
+} from "./chunk-OBDIJ4QS.js";
+
+// plots/dmr/settings/defaults.ts
+function getDefaultDMRSettings(opts) {
+  const overrides = opts.settings || {};
+  const dm = opts?.app?.vocabApi?.termdbConfig?.queries?.dnaMethylation;
+  const chr = opts?.coordinateOverride?.chr;
+  const elementScale = dm?.regionAnalysis == "element" || Array.isArray(dm?.cpgChroms) && !!chr && !dm.cpgChroms.includes(chr);
+  const defaults = {
+    blockWidth: 800,
+    pad: elementScale ? 1e5 : 2e3,
+    lambda: elementScale ? 5e4 : 1e3,
+    C: 2,
+    fdr_cutoff: 0.05,
+    colors: {
+      group1: "#3b5ee6",
+      group2: "#c04e00",
+      hyper: HYPER_COLOR,
+      hypo: HYPO_COLOR
+    },
+    maxLoessRegion: 25e4,
+    minProbesForCi: 10,
+    backend: "rust",
+    maxRegionSize: 5e6
+  };
+  if (overrides.colors) {
+    Object.assign(defaults.colors, overrides.colors);
+    delete overrides.colors;
+  }
+  return Object.assign(defaults, overrides);
+}
+
+// plots/dmr/viewModel/DmrViewModel.ts
+var CCRE_TRACK_NAME = "ENCODE cCREs";
+var DmrViewModel = class {
+  constructor(dmrResult, config, genomeObj, queryChr, queryStart, queryStop) {
+    const { settings } = config;
+    const dmrBedItems = this.makeDmrBedItems(dmrResult, settings);
+    const sigCpgBedItems = this.makeSigCpgBedItems(dmrResult, settings, queryChr, queryStart, queryStop);
+    const xRange = (queryStop ?? 0) - (queryStart ?? 0);
+    const loess = dmrResult.diagnostic?.loess;
+    const showLoess = !!(loess && loess.group1_fitted.length > 0 && loess.group2_fitted.length > 0);
+    const showDots = xRange <= settings.dmr.maxLoessRegion;
+    const betaTrackResult = dmrResult.diagnostic ? this.renderBetaTrack(
+      dmrResult.diagnostic,
+      config,
+      settings.dmr.blockWidth,
+      showLoess,
+      showDots,
+      queryStart,
+      queryStop
+    ) : void 0;
+    this.viewData = {
+      tklst: this.buildTrackList(dmrBedItems, sigCpgBedItems, genomeObj, betaTrackResult?.img),
+      legendRows: this.buildLegendData(
+        config,
+        dmrResult.dmrs,
+        sigCpgBedItems,
+        showLoess,
+        showDots,
+        betaTrackResult?.showCi ?? false
+      ),
+      diagnostic: dmrResult.diagnostic,
+      dmrs: dmrResult.dmrs,
+      dmrBedItems,
+      showLoess,
+      showDots
+    };
+  }
+  buildTrackList(dmrBedItems, sigCpgBedItems, genomeObj, betaTrackImg) {
+    const tklst = [];
+    first_genetrack_tolist(genomeObj, tklst);
+    const ccre = (genomeObj?.tracks || []).find((t) => t.name == CCRE_TRACK_NAME);
+    if (ccre) tklst.push(structuredClone(ccre));
+    tklst.push({ type: "bedj", name: "DMRs", bedItems: dmrBedItems });
+    tklst.push({ type: "bedj", name: "Sig. CpGs", bedItems: sigCpgBedItems });
+    if (betaTrackImg) {
+      tklst.push({
+        type: "bigwig",
+        name: "Per-CpG Means",
+        height: 150,
+        imgData: betaTrackImg
+      });
+    }
+    return tklst;
+  }
+  buildLegendData(config, dmrs, sigCpgBedItems, showLoess, showDots, showCi) {
+    const { colors } = config.settings.dmr;
+    const g1 = config.group1Name || "Group 1";
+    const g2 = config.group2Name || "Group 2";
+    const meansItems = [];
+    if (showDots) {
+      meansItems.push({ text: `${g1} (control)`, color: colors.group1 }, { text: `${g2} (case)`, color: colors.group2 });
+    }
+    if (showLoess) {
+      const ciLabel = showCi ? " + 95% CI" : "";
+      meansItems.push(
+        { text: `${g1} LOESS trend${ciLabel}`, color: colors.group1, style: showCi ? "shaded" : "dashed" },
+        { text: `${g2} LOESS trend${ciLabel}`, color: colors.group2, style: showCi ? "shaded" : "dashed" }
+      );
+    }
+    const rows = [{ label: "Per-CpG Means", items: meansItems }];
+    const hasHyper = dmrs.some((d) => d.direction === "hyper");
+    const hasHypo = dmrs.some((d) => d.direction === "hypo");
+    if (hasHyper || hasHypo) {
+      const items = [];
+      if (hasHyper) items.push({ text: "Hypermethylated", color: colors.hyper });
+      if (hasHypo) items.push({ text: "Hypomethylated", color: colors.hypo });
+      rows.push({ label: "DMR", items });
+    }
+    if (sigCpgBedItems.length) {
+      const items = [];
+      const hasHyperCpg = sigCpgBedItems.some((b) => b.color === colors.hyper);
+      const hasHypoCpg = sigCpgBedItems.some((b) => b.color === colors.hypo);
+      if (hasHyperCpg) items.push({ text: "Hyper (FDR sig.)", color: colors.hyper });
+      if (hasHypoCpg) items.push({ text: "Hypo (FDR sig.)", color: colors.hypo });
+      rows.push({ label: "Sig. CpGs", items });
+    }
+    return rows;
+  }
+  /**
+   * Render the per-CpG means scatter plot to an offscreen canvas and return
+   * a data URI suitable for the bigwig imgData track.
+   */
+  renderBetaTrack(diagnostic, config, blockWidth, showLoess, showDots, queryStart, queryStop) {
+    const { probes } = diagnostic;
+    if (!probes.positions.length) return void 0;
+    const { colors, fdr_cutoff, minProbesForCi } = config.settings.dmr;
+    const dpr = typeof window !== "undefined" && window.devicePixelRatio > 1 ? window.devicePixelRatio : 1;
+    const width = blockWidth;
+    const height = 150;
+    const canvas = document.createElement("canvas");
+    canvas.width = width * dpr;
+    canvas.height = height * dpr;
+    const ctx = canvas.getContext("2d");
+    if (!ctx) return void 0;
+    ctx.scale(dpr, dpr);
+    const xMin = queryStart ?? probes.positions[0];
+    const xMax = queryStop ?? probes.positions[probes.positions.length - 1];
+    const xRange = xMax - xMin || 1;
+    const scaleX = (val) => (val - xMin) / xRange * width;
+    const scaleY = (val) => height - val * height;
+    ctx.clearRect(0, 0, width, height);
+    let showCi = false;
+    if (showLoess && diagnostic.loess) {
+      const { loess } = diagnostic;
+      const firstProbePos = probes.positions[0];
+      const lastProbePos = probes.positions[probes.positions.length - 1];
+      showCi = probes.positions.length >= minProbesForCi;
+      for (const [fitted, ciLower, ciUpper, color] of [
+        [loess.group1_fitted, loess.group1_ci_lower, loess.group1_ci_upper, colors.group1],
+        [loess.group2_fitted, loess.group2_ci_lower, loess.group2_ci_upper, colors.group2]
+      ]) {
+        if (!fitted.length) continue;
+        const lPos = loess.positions;
+        let iStart = 0;
+        let iEnd = lPos.length - 1;
+        while (iStart < lPos.length && lPos[iStart] < firstProbePos) iStart++;
+        while (iEnd >= 0 && lPos[iEnd] > lastProbePos) iEnd--;
+        if (iStart > iEnd) continue;
+        if (showCi) {
+          ctx.globalAlpha = 0.12;
+          ctx.fillStyle = color;
+          ctx.beginPath();
+          for (let i = iStart; i <= iEnd; i++) {
+            ctx.lineTo(scaleX(lPos[i]), scaleY(Math.max(0, Math.min(1, ciUpper[i]))));
+          }
+          for (let i = iEnd; i >= iStart; i--) {
+            ctx.lineTo(scaleX(lPos[i]), scaleY(Math.max(0, Math.min(1, ciLower[i]))));
+          }
+          ctx.closePath();
+          ctx.fill();
+        }
+        ctx.globalAlpha = 0.8;
+        ctx.strokeStyle = color;
+        ctx.lineWidth = 2;
+        ctx.setLineDash(showCi ? [] : [6, 4]);
+        ctx.beginPath();
+        for (let i = iStart; i <= iEnd; i++) {
+          ctx.lineTo(scaleX(lPos[i]), scaleY(Math.max(0, Math.min(1, fitted[i]))));
+        }
+        ctx.stroke();
+        ctx.setLineDash([]);
+      }
+    }
+    if (!showDots) {
+      ctx.globalAlpha = 1;
+      return { img: { minv: 0, maxv: 1, src: canvas.toDataURL("image/png") }, showCi };
+    }
+    for (let i = 0; i < probes.positions.length; i++) {
+      const x = scaleX(probes.positions[i]);
+      const isSig = probes.fdr[i] < fdr_cutoff;
+      const alpha = isSig ? 0.85 : 0.3;
+      ctx.globalAlpha = alpha;
+      ctx.fillStyle = colors.group1;
+      const m1 = probes.mean_group1[i];
+      if (m1 != null) {
+        ctx.beginPath();
+        ctx.arc(x, scaleY(m1), 4, 0, Math.PI * 2);
+        ctx.fill();
+      }
+      ctx.fillStyle = colors.group2;
+      const m2 = probes.mean_group2[i];
+      if (m2 != null) {
+        ctx.beginPath();
+        ctx.arc(x, scaleY(m2), 4, 0, Math.PI * 2);
+        ctx.fill();
+      }
+    }
+    ctx.globalAlpha = 1;
+    return { img: { minv: 0, maxv: 1, src: canvas.toDataURL("image/png") }, showCi };
+  }
+  makeDmrBedItems(dmrResult, settings) {
+    return dmrResult.dmrs.map((dmr) => {
+      const negLog = -Math.log10(Math.max(dmr.min_smoothed_fdr, 1e-300));
+      const alpha = Math.round(Math.min(255, Math.max(50, negLog / 10 * 255)));
+      const hex = alpha.toString(16).padStart(2, "0");
+      const base = dmr.direction === "hyper" ? settings.dmr.colors.hyper : settings.dmr.colors.hypo;
+      return { chr: dmr.chr, start: dmr.start, stop: dmr.stop, color: base + hex };
+    });
+  }
+  makeSigCpgBedItems(dmrResult, settings, chr, queryStart, queryStop) {
+    const diag = dmrResult.diagnostic;
+    if (!diag) return [];
+    const { probes } = diag;
+    const items = [];
+    const minDeltaBeta = 0.05;
+    for (let i = 0; i < probes.positions.length; i++) {
+      if (probes.fdr[i] >= settings.dmr.fdr_cutoff) continue;
+      const pos = probes.positions[i];
+      if (queryStart != null && queryStop != null && (pos < queryStart || pos > queryStop)) continue;
+      const mg1 = probes.mean_group1[i];
+      const mg2 = probes.mean_group2[i];
+      if (mg1 == null || mg2 == null) continue;
+      const deltaBeta = mg2 - mg1;
+      if (Math.abs(deltaBeta) < minDeltaBeta) continue;
+      const color = deltaBeta >= 0 ? settings.dmr.colors.hyper : settings.dmr.colors.hypo;
+      items.push({ chr, start: pos, stop: pos + 1, color });
+    }
+    return items;
+  }
+};
+
+export {
+  getDefaultDMRSettings,
+  CCRE_TRACK_NAME,
+  DmrViewModel
+};
+//# sourceMappingURL=chunk-V3SOBDIT.js.map

package/dist/chunk-V3SOBDIT.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../plots/dmr/settings/defaults.ts", "../plots/dmr/viewModel/DmrViewModel.ts"],
+  "sourcesContent": ["import type { DMRSettings } from './Settings.ts'\n\n// direction colours are shared with the server-rendered scan Manhattan\nexport { HYPER_COLOR, HYPO_COLOR } from '#shared/dmrColors.js'\nimport { HYPER_COLOR, HYPO_COLOR } from '#shared/dmrColors.js'\n\nexport function getDefaultDMRSettings(opts: any): DMRSettings {\n\tconst overrides = opts.settings || {}\n\t/* A dataset with no CpG-level matrix runs the region analysis on its element matrix, where one\n\trow is a cCRE, not a CpG. Those sit ~10kb apart against a CpG's ~100bp, so the CpG-scale window\n\tand kernel below would frame one element and smooth nothing. Scaled to element spacing instead.\n\tponytail: fixed values, not derived from the matrix's actual spacing \u2014 worth deriving only if a\n\tdataset shows up whose element density is far off this one's. The user can pan/zoom either way. */\n\tconst dm = opts?.app?.vocabApi?.termdbConfig?.queries?.dnaMethylation\n\t/* Per chromosome, not per dataset: a cohort with some shards built runs CpG resolution where one\n\texists and elements everywhere else (server resolveMethylationMatrix), so a dataset-wide flag gave\n\tthe fallback chromosomes a CpG-scale window over rows ~10 kb apart. cpgChroms is sent only for a\n\tshard-backed dataset without a genome-wide file. */\n\tconst chr = opts?.coordinateOverride?.chr\n\tconst elementScale =\n\t\tdm?.regionAnalysis == 'element' || (Array.isArray(dm?.cpgChroms) && !!chr && !dm.cpgChroms.includes(chr))\n\tconst defaults = {\n\t\tblockWidth: 800,\n\t\tpad: elementScale ? 100_000 : 2000,\n\t\tlambda: elementScale ? 50_000 : 1000,\n\t\tC: 2,\n\t\tfdr_cutoff: 0.05,\n\t\tcolors: {\n\t\t\tgroup1: '#3b5ee6',\n\t\t\tgroup2: '#c04e00',\n\t\t\thyper: HYPER_COLOR,\n\t\t\thypo: HYPO_COLOR\n\t\t},\n\t\tmaxLoessRegion: 250_000,\n\t\tminProbesForCi: 10,\n\t\tbackend: 'rust' as const,\n\t\tmaxRegionSize: 5_000_000\n\t}\n\n\t// Deep-merge colors so hyper/hypo defaults are preserved\n\t// when only group colors are overridden\n\tif (overrides.colors) {\n\t\tObject.assign(defaults.colors, overrides.colors)\n\t\tdelete overrides.colors\n\t}\n\n\treturn Object.assign(defaults, overrides)\n}\n", "import { first_genetrack_tolist } from '#common/1stGenetk'\n\n/** Name of the regulatory-element track in the genome config, if it declares one. */\nexport const CCRE_TRACK_NAME = 'ENCODE cCREs'\nimport type { TermdbDmrSuccessResponse, DmrDiagnostic } from '#types'\nimport type { DmrConfig, BedItem, LegendRow, DmrViewData } from '../DmrTypes.ts'\n\nexport class DmrViewModel {\n\tviewData: DmrViewData\n\n\tconstructor(\n\t\tdmrResult: TermdbDmrSuccessResponse,\n\t\tconfig: DmrConfig,\n\t\tgenomeObj: any,\n\t\tqueryChr: string,\n\t\tqueryStart?: number,\n\t\tqueryStop?: number\n\t) {\n\t\tconst { settings } = config\n\t\tconst dmrBedItems = this.makeDmrBedItems(dmrResult, settings)\n\t\tconst sigCpgBedItems = this.makeSigCpgBedItems(dmrResult, settings, queryChr, queryStart, queryStop)\n\n\t\tconst xRange = (queryStop ?? 0) - (queryStart ?? 0)\n\t\tconst loess = dmrResult.diagnostic?.loess\n\t\tconst showLoess = !!(loess && loess.group1_fitted.length > 0 && loess.group2_fitted.length > 0)\n\t\tconst showDots = xRange <= settings.dmr.maxLoessRegion\n\n\t\tconst betaTrackResult = dmrResult.diagnostic\n\t\t\t? this.renderBetaTrack(\n\t\t\t\t\tdmrResult.diagnostic,\n\t\t\t\t\tconfig,\n\t\t\t\t\tsettings.dmr.blockWidth,\n\t\t\t\t\tshowLoess,\n\t\t\t\t\tshowDots,\n\t\t\t\t\tqueryStart,\n\t\t\t\t\tqueryStop\n\t\t\t  )\n\t\t\t: undefined\n\n\t\tthis.viewData = {\n\t\t\ttklst: this.buildTrackList(dmrBedItems, sigCpgBedItems, genomeObj, betaTrackResult?.img),\n\t\t\tlegendRows: this.buildLegendData(\n\t\t\t\tconfig,\n\t\t\t\tdmrResult.dmrs,\n\t\t\t\tsigCpgBedItems,\n\t\t\t\tshowLoess,\n\t\t\t\tshowDots,\n\t\t\t\tbetaTrackResult?.showCi ?? false\n\t\t\t),\n\t\t\tdiagnostic: dmrResult.diagnostic,\n\t\t\tdmrs: dmrResult.dmrs,\n\t\t\tdmrBedItems,\n\t\t\tshowLoess,\n\t\t\tshowDots\n\t\t}\n\t}\n\n\tprivate buildTrackList(\n\t\tdmrBedItems: BedItem[],\n\t\tsigCpgBedItems: BedItem[],\n\t\tgenomeObj: any,\n\t\tbetaTrackImg?: { minv: number; maxv: number; src: string }\n\t): any[] {\n\t\tconst tklst: any[] = []\n\t\tfirst_genetrack_tolist(genomeObj, tklst)\n\t\t/* Regulatory context, switched on here rather than left in the Tracks menu. A DMR next to a\n\t\tgene model says where it is; a DMR next to the cCREs says what it is sitting on, which is the\n\t\tquestion the element-level view was answering. Taken from the genome's own declaration by\n\t\tname, so a genome that does not declare it simply renders without the row. */\n\t\tconst ccre = (genomeObj?.tracks || []).find((t: any) => t.name == CCRE_TRACK_NAME)\n\t\tif (ccre) tklst.push(structuredClone(ccre))\n\t\ttklst.push({ type: 'bedj', name: 'DMRs', bedItems: dmrBedItems })\n\t\ttklst.push({ type: 'bedj', name: 'Sig. CpGs', bedItems: sigCpgBedItems })\n\t\tif (betaTrackImg) {\n\t\t\ttklst.push({\n\t\t\t\ttype: 'bigwig',\n\t\t\t\tname: 'Per-CpG Means',\n\t\t\t\theight: 150,\n\t\t\t\timgData: betaTrackImg\n\t\t\t})\n\t\t}\n\t\treturn tklst\n\t}\n\n\tprivate buildLegendData(\n\t\tconfig: DmrConfig,\n\t\tdmrs: TermdbDmrSuccessResponse['dmrs'],\n\t\tsigCpgBedItems: BedItem[],\n\t\tshowLoess: boolean,\n\t\tshowDots: boolean,\n\t\tshowCi: boolean\n\t): LegendRow[] {\n\t\tconst { colors } = config.settings.dmr\n\t\tconst g1 = config.group1Name || 'Group 1'\n\t\tconst g2 = config.group2Name || 'Group 2'\n\t\tconst meansItems: LegendRow['items'] = []\n\t\tif (showDots) {\n\t\t\tmeansItems.push({ text: `${g1} (control)`, color: colors.group1 }, { text: `${g2} (case)`, color: colors.group2 })\n\t\t}\n\t\tif (showLoess) {\n\t\t\tconst ciLabel = showCi ? ' + 95% CI' : ''\n\t\t\tmeansItems.push(\n\t\t\t\t{ text: `${g1} LOESS trend${ciLabel}`, color: colors.group1, style: showCi ? 'shaded' : 'dashed' },\n\t\t\t\t{ text: `${g2} LOESS trend${ciLabel}`, color: colors.group2, style: showCi ? 'shaded' : 'dashed' }\n\t\t\t)\n\t\t}\n\t\tconst rows: LegendRow[] = [{ label: 'Per-CpG Means', items: meansItems }]\n\t\t// Only show DMR legend entries for directions present in the results\n\t\tconst hasHyper = dmrs.some(d => d.direction === 'hyper')\n\t\tconst hasHypo = dmrs.some(d => d.direction === 'hypo')\n\t\tif (hasHyper || hasHypo) {\n\t\t\tconst items: LegendRow['items'] = []\n\t\t\tif (hasHyper) items.push({ text: 'Hypermethylated', color: colors.hyper })\n\t\t\tif (hasHypo) items.push({ text: 'Hypomethylated', color: colors.hypo })\n\t\t\trows.push({ label: 'DMR', items })\n\t\t}\n\t\tif (sigCpgBedItems.length) {\n\t\t\tconst items: LegendRow['items'] = []\n\t\t\tconst hasHyperCpg = sigCpgBedItems.some(b => b.color === colors.hyper)\n\t\t\tconst hasHypoCpg = sigCpgBedItems.some(b => b.color === colors.hypo)\n\t\t\tif (hasHyperCpg) items.push({ text: 'Hyper (FDR sig.)', color: colors.hyper })\n\t\t\tif (hasHypoCpg) items.push({ text: 'Hypo (FDR sig.)', color: colors.hypo })\n\t\t\trows.push({ label: 'Sig. CpGs', items })\n\t\t}\n\t\treturn rows\n\t}\n\n\t/**\n\t * Render the per-CpG means scatter plot to an offscreen canvas and return\n\t * a data URI suitable for the bigwig imgData track.\n\t */\n\tprivate renderBetaTrack(\n\t\tdiagnostic: DmrDiagnostic,\n\t\tconfig: DmrConfig,\n\t\tblockWidth: number,\n\t\tshowLoess: boolean,\n\t\tshowDots: boolean,\n\t\tqueryStart?: number,\n\t\tqueryStop?: number\n\t): { img: { minv: number; maxv: number; src: string }; showCi: boolean } | undefined {\n\t\tconst { probes } = diagnostic\n\t\tif (!probes.positions.length) return undefined\n\n\t\tconst { colors, fdr_cutoff, minProbesForCi } = config.settings.dmr\n\t\tconst dpr = typeof window !== 'undefined' && window.devicePixelRatio > 1 ? window.devicePixelRatio : 1\n\t\tconst width = blockWidth\n\t\tconst height = 150\n\n\t\tconst canvas = document.createElement('canvas')\n\t\tcanvas.width = width * dpr\n\t\tcanvas.height = height * dpr\n\t\tconst ctx = canvas.getContext('2d')\n\t\tif (!ctx) return undefined\n\t\tctx.scale(dpr, dpr)\n\n\t\t// Use the full block view range so dots align with bedj tracks above.\n\t\t// The block stretches the image to fill the view from queryStart to queryStop.\n\t\tconst xMin = queryStart ?? probes.positions[0]\n\t\tconst xMax = queryStop ?? probes.positions[probes.positions.length - 1]\n\t\tconst xRange = xMax - xMin || 1\n\t\tconst scaleX = (val: number) => ((val - xMin) / xRange) * width\n\t\tconst scaleY = (val: number) => height - val * height // beta 0-1\n\n\t\t// Transparent background so block mouse events (yellow line) show through\n\t\tctx.clearRect(0, 0, width, height)\n\n\t\t// DMR region shading omitted \u2014 already shown as a bedj track above\n\n\t\t// Draw LOESS curves with shaded CI regions, clipped to probe data range.\n\t\t// Only show CIs when there are enough probes for a reliable estimate.\n\t\tlet showCi = false\n\t\tif (showLoess && diagnostic.loess) {\n\t\t\tconst { loess } = diagnostic\n\t\t\tconst firstProbePos = probes.positions[0]\n\t\t\tconst lastProbePos = probes.positions[probes.positions.length - 1]\n\t\t\tshowCi = probes.positions.length >= minProbesForCi\n\n\t\t\tfor (const [fitted, ciLower, ciUpper, color] of [\n\t\t\t\t[loess.group1_fitted, loess.group1_ci_lower, loess.group1_ci_upper, colors.group1],\n\t\t\t\t[loess.group2_fitted, loess.group2_ci_lower, loess.group2_ci_upper, colors.group2]\n\t\t\t] as [number[], number[], number[], string][]) {\n\t\t\t\tif (!fitted.length) continue\n\t\t\t\tconst lPos = loess.positions\n\n\t\t\t\t// Find LOESS indices within the range of actual probe positions\n\t\t\t\tlet iStart = 0\n\t\t\t\tlet iEnd = lPos.length - 1\n\t\t\t\twhile (iStart < lPos.length && lPos[iStart] < firstProbePos) iStart++\n\t\t\t\twhile (iEnd >= 0 && lPos[iEnd] > lastProbePos) iEnd--\n\t\t\t\tif (iStart > iEnd) continue\n\n\t\t\t\tif (showCi) {\n\t\t\t\t\t// Draw CI as shaded region\n\t\t\t\t\tctx.globalAlpha = 0.12\n\t\t\t\t\tctx.fillStyle = color\n\t\t\t\t\tctx.beginPath()\n\t\t\t\t\tfor (let i = iStart; i <= iEnd; i++) {\n\t\t\t\t\t\tctx.lineTo(scaleX(lPos[i]), scaleY(Math.max(0, Math.min(1, ciUpper[i]))))\n\t\t\t\t\t}\n\t\t\t\t\tfor (let i = iEnd; i >= iStart; i--) {\n\t\t\t\t\t\tctx.lineTo(scaleX(lPos[i]), scaleY(Math.max(0, Math.min(1, ciLower[i]))))\n\t\t\t\t\t}\n\t\t\t\t\tctx.closePath()\n\t\t\t\t\tctx.fill()\n\t\t\t\t}\n\n\t\t\t\t// Draw LOESS fitted curve (dashed when no CI, solid otherwise)\n\t\t\t\tctx.globalAlpha = 0.8\n\t\t\t\tctx.strokeStyle = color\n\t\t\t\tctx.lineWidth = 2\n\t\t\t\tctx.setLineDash(showCi ? [] : [6, 4])\n\t\t\t\tctx.beginPath()\n\t\t\t\tfor (let i = iStart; i <= iEnd; i++) {\n\t\t\t\t\tctx.lineTo(scaleX(lPos[i]), scaleY(Math.max(0, Math.min(1, fitted[i]))))\n\t\t\t\t}\n\t\t\t\tctx.stroke()\n\t\t\t\tctx.setLineDash([])\n\t\t\t}\n\t\t}\n\n\t\t// Draw dots (hidden for large regions where only LOESS is shown)\n\t\tif (!showDots) {\n\t\t\tctx.globalAlpha = 1\n\t\t\treturn { img: { minv: 0, maxv: 1, src: canvas.toDataURL('image/png') }, showCi }\n\t\t}\n\t\tfor (let i = 0; i < probes.positions.length; i++) {\n\t\t\tconst x = scaleX(probes.positions[i])\n\t\t\tconst isSig = probes.fdr[i] < fdr_cutoff\n\t\t\tconst alpha = isSig ? 0.85 : 0.3\n\n\t\t\t// Group 1 (control)\n\t\t\tctx.globalAlpha = alpha\n\t\t\tctx.fillStyle = colors.group1\n\t\t\tconst m1 = probes.mean_group1[i]\n\t\t\tif (m1 != null) {\n\t\t\t\tctx.beginPath()\n\t\t\t\tctx.arc(x, scaleY(m1), 4, 0, Math.PI * 2)\n\t\t\t\tctx.fill()\n\t\t\t}\n\n\t\t\t// Group 2 (case)\n\t\t\tctx.fillStyle = colors.group2\n\t\t\tconst m2 = probes.mean_group2[i]\n\t\t\tif (m2 != null) {\n\t\t\t\tctx.beginPath()\n\t\t\t\tctx.arc(x, scaleY(m2), 4, 0, Math.PI * 2)\n\t\t\t\tctx.fill()\n\t\t\t}\n\t\t}\n\t\tctx.globalAlpha = 1\n\n\t\treturn { img: { minv: 0, maxv: 1, src: canvas.toDataURL('image/png') }, showCi }\n\t}\n\n\tprivate makeDmrBedItems(dmrResult: TermdbDmrSuccessResponse, settings: DmrConfig['settings']): BedItem[] {\n\t\treturn dmrResult.dmrs.map(dmr => {\n\t\t\t// Map -log10(min_smoothed_fdr) to alpha: more significant = more opaque\n\t\t\tconst negLog = -Math.log10(Math.max(dmr.min_smoothed_fdr, 1e-300))\n\t\t\tconst alpha = Math.round(Math.min(255, Math.max(50, (negLog / 10) * 255)))\n\t\t\tconst hex = alpha.toString(16).padStart(2, '0')\n\t\t\tconst base = dmr.direction === 'hyper' ? settings.dmr.colors.hyper : settings.dmr.colors.hypo\n\t\t\treturn { chr: dmr.chr, start: dmr.start, stop: dmr.stop, color: base + hex }\n\t\t})\n\t}\n\n\tprivate makeSigCpgBedItems(\n\t\tdmrResult: TermdbDmrSuccessResponse,\n\t\tsettings: DmrConfig['settings'],\n\t\tchr: string,\n\t\tqueryStart?: number,\n\t\tqueryStop?: number\n\t): BedItem[] {\n\t\tconst diag = dmrResult.diagnostic\n\t\tif (!diag) return []\n\t\tconst { probes } = diag\n\t\tconst items: BedItem[] = []\n\t\tconst minDeltaBeta = 0.05\n\t\tfor (let i = 0; i < probes.positions.length; i++) {\n\t\t\tif (probes.fdr[i] >= settings.dmr.fdr_cutoff) continue\n\t\t\tconst pos = probes.positions[i]\n\t\t\tif (queryStart != null && queryStop != null && (pos < queryStart || pos > queryStop)) continue\n\t\t\tconst mg1 = probes.mean_group1[i]\n\t\t\tconst mg2 = probes.mean_group2[i]\n\t\t\tif (mg1 == null || mg2 == null) continue\n\t\t\tconst deltaBeta = mg2 - mg1\n\t\t\tif (Math.abs(deltaBeta) < minDeltaBeta) continue\n\t\t\tconst color = deltaBeta >= 0 ? settings.dmr.colors.hyper : settings.dmr.colors.hypo\n\t\t\titems.push({ chr, start: pos, stop: pos + 1, color })\n\t\t}\n\t\treturn items\n\t}\n}\n"],
+  "mappings": ";;;;;;;;;AAMO,SAAS,sBAAsB,MAAwB;AAC7D,QAAM,YAAY,KAAK,YAAY,CAAC;AAMpC,QAAM,KAAK,MAAM,KAAK,UAAU,cAAc,SAAS;AAKvD,QAAM,MAAM,MAAM,oBAAoB;AACtC,QAAM,eACL,IAAI,kBAAkB,aAAc,MAAM,QAAQ,IAAI,SAAS,KAAK,CAAC,CAAC,OAAO,CAAC,GAAG,UAAU,SAAS,GAAG;AACxG,QAAM,WAAW;AAAA,IAChB,YAAY;AAAA,IACZ,KAAK,eAAe,MAAU;AAAA,IAC9B,QAAQ,eAAe,MAAS;AAAA,IAChC,GAAG;AAAA,IACH,YAAY;AAAA,IACZ,QAAQ;AAAA,MACP,QAAQ;AAAA,MACR,QAAQ;AAAA,MACR,OAAO;AAAA,MACP,MAAM;AAAA,IACP;AAAA,IACA,gBAAgB;AAAA,IAChB,gBAAgB;AAAA,IAChB,SAAS;AAAA,IACT,eAAe;AAAA,EAChB;AAIA,MAAI,UAAU,QAAQ;AACrB,WAAO,OAAO,SAAS,QAAQ,UAAU,MAAM;AAC/C,WAAO,UAAU;AAAA,EAClB;AAEA,SAAO,OAAO,OAAO,UAAU,SAAS;AACzC;;;AC5CO,IAAM,kBAAkB;AAIxB,IAAM,eAAN,MAAmB;AAAA,EAGzB,YACC,WACA,QACA,WACA,UACA,YACA,WACC;AACD,UAAM,EAAE,SAAS,IAAI;AACrB,UAAM,cAAc,KAAK,gBAAgB,WAAW,QAAQ;AAC5D,UAAM,iBAAiB,KAAK,mBAAmB,WAAW,UAAU,UAAU,YAAY,SAAS;AAEnG,UAAM,UAAU,aAAa,MAAM,cAAc;AACjD,UAAM,QAAQ,UAAU,YAAY;AACpC,UAAM,YAAY,CAAC,EAAE,SAAS,MAAM,cAAc,SAAS,KAAK,MAAM,cAAc,SAAS;AAC7F,UAAM,WAAW,UAAU,SAAS,IAAI;AAExC,UAAM,kBAAkB,UAAU,aAC/B,KAAK;AAAA,MACL,UAAU;AAAA,MACV;AAAA,MACA,SAAS,IAAI;AAAA,MACb;AAAA,MACA;AAAA,MACA;AAAA,MACA;AAAA,IACA,IACA;AAEH,SAAK,WAAW;AAAA,MACf,OAAO,KAAK,eAAe,aAAa,gBAAgB,WAAW,iBAAiB,GAAG;AAAA,MACvF,YAAY,KAAK;AAAA,QAChB;AAAA,QACA,UAAU;AAAA,QACV;AAAA,QACA;AAAA,QACA;AAAA,QACA,iBAAiB,UAAU;AAAA,MAC5B;AAAA,MACA,YAAY,UAAU;AAAA,MACtB,MAAM,UAAU;AAAA,MAChB;AAAA,MACA;AAAA,MACA;AAAA,IACD;AAAA,EACD;AAAA,EAEQ,eACP,aACA,gBACA,WACA,cACQ;AACR,UAAM,QAAe,CAAC;AACtB,2BAAuB,WAAW,KAAK;AAKvC,UAAM,QAAQ,WAAW,UAAU,CAAC,GAAG,KAAK,CAAC,MAAW,EAAE,QAAQ,eAAe;AACjF,QAAI,KAAM,OAAM,KAAK,gBAAgB,IAAI,CAAC;AAC1C,UAAM,KAAK,EAAE,MAAM,QAAQ,MAAM,QAAQ,UAAU,YAAY,CAAC;AAChE,UAAM,KAAK,EAAE,MAAM,QAAQ,MAAM,aAAa,UAAU,eAAe,CAAC;AACxE,QAAI,cAAc;AACjB,YAAM,KAAK;AAAA,QACV,MAAM;AAAA,QACN,MAAM;AAAA,QACN,QAAQ;AAAA,QACR,SAAS;AAAA,MACV,CAAC;AAAA,IACF;AACA,WAAO;AAAA,EACR;AAAA,EAEQ,gBACP,QACA,MACA,gBACA,WACA,UACA,QACc;AACd,UAAM,EAAE,OAAO,IAAI,OAAO,SAAS;AACnC,UAAM,KAAK,OAAO,cAAc;AAChC,UAAM,KAAK,OAAO,cAAc;AAChC,UAAM,aAAiC,CAAC;AACxC,QAAI,UAAU;AACb,iBAAW,KAAK,EAAE,MAAM,GAAG,EAAE,cAAc,OAAO,OAAO,OAAO,GAAG,EAAE,MAAM,GAAG,EAAE,WAAW,OAAO,OAAO,OAAO,CAAC;AAAA,IAClH;AACA,QAAI,WAAW;AACd,YAAM,UAAU,SAAS,cAAc;AACvC,iBAAW;AAAA,QACV,EAAE,MAAM,GAAG,EAAE,eAAe,OAAO,IAAI,OAAO,OAAO,QAAQ,OAAO,SAAS,WAAW,SAAS;AAAA,QACjG,EAAE,MAAM,GAAG,EAAE,eAAe,OAAO,IAAI,OAAO,OAAO,QAAQ,OAAO,SAAS,WAAW,SAAS;AAAA,MAClG;AAAA,IACD;AACA,UAAM,OAAoB,CAAC,EAAE,OAAO,iBAAiB,OAAO,WAAW,CAAC;AAExE,UAAM,WAAW,KAAK,KAAK,OAAK,EAAE,cAAc,OAAO;AACvD,UAAM,UAAU,KAAK,KAAK,OAAK,EAAE,cAAc,MAAM;AACrD,QAAI,YAAY,SAAS;AACxB,YAAM,QAA4B,CAAC;AACnC,UAAI,SAAU,OAAM,KAAK,EAAE,MAAM,mBAAmB,OAAO,OAAO,MAAM,CAAC;AACzE,UAAI,QAAS,OAAM,KAAK,EAAE,MAAM,kBAAkB,OAAO,OAAO,KAAK,CAAC;AACtE,WAAK,KAAK,EAAE,OAAO,OAAO,MAAM,CAAC;AAAA,IAClC;AACA,QAAI,eAAe,QAAQ;AAC1B,YAAM,QAA4B,CAAC;AACnC,YAAM,cAAc,eAAe,KAAK,OAAK,EAAE,UAAU,OAAO,KAAK;AACrE,YAAM,aAAa,eAAe,KAAK,OAAK,EAAE,UAAU,OAAO,IAAI;AACnE,UAAI,YAAa,OAAM,KAAK,EAAE,MAAM,oBAAoB,OAAO,OAAO,MAAM,CAAC;AAC7E,UAAI,WAAY,OAAM,KAAK,EAAE,MAAM,mBAAmB,OAAO,OAAO,KAAK,CAAC;AAC1E,WAAK,KAAK,EAAE,OAAO,aAAa,MAAM,CAAC;AAAA,IACxC;AACA,WAAO;AAAA,EACR;AAAA;AAAA;AAAA;AAAA;AAAA,EAMQ,gBACP,YACA,QACA,YACA,WACA,UACA,YACA,WACoF;AACpF,UAAM,EAAE,OAAO,IAAI;AACnB,QAAI,CAAC,OAAO,UAAU,OAAQ,QAAO;AAErC,UAAM,EAAE,QAAQ,YAAY,eAAe,IAAI,OAAO,SAAS;AAC/D,UAAM,MAAM,OAAO,WAAW,eAAe,OAAO,mBAAmB,IAAI,OAAO,mBAAmB;AACrG,UAAM,QAAQ;AACd,UAAM,SAAS;AAEf,UAAM,SAAS,SAAS,cAAc,QAAQ;AAC9C,WAAO,QAAQ,QAAQ;AACvB,WAAO,SAAS,SAAS;AACzB,UAAM,MAAM,OAAO,WAAW,IAAI;AAClC,QAAI,CAAC,IAAK,QAAO;AACjB,QAAI,MAAM,KAAK,GAAG;AAIlB,UAAM,OAAO,cAAc,OAAO,UAAU,CAAC;AAC7C,UAAM,OAAO,aAAa,OAAO,UAAU,OAAO,UAAU,SAAS,CAAC;AACtE,UAAM,SAAS,OAAO,QAAQ;AAC9B,UAAM,SAAS,CAAC,SAAkB,MAAM,QAAQ,SAAU;AAC1D,UAAM,SAAS,CAAC,QAAgB,SAAS,MAAM;AAG/C,QAAI,UAAU,GAAG,GAAG,OAAO,MAAM;AAMjC,QAAI,SAAS;AACb,QAAI,aAAa,WAAW,OAAO;AAClC,YAAM,EAAE,MAAM,IAAI;AAClB,YAAM,gBAAgB,OAAO,UAAU,CAAC;AACxC,YAAM,eAAe,OAAO,UAAU,OAAO,UAAU,SAAS,CAAC;AACjE,eAAS,OAAO,UAAU,UAAU;AAEpC,iBAAW,CAAC,QAAQ,SAAS,SAAS,KAAK,KAAK;AAAA,QAC/C,CAAC,MAAM,eAAe,MAAM,iBAAiB,MAAM,iBAAiB,OAAO,MAAM;AAAA,QACjF,CAAC,MAAM,eAAe,MAAM,iBAAiB,MAAM,iBAAiB,OAAO,MAAM;AAAA,MAClF,GAA+C;AAC9C,YAAI,CAAC,OAAO,OAAQ;AACpB,cAAM,OAAO,MAAM;AAGnB,YAAI,SAAS;AACb,YAAI,OAAO,KAAK,SAAS;AACzB,eAAO,SAAS,KAAK,UAAU,KAAK,MAAM,IAAI,cAAe;AAC7D,eAAO,QAAQ,KAAK,KAAK,IAAI,IAAI,aAAc;AAC/C,YAAI,SAAS,KAAM;AAEnB,YAAI,QAAQ;AAEX,cAAI,cAAc;AAClB,cAAI,YAAY;AAChB,cAAI,UAAU;AACd,mBAAS,IAAI,QAAQ,KAAK,MAAM,KAAK;AACpC,gBAAI,OAAO,OAAO,KAAK,CAAC,CAAC,GAAG,OAAO,KAAK,IAAI,GAAG,KAAK,IAAI,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;AAAA,UACzE;AACA,mBAAS,IAAI,MAAM,KAAK,QAAQ,KAAK;AACpC,gBAAI,OAAO,OAAO,KAAK,CAAC,CAAC,GAAG,OAAO,KAAK,IAAI,GAAG,KAAK,IAAI,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;AAAA,UACzE;AACA,cAAI,UAAU;AACd,cAAI,KAAK;AAAA,QACV;AAGA,YAAI,cAAc;AAClB,YAAI,cAAc;AAClB,YAAI,YAAY;AAChB,YAAI,YAAY,SAAS,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;AACpC,YAAI,UAAU;AACd,iBAAS,IAAI,QAAQ,KAAK,MAAM,KAAK;AACpC,cAAI,OAAO,OAAO,KAAK,CAAC,CAAC,GAAG,OAAO,KAAK,IAAI,GAAG,KAAK,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;AAAA,QACxE;AACA,YAAI,OAAO;AACX,YAAI,YAAY,CAAC,CAAC;AAAA,MACnB;AAAA,IACD;AAGA,QAAI,CAAC,UAAU;AACd,UAAI,cAAc;AAClB,aAAO,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,GAAG,KAAK,OAAO,UAAU,WAAW,EAAE,GAAG,OAAO;AAAA,IAChF;AACA,aAAS,IAAI,GAAG,IAAI,OAAO,UAAU,QAAQ,KAAK;AACjD,YAAM,IAAI,OAAO,OAAO,UAAU,CAAC,CAAC;AACpC,YAAM,QAAQ,OAAO,IAAI,CAAC,IAAI;AAC9B,YAAM,QAAQ,QAAQ,OAAO;AAG7B,UAAI,cAAc;AAClB,UAAI,YAAY,OAAO;AACvB,YAAM,KAAK,OAAO,YAAY,CAAC;AAC/B,UAAI,MAAM,MAAM;AACf,YAAI,UAAU;AACd,YAAI,IAAI,GAAG,OAAO,EAAE,GAAG,GAAG,GAAG,KAAK,KAAK,CAAC;AACxC,YAAI,KAAK;AAAA,MACV;AAGA,UAAI,YAAY,OAAO;AACvB,YAAM,KAAK,OAAO,YAAY,CAAC;AAC/B,UAAI,MAAM,MAAM;AACf,YAAI,UAAU;AACd,YAAI,IAAI,GAAG,OAAO,EAAE,GAAG,GAAG,GAAG,KAAK,KAAK,CAAC;AACxC,YAAI,KAAK;AAAA,MACV;AAAA,IACD;AACA,QAAI,cAAc;AAElB,WAAO,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,GAAG,KAAK,OAAO,UAAU,WAAW,EAAE,GAAG,OAAO;AAAA,EAChF;AAAA,EAEQ,gBAAgB,WAAqC,UAA4C;AACxG,WAAO,UAAU,KAAK,IAAI,SAAO;AAEhC,YAAM,SAAS,CAAC,KAAK,MAAM,KAAK,IAAI,IAAI,kBAAkB,MAAM,CAAC;AACjE,YAAM,QAAQ,KAAK,MAAM,KAAK,IAAI,KAAK,KAAK,IAAI,IAAK,SAAS,KAAM,GAAG,CAAC,CAAC;AACzE,YAAM,MAAM,MAAM,SAAS,EAAE,EAAE,SAAS,GAAG,GAAG;AAC9C,YAAM,OAAO,IAAI,cAAc,UAAU,SAAS,IAAI,OAAO,QAAQ,SAAS,IAAI,OAAO;AACzF,aAAO,EAAE,KAAK,IAAI,KAAK,OAAO,IAAI,OAAO,MAAM,IAAI,MAAM,OAAO,OAAO,IAAI;AAAA,IAC5E,CAAC;AAAA,EACF;AAAA,EAEQ,mBACP,WACA,UACA,KACA,YACA,WACY;AACZ,UAAM,OAAO,UAAU;AACvB,QAAI,CAAC,KAAM,QAAO,CAAC;AACnB,UAAM,EAAE,OAAO,IAAI;AACnB,UAAM,QAAmB,CAAC;AAC1B,UAAM,eAAe;AACrB,aAAS,IAAI,GAAG,IAAI,OAAO,UAAU,QAAQ,KAAK;AACjD,UAAI,OAAO,IAAI,CAAC,KAAK,SAAS,IAAI,WAAY;AAC9C,YAAM,MAAM,OAAO,UAAU,CAAC;AAC9B,UAAI,cAAc,QAAQ,aAAa,SAAS,MAAM,cAAc,MAAM,WAAY;AACtF,YAAM,MAAM,OAAO,YAAY,CAAC;AAChC,YAAM,MAAM,OAAO,YAAY,CAAC;AAChC,UAAI,OAAO,QAAQ,OAAO,KAAM;AAChC,YAAM,YAAY,MAAM;AACxB,UAAI,KAAK,IAAI,SAAS,IAAI,aAAc;AACxC,YAAM,QAAQ,aAAa,IAAI,SAAS,IAAI,OAAO,QAAQ,SAAS,IAAI,OAAO;AAC/E,YAAM,KAAK,EAAE,KAAK,OAAO,KAAK,MAAM,MAAM,GAAG,MAAM,CAAC;AAAA,IACrD;AACA,WAAO;AAAA,EACR;AACD;",
+  "names": []
+}

package/dist/chunk-VFUSBU43.js

@@ -0,0 +1,14 @@
+// src/block.lazyload.js
+var Block;
+async function blocklazyload(arg) {
+  if (!Block) {
+    const b = await import("./block-XGK6TEGH.js");
+    Block = b.Block;
+  }
+  return new Block(arg);
+}
+
+export {
+  blocklazyload
+};
+//# sourceMappingURL=chunk-VFUSBU43.js.map

package/dist/chunk-VOF6NWTS.js

@@ -0,0 +1,274 @@
+import {
+  termType2label
+} from "./chunk-3XBG5HIV.js";
+import {
+  TermTypes
+} from "./chunk-SB36AUG7.js";
+import {
+  __export
+} from "./chunk-HS5PO5ZQ.js";
+
+// plots/matrix/hierCluster.renderers.js
+var hierCluster_renderers_exports = {};
+__export(hierCluster_renderers_exports, {
+  maySetSandboxHeader: () => maySetSandboxHeader,
+  plotDendrogramHclust: () => plotDendrogramHclust,
+  renderImage: () => renderImage
+});
+function maySetSandboxHeader(appState) {
+  if (!this.dom.header) return;
+  const dataType = this.config.dataType;
+  const headerText = this.config?.headerText ? `${this.config.headerText} ` : "";
+  let title;
+  if (this.config.preBuiltPlotTitle) {
+    title = this.config.preBuiltPlotTitle;
+  } else if (this.config.appName) {
+    title = `${headerText}${this.config.appName} Clustering`;
+  } else if (dataType == TermTypes.PROTEOME_ABUNDANCE) {
+    title = this.config.assayCohortTitle ? `Protein Abundance Clustering (${this.config.assayCohortTitle})` : "Protein Abundance Clustering";
+  } else {
+    title = `${headerText}${termType2label(dataType)} Clustering`;
+  }
+  this.dom.header.text(title);
+}
+function plotDendrogramHclust(plotOnly) {
+  const d = this.dimensions;
+  const s = this.config.settings.matrix;
+  const xOffset = d.seriesXoffset;
+  const pxr = window.devicePixelRatio <= 1 ? 1 : window.devicePixelRatio;
+  const obj = this.hierClusterData.clustering;
+  const row = obj.row;
+  const col = obj.col;
+  const rowHeight = this.settings.matrix.clusterRowh, { xDendrogramHeight, yDendrogramHeight } = this.settings.hierCluster, colWidth = this.dimensions.dx;
+  if (plotOnly !== "left") {
+    if (!this.settings.hierCluster.clusterSamples) {
+      this.dom.topDendrogram.selectAll("*").remove();
+    } else {
+      const height2px = getHclustHeightScalefactor(col.height, yDendrogramHeight);
+      const height = yDendrogramHeight + 1e-7;
+      const width = Math.min(colWidth * col.inputOrder.length, s.imgWMax);
+      if (width <= 0 || height <= 0) {
+        console.warn(
+          "Skipping top dendrogram render: invalid dimensions.",
+          "This may indicate a zoom feedback loop issue.",
+          {
+            width,
+            height,
+            colWidth,
+            sampleCount: col.inputOrder.length,
+            yDendrogramHeight
+          }
+        );
+        this.dom.topDendrogram.selectAll("*").remove();
+        return;
+      }
+      const canvas = new OffscreenCanvas(width * pxr, height * pxr);
+      const ctx = canvas.getContext("2d");
+      ctx.scale(pxr, pxr);
+      ctx.translate(-d.xMin, 0);
+      ctx.imageSmoothingEnabled = false;
+      ctx.imageSmoothingQuality = "high";
+      ctx.strokeStyle = "black";
+      const mergedClusters = /* @__PURE__ */ new Map();
+      for (const [clusterid0, pair] of col.merge.entries()) {
+        const clusterid = clusterid0 + 1;
+        const children = [];
+        const childrenClusters = [];
+        let x1, x2, y1, y2;
+        if (pair.n1 < 0) {
+          const [name, columnNumber] = getLeafNumber(pair.n1, col.inputOrder, col.order);
+          x1 = colWidth * (columnNumber + 0.5);
+          y1 = yDendrogramHeight;
+          children.push({ name });
+        } else {
+          if (!mergedClusters.has(pair.n1)) throw "pair.n1 is positive but not seen before";
+          const c = mergedClusters.get(pair.n1);
+          x1 = c.x;
+          y1 = c.y;
+          children.push(...c.children);
+          childrenClusters.push(pair.n1);
+        }
+        if (pair.n2 < 0) {
+          const [name, columnNumber] = getLeafNumber(pair.n2, col.inputOrder, col.order);
+          x2 = colWidth * (columnNumber + 0.5);
+          y2 = yDendrogramHeight;
+          children.push({ name });
+        } else {
+          if (!mergedClusters.has(pair.n2)) throw "pair.n1 is positive but not seen before";
+          const c = mergedClusters.get(pair.n2);
+          x2 = c.x;
+          y2 = c.y;
+          children.push(...c.children);
+          childrenClusters.push(pair.n2);
+        }
+        const clusterY = yDendrogramHeight - col.height[clusterid0].height * height2px;
+        const highlight = this.clickedClusterIds?.includes(clusterid);
+        ctx.strokeStyle = highlight ? "red" : "black";
+        ctx.beginPath();
+        ctx.moveTo(x1, y1);
+        ctx.lineTo(x1, clusterY);
+        ctx.lineTo(x2, clusterY);
+        ctx.lineTo(x2, y2);
+        ctx.stroke();
+        ctx.closePath();
+        mergedClusters.set(clusterid, {
+          x: (x1 + x2) / 2,
+          y: clusterY,
+          children,
+          childrenClusters,
+          clusterPosition: {
+            x1,
+            x2,
+            y1,
+            y2,
+            clusterY
+          }
+        });
+      }
+      this.renderImage(
+        this.api,
+        this.dom.topDendrogram,
+        canvas,
+        width,
+        height,
+        xDendrogramHeight + 0.5 * colWidth + d.xMin,
+        s.margin.top + s.scrollHeight
+      );
+      col.mergedClusters = mergedClusters;
+    }
+  }
+  if (plotOnly !== "top") {
+    if (!this.settings.hierCluster.clusterRows) {
+      this.dom.leftDendrogram.selectAll("*").remove();
+    } else {
+      const height2px = getHclustHeightScalefactor(row.height, xDendrogramHeight);
+      const width = xDendrogramHeight + 1e-7;
+      const height = rowHeight * row.inputOrder.length;
+      const canvasWidthPx = Number.isFinite(width) && Number.isFinite(pxr) ? Math.max(0, Math.floor(width * pxr)) : 0;
+      const canvasHeightPx = Number.isFinite(height) && Number.isFinite(pxr) ? Math.max(0, Math.floor(height * pxr)) : 0;
+      if (!Number.isFinite(width) || !Number.isFinite(height) || !Number.isFinite(pxr) || width <= 0 || height <= 0 || pxr <= 0 || canvasWidthPx < 1 || canvasHeightPx < 1) {
+        console.warn(
+          "Skipping left dendrogram render: invalid dimensions.",
+          "This may indicate a zoom feedback loop issue.",
+          {
+            width,
+            height,
+            pxr,
+            canvasWidthPx,
+            canvasHeightPx,
+            rowHeight,
+            termCount: row.inputOrder.length,
+            xDendrogramHeight
+          }
+        );
+        this.dom.leftDendrogram.selectAll("*").remove();
+        return;
+      }
+      const canvas = new OffscreenCanvas(canvasWidthPx, canvasHeightPx);
+      const ctx = canvas.getContext("2d");
+      ctx.scale(pxr, pxr);
+      ctx.imageSmoothingEnabled = false;
+      ctx.imageSmoothingQuality = "high";
+      ctx.strokeStyle = "black";
+      const mergedClusters = /* @__PURE__ */ new Map();
+      for (const [clusterid0, pair] of row.merge.entries()) {
+        const clusterid = clusterid0 + 1;
+        const children = [];
+        const childrenClusters = [];
+        let x1, x2, y1, y2;
+        if (pair.n1 < 0) {
+          const [name, rowNumber] = getLeafNumber(pair.n1, row.inputOrder, row.order);
+          y1 = rowHeight * (rowNumber + 0.5);
+          x1 = xDendrogramHeight;
+          children.push({ name });
+        } else {
+          if (!mergedClusters.has(pair.n1)) throw "pair.n1 is positive but not seen before";
+          const c = mergedClusters.get(pair.n1);
+          x1 = c.x;
+          y1 = c.y;
+          children.push(...c.children);
+          childrenClusters.push(pair.n1);
+        }
+        if (pair.n2 < 0) {
+          const [name, rowNumber] = getLeafNumber(pair.n2, row.inputOrder, row.order);
+          y2 = rowHeight * (rowNumber + 0.5);
+          x2 = xDendrogramHeight;
+          children.push({ name });
+        } else {
+          if (!mergedClusters.has(pair.n2)) throw "pair.n1 is positive but not seen before";
+          const c = mergedClusters.get(pair.n2);
+          x2 = c.x;
+          y2 = c.y;
+          children.push(...c.children);
+          childrenClusters.push(pair.n2);
+        }
+        const clusterX = xDendrogramHeight - row.height[clusterid0].height * height2px;
+        const highlight = this.clickedLeftClusterIds?.includes(clusterid);
+        ctx.strokeStyle = highlight ? "red" : "black";
+        ctx.beginPath();
+        ctx.moveTo(x1, y1);
+        ctx.lineTo(clusterX, y1);
+        ctx.lineTo(clusterX, y2);
+        ctx.lineTo(x2, y2);
+        ctx.stroke();
+        ctx.closePath();
+        mergedClusters.set(clusterid, {
+          x: clusterX,
+          y: (y1 + y2) / 2,
+          children,
+          childrenClusters,
+          clusterPosition: {
+            x1,
+            x2,
+            y1,
+            y2,
+            clusterX
+          }
+        });
+      }
+      const t = this.termOrder.find((t2) => t2.grp.type == "hierCluster" || t2.grp.name == this.hcTermGroup.name);
+      const y = (
+        // t.labelOffset is commented out because it is already handled in adjustSvgDimensions
+        t.grpIndex * s.rowgspace + t.prevGrpTotalIndex * s.rowh + t.totalHtAdjustments + s.margin.top + s.scrollHeight + // left dendrogram image must be lower than the top dendrogram image height
+        yDendrogramHeight
+      );
+      this.renderImage(this.api, this.dom.leftDendrogram, canvas, width, height, 0, y);
+      row.mergedClusters = mergedClusters;
+    }
+  }
+}
+async function renderImage(componentApi, g, canvas, width, height, x, y) {
+  const sequenceId = componentApi.getSequenceId();
+  const reader = new FileReader();
+  reader.addEventListener(
+    "load",
+    () => {
+      if (componentApi.isStaleSequenceId(sequenceId)) return;
+      g.selectAll("*").remove();
+      g.append("image").attr("x", x + 0.033).attr("y", y + 0.033).attr("xlink:href", reader.result).attr("width", width).attr("height", height);
+    },
+    false
+  );
+  const blob = await canvas.convertToBlob({ quality: 1 });
+  reader.readAsDataURL(blob);
+}
+function getHclustHeightScalefactor(lst, ph) {
+  let max = lst[0].height;
+  for (const h of lst) max = Math.max(max, h.height);
+  return ph / max;
+}
+function getLeafNumber(minus, inputOrder, order) {
+  const name = inputOrder[-minus - 1];
+  if (!name) throw "minus not in inputOrder";
+  const i = order.findIndex((j) => j.name == name);
+  if (i == -1) throw "name not found in hc$order";
+  return [name, i];
+}
+
+export {
+  maySetSandboxHeader,
+  plotDendrogramHclust,
+  renderImage,
+  hierCluster_renderers_exports
+};
+//# sourceMappingURL=chunk-VOF6NWTS.js.map