@sjcrh/proteinpaint-client 2.204.0 → 2.205.0

package/dist/brainRegions-NTEAXNZJ.js

@@ -0,0 +1,234 @@
+import {
+  loadBrainAssets,
+  makeDiseaseTabs,
+  renderBrainSvg
+} from "./chunk-2Z4ZSINZ.js";
+import {
+  PlotBase,
+  addGeneSearchbox
+} from "./chunk-J7JDCNLU.js";
+import "./chunk-HJ6L54YS.js";
+import "./chunk-KV4W2ACA.js";
+import "./chunk-N635HDJ4.js";
+import "./chunk-T46FA72N.js";
+import {
+  Menu
+} from "./chunk-ELJX3QIQ.js";
+import "./chunk-EEB5VE2A.js";
+import "./chunk-6RRZRISL.js";
+import "./chunk-2KM4PRQM.js";
+import {
+  dofetch3
+} from "./chunk-JYOIO5UY.js";
+import "./chunk-YN5NY3D3.js";
+import "./chunk-RN4BOWRH.js";
+import "./chunk-IK2BO37K.js";
+import {
+  copyMerge,
+  getCompInit
+} from "./chunk-WINIL2KN.js";
+import "./chunk-PF4DSFDR.js";
+import "./chunk-D6G64XPJ.js";
+import "./chunk-W5J3LTYS.js";
+import "./chunk-YLJOZP4P.js";
+import {
+  linear
+} from "./chunk-HDTFYTEL.js";
+import "./chunk-FXQXCOII.js";
+import "./chunk-TLT4YIG3.js";
+import "./chunk-5R63Q5KH.js";
+import "./chunk-I6Y4O3RR.js";
+import "./chunk-Q5RDQNIT.js";
+import "./chunk-DQC5FFGV.js";
+import "./chunk-HS5PO5ZQ.js";
+
+// plots/brainRegions.ts
+var defaultConfig = {
+  chartType: "brainRegions"
+};
+var P_VALUE_THRESHOLD = 0.05;
+var NONSIG_COLOR = "#ccc";
+var BRAIN_RENDER_W = 520;
+var gradientSeq = 0;
+var BrainRegions = class _BrainRegions extends PlotBase {
+  static {
+    this.type = "brainRegions";
+  }
+  constructor(opts, api) {
+    super(opts, api);
+    this.type = _BrainRegions.type;
+  }
+  async init() {
+    const holder = this.opts.holder.append("div").style("padding", "10px");
+    this.dom = {
+      holder,
+      body: holder.append("div"),
+      tip: new Menu({ padding: "" }),
+      header: this.opts.header
+    };
+    if (this.dom.header) this.dom.header.html("Brain Regional Proteome");
+  }
+  getState(appState) {
+    const config = appState.plots.find((p) => p.id === this.id);
+    if (!config) throw `No plot with id='${this.id}' found`;
+    return { config };
+  }
+  async main() {
+    const gene = this.state.config?.gene;
+    if (!gene) throw new Error("brainRegions: gene is missing");
+    if (this.dom.header) this.dom.header.text(`Brain Regional Proteome: ${gene}`);
+    const body = {
+      genome: this.app.opts.state.vocab.genome,
+      dslabel: this.app.opts.state.vocab.dslabel,
+      gene
+    };
+    const data = await dofetch3("termdb/brainRegions", { body });
+    if (data.error) throw data.error;
+    this.dom.body.selectAll("*").remove();
+    const description = this.app.vocabApi.termdbConfig?.queries?.proteome?.brainRegions?.description;
+    if (description) {
+      this.dom.body.append("div").style("font-size", "0.85em").style("color", "#555").style("margin-bottom", "10px").style("line-height", "1.4").style("max-width", "600px").style("white-space", "normal").style("overflow-wrap", "break-word").text(description);
+    }
+    const isoformIds = Object.keys(data.isoforms);
+    if (isoformIds.length === 0) {
+      this.dom.body.append("div").style("padding", "20px").style("color", "#666").text(`No brain-region data found for gene "${gene}".`);
+      return;
+    }
+    const brainAssets = await loadBrainAssets(data.svgUrl, Object.keys(data.regions));
+    const controlRow = this.dom.body.append("div").style("margin-bottom", "15px");
+    controlRow.append("span").style("font-weight", "bold").text("Isoform: ");
+    let selectedIsoform = isoformIds[0];
+    let selectedDisease = data.diseases[0];
+    const tabsHolder = this.dom.body.append("div");
+    const redraw = () => this.renderBrains(data, selectedIsoform, selectedDisease, brainAssets);
+    if (data.diseases.length > 1) {
+      makeDiseaseTabs(
+        tabsHolder,
+        data.diseases,
+        selectedDisease,
+        (d) => {
+          selectedDisease = d;
+          redraw();
+        },
+        ".9em"
+      );
+    }
+    if (isoformIds.length > 1) {
+      const sel = controlRow.append("select").style("margin-left", "5px").style("padding", "3px 6px").on("change", () => {
+        selectedIsoform = sel.node().value;
+        redraw();
+      });
+      sel.selectAll("option").data(isoformIds).enter().append("option").attr("value", (d) => d).text((d) => `${data.isoforms[d].gene_name} \u2014 ${d}`);
+    } else {
+      controlRow.append("span").style("margin-left", "5px").text(`${data.isoforms[selectedIsoform].gene_name} \u2014 ${selectedIsoform}`);
+    }
+    redraw();
+  }
+  renderBrains(data, selectedIsoform, selectedDisease, brainAssets) {
+    const existing = this.dom.body.select(".sjpp-brain-regions-container");
+    if (!existing.empty()) existing.remove();
+    const container = this.dom.body.append("div").attr("class", "sjpp-brain-regions-container").style("display", "flex").style("gap", "40px").style("flex-wrap", "wrap");
+    const isoformData = data.isoforms[selectedIsoform];
+    if (!isoformData) return;
+    const regionData = isoformData.data[selectedDisease] || {};
+    const allFCs = [];
+    for (const entry of Object.values(regionData)) {
+      if (entry.p_value < P_VALUE_THRESHOLD) allFCs.push(entry.fold_change);
+    }
+    const maxAbsFC = allFCs.length > 0 ? Math.max(...allFCs.map((v) => Math.abs(v))) : 1;
+    const colorScale = linear().domain([-maxAbsFC, 0, maxAbsFC]).range(["#2166ac", "#f7f7f7", "#b2182b"]).clamp(true);
+    renderBrainSvg({
+      holder: container,
+      width: BRAIN_RENDER_W,
+      templateUrl: data.templateUrl,
+      assets: brainAssets,
+      regions: data.regions,
+      title: selectedDisease,
+      tip: this.dom.tip,
+      fillByRegion: (code) => {
+        const entry = regionData[code];
+        if (entry && entry.p_value < P_VALUE_THRESHOLD) return colorScale(entry.fold_change);
+        return NONSIG_COLOR;
+      },
+      tooltipByRegion: (code, label) => {
+        const entry = regionData[code];
+        if (!entry) return `${label} (${code})
+No data`;
+        const fc = entry.fold_change.toFixed(4);
+        const fmt = (v) => v >= 1e-4 ? v.toFixed(4) : v.toExponential(3);
+        const fdr = Number.isFinite(entry.fdr) ? `
+FDR: ${fmt(entry.fdr)}` : "";
+        return `${label} (${code})
+log\u2082 fold change: ${fc}
+p-value: ${fmt(entry.p_value)}${fdr}`;
+      }
+    });
+    this.renderLegend(container, colorScale, maxAbsFC, allFCs.length, selectedDisease);
+  }
+  renderLegend(container, colorScale, maxAbsFC, nSig, disease) {
+    const legendDiv = container.append("div").style("display", "flex").style("flex-direction", "column").style("justify-content", "center").style("padding", "10px");
+    if (!nSig) {
+      legendDiv.append("div").style("font-size", "13px").style("color", "#666").style("max-width", "220px").style("line-height", "1.4").html(
+        `<span style="display:inline-block;width:14px;height:14px;background:${NONSIG_COLOR};border:1px solid #999;vertical-align:middle;margin-right:4px"></span> No region reaches p &lt; ${P_VALUE_THRESHOLD} for this isoform in ${disease}.`
+      );
+      return;
+    }
+    legendDiv.append("div").style("font-weight", "bold").style("font-size", "13px").style("margin-bottom", "8px").text("Fold Change (log\u2082)");
+    const legendWidth = 20;
+    const legendHeight = 200;
+    const svg = legendDiv.append("svg").attr("width", legendWidth + 60).attr("height", legendHeight + 30);
+    const defs = svg.append("defs");
+    const gradientId = `brain-fc-gradient-${gradientSeq++}`;
+    const gradient = defs.append("linearGradient").attr("id", gradientId).attr("x1", "0").attr("y1", "0").attr("x2", "0").attr("y2", "1");
+    const steps = 10;
+    for (let i = 0; i <= steps; i++) {
+      const t = i / steps;
+      const val = maxAbsFC * (1 - 2 * t);
+      gradient.append("stop").attr("offset", `${t * 100}%`).attr("stop-color", colorScale(val));
+    }
+    svg.append("rect").attr("x", 0).attr("y", 10).attr("width", legendWidth).attr("height", legendHeight).style("fill", `url(#${gradientId})`).attr("stroke", "#999");
+    const legendScale = linear().domain([maxAbsFC, -maxAbsFC]).range([10, legendHeight + 10]);
+    const ticks = [-maxAbsFC, -maxAbsFC / 2, 0, maxAbsFC / 2, maxAbsFC];
+    for (const tick of ticks) {
+      const y = legendScale(tick);
+      svg.append("line").attr("x1", legendWidth).attr("y1", y).attr("x2", legendWidth + 5).attr("y2", y).attr("stroke", "#666");
+      svg.append("text").attr("x", legendWidth + 8).attr("y", y).attr("dominant-baseline", "central").attr("font-size", "10px").text(tick.toFixed(2));
+    }
+    legendDiv.append("div").style("margin-top", "10px").style("font-size", "12px").style("color", "#666").html(
+      `<span style="display:inline-block;width:14px;height:14px;background:${NONSIG_COLOR};border:1px solid #999;vertical-align:middle;margin-right:4px"></span> Not significant (p \u2265 0.05)`
+    );
+  }
+};
+var componentInit = getCompInit(BrainRegions);
+async function getPlotConfig(opts) {
+  const config = structuredClone(defaultConfig);
+  if (!opts.gene) throw new Error("brainRegions requires opts.gene");
+  return copyMerge(config, opts);
+}
+function makeChartBtnMenu(holder, chartsInstance) {
+  const row = holder.append("div").style("padding", "5px");
+  row.append("span").style("font-weight", "bold").text("Enter a gene name:");
+  const geneSearch = addGeneSearchbox({
+    row,
+    genome: chartsInstance.app.opts.genome,
+    tip: new Menu({ padding: "0px" }),
+    searchOnly: "gene",
+    callback: async () => {
+      if (!geneSearch.geneSymbol) throw new Error("A valid gene selection is required");
+      chartsInstance.dom.tip.hide();
+      chartsInstance.app.dispatch({
+        type: "plot_create",
+        config: {
+          chartType: "brainRegions",
+          gene: geneSearch.geneSymbol
+        }
+      });
+    }
+  });
+}
+export {
+  componentInit,
+  getPlotConfig,
+  makeChartBtnMenu
+};
+//# sourceMappingURL=brainRegions-NTEAXNZJ.js.map

package/dist/brainRegions-NTEAXNZJ.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../plots/brainRegions.ts"],
+  "sourcesContent": ["import type { MassState, BasePlotConfig } from '#mass/types/mass'\nimport { getCompInit, copyMerge, type RxComponent, type ComponentApi } from '#rx'\nimport { PlotBase } from './PlotBase'\nimport { Menu, addGeneSearchbox } from '#dom'\nimport { dofetch3 } from '#common/dofetch'\nimport { scaleLinear } from 'd3'\nimport { loadBrainAssets, renderBrainSvg, makeDiseaseTabs, type BrainAssets } from './brainRegions.svg'\n\nconst defaultConfig = {\n\tchartType: 'brainRegions'\n}\n\nconst P_VALUE_THRESHOLD = 0.05\nconst NONSIG_COLOR = '#ccc'\n\nconst BRAIN_RENDER_W = 520\n\n// Monotonic counter for unique <linearGradient> ids (Date.now() can collide when\n// two legends render within the same millisecond).\nlet gradientSeq = 0\n\nclass BrainRegions extends PlotBase implements RxComponent {\n\tstatic type = 'brainRegions'\n\ttype: string\n\tdom!: {\n\t\tholder: any\n\t\tbody: any\n\t\ttip: Menu\n\t\theader?: any\n\t}\n\n\tconstructor(opts: any, api: ComponentApi) {\n\t\tsuper(opts, api)\n\t\tthis.type = BrainRegions.type\n\t}\n\n\tasync init() {\n\t\tconst holder = this.opts.holder.append('div').style('padding', '10px')\n\t\tthis.dom = {\n\t\t\tholder,\n\t\t\tbody: holder.append('div'),\n\t\t\ttip: new Menu({ padding: '' }),\n\t\t\theader: this.opts.header\n\t\t}\n\t\tif (this.dom.header) this.dom.header.html('Brain Regional Proteome')\n\t}\n\n\tgetState(appState: MassState) {\n\t\tconst config: any = appState.plots.find((p: BasePlotConfig) => p.id === this.id)\n\t\tif (!config) throw `No plot with id='${this.id}' found`\n\t\treturn { config }\n\t}\n\n\tasync main() {\n\t\tconst gene = this.state.config?.gene\n\t\tif (!gene) throw new Error('brainRegions: gene is missing')\n\n\t\tif (this.dom.header) this.dom.header.text(`Brain Regional Proteome: ${gene}`)\n\n\t\tconst body = {\n\t\t\tgenome: this.app.opts.state.vocab.genome,\n\t\t\tdslabel: this.app.opts.state.vocab.dslabel,\n\t\t\tgene\n\t\t}\n\n\t\tconst data = await dofetch3('termdb/brainRegions', { body })\n\t\tif (data.error) throw data.error\n\n\t\tthis.dom.body.selectAll('*').remove()\n\n\t\t// Intro paragraph (config-driven), styled like the gene-ranking description note.\n\t\tconst description = this.app.vocabApi.termdbConfig?.queries?.proteome?.brainRegions?.description\n\t\tif (description) {\n\t\t\tthis.dom.body\n\t\t\t\t.append('div')\n\t\t\t\t.style('font-size', '0.85em')\n\t\t\t\t.style('color', '#555')\n\t\t\t\t.style('margin-bottom', '10px')\n\t\t\t\t.style('line-height', '1.4')\n\t\t\t\t.style('max-width', '600px')\n\t\t\t\t.style('white-space', 'normal')\n\t\t\t\t.style('overflow-wrap', 'break-word')\n\t\t\t\t.text(description)\n\t\t}\n\n\t\tconst isoformIds = Object.keys(data.isoforms)\n\t\tif (isoformIds.length === 0) {\n\t\t\tthis.dom.body\n\t\t\t\t.append('div')\n\t\t\t\t.style('padding', '20px')\n\t\t\t\t.style('color', '#666')\n\t\t\t\t.text(`No brain-region data found for gene \"${gene}\".`)\n\t\t\treturn\n\t\t}\n\n\t\tconst brainAssets = await loadBrainAssets(data.svgUrl, Object.keys(data.regions))\n\n\t\tconst controlRow = this.dom.body.append('div').style('margin-bottom', '15px')\n\t\tcontrolRow.append('span').style('font-weight', 'bold').text('Isoform: ')\n\n\t\tlet selectedIsoform = isoformIds[0]\n\t\t// one brain at a time, disease chosen via tabs (same UI as the\n\t\t// proteinView brain-region tile)\n\t\tlet selectedDisease: string = data.diseases[0]\n\t\tconst tabsHolder = this.dom.body.append('div')\n\t\tconst redraw = () => this.renderBrains(data, selectedIsoform, selectedDisease, brainAssets)\n\t\tif (data.diseases.length > 1) {\n\t\t\tmakeDiseaseTabs(\n\t\t\t\ttabsHolder,\n\t\t\t\tdata.diseases,\n\t\t\t\tselectedDisease,\n\t\t\t\t(d: string) => {\n\t\t\t\t\tselectedDisease = d\n\t\t\t\t\tredraw()\n\t\t\t\t},\n\t\t\t\t'.9em'\n\t\t\t)\n\t\t}\n\n\t\tif (isoformIds.length > 1) {\n\t\t\tconst sel = controlRow\n\t\t\t\t.append('select')\n\t\t\t\t.style('margin-left', '5px')\n\t\t\t\t.style('padding', '3px 6px')\n\t\t\t\t.on('change', () => {\n\t\t\t\t\tselectedIsoform = sel.node().value\n\t\t\t\t\tredraw()\n\t\t\t\t})\n\n\t\t\tsel\n\t\t\t\t.selectAll('option')\n\t\t\t\t.data(isoformIds)\n\t\t\t\t.enter()\n\t\t\t\t.append('option')\n\t\t\t\t.attr('value', (d: string) => d)\n\t\t\t\t.text((d: string) => `${data.isoforms[d].gene_name} \u2014 ${d}`)\n\t\t} else {\n\t\t\tcontrolRow\n\t\t\t\t.append('span')\n\t\t\t\t.style('margin-left', '5px')\n\t\t\t\t.text(`${data.isoforms[selectedIsoform].gene_name} \u2014 ${selectedIsoform}`)\n\t\t}\n\n\t\tredraw()\n\t}\n\n\trenderBrains(data: any, selectedIsoform: string, selectedDisease: string, brainAssets: BrainAssets) {\n\t\tconst existing = this.dom.body.select('.sjpp-brain-regions-container')\n\t\tif (!existing.empty()) existing.remove()\n\n\t\tconst container = this.dom.body\n\t\t\t.append('div')\n\t\t\t.attr('class', 'sjpp-brain-regions-container')\n\t\t\t.style('display', 'flex')\n\t\t\t.style('gap', '40px')\n\t\t\t.style('flex-wrap', 'wrap')\n\n\t\tconst isoformData = data.isoforms[selectedIsoform]\n\t\tif (!isoformData) return\n\n\t\t// color scale is scoped to the selected disease only\n\t\tconst regionData = isoformData.data[selectedDisease] || {}\n\t\tconst allFCs: number[] = []\n\t\tfor (const entry of Object.values(regionData) as any[]) {\n\t\t\tif (entry.p_value < P_VALUE_THRESHOLD) allFCs.push(entry.fold_change)\n\t\t}\n\n\t\tconst maxAbsFC = allFCs.length > 0 ? Math.max(...allFCs.map(v => Math.abs(v))) : 1\n\t\tconst colorScale = scaleLinear<string>()\n\t\t\t.domain([-maxAbsFC, 0, maxAbsFC])\n\t\t\t.range(['#2166ac', '#f7f7f7', '#b2182b'])\n\t\t\t.clamp(true)\n\n\t\trenderBrainSvg({\n\t\t\tholder: container,\n\t\t\twidth: BRAIN_RENDER_W,\n\t\t\ttemplateUrl: data.templateUrl,\n\t\t\tassets: brainAssets,\n\t\t\tregions: data.regions,\n\t\t\ttitle: selectedDisease,\n\t\t\ttip: this.dom.tip,\n\t\t\tfillByRegion: (code: string) => {\n\t\t\t\tconst entry = regionData[code]\n\t\t\t\tif (entry && entry.p_value < P_VALUE_THRESHOLD) return colorScale(entry.fold_change) as string\n\t\t\t\treturn NONSIG_COLOR\n\t\t\t},\n\t\t\ttooltipByRegion: (code: string, label: string) => {\n\t\t\t\tconst entry = regionData[code]\n\t\t\t\tif (!entry) return `${label} (${code})\\nNo data`\n\t\t\t\tconst fc = entry.fold_change.toFixed(4)\n\t\t\t\tconst fmt = (v: number) => (v >= 0.0001 ? v.toFixed(4) : v.toExponential(3))\n\t\t\t\tconst fdr = Number.isFinite(entry.fdr) ? `\\nFDR: ${fmt(entry.fdr)}` : ''\n\t\t\t\treturn `${label} (${code})\\nlog\u2082 fold change: ${fc}\\np-value: ${fmt(entry.p_value)}${fdr}`\n\t\t\t}\n\t\t})\n\n\t\tthis.renderLegend(container, colorScale, maxAbsFC, allFCs.length, selectedDisease)\n\t}\n\n\trenderLegend(container: any, colorScale: any, maxAbsFC: number, nSig: number, disease: string) {\n\t\tconst legendDiv = container\n\t\t\t.append('div')\n\t\t\t.style('display', 'flex')\n\t\t\t.style('flex-direction', 'column')\n\t\t\t.style('justify-content', 'center')\n\t\t\t.style('padding', '10px')\n\n\t\tif (!nSig) {\n\t\t\t// no region passes the threshold: a gradient would only show the \u00B11 fallback\n\t\t\t// domain, which means nothing. Say so instead.\n\t\t\tlegendDiv\n\t\t\t\t.append('div')\n\t\t\t\t.style('font-size', '13px')\n\t\t\t\t.style('color', '#666')\n\t\t\t\t.style('max-width', '220px')\n\t\t\t\t.style('line-height', '1.4')\n\t\t\t\t.html(\n\t\t\t\t\t`<span style=\"display:inline-block;width:14px;height:14px;background:${NONSIG_COLOR};border:1px solid #999;vertical-align:middle;margin-right:4px\"></span> No region reaches p &lt; ${P_VALUE_THRESHOLD} for this isoform in ${disease}.`\n\t\t\t\t)\n\t\t\treturn\n\t\t}\n\n\t\tlegendDiv\n\t\t\t.append('div')\n\t\t\t.style('font-weight', 'bold')\n\t\t\t.style('font-size', '13px')\n\t\t\t.style('margin-bottom', '8px')\n\t\t\t.text('Fold Change (log\u2082)')\n\n\t\tconst legendWidth = 20\n\t\tconst legendHeight = 200\n\t\tconst svg = legendDiv\n\t\t\t.append('svg')\n\t\t\t.attr('width', legendWidth + 60)\n\t\t\t.attr('height', legendHeight + 30)\n\n\t\tconst defs = svg.append('defs')\n\t\tconst gradientId = `brain-fc-gradient-${gradientSeq++}`\n\t\tconst gradient = defs\n\t\t\t.append('linearGradient')\n\t\t\t.attr('id', gradientId)\n\t\t\t.attr('x1', '0')\n\t\t\t.attr('y1', '0')\n\t\t\t.attr('x2', '0')\n\t\t\t.attr('y2', '1')\n\n\t\tconst steps = 10\n\t\tfor (let i = 0; i <= steps; i++) {\n\t\t\tconst t = i / steps\n\t\t\tconst val = maxAbsFC * (1 - 2 * t)\n\t\t\tgradient\n\t\t\t\t.append('stop')\n\t\t\t\t.attr('offset', `${t * 100}%`)\n\t\t\t\t.attr('stop-color', colorScale(val))\n\t\t}\n\n\t\tsvg\n\t\t\t.append('rect')\n\t\t\t.attr('x', 0)\n\t\t\t.attr('y', 10)\n\t\t\t.attr('width', legendWidth)\n\t\t\t.attr('height', legendHeight)\n\t\t\t.style('fill', `url(#${gradientId})`)\n\t\t\t.attr('stroke', '#999')\n\n\t\tconst legendScale = scaleLinear()\n\t\t\t.domain([maxAbsFC, -maxAbsFC])\n\t\t\t.range([10, legendHeight + 10])\n\n\t\tconst ticks = [-maxAbsFC, -maxAbsFC / 2, 0, maxAbsFC / 2, maxAbsFC]\n\t\tfor (const tick of ticks) {\n\t\t\tconst y = legendScale(tick)\n\t\t\tsvg\n\t\t\t\t.append('line')\n\t\t\t\t.attr('x1', legendWidth)\n\t\t\t\t.attr('y1', y)\n\t\t\t\t.attr('x2', legendWidth + 5)\n\t\t\t\t.attr('y2', y)\n\t\t\t\t.attr('stroke', '#666')\n\t\t\tsvg\n\t\t\t\t.append('text')\n\t\t\t\t.attr('x', legendWidth + 8)\n\t\t\t\t.attr('y', y)\n\t\t\t\t.attr('dominant-baseline', 'central')\n\t\t\t\t.attr('font-size', '10px')\n\t\t\t\t.text(tick.toFixed(2))\n\t\t}\n\n\t\tlegendDiv\n\t\t\t.append('div')\n\t\t\t.style('margin-top', '10px')\n\t\t\t.style('font-size', '12px')\n\t\t\t.style('color', '#666')\n\t\t\t.html(\n\t\t\t\t`<span style=\"display:inline-block;width:14px;height:14px;background:${NONSIG_COLOR};border:1px solid #999;vertical-align:middle;margin-right:4px\"></span> Not significant (p \u2265 0.05)`\n\t\t\t)\n\t}\n}\n\nexport const componentInit = getCompInit(BrainRegions)\n\nexport async function getPlotConfig(opts: any) {\n\tconst config = structuredClone(defaultConfig)\n\tif (!opts.gene) throw new Error('brainRegions requires opts.gene')\n\treturn copyMerge(config, opts)\n}\n\nexport function makeChartBtnMenu(holder: any, chartsInstance: any) {\n\tconst row = holder.append('div').style('padding', '5px')\n\trow.append('span').style('font-weight', 'bold').text('Enter a gene name:')\n\n\tconst geneSearch = addGeneSearchbox({\n\t\trow,\n\t\tgenome: chartsInstance.app.opts.genome,\n\t\ttip: new Menu({ padding: '0px' }),\n\t\tsearchOnly: 'gene',\n\t\tcallback: async () => {\n\t\t\tif (!geneSearch.geneSymbol) throw new Error('A valid gene selection is required')\n\t\t\tchartsInstance.dom.tip.hide()\n\t\t\tchartsInstance.app.dispatch({\n\t\t\t\ttype: 'plot_create',\n\t\t\t\tconfig: {\n\t\t\t\t\tchartType: 'brainRegions',\n\t\t\t\t\tgene: geneSearch.geneSymbol\n\t\t\t\t}\n\t\t\t})\n\t\t}\n\t})\n}\n"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAQA,IAAM,gBAAgB;AAAA,EACrB,WAAW;AACZ;AAEA,IAAM,oBAAoB;AAC1B,IAAM,eAAe;AAErB,IAAM,iBAAiB;AAIvB,IAAI,cAAc;AAElB,IAAM,eAAN,MAAM,sBAAqB,SAAgC;AAAA,EAC1D;AAAA,SAAO,OAAO;AAAA;AAAA,EASd,YAAY,MAAW,KAAmB;AACzC,UAAM,MAAM,GAAG;AACf,SAAK,OAAO,cAAa;AAAA,EAC1B;AAAA,EAEA,MAAM,OAAO;AACZ,UAAM,SAAS,KAAK,KAAK,OAAO,OAAO,KAAK,EAAE,MAAM,WAAW,MAAM;AACrE,SAAK,MAAM;AAAA,MACV;AAAA,MACA,MAAM,OAAO,OAAO,KAAK;AAAA,MACzB,KAAK,IAAI,KAAK,EAAE,SAAS,GAAG,CAAC;AAAA,MAC7B,QAAQ,KAAK,KAAK;AAAA,IACnB;AACA,QAAI,KAAK,IAAI,OAAQ,MAAK,IAAI,OAAO,KAAK,yBAAyB;AAAA,EACpE;AAAA,EAEA,SAAS,UAAqB;AAC7B,UAAM,SAAc,SAAS,MAAM,KAAK,CAAC,MAAsB,EAAE,OAAO,KAAK,EAAE;AAC/E,QAAI,CAAC,OAAQ,OAAM,oBAAoB,KAAK,EAAE;AAC9C,WAAO,EAAE,OAAO;AAAA,EACjB;AAAA,EAEA,MAAM,OAAO;AACZ,UAAM,OAAO,KAAK,MAAM,QAAQ;AAChC,QAAI,CAAC,KAAM,OAAM,IAAI,MAAM,+BAA+B;AAE1D,QAAI,KAAK,IAAI,OAAQ,MAAK,IAAI,OAAO,KAAK,4BAA4B,IAAI,EAAE;AAE5E,UAAM,OAAO;AAAA,MACZ,QAAQ,KAAK,IAAI,KAAK,MAAM,MAAM;AAAA,MAClC,SAAS,KAAK,IAAI,KAAK,MAAM,MAAM;AAAA,MACnC;AAAA,IACD;AAEA,UAAM,OAAO,MAAM,SAAS,uBAAuB,EAAE,KAAK,CAAC;AAC3D,QAAI,KAAK,MAAO,OAAM,KAAK;AAE3B,SAAK,IAAI,KAAK,UAAU,GAAG,EAAE,OAAO;AAGpC,UAAM,cAAc,KAAK,IAAI,SAAS,cAAc,SAAS,UAAU,cAAc;AACrF,QAAI,aAAa;AAChB,WAAK,IAAI,KACP,OAAO,KAAK,EACZ,MAAM,aAAa,QAAQ,EAC3B,MAAM,SAAS,MAAM,EACrB,MAAM,iBAAiB,MAAM,EAC7B,MAAM,eAAe,KAAK,EAC1B,MAAM,aAAa,OAAO,EAC1B,MAAM,eAAe,QAAQ,EAC7B,MAAM,iBAAiB,YAAY,EACnC,KAAK,WAAW;AAAA,IACnB;AAEA,UAAM,aAAa,OAAO,KAAK,KAAK,QAAQ;AAC5C,QAAI,WAAW,WAAW,GAAG;AAC5B,WAAK,IAAI,KACP,OAAO,KAAK,EACZ,MAAM,WAAW,MAAM,EACvB,MAAM,SAAS,MAAM,EACrB,KAAK,wCAAwC,IAAI,IAAI;AACvD;AAAA,IACD;AAEA,UAAM,cAAc,MAAM,gBAAgB,KAAK,QAAQ,OAAO,KAAK,KAAK,OAAO,CAAC;AAEhF,UAAM,aAAa,KAAK,IAAI,KAAK,OAAO,KAAK,EAAE,MAAM,iBAAiB,MAAM;AAC5E,eAAW,OAAO,MAAM,EAAE,MAAM,eAAe,MAAM,EAAE,KAAK,WAAW;AAEvE,QAAI,kBAAkB,WAAW,CAAC;AAGlC,QAAI,kBAA0B,KAAK,SAAS,CAAC;AAC7C,UAAM,aAAa,KAAK,IAAI,KAAK,OAAO,KAAK;AAC7C,UAAM,SAAS,MAAM,KAAK,aAAa,MAAM,iBAAiB,iBAAiB,WAAW;AAC1F,QAAI,KAAK,SAAS,SAAS,GAAG;AAC7B;AAAA,QACC;AAAA,QACA,KAAK;AAAA,QACL;AAAA,QACA,CAAC,MAAc;AACd,4BAAkB;AAClB,iBAAO;AAAA,QACR;AAAA,QACA;AAAA,MACD;AAAA,IACD;AAEA,QAAI,WAAW,SAAS,GAAG;AAC1B,YAAM,MAAM,WACV,OAAO,QAAQ,EACf,MAAM,eAAe,KAAK,EAC1B,MAAM,WAAW,SAAS,EAC1B,GAAG,UAAU,MAAM;AACnB,0BAAkB,IAAI,KAAK,EAAE;AAC7B,eAAO;AAAA,MACR,CAAC;AAEF,UACE,UAAU,QAAQ,EAClB,KAAK,UAAU,EACf,MAAM,EACN,OAAO,QAAQ,EACf,KAAK,SAAS,CAAC,MAAc,CAAC,EAC9B,KAAK,CAAC,MAAc,GAAG,KAAK,SAAS,CAAC,EAAE,SAAS,WAAM,CAAC,EAAE;AAAA,IAC7D,OAAO;AACN,iBACE,OAAO,MAAM,EACb,MAAM,eAAe,KAAK,EAC1B,KAAK,GAAG,KAAK,SAAS,eAAe,EAAE,SAAS,WAAM,eAAe,EAAE;AAAA,IAC1E;AAEA,WAAO;AAAA,EACR;AAAA,EAEA,aAAa,MAAW,iBAAyB,iBAAyB,aAA0B;AACnG,UAAM,WAAW,KAAK,IAAI,KAAK,OAAO,+BAA+B;AACrE,QAAI,CAAC,SAAS,MAAM,EAAG,UAAS,OAAO;AAEvC,UAAM,YAAY,KAAK,IAAI,KACzB,OAAO,KAAK,EACZ,KAAK,SAAS,8BAA8B,EAC5C,MAAM,WAAW,MAAM,EACvB,MAAM,OAAO,MAAM,EACnB,MAAM,aAAa,MAAM;AAE3B,UAAM,cAAc,KAAK,SAAS,eAAe;AACjD,QAAI,CAAC,YAAa;AAGlB,UAAM,aAAa,YAAY,KAAK,eAAe,KAAK,CAAC;AACzD,UAAM,SAAmB,CAAC;AAC1B,eAAW,SAAS,OAAO,OAAO,UAAU,GAAY;AACvD,UAAI,MAAM,UAAU,kBAAmB,QAAO,KAAK,MAAM,WAAW;AAAA,IACrE;AAEA,UAAM,WAAW,OAAO,SAAS,IAAI,KAAK,IAAI,GAAG,OAAO,IAAI,OAAK,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI;AACjF,UAAM,aAAa,OAAoB,EACrC,OAAO,CAAC,CAAC,UAAU,GAAG,QAAQ,CAAC,EAC/B,MAAM,CAAC,WAAW,WAAW,SAAS,CAAC,EACvC,MAAM,IAAI;AAEZ,mBAAe;AAAA,MACd,QAAQ;AAAA,MACR,OAAO;AAAA,MACP,aAAa,KAAK;AAAA,MAClB,QAAQ;AAAA,MACR,SAAS,KAAK;AAAA,MACd,OAAO;AAAA,MACP,KAAK,KAAK,IAAI;AAAA,MACd,cAAc,CAAC,SAAiB;AAC/B,cAAM,QAAQ,WAAW,IAAI;AAC7B,YAAI,SAAS,MAAM,UAAU,kBAAmB,QAAO,WAAW,MAAM,WAAW;AACnF,eAAO;AAAA,MACR;AAAA,MACA,iBAAiB,CAAC,MAAc,UAAkB;AACjD,cAAM,QAAQ,WAAW,IAAI;AAC7B,YAAI,CAAC,MAAO,QAAO,GAAG,KAAK,KAAK,IAAI;AAAA;AACpC,cAAM,KAAK,MAAM,YAAY,QAAQ,CAAC;AACtC,cAAM,MAAM,CAAC,MAAe,KAAK,OAAS,EAAE,QAAQ,CAAC,IAAI,EAAE,cAAc,CAAC;AAC1E,cAAM,MAAM,OAAO,SAAS,MAAM,GAAG,IAAI;AAAA,OAAU,IAAI,MAAM,GAAG,CAAC,KAAK;AACtE,eAAO,GAAG,KAAK,KAAK,IAAI;AAAA,yBAAwB,EAAE;AAAA,WAAc,IAAI,MAAM,OAAO,CAAC,GAAG,GAAG;AAAA,MACzF;AAAA,IACD,CAAC;AAED,SAAK,aAAa,WAAW,YAAY,UAAU,OAAO,QAAQ,eAAe;AAAA,EAClF;AAAA,EAEA,aAAa,WAAgB,YAAiB,UAAkB,MAAc,SAAiB;AAC9F,UAAM,YAAY,UAChB,OAAO,KAAK,EACZ,MAAM,WAAW,MAAM,EACvB,MAAM,kBAAkB,QAAQ,EAChC,MAAM,mBAAmB,QAAQ,EACjC,MAAM,WAAW,MAAM;AAEzB,QAAI,CAAC,MAAM;AAGV,gBACE,OAAO,KAAK,EACZ,MAAM,aAAa,MAAM,EACzB,MAAM,SAAS,MAAM,EACrB,MAAM,aAAa,OAAO,EAC1B,MAAM,eAAe,KAAK,EAC1B;AAAA,QACA,uEAAuE,YAAY,mGAAmG,iBAAiB,wBAAwB,OAAO;AAAA,MACvO;AACD;AAAA,IACD;AAEA,cACE,OAAO,KAAK,EACZ,MAAM,eAAe,MAAM,EAC3B,MAAM,aAAa,MAAM,EACzB,MAAM,iBAAiB,KAAK,EAC5B,KAAK,yBAAoB;AAE3B,UAAM,cAAc;AACpB,UAAM,eAAe;AACrB,UAAM,MAAM,UACV,OAAO,KAAK,EACZ,KAAK,SAAS,cAAc,EAAE,EAC9B,KAAK,UAAU,eAAe,EAAE;AAElC,UAAM,OAAO,IAAI,OAAO,MAAM;AAC9B,UAAM,aAAa,qBAAqB,aAAa;AACrD,UAAM,WAAW,KACf,OAAO,gBAAgB,EACvB,KAAK,MAAM,UAAU,EACrB,KAAK,MAAM,GAAG,EACd,KAAK,MAAM,GAAG,EACd,KAAK,MAAM,GAAG,EACd,KAAK,MAAM,GAAG;AAEhB,UAAM,QAAQ;AACd,aAAS,IAAI,GAAG,KAAK,OAAO,KAAK;AAChC,YAAM,IAAI,IAAI;AACd,YAAM,MAAM,YAAY,IAAI,IAAI;AAChC,eACE,OAAO,MAAM,EACb,KAAK,UAAU,GAAG,IAAI,GAAG,GAAG,EAC5B,KAAK,cAAc,WAAW,GAAG,CAAC;AAAA,IACrC;AAEA,QACE,OAAO,MAAM,EACb,KAAK,KAAK,CAAC,EACX,KAAK,KAAK,EAAE,EACZ,KAAK,SAAS,WAAW,EACzB,KAAK,UAAU,YAAY,EAC3B,MAAM,QAAQ,QAAQ,UAAU,GAAG,EACnC,KAAK,UAAU,MAAM;AAEvB,UAAM,cAAc,OAAY,EAC9B,OAAO,CAAC,UAAU,CAAC,QAAQ,CAAC,EAC5B,MAAM,CAAC,IAAI,eAAe,EAAE,CAAC;AAE/B,UAAM,QAAQ,CAAC,CAAC,UAAU,CAAC,WAAW,GAAG,GAAG,WAAW,GAAG,QAAQ;AAClE,eAAW,QAAQ,OAAO;AACzB,YAAM,IAAI,YAAY,IAAI;AAC1B,UACE,OAAO,MAAM,EACb,KAAK,MAAM,WAAW,EACtB,KAAK,MAAM,CAAC,EACZ,KAAK,MAAM,cAAc,CAAC,EAC1B,KAAK,MAAM,CAAC,EACZ,KAAK,UAAU,MAAM;AACvB,UACE,OAAO,MAAM,EACb,KAAK,KAAK,cAAc,CAAC,EACzB,KAAK,KAAK,CAAC,EACX,KAAK,qBAAqB,SAAS,EACnC,KAAK,aAAa,MAAM,EACxB,KAAK,KAAK,QAAQ,CAAC,CAAC;AAAA,IACvB;AAEA,cACE,OAAO,KAAK,EACZ,MAAM,cAAc,MAAM,EAC1B,MAAM,aAAa,MAAM,EACzB,MAAM,SAAS,MAAM,EACrB;AAAA,MACA,uEAAuE,YAAY;AAAA,IACpF;AAAA,EACF;AACD;AAEO,IAAM,gBAAgB,YAAY,YAAY;AAErD,eAAsB,cAAc,MAAW;AAC9C,QAAM,SAAS,gBAAgB,aAAa;AAC5C,MAAI,CAAC,KAAK,KAAM,OAAM,IAAI,MAAM,iCAAiC;AACjE,SAAO,UAAU,QAAQ,IAAI;AAC9B;AAEO,SAAS,iBAAiB,QAAa,gBAAqB;AAClE,QAAM,MAAM,OAAO,OAAO,KAAK,EAAE,MAAM,WAAW,KAAK;AACvD,MAAI,OAAO,MAAM,EAAE,MAAM,eAAe,MAAM,EAAE,KAAK,oBAAoB;AAEzE,QAAM,aAAa,iBAAiB;AAAA,IACnC;AAAA,IACA,QAAQ,eAAe,IAAI,KAAK;AAAA,IAChC,KAAK,IAAI,KAAK,EAAE,SAAS,MAAM,CAAC;AAAA,IAChC,YAAY;AAAA,IACZ,UAAU,YAAY;AACrB,UAAI,CAAC,WAAW,WAAY,OAAM,IAAI,MAAM,oCAAoC;AAChF,qBAAe,IAAI,IAAI,KAAK;AAC5B,qBAAe,IAAI,SAAS;AAAA,QAC3B,MAAM;AAAA,QACN,QAAQ;AAAA,UACP,WAAW;AAAA,UACX,MAAM,WAAW;AAAA,QAClB;AAAA,MACD,CAAC;AAAA,IACF;AAAA,EACD,CAAC;AACF;",
+  "names": []
+}

package/dist/bubbleHeatmap-7DQNWBQ2.js

@@ -0,0 +1,378 @@
+import {
+  LegendCircleReference,
+  PlotBase,
+  addGeneSearchbox
+} from "./chunk-J7JDCNLU.js";
+import "./chunk-HJ6L54YS.js";
+import "./chunk-KV4W2ACA.js";
+import "./chunk-N635HDJ4.js";
+import "./chunk-T46FA72N.js";
+import {
+  Menu
+} from "./chunk-ELJX3QIQ.js";
+import "./chunk-EEB5VE2A.js";
+import "./chunk-6RRZRISL.js";
+import "./chunk-2KM4PRQM.js";
+import {
+  dofetch3
+} from "./chunk-JYOIO5UY.js";
+import "./chunk-YN5NY3D3.js";
+import "./chunk-RN4BOWRH.js";
+import "./chunk-IK2BO37K.js";
+import {
+  copyMerge,
+  getCompInit
+} from "./chunk-WINIL2KN.js";
+import "./chunk-PF4DSFDR.js";
+import "./chunk-D6G64XPJ.js";
+import "./chunk-W5J3LTYS.js";
+import "./chunk-YLJOZP4P.js";
+import {
+  linear,
+  sqrt
+} from "./chunk-HDTFYTEL.js";
+import "./chunk-FXQXCOII.js";
+import "./chunk-TLT4YIG3.js";
+import "./chunk-5R63Q5KH.js";
+import "./chunk-I6Y4O3RR.js";
+import "./chunk-Q5RDQNIT.js";
+import "./chunk-DQC5FFGV.js";
+import "./chunk-HS5PO5ZQ.js";
+
+// plots/bubbleHeatmap.ts
+var defaultConfig = { chartType: "bubbleHeatmap" };
+var CELL_W = 92;
+var CELL_H = 64;
+var ROW_LABEL_W = 170;
+var COL_LABEL_H = 92;
+var SITE_DOT_R = 5;
+var SITE_DOT_SP = 13;
+var CELL_PAD = 8;
+var MIN_DOT_R = 8;
+var MAX_DOT_R = 20;
+var NEG_LOG_FDR_CAP = 10;
+var BubbleHeatmap = class _BubbleHeatmap extends PlotBase {
+  constructor(opts, api) {
+    super(opts, api);
+    this.currentIsoform = "";
+    this.useAdjusted = false;
+    this.type = _BubbleHeatmap.type;
+    this.components = {};
+  }
+  static {
+    this.type = "bubbleHeatmap";
+  }
+  async init() {
+    const holder = this.opts.holder.append("div").style("padding", "10px");
+    this.dom = {
+      holder,
+      body: holder.append("div"),
+      tip: new Menu({ padding: "" }),
+      header: this.opts.header
+    };
+    if (this.dom.header) this.dom.header.html("Bubble Heatmap");
+  }
+  getState(appState) {
+    const config = appState.plots.find((p) => p.id === this.id);
+    if (!config) throw `No plot with id='${this.id}' found`;
+    return { config };
+  }
+  async main() {
+    const gene = this.state.config?.gene;
+    if (!gene) throw new Error("bubbleHeatmap: gene is missing");
+    if (this.dom.header) this.dom.header.text(`Bubble Heatmap: ${gene}`);
+    const body = {
+      genome: this.app.opts.state.vocab.genome,
+      dslabel: this.app.opts.state.vocab.dslabel,
+      gene
+    };
+    const data = await dofetch3("termdb/bubbleHeatmap", { body });
+    if (data.error) throw data.error;
+    this.data = data;
+    this.dom.body.selectAll("*").remove();
+    const isoformIds = Object.keys(data.isoforms);
+    if (isoformIds.length === 0) {
+      this.dom.body.append("div").style("padding", "20px").style("color", "#666").text(`No data found for gene "${gene}" in any (assay, cohort) DAPfile.`);
+      return;
+    }
+    this.useAdjusted = !!data.proteinReferenceAssay;
+    this.currentIsoform = isoformIds[0];
+    const isoBlock = this.dom.body.append("div").style("margin-bottom", "12px");
+    isoBlock.append("span").style("font-weight", "bold").text("Isoform: ");
+    if (isoformIds.length > 1) {
+      const sel = isoBlock.append("select").style("margin-left", "5px").style("padding", "3px 6px").on("change", () => {
+        this.currentIsoform = sel.node().value;
+        this.renderGrid();
+      });
+      sel.selectAll("option").data(isoformIds).enter().append("option").attr("value", (d) => d).text((d) => `${data.isoforms[d].gene_name} \u2014 ${d}`);
+    } else {
+      isoBlock.append("span").style("margin-left", "5px").text(`${data.isoforms[this.currentIsoform].gene_name} \u2014 ${this.currentIsoform}`);
+    }
+    this.gridHolder = this.dom.body.append("div");
+    this.renderGrid();
+  }
+  renderGrid() {
+    const data = this.data;
+    const selectedIsoform = this.currentIsoform;
+    const useAdjusted = this.useAdjusted;
+    const refAssay = data.proteinReferenceAssay;
+    const threshold = data.fdrThreshold;
+    this.gridHolder.selectAll("*").remove();
+    const container = this.gridHolder.append("div").style("display", "flex").style("gap", "24px").style("align-items", "flex-start").style("flex-wrap", "wrap");
+    const isoformData = data.isoforms[selectedIsoform];
+    if (!isoformData) return;
+    const assays = data.assays;
+    const cohorts = data.cohorts;
+    const nRows = assays.length;
+    const nCols = cohorts.length;
+    const ptmAssays = new Set(data.ptmAssays || []);
+    const isPTMassay = (assay) => ptmAssays.has(assay);
+    const valueOf = (s) => this.valueFor(s, useAdjusted);
+    const negLogFdr = (fdr) => fdr > 0 ? Math.min(-Math.log10(fdr), NEG_LOG_FDR_CAP) : NEG_LOG_FDR_CAP;
+    const slotIndex = /* @__PURE__ */ new Map();
+    const assaySlotCount = /* @__PURE__ */ new Map();
+    let maxAbs = 0;
+    const thresholdNegLog = negLogFdr(threshold);
+    let maxNegLog = thresholdNegLog;
+    for (const assay of assays) {
+      const ptm = isPTMassay(assay);
+      const rawSum = /* @__PURE__ */ new Map();
+      const rawN = /* @__PURE__ */ new Map();
+      const significantSomewhere = /* @__PURE__ */ new Set();
+      for (const cohort of cohorts) {
+        const cell = isoformData.data[assay]?.[cohort];
+        if (!cell) continue;
+        if (ptm) {
+          for (const s of cell.sites) {
+            if (s.significant) {
+              const v = Math.abs(valueOf(s));
+              if (v > maxAbs) maxAbs = v;
+            }
+            rawSum.set(s.id, (rawSum.get(s.id) ?? 0) + s.log2FC);
+            rawN.set(s.id, (rawN.get(s.id) ?? 0) + 1);
+            if (s.significant) significantSomewhere.add(s.id);
+          }
+        } else {
+          const s = cell.sites[0];
+          if (!s) continue;
+          const v = Math.abs(valueOf(s));
+          if (v > maxAbs) maxAbs = v;
+          const nl = negLogFdr(s.fdr);
+          if (nl > maxNegLog) maxNegLog = nl;
+        }
+      }
+      if (ptm) {
+        const meanRaw = (id) => rawSum.get(id) / rawN.get(id);
+        const ordered = [...significantSomewhere].sort((a, b) => meanRaw(b) - meanRaw(a));
+        ordered.forEach((id, i) => slotIndex.set(`${assay}|${id}`, i));
+        assaySlotCount.set(assay, ordered.length);
+      } else {
+        assaySlotCount.set(assay, 1);
+      }
+    }
+    if (maxAbs === 0) maxAbs = 1;
+    if (maxNegLog <= thresholdNegLog) maxNegLog = thresholdNegLog + 1;
+    const colorScale = linear().domain([-maxAbs, 0, maxAbs]).range(["#2166ac", "#f7f7f7", "#b2182b"]).clamp(true);
+    const sizeScale = sqrt().domain([thresholdNegLog, maxNegLog]).range([MIN_DOT_R, MAX_DOT_R]).clamp(true);
+    const layout = assays.map((assay) => {
+      const m = assaySlotCount.get(assay);
+      const subCols = Math.max(1, Math.min(m, Math.floor((CELL_W - 2 * CELL_PAD) / SITE_DOT_SP)));
+      const rows = Math.ceil(m / subCols);
+      return { subCols, rows, height: Math.max(CELL_H, rows * SITE_DOT_SP + 2 * CELL_PAD) };
+    });
+    const rowY = [];
+    let yAcc = COL_LABEL_H;
+    for (let r = 0; r < nRows; r++) {
+      rowY[r] = yAcc;
+      yAcc += layout[r].height;
+    }
+    const gridW = ROW_LABEL_W + nCols * CELL_W + 20;
+    const gridH = yAcc + 20;
+    const svg = container.append("svg").attr("width", gridW).attr("height", gridH).style("flex", "0 0 auto");
+    const grid = svg.append("g");
+    for (let c = 0; c < nCols; c++) {
+      const cx = ROW_LABEL_W + c * CELL_W + CELL_W / 2;
+      grid.append("text").attr("x", cx).attr("y", COL_LABEL_H - 10).attr("text-anchor", "start").attr("font-size", "12px").attr("font-weight", "bold").attr("transform", `rotate(-35 ${cx} ${COL_LABEL_H - 10})`).text(cohorts[c]);
+    }
+    for (let r = 0; r < nRows; r++) {
+      const cy = rowY[r] + layout[r].height / 2;
+      const m = assaySlotCount.get(assays[r]);
+      const lbl = grid.append("text").attr("x", ROW_LABEL_W - 10).attr("y", cy).attr("text-anchor", "end").attr("dominant-baseline", "central").attr("font-size", "12px").attr("font-weight", "bold");
+      lbl.append("tspan").text(assays[r]);
+      lbl.append("tspan").attr("x", ROW_LABEL_W - 10).attr("dy", "1.3em").attr("font-weight", "normal").attr("font-size", "10px").attr("fill", "#888").text(m > 1 ? `${m} sites` : "");
+    }
+    for (let r = 0; r < nRows; r++) {
+      const assay = assays[r];
+      const ptm = isPTMassay(assay);
+      const { subCols, height } = layout[r];
+      for (let c = 0; c < nCols; c++) {
+        const x0 = ROW_LABEL_W + c * CELL_W;
+        const y0 = rowY[r];
+        grid.append("rect").attr("x", x0).attr("y", y0).attr("width", CELL_W).attr("height", height).attr("fill", "none").attr("stroke", "#eee").attr("stroke-width", 1);
+        const cell = isoformData.data[assay]?.[cohorts[c]];
+        if (!cell || !cell.sites.length) continue;
+        const addDot = (s, cx, cy, radius) => {
+          return grid.append("circle").attr("cx", cx).attr("cy", cy).attr("r", radius).attr("fill", colorScale(valueOf(s))).attr("stroke", "#888").attr("stroke-width", 0.8).style("opacity", s.significant ? 1 : 0.35).on(
+            "mouseover",
+            (event) => this.showSiteTip(
+              event,
+              isoformData.gene_name,
+              selectedIsoform,
+              assay,
+              cohorts[c],
+              s,
+              useAdjusted,
+              refAssay
+            )
+          ).on("mouseout", () => this.dom.tip.hide());
+        };
+        if (!ptm) {
+          const s = cell.sites[0];
+          const cx = x0 + CELL_W / 2;
+          const cy = y0 + height / 2;
+          addDot(s, cx, cy, sizeScale(negLogFdr(s.fdr)));
+          continue;
+        }
+        const blockW = subCols * SITE_DOT_SP;
+        const blockH = layout[r].rows * SITE_DOT_SP;
+        const startX = x0 + (CELL_W - blockW) / 2 + SITE_DOT_SP / 2;
+        const startY = y0 + (height - blockH) / 2 + SITE_DOT_SP / 2;
+        for (const s of cell.sites) {
+          if (!s.significant) continue;
+          const slot = slotIndex.get(`${assay}|${s.id}`);
+          const cx = startX + slot % subCols * SITE_DOT_SP;
+          const cy = startY + Math.floor(slot / subCols) * SITE_DOT_SP;
+          addDot(s, cx, cy, SITE_DOT_R);
+        }
+      }
+    }
+    this.renderLegend(container, colorScale, maxAbs, threshold, useAdjusted, refAssay, maxNegLog);
+  }
+  fmtFdr(v) {
+    return v >= 1e-4 ? v.toFixed(4) : v.toExponential(2);
+  }
+  /** true when the protein-adjusted value should be shown instead of raw log2FC */
+  showsAdjusted(s, useAdjusted) {
+    return !!(useAdjusted && s.adjustedAvailable && s.adjustedLog2FC != null);
+  }
+  /** value encoded by color: protein-adjusted when requested & available, else raw */
+  valueFor(s, useAdjusted) {
+    return this.showsAdjusted(s, useAdjusted) ? s.adjustedLog2FC : s.log2FC;
+  }
+  showSiteTip(event, geneName, isoform, assay, cohort, s, useAdjusted, refAssay) {
+    this.dom.tip.clear().show(event.clientX, event.clientY);
+    const t = this.dom.tip.d.append("div").style("padding", "8px").style("font-size", "13px");
+    t.append("div").style("font-weight", "bold").style("margin-bottom", "4px").text(`${geneName} \u2014 ${isoform}`);
+    t.append("div").text(`Assay: ${assay}`);
+    t.append("div").text(`Sample set: ${cohort}`);
+    const isPTM = (this.data.ptmAssays || []).includes(assay);
+    t.append("div").text(`${isPTM ? "Site" : "Protein"}: ${s.id}`);
+    t.append("div").text(`raw log\u2082FC: ${s.log2FC.toFixed(3)}`);
+    if (s.adjustedAvailable) {
+      t.append("div").text(`protein log\u2082FC: ${s.proteinLog2FC.toFixed(3)}`);
+      t.append("div").text(`adjusted log\u2082FC: ${s.adjustedLog2FC.toFixed(3)}`);
+    } else if (refAssay && isPTM) {
+      t.append("div").style("color", "#999").text("adjusted: n/a (protein not measured)");
+    }
+    t.append("div").text(`FDR: ${this.fmtFdr(s.fdr)}`);
+    const shown = this.showsAdjusted(s, useAdjusted) ? "adjusted" : "raw";
+    t.append("div").style("color", "#666").style("margin-top", "4px").text(`Color = ${shown} log\u2082FC.`);
+  }
+  renderLegend(container, colorScale, maxAbs, threshold, useAdjusted, refAssay, maxNegLog) {
+    const legend = container.append("div").style("display", "flex").style("flex-direction", "column").style("gap", "16px").style("padding", "8px 0").style("min-width", "180px").style("max-width", "260px");
+    const colorBlock = legend.append("div");
+    colorBlock.append("div").style("font-weight", "bold").style("font-size", "13px").style("margin-bottom", "6px").text(useAdjusted && refAssay ? "log\u2082FC (PTM-adjusted)" : "log\u2082FC");
+    const cW = 22;
+    const cH = 130;
+    const cSvg = colorBlock.append("svg").attr("width", cW + 60).attr("height", cH + 16);
+    const gid = `bh-grad-${this.id}`;
+    const grad = cSvg.append("defs").append("linearGradient").attr("id", gid).attr("x1", "0").attr("y1", "0").attr("x2", "0").attr("y2", "1");
+    const steps = 10;
+    for (let i = 0; i <= steps; i++) {
+      const t = i / steps;
+      grad.append("stop").attr("offset", `${t * 100}%`).attr("stop-color", colorScale(maxAbs * (1 - 2 * t)));
+    }
+    cSvg.append("rect").attr("x", 0).attr("y", 8).attr("width", cW).attr("height", cH).style("fill", `url(#${gid})`).attr("stroke", "#999");
+    const cScale = linear().domain([maxAbs, -maxAbs]).range([8, cH + 8]);
+    for (const tick of [maxAbs, maxAbs / 2, 0, -maxAbs / 2, -maxAbs]) {
+      const y = cScale(tick);
+      cSvg.append("line").attr("x1", cW).attr("y1", y).attr("x2", cW + 5).attr("y2", y).attr("stroke", "#666");
+      cSvg.append("text").attr("x", cW + 8).attr("y", y).attr("dominant-baseline", "central").attr("font-size", "10px").text(tick.toFixed(2));
+    }
+    const sizeBlock = legend.append("div");
+    sizeBlock.append("div").style("font-weight", "bold").style("font-size", "13px").style("margin-bottom", "6px").text("Non-PTM dot size: significance (\u2212log\u2081\u2080 FDR)");
+    const sSvg = sizeBlock.append("svg");
+    const sG = sSvg.append("g");
+    new LegendCircleReference({
+      g: sG,
+      inputMin: 0,
+      inputMax: MAX_DOT_R * 2,
+      minRadius: MIN_DOT_R,
+      maxRadius: MAX_DOT_R,
+      // capped to match the size scale's domain min (thresholdNegLog in renderGrid)
+      minLabel: Number(Math.min(-Math.log10(threshold), NEG_LOG_FDR_CAP).toFixed(1)),
+      maxLabel: Number(maxNegLog.toFixed(1))
+    });
+    const sPad = 4;
+    const sBox = sG.node().getBBox();
+    sG.attr("transform", `translate(${sPad - sBox.x}, ${sPad - sBox.y})`);
+    sSvg.attr("width", Math.ceil(sBox.width + 2 * sPad)).attr("height", Math.ceil(sBox.height + 2 * sPad));
+    if (refAssay) {
+      const adjLabel = legend.append("div").append("label").style("display", "flex").style("align-items", "center").style("gap", "6px").style("cursor", "pointer").style("font-size", "13px").style("font-weight", "bold").attr(
+        "title",
+        `When checked, the PTM assays have the ${refAssay} log\u2082FC subtracted; other assays are shown unchanged.`
+      );
+      const adjCb = adjLabel.append("input").attr("type", "checkbox").property("checked", this.useAdjusted).on("change", () => {
+        this.useAdjusted = adjCb.property("checked");
+        this.renderGrid();
+      });
+      adjLabel.append("span").style("font-weight", "normal").text("Adjust PTM for total protein abundance");
+    }
+    const notes = legend.append("div").style("font-size", "11px").style("color", "#666").style("line-height", "1.5").style("max-width", "240px").style("overflow-wrap", "break-word");
+    notes.append("div").text(
+      `Color = log\u2082FC. Dot size = significance, \u2212log\u2081\u2080 FDR (non-PTM rows); the smallest size marks the FDR < ${threshold} cutoff. Non-significant dots are faded.`
+    );
+    notes.append("div").style("margin-top", "4px").text(
+      "PTM rows: one fixed-size dot per site significant in that cohort, positions stable across cohorts; non-significant sites are not shown."
+    );
+    notes.append("div").style("margin-top", "4px").text(
+      "A slot stays empty where the site is not significant in that cohort, the assay was not performed, or the protein was not detected."
+    );
+    if (refAssay) {
+      notes.append("div").style("margin-top", "4px").text(`Adjusted log\u2082FC = a PTM site's log\u2082FC \u2212 ${refAssay} log\u2082FC (PTM assays only).`);
+    }
+  }
+};
+var componentInit = getCompInit(BubbleHeatmap);
+async function getPlotConfig(opts) {
+  const config = structuredClone(defaultConfig);
+  if (!opts.gene) throw new Error("bubbleHeatmap requires opts.gene");
+  return copyMerge(config, opts);
+}
+function makeChartBtnMenu(holder, chartsInstance) {
+  const row = holder.append("div").style("padding", "5px");
+  row.append("span").style("font-weight", "bold").text("Enter a gene name:");
+  const geneSearch = addGeneSearchbox({
+    row,
+    genome: chartsInstance.app.opts.genome,
+    tip: new Menu({ padding: "0px" }),
+    searchOnly: "gene",
+    callback: async () => {
+      if (!geneSearch.geneSymbol) throw new Error("A valid gene selection is required");
+      chartsInstance.dom.tip.hide();
+      chartsInstance.app.dispatch({
+        type: "plot_create",
+        config: {
+          chartType: "bubbleHeatmap",
+          gene: geneSearch.geneSymbol
+        }
+      });
+    }
+  });
+}
+export {
+  componentInit,
+  getPlotConfig,
+  makeChartBtnMenu
+};
+//# sourceMappingURL=bubbleHeatmap-7DQNWBQ2.js.map