@sjcrh/proteinpaint-client 2.187.0 → 2.188.0

package/dist/chunk-FUSVXQ6C.js

@@ -0,0 +1,368 @@
+import {
+  fillTermWrapper,
+  termsettingInit
+} from "./chunk-2XPDA6NL.js";
+import {
+  isNumericTerm
+} from "./chunk-FDK4UJMN.js";
+import {
+  select_default
+} from "./chunk-I6Y4O3RR.js";
+
+// plots/matrix/matrix.renderers.js
+function setRenderers(self) {
+  self.render = function() {
+    const s = self.settings.matrix;
+    const l = self.layout;
+    const d = self.dimensions;
+    const duration = self.dom.svg.attr("width") ? s.duration : 0;
+    self.dom.clipRect.attr("x", d.xOffset - 1).attr("y", 0).attr("width", d.mainw + 3).attr("height", d.mainh + 500);
+    self.renderSerieses(s, l, d, duration);
+    self.renderLabels(s, l, d, duration);
+    self.renderDivideByLabel(s, l, d, duration);
+    self.dom.colBeam.attr("width", d.dx).attr("height", d.mainh).style("stroke", s.beamStroke);
+    self.dom.rowBeam.attr("width", d.zoomedMainW).attr("height", s.rowh).style("stroke", s.beamStroke);
+  };
+  self.renderSerieses = function(s, l, d, duration) {
+    if (self.prevUseCanvas != s.useCanvas) {
+      self.dom.seriesesG.selectAll("g").remove();
+    }
+    if (s.useCanvas) {
+      const _g = self.dom.seriesesG.selectAll("g");
+      const g = (
+        /*(_g.size() && _g) ||*/
+        self.dom.seriesesG.append("g").datum(this.serieses)
+      );
+      self.renderCanvas(this.serieses, g, d, s, _g, duration);
+    } else {
+      self.dom.seriesesG.attr("transform", `translate(${d.xOffset + d.seriesXoffset},${d.yOffset})`);
+      const sg = self.dom.seriesesG.selectAll(".sjpp-mass-series-g").data(this.serieses, (series) => series.tw.$id);
+      sg.exit().remove();
+      sg.each(self.renderSeries);
+      sg.enter().append("g").attr("class", "sjpp-mass-series-g").style("opacity", 1e-3).each(self.renderSeries);
+      self.mouseout();
+    }
+    self.prevUseCanvas = s.useCanvas;
+  };
+  self.renderSeries = async function(series) {
+    const s = self.settings.matrix;
+    const d = self.dimensions;
+    const g = select_default(this);
+    const duration = g.attr("transform") ? s.duration : 0;
+    g.attr("transform", `translate(${series.x},${series.y})`).style("opacity", 1);
+    const last = series.cells[series.cells.length - 1];
+    const height = series.y + last?.y + s.rowh;
+    const rects = g.selectAll("rect").data(series.cells, (cell) => cell.sample + ";;" + cell.tw.$id + ";;" + cell.valueIndex);
+    rects.exit().remove();
+    rects.each(self.renderCell);
+    rects.enter().append("rect").each(self.renderCell);
+  };
+  self.renderCanvas = async function(serieses, g, d, s, _g, duration) {
+    const pxr = window.devicePixelRatio <= 1 ? 1 : window.devicePixelRatio;
+    g.selectAll("*").remove();
+    const width = d.imgW;
+    const height = self.dimensions.mainh;
+    const canvas = window.OffscreenCanvas ? new OffscreenCanvas(width * pxr, height * pxr) : (
+      // TODO: no need to support older browser versions???
+      self.dom.holder.append("canvas").attr("width", pxr * width).attr("height", pxr * height).style("opacity", 0).node()
+    );
+    const ctx = canvas.getContext("2d");
+    ctx.imageSmoothingEnabled = false;
+    ctx.imageSmoothingQuality = "high";
+    ctx.scale(pxr, pxr);
+    for (const series of serieses) {
+      for (const cell of series.cells) {
+        self.renderCellWithCanvas(ctx, cell, series, s, d, series.y);
+      }
+    }
+    if (window.OffscreenCanvas) {
+      const reader = new FileReader();
+      reader.addEventListener(
+        "load",
+        () => {
+          _g?.remove();
+          self.dom.seriesesG.attr("transform", `translate(${d.xOffset + d.seriesXoffset},${d.yOffset})`);
+          g.selectAll("image").remove();
+          g.append("image").attr("xlink:href", reader.result).attr("x", d.xMin).attr("width", width).attr("height", height);
+        },
+        false
+      );
+      const blob = await canvas.convertToBlob({ quality: 1 });
+      const dataURL = reader.readAsDataURL(blob);
+    } else {
+      _g?.remove();
+      self.dom.seriesesG.attr("transform", `translate(${d.xOffset + d.seriesXoffset},${d.yOffset})`);
+      const dataURL = canvas.toDataURL();
+      const ratio = window.devicePixelRatio * window.devicePixelRatio;
+      g.append("image").attr("width", width).attr("height", height).attr("xlink:href", dataURL);
+      if (!window.OffscreenCanvas) canvas.remove();
+    }
+    self.mouseout();
+  };
+  self.renderCellWithCanvas = function(ctx, cell, series, s, d, _y) {
+    if (!cell.fill)
+      cell.fill = cell.$id in self.colorScaleByTermId ? self.colorScaleByTermId[cell.$id](cell.key) : getRectFill(cell);
+    const x = cell.x ? cell.x - d.xMin : 0;
+    const y = _y ? _y + cell.y : cell.y || 0;
+    const width = s.useMinPixelWidth ? Math.max(cell.width || d.colw, d.pxw) : cell.width || d.colw;
+    const height = "height" in cell ? cell.height : s.rowh;
+    ctx.fillStyle = cell.fill;
+    ctx.fillRect(x, y, width, height);
+    const borderWidth = Math.min(width, height) * 0.1;
+    if (cell.border) {
+      ctx.lineWidth = borderWidth;
+      ctx.strokeStyle = "white";
+      ctx.strokeRect(x, y, width, height);
+    }
+  };
+  self.renderCell = function(cell) {
+    if (!cell.fill)
+      cell.fill = cell.$id in self.colorScaleByTermId ? self.colorScaleByTermId[cell.$id](cell.key) : getRectFill(cell);
+    const s = self.settings.matrix;
+    const rect = select_default(this).attr("x", cell.x || 0).attr("y", cell.y || 0).attr("width", cell.width || self.dimensions.colw).attr("height", "height" in cell ? Math.max(0, cell.height) : s.rowh).attr("shape-rendering", "crispEdges").attr("fill", cell.fill);
+    if (cell.border) {
+      rect.attr("stroke", "white").attr("stroke-width", 0.8);
+    }
+  };
+  self.renderLabels = function(s, l, d, duration) {
+    for (const direction of ["top", "btm", "left", "right"]) {
+      let renderLabel2 = function(lab) {
+        const g = select_default(this);
+        g.attr("transform", side.attr.labelGTransform);
+        if (!g.select(":scope>text").size()) g.append("text");
+        const showContAxis = !side.isGroup && lab.grp?.type !== "hierCluster" && lab.tw?.q?.mode == "continuous";
+        const labelText = side.label(lab);
+        const text = g.select(":scope>text").attr("fill", "#000");
+        let continuousBarHAdjust;
+        const twSpecificSettings = self.config.settings.matrix.twSpecificSettings;
+        const twSettingsBarH = twSpecificSettings[lab.tw?.$id]?.contBarH;
+        if (twSettingsBarH && s.barh) continuousBarHAdjust = (twSettingsBarH - s.barh) * 0.5;
+        text.attr(
+          "display",
+          lab.grp?.type === "hierCluster" && s.clusterRowh < 6 ? "none" : side.attr.fontSize < 6 || labelText === "configure" ? "none" : ""
+        ).attr("font-size", lab.grp?.type === "hierCluster" ? Math.max(4, s.clusterRowh - 4) : side.attr.fontSize).attr("text-anchor", side.attr.labelAnchor).attr(
+          "transform",
+          side.attr.labelTransform + (continuousBarHAdjust ? ` translate(0,${continuousBarHAdjust})` : "")
+        ).attr("cursor", "pointer").attr(side.attr.textpos.coord, side.attr.textpos.factor * (showContAxis ? 30 : 0));
+        if (!Array.isArray(labelText)) {
+          text.text(labelText);
+          text.attr(
+            "y",
+            lab.grp?.type !== "hierCluster" && lab.tw?.q?.mode == "continuous" ? 10 : lab.grp?.type === "hierCluster" ? 0.1 * s.clusterRowh : 0
+          );
+          if (lab.grp?.type !== "hierCluster" && lab.tw?.q?.mode == "continuous") text.attr("x", -20);
+        } else {
+          text.text("");
+          const tspan = text.selectAll("tspan").data(labelText);
+          tspan.enter().append("tspan").attr("class", getTspanCls2).attr("dx", getTspanDx2).attr("font-size", getTspanFontSize2).text(getTspanText2);
+        }
+        text.on("mouseover", labelText === "configure" ? () => text.attr("opacity", 0.5) : null).on("mouseout", labelText === "configure" ? () => text.attr("opacity", 0) : null);
+        const hasAxis = g.select(".sjpp-matrix-cell-axis").size() && true;
+        if (showContAxis && labelText) {
+          if (!hasAxis) {
+            g.append("g").attr("class", "sjpp-matrix-cell-axis").attr("shape-rendering", "crispEdges");
+          }
+          const axisg = g.select(".sjpp-matrix-cell-axis");
+          axisg.selectAll("*").remove();
+          const domain = [lab.counts.maxval, lab.counts.minval];
+          if (s.transpose) domain.reverse();
+          const twSpecificSettings2 = self.config.settings.matrix.twSpecificSettings;
+          const twSettings = twSpecificSettings2[lab.tw.$id];
+          const x = !s.transpose ? 0 : twSettings.contBarGap - 1 - lab.labelOffset;
+          const y = !s.transpose ? twSettings.contBarGap - 1 - lab.labelOffset : 0;
+          axisg.attr("shape-rendering", "crispEdges").attr("transform", `translate(${x},${y})`).call(side.attr.axisFxn(lab.scales.full.domain(lab.scales.tickValues)).tickValues(lab.scales.tickValues));
+        } else if (hasAxis) {
+          g.select(".sjpp-matrix-cell-axis").remove();
+        }
+      }, getTspanCls2 = function(d2) {
+        return d2.cls;
+      }, getTspanDx2 = function(d2) {
+        return d2.dx;
+      }, getTspanFontSize2 = function(d2) {
+        return d2.fontSize || side.attr.fontSize;
+      }, getTspanText2 = function(d2) {
+        return d2.text;
+      };
+      var renderLabel = renderLabel2, getTspanCls = getTspanCls2, getTspanDx = getTspanDx2, getTspanFontSize = getTspanFontSize2, getTspanText = getTspanText2;
+      const side = l[direction];
+      side.box.style("display", side.display || "").attr("transform", side.attr.boxTransform);
+      const labels = side.box.selectAll(".sjpp-matrix-label").data(side.data, side.key);
+      labels.exit().remove();
+      labels.each(renderLabel2);
+      labels.enter().append("g").attr("class", "sjpp-matrix-label").each(renderLabel2);
+    }
+  };
+  self.colLabelGTransform = (lab, grpIndex) => {
+    const s = self.settings.matrix;
+    const d = self.dimensions;
+    lab.labelOffset = 0.8 * d.colw;
+    const x = lab.grpIndex * s.colgspace + lab.totalIndex * d.dx + lab.labelOffset + lab.totalHtAdjustments;
+    const y = 0;
+    return `translate(${x + d.seriesXoffset},${y})`;
+  };
+  self.colGrpLabelGTransform = (lab, grpIndex) => {
+    const s = self.settings.matrix;
+    const d = self.dimensions;
+    const len = (lab.processedLst || lab.grp.lst).length;
+    const x = lab.grpIndex * s.colgspace + lab.prevGrpTotalIndex * d.dx + len * d.dx / 2 + s.grpLabelFontSize / 2 + lab.totalHtAdjustments;
+    return `translate(${x + d.seriesXoffset},0)`;
+  };
+  self.rowLabelGTransform = (lab, grpIndex) => {
+    const s = self.settings.matrix;
+    const d = self.dimensions;
+    const x = 0;
+    lab.labelOffset = 0.7 * (lab.grp.type == "hierCluster" ? s.clusterRowh : s.rowh);
+    const y = lab.grpIndex * s.rowgspace + lab.totalIndex * d.dy + lab.labelOffset + lab.totalHtAdjustments;
+    return `translate(${x},${y})`;
+  };
+  self.rowGrpLabelGTransform = (lab, grpIndex) => {
+    const s = self.settings.matrix;
+    const d = self.dimensions;
+    const len = (lab.processedLst || lab.grp.lst).length;
+    const x = lab.tw?.q?.mode == "continuous" ? 20 : 0;
+    const y = lab.grpIndex * s.rowgspace + lab.prevGrpTotalIndex * d.dy + len * d.dy / 2 + s.grpLabelFontSize / 2 + lab.totalHtAdjustments;
+    return `translate(${x},${y})`;
+  };
+  self.rowAxisGTransform = (lab, grpIndex) => {
+    const s = self.settings.matrix;
+    const d = self.dimensions;
+    const x = 0;
+    const y = lab.grpIndex * s.rowgspace + lab.totalIndex * d.dy + 0.7 * s.rowh + lab.totalHtAdjustments;
+    return `translate(${x},${y})`;
+  };
+  self.renderDivideByLabel = async (s, l, d) => {
+    self.dom.mainG.selectAll(".sjpp-matrix-divide-by-label").remove();
+    if (!self.config.divideBy) return;
+    const name = self.config.divideBy?.term.name || "";
+    const text = name.length <= s.rowlabelmaxchars ? name : name.slice(0, s.rowlabelmaxchars) + "\u2026";
+    const sides = !s.transpose ? [l.left, l.right] : [l.top, l.bottom];
+    const box = sides.find((d2) => !d2.isGroup)?.box;
+    const y = (s.collabelpos == "top" ? d.mainh + s.collabelmaxchars : -s.collabelmaxchars) + 8;
+    const anchor = s.rowlabelpos == "left" ? "end" : "start";
+    const cl = s.controlLabels;
+    const gNote = box.append("g").attr("class", "sjpp-matrix-divide-by-label").attr("transform", `translate(0, ${y})`);
+    gNote.append("text").attr("text-anchor", anchor).attr("font-style", "italic").attr("y", -20).text(`${cl.Samples} grouped by`);
+    const g = box.datum({ tw: self.config.divideBy }).append("g").attr("class", "sjpp-matrix-divide-by-label").attr("transform", `translate(0, ${y})`).on("click", (event, d2) => {
+      pill.showMenu(event, textElem.node());
+    });
+    const textElem = g.append("text").attr("text-anchor", anchor).attr("font-weight", 600).text(text);
+    g.append("title").text(`${cl.Samples} are grouped by this gene or variable. Click to edit.`);
+    const customMenuOptions = [];
+    const tvsKey = isNumericTerm(self.config.divideBy.term) ? "ranges" : "values";
+    if (self.config.legendValueFilter.lst?.find(
+      (l2) => l2.legendGrpName == self.config.divideBy.term.id || l2.legendGrpName == self.config.divideBy.term.name
+    )?.tvs[tvsKey]?.length) {
+      customMenuOptions.push({ label: `Show filtered ${cl.samples}`, callback: self.showDeletedSampleGroups });
+    }
+    const pill = await termsettingInit({
+      menuOptions: "{edit,replace,remove}",
+      //numericEditMenuVersion: opts.numericEditMenuVersion,
+      customMenuOptions,
+      //custom menu options other than menuOptions
+      vocabApi: self.app.vocabApi,
+      vocab: self.state.vocab,
+      //activeCohort: opts.state?.activeCohort,
+      holder: g,
+      debug: self.opts.debug,
+      usecase: { target: "matrix" },
+      getBodyParams: () => {
+        const currentGeneNames = self.termOrder.filter((t) => t.tw.term.type === "geneVariant").map(
+          (t) => t.tw.term.chr ? `${t.tw.term.chr}:${t.tw.term.start}-${t.tw.term.stop}` : t.tw.term.gene || t.tw.term.name
+        );
+        if (currentGeneNames.length) return { currentGeneNames };
+        return {};
+      },
+      callback: async (tw) => {
+        if (self.dom.loadingDiv && self.dom.svg) {
+          self.dom.loadingDiv.selectAll("*").remove();
+          self.dom.loadingDiv.html("").style("display", "").style("position", "relative").style("left", "45%");
+          self.dom.loadingDiv.html("Processing data ...");
+          self.dom.svg.style("opacity", 0.1).style("pointer-events", "none");
+        }
+        if (tw && !tw.q) throw "data.q{} missing from pill callback";
+        if (tw?.term && isNumericTerm(tw.term)) {
+          tw.q = { ...tw.q, mode: "discrete" };
+        }
+        if (tw) await fillTermWrapper(tw, self.app.vocabApi);
+        await pill.main(tw ? tw : { term: null, q: null });
+        box.datum({ tw });
+        self.app.dispatch({
+          type: "plot_edit",
+          id: self.id,
+          config: {
+            divideBy: tw,
+            legendValueFilter: self.mayRemoveTvsEntry(self.config.divideBy)
+          }
+        });
+      }
+    });
+    const arg = {
+      term: self.config.divideBy.term,
+      q: self.config.divideBy.q
+    };
+    if (self.config.divideBy.$id) arg.$id = self.config.divideBy.$id;
+    pill.main(arg);
+  };
+  self.adjustSvgDimensions = async function(prevTranspose) {
+    const s = self.settings.matrix;
+    const hc = self.settings.hierCluster || {};
+    const l = self.layout;
+    const hcHeight = !hc.yDendrogramHeight ? 0 : hc.yDendrogramHeight + (l.top.display === "none" ? 0 : 10);
+    const hcWidth = hc.xDendrogramHeight || 0;
+    const d = self.dimensions;
+    const duration = self.dom.svg.attr("width") ? s.duration : 0;
+    await sleep(prevTranspose == s.transpose ? duration : s.duration);
+    const topBox = l.top.box.node().getBBox();
+    const btmBox = l.btm.box.node().getBBox();
+    const leftBox = l.left.box.node().getBBox();
+    const rtBox = l.right.box.node().getBBox();
+    const legendBox = self.dom.legendG.node().getBBox();
+    const seriesBox = self.dom.seriesesG.node().getBBox();
+    d.extraWidth = leftBox.width + rtBox.width + s.margin.left + s.margin.right + s.rowlabelgap * 2;
+    d.extraHeight = topBox.height + btmBox.height + s.margin.top + s.margin.bottom + s.collabelgap * 2;
+    d.svgw = d.mainw + d.extraWidth + hcWidth;
+    d.svgh = d.mainh + d.extraHeight + legendBox.height + 20 + s.scrollHeight + hcHeight;
+    self.dom.svg.attr("width", d.svgw).attr("height", d.svgh);
+    let maxLabelWidth = self.type == "hierCluster" ? 0 : leftBox.width, maxLabelNumChars = 0;
+    if (hc.xDendrogramHeight) {
+      self.dom.termLabelG.selectAll(".sjpp-matrix-label").each(function(d2) {
+        if (d2.grp.type !== "hierCluster") return;
+        const box = this.getBBox();
+        if (box.width > maxLabelWidth) {
+          maxLabelWidth = box.width;
+          maxLabelNumChars = d2.label.length;
+        }
+      });
+    }
+    const x = -l.left.offset + hcWidth + maxLabelWidth;
+    const xAdjust = !hc.xDendrogramHeight ? 0 : Math.max(leftBox.width - (hc.xDendrogramHeight + maxLabelWidth), 0);
+    const y = (l.top.display == "none" ? 0 : topBox.height) - l.top.offset + hcHeight;
+    self.dom.mainG.attr("transform", `translate(${x + xAdjust},${y})`);
+    self.dom.clipRect.attr("y", -y).attr("height", d.mainh + 500 + y);
+    const legendX = d.xOffset + (s.transpose ? 20 : 0);
+    const legendY = d.yOffset + d.mainh + s.collabelgap + (l.btm.display == "none" ? 0 : btmBox.height) + 20;
+    self.dom.legendG.attr("transform", `translate(${legendX},${legendY})`);
+    if (hc.xDendrogramHeight) {
+      const dendroX = maxLabelWidth + xAdjust - l.left.offset + d.xOffset - d.dx / 2;
+      self.dom.hcClipRect.attr("x", dendroX + hcWidth + d.dx / 2).attr("y", 0).attr("width", d.mainw + 3).attr("height", d.mainh + hc.yDendrogramHeight + 500);
+      self.topDendroX = dendroX + d.seriesXoffset;
+      self.dom.topDendrogram.attr("transform", `translate(${self.topDendroX}, 0)`);
+      const y2 = l.top.display == "none" ? 0 : topBox.height + s.collabelgap;
+      self.dom.leftDendrogram.attr("transform", `translate(${dendroX - maxLabelWidth - 10}, ${y2})`);
+    }
+  };
+}
+function getRectFill(d) {
+  if (d.fill) return d.fill;
+  const cls = d.class || Array.isArray(d.values) && d.values[0].class;
+  if (!cls) console.log;
+  return cls ? mclass[cls].color : "#555";
+}
+function sleep(ms) {
+  return new Promise((resolve) => setTimeout(resolve, ms));
+}
+
+export {
+  setRenderers
+};
+//# sourceMappingURL=chunk-FUSVXQ6C.js.map

package/dist/chunk-G7JRRTS5.js

@@ -0,0 +1,170 @@
+import {
+  dofetch3
+} from "./chunk-MYWJQRTD.js";
+import {
+  DNA_METHYLATION,
+  GENE_EXPRESSION,
+  PROTEOME_DAP,
+  SINGLECELL_CELLTYPE
+} from "./chunk-FDK4UJMN.js";
+import {
+  rgb
+} from "./chunk-OMR2DT66.js";
+
+// plots/volcano/colors.ts
+function getGroupColors(config) {
+  const groups = config?.samplelst?.groups;
+  const termValues = config?.tw?.term?.values;
+  const rawDown = termValues?.[groups?.[0]?.name]?.color || "red";
+  const rawUp = termValues?.[groups?.[1]?.name]?.color || "blue";
+  return {
+    controlColor: toHex(rawDown, "red"),
+    caseColor: toHex(rawUp, "blue")
+  };
+}
+function toHex(color, fallback) {
+  const c = rgb(color || fallback);
+  return c.displayable() ? c.formatHex() : rgb(fallback).formatHex();
+}
+
+// plots/volcano/model/VolcanoModel.ts
+var VolcanoModel = class {
+  constructor(app, termType) {
+    this.app = app;
+    this.termType = termType;
+  }
+  /** May use mapper instead as more termTypes are added */
+  async getData(config, settings) {
+    this.config = config;
+    this.settings = settings;
+    if (this.termType === GENE_EXPRESSION) {
+      const body = await this.getGERequestBody();
+      const response = await dofetch3("termdb/DE", { body });
+      if (response && !response.error) response.daRequest = body;
+      return response;
+    }
+    if (this.termType === DNA_METHYLATION) {
+      const body = await this.getDMRequestBody();
+      const response = await dofetch3("termdb/diffMeth", { body });
+      if (response && !response.error) response.daRequest = body;
+      return response;
+    }
+    if (this.termType === SINGLECELL_CELLTYPE) {
+      const body = await this.getSCCTRequestBody();
+      return await dofetch3("termdb/singlecellDEgenes", { body });
+    }
+    if (this.termType === PROTEOME_DAP) {
+      const body = this.getDapRequestBody();
+      return await dofetch3("termdb/dapVolcano", { body });
+    }
+    throw new Error(`Volcano plot does not support route for termType='${this.termType}'`);
+  }
+  //Gene expression
+  async getGERequestBody() {
+    await this.getOtherSamples(this.config.samplelst);
+    const state = this.app.getState();
+    const body = {
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      method: this.settings.method,
+      min_count: this.settings.minCount,
+      min_total_count: this.settings.minTotalCount,
+      samplelst: this.config.samplelst,
+      filter: state.termfilter.filter,
+      filter0: state.termfilter.filter0,
+      cpm_cutoff: this.settings.cpmCutoff,
+      volcanoRender: this.getVolcanoRender()
+    };
+    this.addConfounderTw(body);
+    return body;
+  }
+  //DNA methylation
+  async getDMRequestBody() {
+    await this.getOtherSamples(this.config.samplelst);
+    const state = this.app.getState();
+    const body = {
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      samplelst: this.config.samplelst,
+      filter: state.termfilter.filter,
+      filter0: state.termfilter.filter0,
+      min_samples_per_group: this.settings.minSamplesPerGroup,
+      volcanoRender: this.getVolcanoRender()
+    };
+    this.addConfounderTw(body);
+    return body;
+  }
+  /** Parameters telling the server to run the `volcano` Rust renderer and return a
+   * volcano PNG + top-significant rows instead of the full dot list. */
+  getVolcanoRender() {
+    const dotRadius = Math.max(this.settings.width, this.settings.height) / 80;
+    const { caseColor, controlColor } = getGroupColors(this.config);
+    return {
+      significanceThresholds: {
+        pValueCutoff: this.settings.pValue,
+        pValueType: this.settings.pValueType,
+        foldChangeCutoff: this.settings.foldChangeCutoff
+      },
+      pixelWidth: this.settings.width,
+      pixelHeight: this.settings.height,
+      colorSignificant: toHex(this.settings.defaultSignColor, "red"),
+      colorSignificantUp: caseColor,
+      colorSignificantDown: controlColor,
+      colorNonsignificant: toHex(this.settings.defaultNonSignColor, "black"),
+      dotRadius,
+      maxInteractiveDots: this.settings.maxInteractiveDots
+    };
+  }
+  //This is a workaround until the server can accept an arr of confounder tws
+  addConfounderTw(body) {
+    const confounders = this.config?.confounderTws;
+    if (confounders?.length) {
+      body.tw = this.config.confounderTws[0];
+      if (confounders.length > 1) body.tw2 = this.config.confounderTws[1];
+    }
+  }
+  //Single cell cell type
+  getSCCTRequestBody() {
+    const body = {
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      sample: this.config.sample,
+      termId: this.config.termId,
+      categoryName: this.config.categoryName,
+      volcanoRender: this.getVolcanoRender()
+    };
+    return body;
+  }
+  getDapRequestBody() {
+    const { organism, assay, cohort } = this.config.proteomeDetails;
+    return {
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      organism,
+      assay,
+      cohort,
+      volcanoRender: this.getVolcanoRender()
+    };
+  }
+  /** retrieve the sampleId/sampleName for samples in
+   * the "others" group instead of using {in: false} */
+  async getOtherSamples(samplelst) {
+    const othersSamplesGroup = samplelst.groups.find((g) => !g.in);
+    if (!othersSamplesGroup) return;
+    const state = this.app.getState();
+    const samplesGroup = samplelst.groups.find((g) => g.in);
+    othersSamplesGroup.values = [];
+    for (const s of await this.app.vocabApi.getFilteredSampleList(state.termfilter.filter)) {
+      if (samplesGroup.values.indexOf((i) => i.sampleId == s.id) == -1) {
+        othersSamplesGroup.values.push({ sampleId: s.id, sample: s.name });
+      }
+    }
+    othersSamplesGroup.in = true;
+  }
+};
+
+export {
+  getGroupColors,
+  VolcanoModel
+};
+//# sourceMappingURL=chunk-G7JRRTS5.js.map

package/dist/chunk-G7JRRTS5.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../plots/volcano/colors.ts", "../plots/volcano/model/VolcanoModel.ts"],
+  "sourcesContent": ["import { rgb } from 'd3-color'\n\n/** Resolve the case/control dot colors for a volcano plot in one place so the\n * interactive SVG overlay (VolcanoViewModel) and the server-rendered PNG\n * (VolcanoModel \u2192 Rust) paint each dot the same color.\n *\n * `caseColor` maps to points with `fold_change > 0` (group 2 in samplelst),\n * `controlColor` to `fold_change < 0` (group 1). Every returned color is a\n * `#rrggbb` hex string \u2014 CSS names like `'red'` are normalized via d3-color\n * so the Rust renderer's hex-only parser doesn't fall back to a muted tuple.\n */\nexport function getGroupColors(config: any): { caseColor: string; controlColor: string } {\n\tconst groups = config?.samplelst?.groups\n\tconst termValues = config?.tw?.term?.values\n\tconst rawDown = termValues?.[groups?.[0]?.name]?.color || 'red'\n\tconst rawUp = termValues?.[groups?.[1]?.name]?.color || 'blue'\n\treturn {\n\t\tcontrolColor: toHex(rawDown, 'red'),\n\t\tcaseColor: toHex(rawUp, 'blue')\n\t}\n}\n\n/** Normalize any CSS-accepted color string into `#rrggbb`. */\nexport function toHex(color: string | undefined, fallback: string): string {\n\tconst c = rgb(color || fallback)\n\treturn c.displayable() ? c.formatHex() : rgb(fallback).formatHex()\n}\n", "import type { MassAppApi } from '#mass/types/mass'\nimport { dofetch3 } from '#common/dofetch'\nimport type { DERequest, DiffMethRequest, TermdbSingleCellDEgenesRequest, VolcanoRenderRequest } from '#types'\nimport { DNA_METHYLATION, GENE_EXPRESSION, SINGLECELL_CELLTYPE, PROTEOME_DAP } from '#shared/terms.js'\nimport { getGroupColors, toHex } from '../colors'\n\nexport class VolcanoModel {\n\tapp: MassAppApi\n\tconfig!: any\n\tsettings!: any\n\ttermType: string\n\n\tconstructor(app: MassAppApi, termType: string) {\n\t\tthis.app = app\n\t\tthis.termType = termType\n\t}\n\n\t/** May use mapper instead as more termTypes are added */\n\tasync getData(config: any, settings: any) {\n\t\tthis.config = config\n\t\tthis.settings = settings\n\n\t\tif (this.termType === GENE_EXPRESSION) {\n\t\t\tconst body = await this.getGERequestBody()\n\t\t\tconst response = await dofetch3('termdb/DE', { body })\n\t\t\t// Surface the DE request so downstream plots (GSEA) can snapshot\n\t\t\t// it and later ask the server to recompute the DA cache if the\n\t\t\t// file is missing on a peer node or after TTL eviction.\n\t\t\tif (response && !response.error) response.daRequest = body\n\t\t\treturn response\n\t\t}\n\t\tif (this.termType === DNA_METHYLATION) {\n\t\t\tconst body = await this.getDMRequestBody()\n\t\t\tconst response = await dofetch3('termdb/diffMeth', { body })\n\t\t\t// Surface the DM request the same way the GE branch above does so\n\t\t\t// the GSEA tab can snapshot it and the server can recompute the DM\n\t\t\t// cache if the file is missing on a peer node or after TTL.\n\t\t\tif (response && !response.error) response.daRequest = body\n\t\t\treturn response\n\t\t}\n\t\tif (this.termType === SINGLECELL_CELLTYPE) {\n\t\t\tconst body = await this.getSCCTRequestBody()\n\t\t\treturn await dofetch3('termdb/singlecellDEgenes', { body })\n\t\t}\n\t\tif (this.termType === PROTEOME_DAP) {\n\t\t\tconst body = this.getDapRequestBody()\n\t\t\treturn await dofetch3('termdb/dapVolcano', { body })\n\t\t}\n\t\tthrow new Error(`Volcano plot does not support route for termType='${this.termType}'`)\n\t}\n\n\t//Gene expression\n\tasync getGERequestBody() {\n\t\tawait this.getOtherSamples(this.config.samplelst)\n\t\tconst state = this.app.getState()\n\t\tconst body = {\n\t\t\tgenome: this.app.vocabApi.vocab.genome,\n\t\t\tdslabel: this.app.vocabApi.vocab.dslabel,\n\t\t\tmethod: this.settings.method,\n\t\t\tmin_count: this.settings.minCount,\n\t\t\tmin_total_count: this.settings.minTotalCount,\n\t\t\tsamplelst: this.config.samplelst,\n\t\t\tfilter: state.termfilter.filter,\n\t\t\tfilter0: state.termfilter.filter0,\n\t\t\tcpm_cutoff: this.settings.cpmCutoff,\n\t\t\tvolcanoRender: this.getVolcanoRender()\n\t\t} as Partial<DERequest> //remove Partial when storage_type is removed from DERequest\n\n\t\tthis.addConfounderTw(body)\n\n\t\treturn body\n\t}\n\n\t//DNA methylation\n\tasync getDMRequestBody() {\n\t\tawait this.getOtherSamples(this.config.samplelst)\n\t\tconst state = this.app.getState()\n\t\tconst body = {\n\t\t\tgenome: this.app.vocabApi.vocab.genome,\n\t\t\tdslabel: this.app.vocabApi.vocab.dslabel,\n\t\t\tsamplelst: this.config.samplelst,\n\t\t\tfilter: state.termfilter.filter,\n\t\t\tfilter0: state.termfilter.filter0,\n\t\t\tmin_samples_per_group: this.settings.minSamplesPerGroup,\n\t\t\tvolcanoRender: this.getVolcanoRender()\n\t\t} as Partial<DiffMethRequest>\n\n\t\tthis.addConfounderTw(body)\n\n\t\treturn body\n\t}\n\n\t/** Parameters telling the server to run the `volcano` Rust renderer and return a\n\t * volcano PNG + top-significant rows instead of the full dot list. */\n\tgetVolcanoRender(): VolcanoRenderRequest {\n\t\t// Match the client overlay's radius (see VolcanoViewModel.setPointData)\n\t\t// so the PNG rings and the interactive overlay rings line up; otherwise\n\t\t// a smaller PNG ring sits inside the larger overlay ring and looks like\n\t\t// a stray dot at the center.\n\t\tconst dotRadius = Math.max(this.settings.width, this.settings.height) / 80\n\t\t// Resolve case/control colors via the shared helper (see colors.ts) so the\n\t\t// PNG and the SVG overlay paint each side with the exact same hex string.\n\t\tconst { caseColor, controlColor } = getGroupColors(this.config)\n\t\treturn {\n\t\t\tsignificanceThresholds: {\n\t\t\t\tpValueCutoff: this.settings.pValue,\n\t\t\t\tpValueType: this.settings.pValueType,\n\t\t\t\tfoldChangeCutoff: this.settings.foldChangeCutoff\n\t\t\t},\n\t\t\tpixelWidth: this.settings.width,\n\t\t\tpixelHeight: this.settings.height,\n\t\t\tcolorSignificant: toHex(this.settings.defaultSignColor, 'red'),\n\t\t\tcolorSignificantUp: caseColor,\n\t\t\tcolorSignificantDown: controlColor,\n\t\t\tcolorNonsignificant: toHex(this.settings.defaultNonSignColor, 'black'),\n\t\t\tdotRadius,\n\t\t\tmaxInteractiveDots: this.settings.maxInteractiveDots\n\t\t}\n\t}\n\n\t//This is a workaround until the server can accept an arr of confounder tws\n\taddConfounderTw(body) {\n\t\tconst confounders = this.config?.confounderTws\n\t\tif (confounders?.length) {\n\t\t\tbody.tw = this.config.confounderTws[0]\n\t\t\tif (confounders.length > 1) body.tw2 = this.config.confounderTws[1]\n\t\t}\n\t}\n\n\t//Single cell cell type\n\tgetSCCTRequestBody(): TermdbSingleCellDEgenesRequest {\n\t\tconst body = {\n\t\t\tgenome: this.app.vocabApi.vocab.genome,\n\t\t\tdslabel: this.app.vocabApi.vocab.dslabel,\n\t\t\tsample: this.config.sample,\n\t\t\ttermId: this.config.termId,\n\t\t\tcategoryName: this.config.categoryName,\n\t\t\tvolcanoRender: this.getVolcanoRender()\n\t\t}\n\t\treturn body\n\t}\n\n\tgetDapRequestBody() {\n\t\tconst { organism, assay, cohort } = this.config.proteomeDetails\n\t\treturn {\n\t\t\tgenome: this.app.vocabApi.vocab.genome,\n\t\t\tdslabel: this.app.vocabApi.vocab.dslabel,\n\t\t\torganism,\n\t\t\tassay,\n\t\t\tcohort,\n\t\t\tvolcanoRender: this.getVolcanoRender()\n\t\t}\n\t}\n\n\t/** retrieve the sampleId/sampleName for samples in\n\t * the \"others\" group instead of using {in: false} */\n\tasync getOtherSamples(samplelst) {\n\t\tconst othersSamplesGroup = samplelst.groups.find(g => !g.in)\n\t\tif (!othersSamplesGroup) return\n\n\t\tconst state = this.app.getState()\n\t\tconst samplesGroup = samplelst.groups.find(g => g.in)\n\t\tothersSamplesGroup.values = []\n\t\t// retrieve full list of samples based on current filter. put samples not in samplesGroup in \"others\" group\n\t\tfor (const s of await this.app.vocabApi.getFilteredSampleList(state.termfilter.filter)) {\n\t\t\t// s={id,name}, samplelst.groups[].values[]={sampleId,sample}\n\t\t\tif (samplesGroup.values.indexOf(i => i.sampleId == s.id) == -1) {\n\t\t\t\tothersSamplesGroup.values.push({ sampleId: s.id, sample: s.name })\n\t\t\t}\n\t\t}\n\t\tothersSamplesGroup.in = true\n\t}\n}\n"],
+  "mappings": ";;;;;;;;;;;;;;AAWO,SAAS,eAAe,QAA0D;AACxF,QAAM,SAAS,QAAQ,WAAW;AAClC,QAAM,aAAa,QAAQ,IAAI,MAAM;AACrC,QAAM,UAAU,aAAa,SAAS,CAAC,GAAG,IAAI,GAAG,SAAS;AAC1D,QAAM,QAAQ,aAAa,SAAS,CAAC,GAAG,IAAI,GAAG,SAAS;AACxD,SAAO;AAAA,IACN,cAAc,MAAM,SAAS,KAAK;AAAA,IAClC,WAAW,MAAM,OAAO,MAAM;AAAA,EAC/B;AACD;AAGO,SAAS,MAAM,OAA2B,UAA0B;AAC1E,QAAM,IAAI,IAAI,SAAS,QAAQ;AAC/B,SAAO,EAAE,YAAY,IAAI,EAAE,UAAU,IAAI,IAAI,QAAQ,EAAE,UAAU;AAClE;;;ACpBO,IAAM,eAAN,MAAmB;AAAA,EAMzB,YAAY,KAAiB,UAAkB;AAC9C,SAAK,MAAM;AACX,SAAK,WAAW;AAAA,EACjB;AAAA;AAAA,EAGA,MAAM,QAAQ,QAAa,UAAe;AACzC,SAAK,SAAS;AACd,SAAK,WAAW;AAEhB,QAAI,KAAK,aAAa,iBAAiB;AACtC,YAAM,OAAO,MAAM,KAAK,iBAAiB;AACzC,YAAM,WAAW,MAAM,SAAS,aAAa,EAAE,KAAK,CAAC;AAIrD,UAAI,YAAY,CAAC,SAAS,MAAO,UAAS,YAAY;AACtD,aAAO;AAAA,IACR;AACA,QAAI,KAAK,aAAa,iBAAiB;AACtC,YAAM,OAAO,MAAM,KAAK,iBAAiB;AACzC,YAAM,WAAW,MAAM,SAAS,mBAAmB,EAAE,KAAK,CAAC;AAI3D,UAAI,YAAY,CAAC,SAAS,MAAO,UAAS,YAAY;AACtD,aAAO;AAAA,IACR;AACA,QAAI,KAAK,aAAa,qBAAqB;AAC1C,YAAM,OAAO,MAAM,KAAK,mBAAmB;AAC3C,aAAO,MAAM,SAAS,4BAA4B,EAAE,KAAK,CAAC;AAAA,IAC3D;AACA,QAAI,KAAK,aAAa,cAAc;AACnC,YAAM,OAAO,KAAK,kBAAkB;AACpC,aAAO,MAAM,SAAS,qBAAqB,EAAE,KAAK,CAAC;AAAA,IACpD;AACA,UAAM,IAAI,MAAM,qDAAqD,KAAK,QAAQ,GAAG;AAAA,EACtF;AAAA;AAAA,EAGA,MAAM,mBAAmB;AACxB,UAAM,KAAK,gBAAgB,KAAK,OAAO,SAAS;AAChD,UAAM,QAAQ,KAAK,IAAI,SAAS;AAChC,UAAM,OAAO;AAAA,MACZ,QAAQ,KAAK,IAAI,SAAS,MAAM;AAAA,MAChC,SAAS,KAAK,IAAI,SAAS,MAAM;AAAA,MACjC,QAAQ,KAAK,SAAS;AAAA,MACtB,WAAW,KAAK,SAAS;AAAA,MACzB,iBAAiB,KAAK,SAAS;AAAA,MAC/B,WAAW,KAAK,OAAO;AAAA,MACvB,QAAQ,MAAM,WAAW;AAAA,MACzB,SAAS,MAAM,WAAW;AAAA,MAC1B,YAAY,KAAK,SAAS;AAAA,MAC1B,eAAe,KAAK,iBAAiB;AAAA,IACtC;AAEA,SAAK,gBAAgB,IAAI;AAEzB,WAAO;AAAA,EACR;AAAA;AAAA,EAGA,MAAM,mBAAmB;AACxB,UAAM,KAAK,gBAAgB,KAAK,OAAO,SAAS;AAChD,UAAM,QAAQ,KAAK,IAAI,SAAS;AAChC,UAAM,OAAO;AAAA,MACZ,QAAQ,KAAK,IAAI,SAAS,MAAM;AAAA,MAChC,SAAS,KAAK,IAAI,SAAS,MAAM;AAAA,MACjC,WAAW,KAAK,OAAO;AAAA,MACvB,QAAQ,MAAM,WAAW;AAAA,MACzB,SAAS,MAAM,WAAW;AAAA,MAC1B,uBAAuB,KAAK,SAAS;AAAA,MACrC,eAAe,KAAK,iBAAiB;AAAA,IACtC;AAEA,SAAK,gBAAgB,IAAI;AAEzB,WAAO;AAAA,EACR;AAAA;AAAA;AAAA,EAIA,mBAAyC;AAKxC,UAAM,YAAY,KAAK,IAAI,KAAK,SAAS,OAAO,KAAK,SAAS,MAAM,IAAI;AAGxE,UAAM,EAAE,WAAW,aAAa,IAAI,eAAe,KAAK,MAAM;AAC9D,WAAO;AAAA,MACN,wBAAwB;AAAA,QACvB,cAAc,KAAK,SAAS;AAAA,QAC5B,YAAY,KAAK,SAAS;AAAA,QAC1B,kBAAkB,KAAK,SAAS;AAAA,MACjC;AAAA,MACA,YAAY,KAAK,SAAS;AAAA,MAC1B,aAAa,KAAK,SAAS;AAAA,MAC3B,kBAAkB,MAAM,KAAK,SAAS,kBAAkB,KAAK;AAAA,MAC7D,oBAAoB;AAAA,MACpB,sBAAsB;AAAA,MACtB,qBAAqB,MAAM,KAAK,SAAS,qBAAqB,OAAO;AAAA,MACrE;AAAA,MACA,oBAAoB,KAAK,SAAS;AAAA,IACnC;AAAA,EACD;AAAA;AAAA,EAGA,gBAAgB,MAAM;AACrB,UAAM,cAAc,KAAK,QAAQ;AACjC,QAAI,aAAa,QAAQ;AACxB,WAAK,KAAK,KAAK,OAAO,cAAc,CAAC;AACrC,UAAI,YAAY,SAAS,EAAG,MAAK,MAAM,KAAK,OAAO,cAAc,CAAC;AAAA,IACnE;AAAA,EACD;AAAA;AAAA,EAGA,qBAAqD;AACpD,UAAM,OAAO;AAAA,MACZ,QAAQ,KAAK,IAAI,SAAS,MAAM;AAAA,MAChC,SAAS,KAAK,IAAI,SAAS,MAAM;AAAA,MACjC,QAAQ,KAAK,OAAO;AAAA,MACpB,QAAQ,KAAK,OAAO;AAAA,MACpB,cAAc,KAAK,OAAO;AAAA,MAC1B,eAAe,KAAK,iBAAiB;AAAA,IACtC;AACA,WAAO;AAAA,EACR;AAAA,EAEA,oBAAoB;AACnB,UAAM,EAAE,UAAU,OAAO,OAAO,IAAI,KAAK,OAAO;AAChD,WAAO;AAAA,MACN,QAAQ,KAAK,IAAI,SAAS,MAAM;AAAA,MAChC,SAAS,KAAK,IAAI,SAAS,MAAM;AAAA,MACjC;AAAA,MACA;AAAA,MACA;AAAA,MACA,eAAe,KAAK,iBAAiB;AAAA,IACtC;AAAA,EACD;AAAA;AAAA;AAAA,EAIA,MAAM,gBAAgB,WAAW;AAChC,UAAM,qBAAqB,UAAU,OAAO,KAAK,OAAK,CAAC,EAAE,EAAE;AAC3D,QAAI,CAAC,mBAAoB;AAEzB,UAAM,QAAQ,KAAK,IAAI,SAAS;AAChC,UAAM,eAAe,UAAU,OAAO,KAAK,OAAK,EAAE,EAAE;AACpD,uBAAmB,SAAS,CAAC;AAE7B,eAAW,KAAK,MAAM,KAAK,IAAI,SAAS,sBAAsB,MAAM,WAAW,MAAM,GAAG;AAEvF,UAAI,aAAa,OAAO,QAAQ,OAAK,EAAE,YAAY,EAAE,EAAE,KAAK,IAAI;AAC/D,2BAAmB,OAAO,KAAK,EAAE,UAAU,EAAE,IAAI,QAAQ,EAAE,KAAK,CAAC;AAAA,MAClE;AAAA,IACD;AACA,uBAAmB,KAAK;AAAA,EACzB;AACD;",
+  "names": []
+}