@sjcrh/proteinpaint-client 2.207.1 → 2.209.0

package/dist/cellTypeBubbleHeatmap-ZPTAGEWW.js

@@ -1,278 +0,0 @@
-import {
-  LegendCircleReference,
-  PlotBase,
-  addGeneSearchbox
-} from "./chunk-PC4MFDHP.js";
-import "./chunk-HJ6L54YS.js";
-import "./chunk-KV4W2ACA.js";
-import "./chunk-HPAW7XDM.js";
-import {
-  Menu
-} from "./chunk-ELJX3QIQ.js";
-import "./chunk-BZN2O76M.js";
-import "./chunk-EEB5VE2A.js";
-import "./chunk-6RRZRISL.js";
-import "./chunk-2KM4PRQM.js";
-import {
-  dofetch3
-} from "./chunk-52QHIKH2.js";
-import "./chunk-A2ORIMUJ.js";
-import "./chunk-PPSWNLMG.js";
-import "./chunk-RUBZCKIX.js";
-import {
-  copyMerge,
-  getCompInit
-} from "./chunk-WINIL2KN.js";
-import "./chunk-PF4DSFDR.js";
-import "./chunk-7X6NF7NI.js";
-import "./chunk-W5J3LTYS.js";
-import "./chunk-Z2ZITHT4.js";
-import {
-  linear,
-  sqrt
-} from "./chunk-4OLM3KSB.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/cellTypeBubbleHeatmap.ts
-var defaultConfig = { chartType: "cellTypeBubbleHeatmap" };
-var CELL_W = 84;
-var CELL_H = 60;
-var ROW_LABEL_W = 74;
-var GROUP_LABEL_H = 22;
-var GENO_LABEL_H = 40;
-var COL_LABEL_H = GROUP_LABEL_H + GENO_LABEL_H;
-var MIN_DOT_R = 8;
-var MAX_DOT_R = 22;
-var NEG_LOG_FDR_CAP = 10;
-var COLOR_NEG = "#762a83";
-var COLOR_ZERO = "#f7f7f7";
-var COLOR_POS = "#2166ac";
-var CellTypeBubbleHeatmap = class _CellTypeBubbleHeatmap extends PlotBase {
-  constructor(opts, api) {
-    super(opts, api);
-    this.currentIsoform = "";
-    this.type = _CellTypeBubbleHeatmap.type;
-  }
-  static {
-    this.type = "cellTypeBubbleHeatmap";
-  }
-  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("Cell-type 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("cellTypeBubbleHeatmap: gene is missing");
-    if (this.dom.header) this.dom.header.text(`Cell-type 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/cellTypeBubbleHeatmap", { 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 cohort DAPfile.`);
-      return;
-    }
-    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 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 columns = data.columns;
-    const rows = data.rows;
-    const nCols = columns.length;
-    const nRows = rows.length;
-    const negLogFdr = (fdr) => fdr > 0 ? Math.min(-Math.log10(fdr), NEG_LOG_FDR_CAP) : NEG_LOG_FDR_CAP;
-    const cellOf = (colKey, rowKey) => isoformData.data[colKey]?.[rowKey];
-    let maxAbs = 0;
-    const thresholdNegLog = negLogFdr(threshold);
-    let maxNegLog = thresholdNegLog;
-    for (const col of columns) {
-      for (const row of rows) {
-        const s = cellOf(col.key, row.key);
-        if (!s) continue;
-        const v = Math.abs(s.log2FC);
-        if (v > maxAbs) maxAbs = v;
-        const nl = negLogFdr(s.fdr);
-        if (nl > maxNegLog) maxNegLog = nl;
-      }
-    }
-    if (maxAbs === 0) maxAbs = 1;
-    if (maxNegLog <= thresholdNegLog) maxNegLog = thresholdNegLog + 1;
-    const colorScale = linear().domain([-maxAbs, 0, maxAbs]).range([COLOR_NEG, COLOR_ZERO, COLOR_POS]).clamp(true);
-    const sizeScale = sqrt().domain([thresholdNegLog, maxNegLog]).range([MIN_DOT_R, MAX_DOT_R]).clamp(true);
-    const gridW = ROW_LABEL_W + nCols * CELL_W + 20;
-    const gridH = COL_LABEL_H + nRows * CELL_H + 20;
-    const svg = container.append("svg").attr("width", gridW).attr("height", gridH).style("flex", "0 0 auto");
-    const grid = svg.append("g");
-    let c = 0;
-    while (c < nCols) {
-      const cellType = columns[c].cellType;
-      let end = c;
-      while (end + 1 < nCols && columns[end + 1].cellType === cellType) end++;
-      const xStart = ROW_LABEL_W + c * CELL_W;
-      const xEnd = ROW_LABEL_W + (end + 1) * CELL_W;
-      const xMid = (xStart + xEnd) / 2;
-      grid.append("text").attr("x", xMid).attr("y", GROUP_LABEL_H - 7).attr("text-anchor", "middle").attr("font-size", "13px").attr("font-weight", "bold").text(cellType);
-      grid.append("line").attr("x1", xStart + 4).attr("y1", GROUP_LABEL_H - 3).attr("x2", xEnd - 4).attr("y2", GROUP_LABEL_H - 3).attr("stroke", "#bbb").attr("stroke-width", 1);
-      c = end + 1;
-    }
-    for (let col = 0; col < nCols; col++) {
-      const cx = ROW_LABEL_W + col * CELL_W + CELL_W / 2;
-      grid.append("text").attr("x", cx).attr("y", COL_LABEL_H - 14).attr("text-anchor", "middle").attr("font-size", "12px").attr("font-weight", "600").text(columns[col].genotype);
-    }
-    for (let r = 0; r < nRows; r++) {
-      const cy = COL_LABEL_H + r * CELL_H + CELL_H / 2;
-      grid.append("text").attr("x", ROW_LABEL_W - 12).attr("y", cy).attr("text-anchor", "end").attr("dominant-baseline", "central").attr("font-size", "13px").attr("font-weight", "bold").text(rows[r].label);
-    }
-    for (let r = 0; r < nRows; r++) {
-      for (let col = 0; col < nCols; col++) {
-        const x0 = ROW_LABEL_W + col * CELL_W;
-        const y0 = COL_LABEL_H + r * CELL_H;
-        grid.append("rect").attr("x", x0).attr("y", y0).attr("width", CELL_W).attr("height", CELL_H).attr("fill", "none").attr("stroke", "#eee").attr("stroke-width", 1);
-        const s = cellOf(columns[col].key, rows[r].key);
-        if (!s) continue;
-        const cx = x0 + CELL_W / 2;
-        const cy = y0 + CELL_H / 2;
-        grid.append("circle").attr("cx", cx).attr("cy", cy).attr("r", sizeScale(negLogFdr(s.fdr))).attr("fill", colorScale(s.log2FC)).attr("stroke", "#888").attr("stroke-width", 0.8).style("opacity", s.significant ? 1 : 0.35).on(
-          "mouseover",
-          (event) => this.showCellTip(event, isoformData.gene_name, selectedIsoform, columns[col], rows[r], s)
-        ).on("mouseout", () => this.dom.tip.hide());
-      }
-    }
-    this.renderLegend(container, colorScale, maxAbs, threshold, maxNegLog);
-  }
-  fmtFdr(v) {
-    return v >= 1e-4 ? v.toFixed(4) : v.toExponential(2);
-  }
-  showCellTip(event, geneName, isoform, col, row, s) {
-    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(`Cell type: ${col.cellType}`);
-    t.append("div").text(`Genotype: ${col.genotype}`);
-    t.append("div").text(`Timepoint: ${row.label}`);
-    t.append("div").text(`Protein: ${s.id}`);
-    t.append("div").text(`log\u2082FC: ${s.log2FC.toFixed(3)}`);
-    t.append("div").text(`FDR: ${this.fmtFdr(s.fdr)}${s.significant ? "" : " (n.s.)"}`);
-    t.append("div").style("color", "#666").style("margin-top", "4px").text("Color = log\u2082FC (blue up / purple down). Size = \u2212log\u2081\u2080 FDR.");
-  }
-  renderLegend(container, colorScale, maxAbs, threshold, 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("log\u2082FC");
-    const cW = 22;
-    const cH = 130;
-    const cSvg = colorBlock.append("svg").attr("width", cW + 80).attr("height", cH + 16);
-    const gid = `ctbh-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 > 0 ? "+" : ""}${tick.toFixed(2)}`);
-    }
-    colorBlock.append("div").style("font-size", "11px").style("color", "#666").style("margin-top", "2px").text("blue = up (+), purple = down (\u2212)");
-    const sizeBlock = legend.append("div");
-    sizeBlock.append("div").style("font-weight", "bold").style("font-size", "13px").style("margin-bottom", "6px").text("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,
-      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));
-    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 (blue up, purple down). Dot size = significance, \u2212log\u2081\u2080 FDR; the smallest size marks the FDR < ${threshold} cutoff. Non-significant dots (FDR \u2265 ${threshold}) are faded.`
-    );
-    notes.append("div").style("margin-top", "4px").text("An empty cell means the cohort was not assayed (e.g. OPC has no 4m) or the protein was not detected.");
-  }
-};
-var componentInit = getCompInit(CellTypeBubbleHeatmap);
-async function getPlotConfig(opts) {
-  const config = structuredClone(defaultConfig);
-  if (!opts.gene) throw new Error("cellTypeBubbleHeatmap 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: "cellTypeBubbleHeatmap",
-          gene: geneSearch.geneSymbol
-        }
-      });
-    }
-  });
-}
-export {
-  componentInit,
-  getPlotConfig,
-  makeChartBtnMenu
-};
-//# sourceMappingURL=cellTypeBubbleHeatmap-ZPTAGEWW.js.map

package/dist/chunk-2GBG3KA7.js

@@ -1,129 +0,0 @@
-import {
-  keyupEnter
-} from "./chunk-PC4MFDHP.js";
-
-// src/block.mds.svcnv.share.js
-function rnabamtk_initparam(c) {
-  if (!c.dna_mintotalreads) c.dna_mintotalreads = 8;
-  if (!c.rna_mintotalreads) c.rna_mintotalreads = 8;
-  if (!c.hetsnp_minbaf) c.hetsnp_minbaf = 0.3;
-  if (!c.hetsnp_maxbaf) c.hetsnp_maxbaf = 0.7;
-  if (c.rnapileup_q == void 0) c.rnapileup_q = 0;
-  if (!c.rnapileup_Q) c.rnapileup_Q = 13;
-  if (!c.binompvaluecutoff) c.binompvaluecutoff = 0.05;
-  if (!c.clientcolor_snpinuse) c.clientcolor_snpinuse = "blue";
-  if (!c.clientcolor_markernotinuse) c.clientcolor_markernotinuse = "#bbb";
-}
-function configPanel_rnabam(tk, block, loadTk) {
-  const c = tk.checkrnabam;
-  if (!c) return;
-  tk.tkconfigtip.d.append("hr");
-  const d = tk.tkconfigtip.d.append("div").style("margin", "15px 0px");
-  d.append("div").style("opacity", 0.5).style("font-size", ".9em").text("Finding heterozygous SNPs in DNA");
-  {
-    const row = d.append("div").style("margin-top", "5px");
-    row.append("span").html("DNA minimum total read count ");
-    row.append("input").attr("type", "number").style("width", "50px").property("value", c.dna_mintotalreads).on("keyup", (event) => {
-      if (!keyupEnter(event)) return;
-      let v = Number.parseInt(event.target.value);
-      if (!v || v <= 0) return;
-      if (c.dna_mintotalreads == v) {
-        return;
-      }
-      c.dna_mintotalreads = v;
-      loadTk(tk, block);
-    });
-    row.append("div").style("opacity", ".5").style("font-size", ".8em").text("If a SNP's total coverage is below cutoff, it will be skipped.");
-  }
-  {
-    const row = d.append("div").style("margin-top", "5px");
-    row.append("span").html("Heterozygous SNP BAF range&nbsp;&nbsp;");
-    row.append("input").attr("type", "number").style("width", "50px").property("value", c.hetsnp_minbaf).on("keyup", (event) => {
-      if (!keyupEnter(event)) return;
-      let v = Number.parseFloat(event.target.value);
-      if (!v || v <= 0) return;
-      if (c.hetsnp_minbaf == v) {
-        return;
-      }
-      c.hetsnp_minbaf = v;
-      loadTk(tk, block);
-    });
-    row.append("span").style("opacity", ".5").style("font-size", ".8em").html("&nbsp;&leq; BAF &leq;&nbsp;");
-    row.append("input").attr("type", "number").style("width", "50px").property("value", c.hetsnp_maxbaf).on("keyup", (event) => {
-      if (!keyupEnter(event)) return;
-      let v = Number.parseFloat(event.target.value);
-      if (!v || v <= 0) return;
-      if (c.hetsnp_maxbaf == v) {
-        return;
-      }
-      c.hetsnp_maxbaf = v;
-      loadTk(tk, block);
-    });
-    row.append("div").style("opacity", ".5").style("font-size", ".8em").text("If a SNP's BAF (B-allele fraction) is within this range, it is heterozygous.");
-  }
-  d.append("div").style("margin-top", "20px").style("opacity", 0.5).style("font-size", ".9em").text("Counting alleles in RNA-seq BAM file");
-  {
-    const row = d.append("div").style("margin-top", "5px");
-    row.append("span").html("Skip alignments with mapQ smaller than&nbsp;");
-    row.append("input").attr("type", "number").style("width", "50px").property("value", c.rnapileup_q).on("keyup", (event) => {
-      if (!keyupEnter(event)) return;
-      let v = Number.parseInt(event.target.value);
-      if (!v || v < 0) return;
-      if (c.rnapileup_q == v) {
-        return;
-      }
-      c.rnapileup_q = v;
-      loadTk(tk, block);
-    });
-  }
-  {
-    const row = d.append("div").style("margin-top", "5px");
-    row.append("span").html("Skip bases with baseQ/BAQ smaller than&nbsp;");
-    row.append("input").attr("type", "number").style("width", "50px").property("value", c.rnapileup_Q).on("keyup", (event) => {
-      if (!keyupEnter(event)) return;
-      let v = Number.parseInt(event.target.value);
-      if (!v || v <= 0) return;
-      if (c.rnapileup_Q == v) {
-        return;
-      }
-      c.rnapileup_Q = v;
-      loadTk(tk, block);
-    });
-  }
-  d.append("div").style("margin-top", "20px").style("opacity", 0.5).style("font-size", ".9em").text("Binomial test on whether a heterozygous SNP shows allelic bias in RNA");
-  {
-    const row = d.append("div").style("margin-top", "5px");
-    row.append("span").html("P-value cutoff&nbsp;");
-    row.append("input").attr("type", "number").style("width", "50px").property("value", c.binompvaluecutoff).on("keyup", (event) => {
-      if (!keyupEnter(event)) return;
-      let v = Number.parseFloat(event.target.value);
-      if (!v || v <= 0 || v >= 1) return;
-      if (c.binompvaluecutoff == v) {
-        return;
-      }
-      c.binompvaluecutoff = v;
-      loadTk(tk, block);
-    });
-  }
-  {
-    const row = d.append("div").style("margin-top", "5px");
-    row.append("span").html("RNA minimum total read count&nbsp;");
-    row.append("input").attr("type", "number").style("width", "50px").property("value", c.rna_mintotalreads).on("keyup", (event) => {
-      if (!keyupEnter(event)) return;
-      let v = Number.parseInt(event.target.value);
-      if (!v || v <= 0) return;
-      if (c.rna_mintotalreads == v) {
-        return;
-      }
-      c.rna_mintotalreads = v;
-      loadTk(tk, block);
-    });
-    row.append("div").style("opacity", ".5").style("font-size", ".8em").text("If a SNP's total read count from RNA is below cutoff, it won't do binomial test.");
-  }
-}
-
-export {
-  rnabamtk_initparam,
-  configPanel_rnabam
-};
-//# sourceMappingURL=chunk-2GBG3KA7.js.map