@sjcrh/proteinpaint-client 2.220.0 → 2.221.0

package/dist/chunk-AENUJ7UO.js

@@ -1,299 +0,0 @@
-import {
-  first_genetrack_tolist,
-  gmmode,
-  sayerror
-} from "./chunk-PEKCVKJK.js";
-import {
-  dofetch3
-} from "./chunk-4TJO2DGU.js";
-import {
-  codon_stop,
-  nt2aa,
-  proteinDomainColorScale
-} from "./chunk-WJLRU5OF.js";
-import {
-  select_default
-} from "./chunk-I6Y4O3RR.js";
-
-// common/snp.js
-async function string2snp(genome, str) {
-  const data = await dofetch3("snp", {
-    method: "POST",
-    body: JSON.stringify({ byName: true, genome: genome.name, lst: [str] })
-  });
-  if (data.error) throw data.error;
-  if (!data.results || data.results.length == 0) throw str + ": not a SNP";
-  for (const i of data.results) {
-    const chr = genome.chrlookup[i.chrom.toUpperCase()];
-    if (chr && chr.major) {
-      return {
-        chr: i.chrom,
-        start: i.chromStart,
-        stop: i.chromEnd
-      };
-    }
-  }
-  const r = data.results[0];
-  return {
-    chr: r.chrom,
-    start: r.chromStart,
-    stop: r.chromEnd
-  };
-}
-
-// src/block.init.js
-async function block_init_default(arg) {
-  if (!arg.holder) throw "No holder for block.init";
-  if (!arg.genome) throw "no genome";
-  if (arg.holder instanceof Element) arg.holder = select_default(arg.holder);
-  if (!arg.tklst) arg.tklst = [];
-  if (arg.query) {
-    await step1_findgm(arg);
-    return;
-  }
-  if (arg.model && arg.allmodels) {
-    await step2_getseq(arg);
-    return;
-  }
-}
-async function step1_findgm(arg) {
-  const wait = arg.holder.append("p").style("font-size", "2em").style("color", "#858585").text("Searching for " + arg.query + " ...");
-  const data = await dofetch3("genelookup", {
-    body: { deep: 1, input: arg.query, genome: arg.genome.name }
-  });
-  if (!data) throw "querying genes: server error";
-  if (data.error) throw "error querying genes: " + data.error;
-  if (!data.gmlst || data.gmlst.length == 0) {
-    if (arg.genome.hasSNP) {
-      try {
-        const r = await string2snp(arg.genome, arg.query);
-        wait.remove();
-        const par = {
-          genome: arg.genome,
-          holder: arg.holder,
-          chr: r.chr,
-          start: Math.max(0, r.start - 300),
-          stop: r.start + 300,
-          nobox: true,
-          tklst: arg.tklst,
-          debugmode: arg.debugmode
-        };
-        first_genetrack_tolist(arg.genome, par.tklst);
-        const b = await import("./block-VMVJH77F.js");
-        const block = new b.Block(par);
-        block.addhlregion(r.chr, r.start, r.stop - 1);
-      } catch (e) {
-        wait.text("Not a gene or SNP: " + arg.query);
-      }
-    } else {
-      wait.text("No match to gene: " + arg.query);
-    }
-    return;
-  }
-  wait.remove();
-  arg.allmodels = data.gmlst;
-  for (const m of arg.allmodels) {
-    if (m.isoform.toUpperCase() == (data.found_isoform ? data.found_isoform.toUpperCase() : arg.query.toUpperCase())) {
-      arg.model = m;
-      await step2_getseq(arg);
-      return;
-    }
-  }
-  const defaultisoforms = [];
-  for (const m of arg.allmodels) {
-    if (!m.isoform) throw "isoform missing from one gene model: " + JSON.stringify(m);
-    const n = m.isoform.toUpperCase();
-    if (arg.genome.isoformcache.has(n)) {
-      let nothas = true;
-      for (const m2 of arg.genome.isoformcache.get(n)) {
-        if (m2.chr == m.chr && m2.start == m.start && m2.stop == m.stop && m2.strand == m.strand) {
-          nothas = false;
-          break;
-        }
-      }
-      if (nothas) {
-        arg.genome.isoformcache.get(n).push(m);
-      }
-    } else {
-      arg.genome.isoformcache.set(n, [m]);
-    }
-    if (m.isoform.toUpperCase() == arg.query.toUpperCase()) {
-      defaultisoforms.push(m);
-      break;
-    }
-    if (m.isdefault) {
-      defaultisoforms.push(m);
-    }
-  }
-  if (defaultisoforms.length == 1) {
-    arg.model = defaultisoforms[0];
-  } else if (defaultisoforms.length > 1) {
-    for (const m of defaultisoforms) {
-      if (m.chr == "chrY") {
-        continue;
-      }
-      const chr = arg.genome.chrlookup[m.chr.toUpperCase()];
-      if (!chr) {
-        continue;
-      }
-      if (!chr.major) {
-        continue;
-      }
-      arg.model = m;
-      break;
-    }
-    if (!arg.model) {
-      arg.model = defaultisoforms[0];
-    }
-  }
-  if (!arg.model) {
-    arg.model = arg.allmodels[0];
-  }
-  await step2_getseq(arg);
-}
-async function step2_getseq(arg) {
-  if (arg.model.genomicseq) {
-    checker();
-    step2_getpdomain(arg);
-    return;
-  }
-  const par = {
-    genome: arg.genome.name,
-    coord: arg.model.chr + ":" + (arg.model.start + 1) + "-" + arg.model.stop
-  };
-  const data = await dofetch3("ntseq", { method: "POST", body: JSON.stringify(par) });
-  if (!data) throw "getting sequence: server error";
-  if (data.error) throw "getting sequence: " + data.error;
-  if (!data.seq) throw "no nt seq???";
-  arg.model.genomicseq = data.seq.toUpperCase();
-  arg.model.aaseq = nt2aa(arg.model);
-  checker();
-  await step2_getpdomain(arg);
-  function checker() {
-    if (arg.model.aaseq) {
-      const stop = arg.model.aaseq.indexOf(codon_stop);
-      const cdslen = arg.model.cdslen - (arg.model.startCodonFrame ? 3 - arg.model.startCodonFrame : 0);
-      if (stop != -1 && stop < cdslen / 3 - 1) {
-        sayerror(arg.holder, "Translating " + arg.model.isoform + " ends at " + stop + " AA, expecting " + cdslen / 3);
-      }
-    }
-  }
-}
-async function step2_getpdomain(arg) {
-  const isoform2gm = /* @__PURE__ */ new Map();
-  for (const m of arg.allmodels) {
-    if (!m.pdomains) {
-      m.pdomains = [];
-      m.domain_hidden = {};
-      if (!isoform2gm.has(m.isoform)) isoform2gm.set(m.isoform, []);
-      isoform2gm.get(m.isoform).push(m);
-    }
-  }
-  if (isoform2gm.size == 0) {
-    await step3(arg);
-    return;
-  }
-  const data = await dofetch3("pdomain", {
-    method: "POST",
-    body: JSON.stringify({ genome: arg.genome.name, isoforms: [...isoform2gm.keys()] })
-  });
-  if (data.error) throw "error getting protein domain: " + data.error;
-  if (!Array.isArray(data.lst)) throw ".lst[] not array";
-  for (const a of data.lst) {
-    for (const m of isoform2gm.get(a.name)) {
-      m.pdomains = a.pdomains;
-      if (arg.hidePdomain) {
-        for (const i of a.pdomains) {
-          m.domain_hidden[i.name + i.description] = 1;
-        }
-      }
-    }
-  }
-  if (arg.geneDomains) {
-    if (typeof arg.geneDomains != "object") throw "geneDomains not object";
-    for (const isoform in arg.geneDomains) {
-      const lst = isoform2gm.get(isoform);
-      if (!lst) throw `unknown isoform ${isoform} from geneDomains{}`;
-      for (const g of lst) {
-        if (!g.pdomains) g.pdomains = [];
-        if (!Array.isArray(arg.geneDomains[isoform])) throw `geneDomains[${isoform}] not array`;
-        for (const b of arg.geneDomains[isoform]) {
-          if (typeof b != "object") throw "element from geneDomains[] not object";
-          if (!Number.isInteger(b.start)) throw "start not integer from geneDomains[]";
-          if (!Number.isInteger(b.stop)) throw "stop not integer from geneDomains[]";
-          if (b.start > b.stop) throw "start>stop from geneDomains[]";
-          if (!b.name) b.name = "Custom domain";
-          if (!g.pdomains.find((a) => a.start == b.start && a.stop == b.stop && a.name == b.name)) g.pdomains.push(b);
-        }
-      }
-    }
-  }
-  const s = proteinDomainColorScale();
-  for (const lst of isoform2gm.values()) {
-    for (const g of lst) {
-      for (const d of g.pdomains || []) {
-        if (!d.color) d.color = s(d.name + d.description);
-      }
-    }
-  }
-  await step3(arg);
-}
-async function step3(arg) {
-  let mode = arg.gmmode;
-  if (!mode) {
-    if (arg.model.cdslen) {
-      mode = gmmode.protein;
-    } else {
-      mode = gmmode.exononly;
-    }
-  }
-  if (arg.dataset) {
-    if (!Array.isArray(arg.dataset)) throw "dataset is not array";
-    for (const dsname of arg.dataset) {
-      if (arg.genome.datasets[dsname] && !arg.genome.datasets[dsname].legacyDsIsUninitiated) continue;
-      const d = await dofetch3(`getDataset?genome=${arg.genome.name}&dsname=${dsname}`);
-      if (d.error) throw `invalid name from dataset[]: ${d.error}`;
-      if (!d.ds) throw ".ds missing";
-      const ds = arg.genome.datasets[d.ds.label];
-      Object.assign(ds, d.ds);
-      const _ = await import("./legacyDataset-IEFWFVS6.js");
-      _.validate_oldds(ds);
-      delete ds.legacyDsIsUninitiated;
-    }
-  }
-  const b = await import("./block-VMVJH77F.js");
-  arg.__blockInstance = new b.Block({
-    genome: arg.genome,
-    holder: arg.holder,
-    nobox: true,
-    usegm: arg.model,
-    gmstackheight: 37,
-    allgm: arg.allmodels,
-    datasetlst: arg.dataset,
-    legacyDsFilter: arg.legacyDsFilter,
-    mset: arg.mset,
-    hlaachange: arg.hlaachange,
-    hlvariants: arg.hlvariants,
-    hlregions: arg.hlregions,
-    aarange: arg.aarange,
-    gmmode: mode,
-    hidedatasetexpression: arg.hidedatasetexpression,
-    hidegenecontrol: arg.hidegenecontrol,
-    hidegenelegend: arg.hidegenelegend,
-    variantPageCall_snv: arg.variantPageCall_snv,
-    datasetqueries: arg.datasetqueries,
-    samplecart: arg.samplecart,
-    debugmode: arg.debugmode,
-    tklst: arg.tklst,
-    mclassOverride: arg.mclassOverride,
-    hide_dsHandles: arg.hide_dsHandles,
-    onloadalltk_always: arg.onloadalltk_always,
-    onAddRemoveTk: arg.onAddRemoveTk
-  });
-}
-
-export {
-  string2snp,
-  block_init_default
-};
-//# sourceMappingURL=chunk-AENUJ7UO.js.map

package/dist/chunk-AZSWQWKL.js

@@ -1,263 +0,0 @@
-import {
-  hierCluster_renderers_exports
-} from "./chunk-GBKQHVKJ.js";
-import {
-  hierCluster_interactivity_exports
-} from "./chunk-4GKC43U3.js";
-import {
-  Matrix
-} from "./chunk-HN67TOA7.js";
-import {
-  filterJoin,
-  getNormalRoot
-} from "./chunk-PEKCVKJK.js";
-import {
-  clusterMethodLst,
-  distanceMethodLst,
-  dofetch3
-} from "./chunk-4TJO2DGU.js";
-import {
-  TermTypes2Dt,
-  dictionaryNumericTypes
-} from "./chunk-JHLF2ZC2.js";
-import {
-  colorScaleMap
-} from "./chunk-WJLRU5OF.js";
-import {
-  deepEqual,
-  getCompInit
-} from "./chunk-HBNB5TRH.js";
-import {
-  extent,
-  linear
-} from "./chunk-4OLM3KSB.js";
-
-// plots/matrix/hierCluster.js
-var HierCluster = class _HierCluster extends Matrix {
-  static type = "hierCluster";
-  constructor(opts) {
-    super(opts);
-    this.type = _HierCluster.type;
-    this.chartType = _HierCluster.type;
-  }
-  async init(appState) {
-    await super.init(appState);
-    this.maySetSandboxHeader(appState);
-    this.hcClipId = this.seriesClipId + "-hc";
-    this.dom.hcClipRect = this.dom.svg.select("defs").append("clipPath").attr("id", this.hcClipId).attr("clipPathUnits", "userSpaceOnUse").append("rect").attr("display", "block");
-    this.dom.topDendrogram = this.dom.svg.insert("g", "g").attr("clip-path", `url(#${this.hcClipId})`).append("g").attr("class", "sjpp-matrix-dendrogram").attr("data-testid", "hierCluster_top_dendrogram").on("click", (event) => {
-      const clickedClusterId = this.getClusterFromTopDendrogram(event);
-      if (clickedClusterId) {
-        this.clickedClusterIds = this.getAllChildrenClusterIds(clickedClusterId);
-        this.clickedClusterIds.push(clickedClusterId);
-        const clickedCluster = this.hierClusterData.clustering.col.mergedClusters.get(clickedClusterId);
-        const clickedClusterSampleNames = clickedCluster.children.map((c) => c.name);
-        this.addSelectedSamplesOptions(clickedClusterSampleNames, event);
-      } else {
-        delete this.clickedClusterIds;
-      }
-      if (this.clickedLeftClusterIds) {
-        delete this.clickedLeftClusterIds;
-        this.plotDendrogramHclust();
-      } else this.plotDendrogramHclust("top");
-    });
-    this.dom.leftDendrogram = this.dom.svg.insert("g", "g").attr("class", "sjpp-matrix-dendrogram").attr("data-testid", "hierCluster_left_dendrogram").on("click", (event) => {
-      const clickedLeftClusterId = this.getClusterFromLeftDendrogram(event);
-      if (clickedLeftClusterId) {
-        this.clickedLeftClusterIds = this.getAllChildrenClusterIds(clickedLeftClusterId, true);
-        this.clickedLeftClusterIds.push(clickedLeftClusterId);
-        const clickedLeftCluster = this.hierClusterData.clustering.row.mergedClusters.get(clickedLeftClusterId);
-        const clickedLeftClusterRowsNames = clickedLeftCluster.children.map((c) => c.name);
-        this.addSelectedRowsOptions(clickedLeftClusterRowsNames, event);
-      } else {
-        delete this.clickedLeftClusterIds;
-      }
-      if (this.clickedClusterIds) {
-        delete this.clickedClusterIds;
-        this.plotDendrogramHclust();
-      } else this.plotDendrogramHclust("left");
-    });
-  }
-  async setHierClusterData(_data = {}) {
-    this.prevServerData = this.currServerData;
-    const [d, twlst] = await this.requestData({});
-    if (d.error) throw d.error;
-    this.currServerData = structuredClone(d);
-    if (!deepEqual(this.prevServerData, this.currServerData)) {
-      delete this.clickedClusterIds;
-      delete this.clickedLeftClusterIds;
-    }
-    const s = this.settings.hierCluster;
-    if (!d.clustering) {
-      if (d.gene) {
-        throw `Cannot do clustering: data is only available for 1 gene (${d.gene}). Try again by adding more genes.`;
-      }
-    }
-    this.hierClusterData = d;
-    const c = this.hierClusterData.clustering;
-    this.setHierColorScale(c);
-    const samples = {};
-    for (const [i, column] of c.col.order.entries()) {
-      samples[column.name] = { sample: column.name };
-      for (const [j, row] of c.row.order.entries()) {
-        const tw = twlst.find((tw2) => tw2.$id === row.name || tw2.id === row.name);
-        const value = c.matrix[j][i];
-        samples[column.name][tw.$id] = {
-          key: tw.term.name,
-          values: [
-            {
-              sample: column.name,
-              dt: TermTypes2Dt[this.state.config.dataType],
-              label: s.termGroupName,
-              // gene: tw.term.name,
-              // chr: tw.term.chr,
-              // pos: `${tw.term.start}-${tw.term.stop}`,
-              value
-              // the color will be computed in matrix.cells, so that
-              // it can get updated even when there are no nonsetting state diff
-            }
-          ]
-        };
-      }
-    }
-    this.hcTermNameOrder = this.settings.hierCluster.sortClusterRows == "asListed" ? twlst.map((t) => t.term.name) : this.settings.hierCluster.sortClusterRows == "byName" ? twlst.map((t) => t.term.name).sort() : dictionaryNumericTypes.has(this.config.dataType) ? c.row.order.map((row) => twlst.find((t) => t.$id == row.name || t.term.id == row.name)?.term.name) : c.row.order.map((row) => twlst.find((t) => t.$id == row.name)?.term.name);
-    if (this.hcTermNameOrder.includes(void 0)) throw `unable to map row.name to term.name`;
-    this.hcTermSorter = (a, b) => {
-      const i = this.hcTermNameOrder.indexOf(a.tw.term.name);
-      const j = this.hcTermNameOrder.indexOf(b.tw.term.name);
-      if (i == -1 && j == -1) return 0;
-      if (i == -1) return 1;
-      if (j == -1) return -1;
-      return i - j;
-    };
-    this.hcSampleNameOrder = c.col.order.map((col) => col.name);
-    this.hcSampleSorter = (a, b) => {
-      const i = this.hcSampleNameOrder.indexOf(a.sample);
-      const j = this.hcSampleNameOrder.indexOf(b.sample);
-      if (i == -1 && j == -1) return 0;
-      if (i == -1) return 1;
-      if (j == -1) return -1;
-      return i - j;
-    };
-    const byTermId = {};
-    for (const tw of twlst) {
-      if (d.byTermId?.[tw.term.name]) byTermId[tw.$id] = d.byTermId[tw.term.name];
-    }
-    this.hierClusterSamples = {
-      refs: { byTermId, bySampleId: d.bySampleId },
-      lst: c.col.order.map((c2) => samples[c2.name]),
-      samples,
-      removedHierClusterTerms: d.removedHierClusterTerms
-    };
-  }
-  async requestData() {
-    const body = this.currRequestOpts?.hierCluster || this.getHCRequestBody(this.state);
-    const twlst = this.hcTermGroup.lst;
-    const data = await dofetch3("termdb/cluster", { body, signal: this.api.getAbortSignal?.() });
-    return [data, twlst];
-  }
-  getHCRequestBody(state) {
-    this.hcTermGroup = this.config.termgroups.find((grp) => grp.type == "hierCluster") || this.termOrder?.find((t) => t.grp.type == "hierCluster")?.grp;
-    const s = state.config.settings.hierCluster;
-    const dictionaryLegendFilter = {
-      type: "tvslst",
-      in: true,
-      join: "and",
-      lst: state.config.legendValueFilter.lst.filter((f) => !f.tvs.legendFilterType)
-    };
-    const terms = this.getClusterRowTermsAsParameter();
-    if (!terms.length) throw "no data";
-    if (!clusterMethodLst.find((i) => i.value == s.clusterMethod)) throw "Invalid cluster method";
-    if (!distanceMethodLst.find((i) => i.value == s.distanceMethod)) throw "Invalid distance method";
-    const body = {
-      genome: state.vocab.genome,
-      dslabel: state.vocab.dslabel,
-      dataType: state.config.dataType,
-      clusterMethod: s.clusterMethod,
-      distanceMethod: s.distanceMethod,
-      zScoreTransformation: s.zScoreTransformation,
-      terms,
-      filter: getNormalRoot(filterJoin([state.filter, dictionaryLegendFilter])),
-      filter0: state.filter0
-    };
-    if (state.config.dataType == "proteomeAbundance") {
-      body.proteomeDetails = {
-        organism: state.config.proteomeDetails?.organism,
-        assay: state.config.proteomeDetails?.assay,
-        cohort: state.config.proteomeDetails?.cohort
-      };
-    }
-    return body;
-  }
-  combineData() {
-    if (!this.hierClusterSamples) return;
-    const d = this.data;
-    const removedHierClusterTerms = this.hierClusterSamples.removedHierClusterTerms;
-    const samples = {};
-    const lst = [];
-    for (const sampleId in this.hierClusterSamples.samples) {
-      const s = this.hierClusterSamples.samples[sampleId];
-      samples[sampleId] = s;
-      lst.push(s);
-      if (sampleId in d.samples) Object.assign(s, d.samples[sampleId]);
-      const _ref_ = this.hierClusterSamples.refs.bySampleId[sampleId] || {};
-      if (!s._ref_) s._ref_ = _ref_;
-      else Object.assign(s._ref_, _ref_);
-    }
-    const t = this.hierClusterSamples.refs.byTermId;
-    for (const $id of Object.keys(t)) {
-      d.refs.byTermId[$id] = Object.assign({}, d.refs.byTermId[$id] || {}, t[$id]);
-    }
-    this.data = { samples, lst, refs: d.refs, removedHierClusterTerms };
-  }
-  setHierColorScale(c) {
-    const hc = this.settings.hierCluster;
-    const scale = linear(colorScaleMap[hc.colorScale].domain, colorScaleMap[hc.colorScale].range).clamp(true);
-    const globalMinMaxes = [];
-    for (const row of c.matrix) {
-      globalMinMaxes.push(...extent(row));
-    }
-    const absMax = Math.min(hc.zScoreCap, Math.max(...extent(globalMinMaxes).map(Math.abs)));
-    const [min, max] = hc.zScoreTransformation ? [-absMax, absMax] : [Math.min(...globalMinMaxes), Math.max(...globalMinMaxes)];
-    this.hierClusterValues = { scale, min, max };
-  }
-  getValueColor(value) {
-    const hc = this.settings.hierCluster;
-    if (hc.zScoreTransformation) {
-      const zScoreCap = this.settings.hierCluster.zScoreCap;
-      return this.hierClusterValues.scale((value - -zScoreCap) / (zScoreCap * 2));
-    } else {
-      return this.hierClusterValues.scale(value / this.hierClusterValues.max);
-    }
-  }
-  /* returns list of gene terms as request parameter, e.g. {gene,chr,start,stop}
-  	request parameter only need term but not tw, as it will simply fetch continuous sample values on terms without transform
-  
-  	use of this function is unfortunate because:
-  		the incomplete migration of {name} to {gene} for gene-based term
-  		geneset edit ui is hardcoded to return {name}
-  		existing plot states contain {name}
-  
-  	!!! migration instruction !!!
-  	- term.name is for display only, if a term is gene-based, it has term.gene=str
-  	- a geneVariant term can be based on a genomic range (and not a gene), in that case it won't have term.gene and cannot be used where gene is expected, e.g. gene-based clustering analysis
-  
-  	*/
-  getClusterRowTermsAsParameter() {
-    const lst = this.hcTermGroup.lst.map(this.opts.app.vocabApi.getTwMinCopy);
-    lst.sort((a, b) => a.term.name < b.term.name ? -1 : 1);
-    return lst;
-  }
-};
-for (const methods of [hierCluster_renderers_exports, hierCluster_interactivity_exports]) {
-  for (const methodName in methods) HierCluster.prototype[methodName] = methods[methodName];
-}
-var hierClusterInit = getCompInit(HierCluster);
-var componentInit = hierClusterInit;
-
-export {
-  HierCluster,
-  hierClusterInit,
-  componentInit
-};
-//# sourceMappingURL=chunk-AZSWQWKL.js.map