@sjcrh/proteinpaint-client 2.224.0 → 2.224.1

package/dist/chunk-JDRUY7IR.js

@@ -0,0 +1,266 @@
+import {
+  DATermTypes,
+  DMRCATE_DEFAULTS,
+  getOthersGroup
+} from "./chunk-3YQYWX4M.js";
+import {
+  dofetch3
+} from "./chunk-KBTFT76L.js";
+import {
+  DMR_SCAN_ELEMENT_TYPE
+} from "./chunk-WJLRU5OF.js";
+import {
+  rgb
+} from "./chunk-Q5RDQNIT.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 {
+  /** TODO: This model is used in both the volcano and gsea.
+   * In the future, create base model in DA and use specific
+   * classes for the volcano and gsea. */
+  constructor(plot, termType) {
+    this.plot = plot;
+    this.app = plot.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 === DATermTypes.GENE_EXPRESSION) {
+      const body = await this.getGERequestBody();
+      const response = await dofetch3("termdb/DE", { body, signal: this.plot.api?.getAbortSignal() });
+      if (response && !response.error) {
+        response.daRequest = body;
+        this.setGroupSizes(response);
+      }
+      return response;
+    }
+    if (this.termType === DATermTypes.DNA_METHYLATION) {
+      const body = await this.getDMRequestBody();
+      const response = await dofetch3("termdb/diffMeth", { body, signal: this.plot.api?.getAbortSignal() });
+      if (response && !response.error) {
+        response.daRequest = body;
+        this.setGroupSizes(response);
+      }
+      return response;
+    }
+    if (this.termType === DATermTypes.JUNCTION) {
+      const body = await this.getDSRequestBody();
+      const response = await dofetch3("termdb/diffSplice", { body, signal: this.plot.api?.getAbortSignal() });
+      if (response && !response.error) {
+        response.daRequest = body;
+        this.setGroupSizes(response);
+      }
+      return response;
+    }
+    if (this.termType === DATermTypes.SINGLECELL_CELLTYPE) {
+      const body = await this.getSCCTRequestBody();
+      return await dofetch3("termdb/singlecellDEgenes", { body, signal: this.plot.api?.getAbortSignal() });
+    }
+    if (this.termType === DATermTypes.PROTEOME_DAP) {
+      const body = this.getDapRequestBody();
+      return await dofetch3("termdb/dapVolcano", { body, signal: this.plot.api?.getAbortSignal() });
+    }
+    if (this.termType === DATermTypes.SINGLECELL_GENE_EXPRESSION) {
+    }
+    throw new Error(`Volcano plot does not support route for termType='${this.termType}'`);
+  }
+  //Gene expression
+  async getGERequestBody() {
+    this.setOthersGroup(this.config.samplelst);
+    const state = this.app.getState();
+    const body = {
+      kind: "DE",
+      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()
+    };
+    const pseudobulk = this.config.tw?.pseudobulk;
+    if (pseudobulk) body.pseudobulk = pseudobulk;
+    this.addConfounderTw(body);
+    return body;
+  }
+  //DNA methylation
+  async getDMRequestBody() {
+    this.setOthersGroup(this.config.samplelst);
+    const state = this.app.getState();
+    const body = {
+      kind: "DM",
+      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,
+      exclude_sex_chr: this.settings.excludeSexChr,
+      /* Omitted rather than sent as 'promoter' when it is the default, so a request
+      from a promoter-only dataset is byte-identical to what this client sent before
+      the element picker existed. The server resolves an absent element_type to
+      'promoter'. This does NOT preserve cache keys -- the key object gained the
+      field server-side, so every pre-existing dm/ entry is orphaned on deploy
+      regardless of what the client sends. */
+      ...this.settings.elementType && this.settings.elementType != "promoter" ? { element_type: this.settings.elementType } : {},
+      // the scan's own knobs; the server ignores them for any other element type
+      ...this.settings.elementType == DMR_SCAN_ELEMENT_TYPE ? {
+        scan: {
+          ...this.settings.scanChromosome ? { chromosome: this.settings.scanChromosome } : {},
+          backgroundCorrection: !!this.settings.backgroundCorrection,
+          minCpgs: this.settings.minCpgs,
+          profileBinBp: this.settings.profileBinBp,
+          // only when changed: an explicit default would orphan every cached scan
+          ...this.settings.lambda != DMRCATE_DEFAULTS.lambda ? { lambda: this.settings.lambda } : {},
+          ...this.settings.C != DMRCATE_DEFAULTS.C ? { C: this.settings.C } : {},
+          ...this.settings.fdrCutoff != DMRCATE_DEFAULTS.fdrCutoff ? { fdrCutoff: this.settings.fdrCutoff } : {}
+        }
+      } : {},
+      volcanoRender: this.getVolcanoRender()
+    };
+    if (this.settings.elementType != DMR_SCAN_ELEMENT_TYPE) this.addConfounderTw(body);
+    return body;
+  }
+  //Splicing
+  async getDSRequestBody() {
+    this.setOthersGroup(this.config.samplelst);
+    const state = this.app.getState();
+    const body = {
+      kind: "DS",
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      samplelst: this.config.samplelst,
+      filter: state.termfilter.filter,
+      filter0: state.termfilter.filter0,
+      method: this.settings.method,
+      /* leafcutter's own cluster filters, applied by the edgeR screen too so both engines
+      rank the same cluster set. Sent already scaled to cohort size by validateDSSettings --
+      as absolute counts they silently loosen as N grows. */
+      minSamplesPerIntron: this.settings.minSamplesPerIntron,
+      minSamplesPerGroup: this.settings.minSamplesPerGroup,
+      minCountsPerCluster: this.settings.minCountsPerCluster,
+      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);
+    const useDeltaBeta = this.termType === DATermTypes.DNA_METHYLATION && this.settings.xAxis == "delta_beta";
+    const render = {
+      significanceThresholds: {
+        pValueCutoff: this.settings.pValue,
+        pValueType: this.settings.pValueType,
+        foldChangeCutoff: useDeltaBeta ? this.settings.deltaBetaCutoff : this.settings.foldChangeCutoff
+      },
+      ...useDeltaBeta ? { xField: "delta_beta" } : {},
+      // Tied to the delta-beta axis: the control that sets it is only offered for
+      // methylation, and centering a log2 fold-change axis is a different conversation.
+      ...useDeltaBeta && this.settings.centerDeltaBeta ? { centerX: true } : {},
+      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,
+      // Render the PNG at device-pixel resolution so it stays sharp on
+      // retina screens. The server reports the plot extent in CSS-space,
+      // so SVG overlay coords are unaffected.
+      //
+      // Oversample by 2× so the PNG also stays sharp when the user
+      // *zooms in after* the initial render (the captured DPR is frozen
+      // at fetch time — bigger headroom = more tolerable post-render
+      // zoom before pixelation appears). The server clamp keeps the
+      // bitmap memory bounded.
+      devicePixelRatio: (typeof window !== "undefined" ? window.devicePixelRatio : 1) * 2
+    };
+    if (this.termType === DATermTypes.JUNCTION) {
+      render.xField = "delta_psi";
+      render.significanceThresholds.foldChangeCutoff = this.settings.deltaPsiCutoff;
+    }
+    return render;
+  }
+  //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()
+    };
+  }
+  /** The "others" group arrives as a placeholder with {in: false}. Replace it with the group that
+   * it stands for: the filter of the app without the samples of the other group, which the server
+   * resolves. */
+  setOthersGroup(samplelst) {
+    const i = samplelst.groups.findIndex((g) => !g.in);
+    if (i == -1) return;
+    const samplesGroup = samplelst.groups.find((g) => g.in);
+    samplelst.groups[i] = getOthersGroup(
+      { name: samplelst.groups[i].name, values: samplesGroup.values },
+      this.app.getState().termfilter
+    );
+  }
+  /** A group that the server resolved from a filter lists no sample here. Keep the number of its
+   * samples that the analysis ran on, where the group does not carry one, for getSampleNum(). */
+  setGroupSizes(response) {
+    for (const [i, g] of this.config.samplelst.groups.entries()) {
+      const n = response[`sample_size${i + 1}`];
+      if (!Array.isArray(g.values) && g.sampleCount === void 0 && Number.isFinite(n)) g.sampleCount = n;
+    }
+  }
+};
+
+export {
+  getGroupColors,
+  VolcanoModel
+};
+//# sourceMappingURL=chunk-JDRUY7IR.js.map

package/dist/chunk-JUPLOOCO.js

@@ -0,0 +1,141 @@
+import {
+  skipPrevActionAbort,
+  storeInit
+} from "./chunk-ZV7LZY67.js";
+import {
+  recoverInit
+} from "./chunk-4E2F26CC.js";
+import {
+  AppBase,
+  sayerror,
+  vocabInit
+} from "./chunk-3YQYWX4M.js";
+import {
+  importPlot
+} from "./chunk-R3RD3VEX.js";
+import {
+  Menu
+} from "./chunk-CNFTQSNE.js";
+import {
+  AppApi
+} from "./chunk-HBNB5TRH.js";
+
+// plots/plot.app.ts
+var PlotApp = class _PlotApp extends AppBase {
+  constructor(opts, api) {
+    super(opts);
+    this.components = {};
+    this.wasDestroyed = false;
+    this.api = api;
+    this.type = _PlotApp.type;
+    this.dom = this.getDom(opts);
+  }
+  static {
+    this.type = "app";
+  }
+  getDom(opts) {
+    const dom = {
+      holder: opts.holder,
+      errdiv: opts.holder.append("div"),
+      plotDiv: opts.holder.append("div")
+    };
+    const controls = opts.violin?.mode == "minimal" ? null : opts.holder.append("div").style("white-space", "nowrap");
+    if (controls) {
+      dom.plotControls = controls.append("div").style("display", "inline-block");
+      dom.recoverControls = controls.append("div").style("display", "inline-block");
+    }
+    return dom;
+  }
+  async preApiFreeze(api) {
+    try {
+      api.tip = new Menu({ padding: "5px" });
+      api.printError = (e) => this.printError(e);
+      const vocab = this.opts.state.vocab;
+      api.vocabApi = this.opts.vocabApi ? this.opts.vocabApi : await vocabInit({
+        app: api,
+        state: {
+          vocab: {
+            // either (genome + dslabel) XOR (terms) can be undefined
+            genome: vocab?.genome || this.opts.state.genome,
+            dslabel: vocab?.dslabel || this.opts.state.dslabel,
+            terms: vocab?.terms
+          }
+        },
+        fetchOpts: this.opts.fetchOpts
+      });
+      this.opts.state.vocab = api.vocabApi.vocab;
+    } catch (e) {
+      console.log(`preApiFreeze error`, e);
+      throw e;
+    }
+  }
+  async init() {
+    try {
+      this.opts.state.nav = { header_mode: "hidden" };
+      this.store = await storeInit({ app: this.api, state: this.opts.state });
+      this.state = await this.store.copyState();
+      this.components = {
+        plots: {}
+      };
+      if (this.opts.app?.features?.includes("recover"))
+        this.components.recover = await recoverInit({
+          app: this.api,
+          holder: this.dom.recoverControls,
+          // TODO: ???? may limit the tracked state to only the filter, activeCohort ???
+          getState: (appState) => appState,
+          //reactsTo: action => true, //action.type != 'plot_edit' || action.type == 'app_refresh',
+          maxHistoryLen: 10
+        });
+      if (this.opts.app?.doNotAwaitInitRender) {
+        this.api.dispatch();
+      } else {
+        await this.api.dispatch();
+      }
+    } catch (e) {
+      this.printError(e);
+      throw e;
+    }
+  }
+  async main() {
+    this.api.vocabApi.main();
+    for (const id in this.components.plots) {
+      const plot = this.components.plots[id];
+      if (!this.state.plots.find((p) => p.id === plot.id)) {
+        plot.destroy();
+        delete this.components.plots[id];
+      }
+    }
+    for (const plot of this.state.plots.values()) {
+      if (plot.parentId) continue;
+      if (!this.components.plots[plot.id]) {
+        const holder = this.opts?.app?.getPlotHolder ? this.opts.app.getPlotHolder(plot, this.dom.holder) : this.dom.holder.append("div");
+        if (!this.dom.plotDiv) this.dom.plotDiv = holder;
+        const { componentInit } = await importPlot(plot.chartType);
+        const plotApi = await componentInit({
+          id: plot.id,
+          app: this.api,
+          holder,
+          controls: this.dom.plotControls
+        });
+        this.components.plots[plot.id] = plotApi;
+      }
+    }
+  }
+  printError(e) {
+    sayerror(this.dom.errdiv, "Error: " + (e.message || e));
+    if (e.stack) console.log(e.stack);
+    this.bus.emit("error");
+  }
+  skipPrevActionAbort(action) {
+    return skipPrevActionAbort(action);
+  }
+  destroy() {
+    if (this.dom?.holder) this.dom.holder.selectAll("*").remove();
+  }
+};
+var appInit = AppApi.getInitFxn(PlotApp);
+
+export {
+  appInit
+};
+//# sourceMappingURL=chunk-JUPLOOCO.js.map

package/dist/chunk-JXVRNCAZ.js

@@ -0,0 +1,96 @@
+import {
+  keyupEnter,
+  require_debounce
+} from "./chunk-3YQYWX4M.js";
+import {
+  Menu
+} from "./chunk-CNFTQSNE.js";
+import {
+  dofetch,
+  dofetch3
+} from "./chunk-KBTFT76L.js";
+import {
+  __toESM
+} from "./chunk-HS5PO5ZQ.js";
+
+// src/gene.js
+var import_debounce = __toESM(require_debounce(), 1);
+var tip;
+function gene_searchbox(p) {
+  if (!tip) {
+    tip = new Menu({ padding: "" });
+    tip.d.style("z-index", 1e3);
+  }
+  const input = p.div.append("input").attr("placeholder", "Search gene").style("width", p.width || "100px");
+  const printdiv = p.resultdiv || (p.tip ? p.tip.d : tip.d);
+  function fold() {
+    if (p.resultdiv) {
+      p.resultdiv.selectAll("*").remove();
+    } else if (p.tip) {
+      p.tip.hide();
+    } else {
+      tip.hide();
+    }
+  }
+  input.on("keyup", (event) => {
+    const str = event.target.value;
+    if (str.length <= 1) {
+      fold();
+      return;
+    }
+    if (keyupEnter(event)) {
+      const hitgene = printdiv.select(".sja_menuoption");
+      if (hitgene.size() > 0) {
+        p.callback(hitgene.text());
+        fold();
+      }
+      return;
+    }
+    debouncer();
+  });
+  input.node().focus();
+  function genesearch() {
+    dofetch("genelookup", { genome: p.genome, input: input.property("value") }).then((data) => {
+      if (data.error) throw data.error;
+      if (!data.hits) throw ".hits[] missing";
+      if (p.resultdiv) {
+        p.resultdiv.selectAll("*").remove();
+      } else if (p.tip) {
+        p.tip.clear().showunder(input.node());
+      } else {
+        tip.clear().showunder(input.node());
+      }
+      for (const name of data.hits) {
+        printdiv.append("div").attr("class", "sja_menuoption").text(name).on("click", () => {
+          p.callback(name);
+          fold();
+        });
+      }
+    }).catch((err) => {
+      printdiv.append("div").text(err.message || err);
+      if (err.stack) console.log(err.stack);
+    });
+  }
+  const debouncer = (0, import_debounce.debounce)(genesearch, 300);
+}
+function findgenemodel_bysymbol(genome, str) {
+  return dofetch3("genelookup", {
+    body: {
+      deep: 1,
+      input: str,
+      genome
+    }
+  }).then((data) => {
+    if (data.error) throw data.error;
+    if (!data.gmlst || data.gmlst.length == 0) return null;
+    return data.gmlst;
+  }).catch((e) => {
+    throw e;
+  });
+}
+
+export {
+  gene_searchbox,
+  findgenemodel_bysymbol
+};
+//# sourceMappingURL=chunk-JXVRNCAZ.js.map