@sjcrh/proteinpaint-client 2.166.0 → 2.167.0

package/dist/chunk-5GQEXGQC.js

@@ -0,0 +1,256 @@
+import {
+  Matrix
+} from "./chunk-TRS2KOXE.js";
+import {
+  hierCluster_renderers_exports
+} from "./chunk-D53HBTYQ.js";
+import {
+  hierCluster_interactivity_exports
+} from "./chunk-WVZV43PF.js";
+import {
+  filterJoin,
+  getNormalRoot
+} from "./chunk-PX7A6LBY.js";
+import {
+  clusterMethodLst,
+  distanceMethodLst,
+  dofetch3
+} from "./chunk-6ADTREVX.js";
+import {
+  NUMERIC_DICTIONARY_TERM,
+  TermTypes2Dt
+} from "./chunk-DKI7YOTJ.js";
+import {
+  colorScaleMap
+} from "./chunk-6VQ4C735.js";
+import {
+  deepEqual,
+  getCompInit
+} from "./chunk-DDOM4XYV.js";
+import {
+  extent,
+  linear
+} from "./chunk-W2IWHXLL.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() : this.config.dataType === NUMERIC_DICTIONARY_TERM ? 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
+    };
+    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-5GQEXGQC.js.map

package/dist/chunk-5GQEXGQC.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../plots/matrix/hierCluster.js"],
+  "sourcesContent": ["import { Matrix } from './matrix'\nimport { getCompInit, deepEqual } from '#rx'\nimport * as renderers from './hierCluster.renderers'\nimport * as interactivity from './hierCluster.interactivity'\nimport { dofetch3 } from '#common/dofetch'\nimport { extent } from 'd3-array'\nimport { scaleLinear } from 'd3-scale'\nimport { filterJoin, getNormalRoot } from '#filter'\nimport { clusterMethodLst, distanceMethodLst } from '#shared/clustering.js'\nimport { TermTypes, TermTypes2Dt, NUMERIC_DICTIONARY_TERM } from '#shared/terms.js'\nimport { colorScaleMap } from '#shared/common.js'\n\nexport class HierCluster extends Matrix {\n\tstatic type = 'hierCluster'\n\n\tconstructor(opts) {\n\t\tsuper(opts)\n\t\tthis.type = HierCluster.type\n\t\tthis.chartType = HierCluster.type\n\t}\n\n\tasync init(appState) {\n\t\tawait super.init(appState)\n\n\t\tthis.maySetSandboxHeader(appState)\n\n\t\tthis.hcClipId = this.seriesClipId + '-hc'\n\t\tthis.dom.hcClipRect = this.dom.svg\n\t\t\t.select('defs')\n\t\t\t.append('clipPath')\n\t\t\t.attr('id', this.hcClipId)\n\t\t\t//.attr('clipPathUnits', 'objectBoundingBox')\n\t\t\t.attr('clipPathUnits', 'userSpaceOnUse')\n\t\t\t.append('rect')\n\t\t\t.attr('display', 'block')\n\n\t\tthis.dom.topDendrogram = this.dom.svg\n\t\t\t.insert('g', 'g')\n\t\t\t.attr('clip-path', `url(#${this.hcClipId})`)\n\t\t\t.append('g')\n\t\t\t.attr('class', 'sjpp-matrix-dendrogram')\n\t\t\t.attr('data-testid', 'hierCluster_top_dendrogram')\n\t\t\t.on('click', event => {\n\t\t\t\tconst clickedClusterId = this.getClusterFromTopDendrogram(event)\n\t\t\t\tif (clickedClusterId) {\n\t\t\t\t\tthis.clickedClusterIds = this.getAllChildrenClusterIds(clickedClusterId)\n\t\t\t\t\tthis.clickedClusterIds.push(clickedClusterId)\n\n\t\t\t\t\tconst clickedCluster = this.hierClusterData.clustering.col.mergedClusters.get(clickedClusterId)\n\t\t\t\t\tconst clickedClusterSampleNames = clickedCluster.children.map(c => c.name)\n\t\t\t\t\tthis.addSelectedSamplesOptions(clickedClusterSampleNames, event)\n\t\t\t\t} else {\n\t\t\t\t\t// if not clicking on a cluster, change highlighted cluster color from red back to black\n\t\t\t\t\tdelete this.clickedClusterIds\n\t\t\t\t}\n\n\t\t\t\t// rerender the row Dendrogram\n\t\t\t\tif (this.clickedLeftClusterIds) {\n\t\t\t\t\t// when left dendrogram has highlight, clicking top dendro should cancel it\n\t\t\t\t\tdelete this.clickedLeftClusterIds\n\t\t\t\t\tthis.plotDendrogramHclust()\n\t\t\t\t} else this.plotDendrogramHclust('top')\n\t\t\t})\n\n\t\tthis.dom.leftDendrogram = this.dom.svg\n\t\t\t.insert('g', 'g')\n\t\t\t.attr('class', 'sjpp-matrix-dendrogram')\n\t\t\t.attr('data-testid', 'hierCluster_left_dendrogram')\n\t\t\t.on('click', event => {\n\t\t\t\tconst clickedLeftClusterId = this.getClusterFromLeftDendrogram(event)\n\t\t\t\tif (clickedLeftClusterId) {\n\t\t\t\t\tthis.clickedLeftClusterIds = this.getAllChildrenClusterIds(clickedLeftClusterId, true)\n\t\t\t\t\tthis.clickedLeftClusterIds.push(clickedLeftClusterId)\n\n\t\t\t\t\tconst clickedLeftCluster = this.hierClusterData.clustering.row.mergedClusters.get(clickedLeftClusterId)\n\n\t\t\t\t\tconst clickedLeftClusterRowsNames = clickedLeftCluster.children.map(c => c.name)\n\n\t\t\t\t\tthis.addSelectedRowsOptions(clickedLeftClusterRowsNames, event)\n\t\t\t\t} else {\n\t\t\t\t\t// if not clicking on a cluster, change highlighted cluster color from red back to black\n\t\t\t\t\tdelete this.clickedLeftClusterIds\n\t\t\t\t}\n\n\t\t\t\t// rerender the row Dendrogram\n\t\t\t\tif (this.clickedClusterIds) {\n\t\t\t\t\t// when top dendrogram has highlight, clicking left dendro should cancel it\n\t\t\t\t\tdelete this.clickedClusterIds\n\t\t\t\t\tthis.plotDendrogramHclust()\n\t\t\t\t} else this.plotDendrogramHclust('left')\n\t\t\t})\n\t\t//.attr('clip-path', `url(#${this.seriesClipId})`)\n\t}\n\n\tasync setHierClusterData(_data = {}) {\n\t\tthis.prevServerData = this.currServerData\n\t\tconst [d, twlst] = await this.requestData({})\n\t\tif (d.error) throw d.error\n\t\tthis.currServerData = structuredClone(d)\n\t\tif (!deepEqual(this.prevServerData, this.currServerData)) {\n\t\t\t// do not persist highlighted top/left dendrogram branch selection\n\t\t\t// when the cohort, clustering method, or other config changes the server data\n\t\t\tdelete this.clickedClusterIds\n\t\t\tdelete this.clickedLeftClusterIds\n\t\t}\n\t\tconst s = this.settings.hierCluster\n\n\t\tif (!d.clustering) {\n\t\t\t// stop-gap data validation, lacks essential data part\n\t\t\tif (d.gene) {\n\t\t\t\t// for now backend returns {gene:str, data:{}} if there's only 1 eligible gene\n\t\t\t\tthrow `Cannot do clustering: data is only available for 1 gene (${d.gene}). Try again by adding more genes.`\n\t\t\t}\n\t\t\t//throw 'Cannot do clustering: invalid server response (lacks .clustering{})'\n\t\t}\n\t\tthis.hierClusterData = d\n\n\t\tconst c = this.hierClusterData.clustering\n\t\tthis.setHierColorScale(c)\n\n\t\tconst samples = {}\n\n\t\t/* see comments inside plotDendrogramHclust() on structure of d.clustering.row{} and col{}\n\t\tassumes c.col is samples and c.row is non-sample things (genes for now); later may flip to c.row be samples instead!!\n\t\t*/\n\t\tfor (const [i, column] of c.col.order.entries()) {\n\t\t\tsamples[column.name] = { sample: column.name }\n\t\t\tfor (const [j, row] of c.row.order.entries()) {\n\t\t\t\tconst tw = twlst.find(tw => tw.$id === row.name || tw.id === row.name)\n\t\t\t\t// if (!tw) {\n\t\t\t\t// \tconsole.warn(`no matching tw for row.name='${row.name}'`)\n\t\t\t\t// \tcontinue\n\t\t\t\t// }\n\t\t\t\tconst value = c.matrix[j][i]\n\t\t\t\tsamples[column.name][tw.$id] = {\n\t\t\t\t\tkey: tw.term.name,\n\t\t\t\t\tvalues: [\n\t\t\t\t\t\t{\n\t\t\t\t\t\t\tsample: column.name,\n\t\t\t\t\t\t\tdt: TermTypes2Dt[this.state.config.dataType],\n\t\t\t\t\t\t\tlabel: s.termGroupName,\n\t\t\t\t\t\t\t// gene: tw.term.name,\n\t\t\t\t\t\t\t// chr: tw.term.chr,\n\t\t\t\t\t\t\t// pos: `${tw.term.start}-${tw.term.stop}`,\n\t\t\t\t\t\t\tvalue\n\t\t\t\t\t\t\t// the color will be computed in matrix.cells, so that\n\t\t\t\t\t\t\t// it can get updated even when there are no nonsetting state diff\n\t\t\t\t\t\t}\n\t\t\t\t\t]\n\t\t\t\t}\n\t\t\t}\n\t\t}\n\n\t\tthis.hcTermNameOrder =\n\t\t\tthis.settings.hierCluster.sortClusterRows == 'asListed'\n\t\t\t\t? twlst.map(t => t.term.name)\n\t\t\t\t: this.settings.hierCluster.sortClusterRows == 'byName'\n\t\t\t\t? twlst.map(t => t.term.name).sort()\n\t\t\t\t: this.config.dataType === NUMERIC_DICTIONARY_TERM\n\t\t\t\t? c.row.order.map(row => twlst.find(t => t.$id == row.name || t.term.id == row.name)?.term.name)\n\t\t\t\t: c.row.order.map(row => twlst.find(t => t.$id == row.name)?.term.name)\n\n\t\tif (this.hcTermNameOrder.includes(undefined)) throw `unable to map row.name to term.name`\n\n\t\tthis.hcTermSorter = (a, b) => {\n\t\t\tconst i = this.hcTermNameOrder.indexOf(a.tw.term.name)\n\t\t\tconst j = this.hcTermNameOrder.indexOf(b.tw.term.name)\n\t\t\tif (i == -1 && j == -1) return 0\n\t\t\tif (i == -1) return 1\n\t\t\tif (j == -1) return -1\n\t\t\treturn i - j\n\t\t}\n\n\t\tthis.hcSampleNameOrder = c.col.order.map(col => col.name)\n\t\tthis.hcSampleSorter = (a, b) => {\n\t\t\tconst i = this.hcSampleNameOrder.indexOf(a.sample)\n\t\t\tconst j = this.hcSampleNameOrder.indexOf(b.sample)\n\t\t\tif (i == -1 && j == -1) return 0\n\t\t\tif (i == -1) return 1\n\t\t\tif (j == -1) return -1\n\t\t\treturn i - j\n\t\t}\n\n\t\t// from d.byTermId to byTermId: change byTermId keys from gene names to $ids\n\t\tconst byTermId = {}\n\t\tfor (const tw of twlst) {\n\t\t\tif (d.byTermId[tw.term.name]) byTermId[tw.$id] = d.byTermId[tw.term.name]\n\t\t}\n\t\tthis.hierClusterSamples = {\n\t\t\trefs: { byTermId, bySampleId: d.bySampleId },\n\t\t\tlst: c.col.order.map(c => samples[c.name]),\n\t\t\tsamples,\n\t\t\tremovedHierClusterTerms: d.removedHierClusterTerms\n\t\t}\n\t}\n\n\tasync requestData() {\n\t\t// may revert to using {signal} argument if detectStale() is used to wrap this function\n\t\tconst body = this.currRequestOpts?.hierCluster || this.getHCRequestBody(this.state)\n\t\tconst twlst = this.hcTermGroup.lst\n\t\tconst data = await dofetch3('termdb/cluster', { body, signal: this.api.getAbortSignal?.() })\n\t\t// return the twlst that was submitted in the data request, this data-to-twlst reference will not be affected by race condition\n\t\treturn [data, twlst]\n\t}\n\n\tgetHCRequestBody(state) {\n\t\tthis.hcTermGroup =\n\t\t\tthis.config.termgroups.find(grp => grp.type == 'hierCluster') ||\n\t\t\tthis.termOrder?.find(t => t.grp.type == 'hierCluster')?.grp\n\n\t\tconst s = state.config.settings.hierCluster\n\t\t// temporary fix to get rid of hard/soft filter and only keep dictionary legend filter,\n\t\t// soft filter shouldn't be used to filter out any samples for hierCluster\n\t\t// TODO: add hard filter back to filter out samples\n\t\tconst dictionaryLegendFilter = {\n\t\t\ttype: 'tvslst',\n\t\t\tin: true,\n\t\t\tjoin: 'and',\n\t\t\tlst: state.config.legendValueFilter.lst.filter(f => !f.tvs.legendFilterType)\n\t\t}\n\t\tconst terms = this.getClusterRowTermsAsParameter()\n\t\tif (!terms.length) throw 'no data'\n\t\t// !!! NOTE !!!\n\t\t// all parameters here must remove payload properties that are\n\t\t// not relevant to the data request, so that the dofetch and/or\n\t\t// browser caching would work\n\n\t\t// Checking if cluster and distance method for hierarchial clustering is valid\n\t\tif (!clusterMethodLst.find(i => i.value == s.clusterMethod)) throw 'Invalid cluster method'\n\t\tif (!distanceMethodLst.find(i => i.value == s.distanceMethod)) throw 'Invalid distance method'\n\t\tconst body = {\n\t\t\tgenome: state.vocab.genome,\n\t\t\tdslabel: state.vocab.dslabel,\n\t\t\tdataType: state.config.dataType,\n\t\t\tclusterMethod: s.clusterMethod,\n\t\t\tdistanceMethod: s.distanceMethod,\n\t\t\tzScoreTransformation: s.zScoreTransformation,\n\t\t\tterms,\n\t\t\tfilter: getNormalRoot(filterJoin([state.filter, dictionaryLegendFilter])),\n\t\t\tfilter0: state.filter0\n\t\t}\n\t\treturn body\n\t}\n\n\tcombineData() {\n\t\tif (!this.hierClusterSamples) return\n\t\tconst d = this.data // matrix data\n\t\tconst removedHierClusterTerms = this.hierClusterSamples.removedHierClusterTerms\n\t\tconst samples = {}\n\t\tconst lst = []\n\t\t// the gene expression samples will be used as a filter for the matrix samples\n\t\tfor (const sampleId in this.hierClusterSamples.samples) {\n\t\t\tconst s = this.hierClusterSamples.samples[sampleId]\n\t\t\tsamples[sampleId] = s\n\t\t\tlst.push(s)\n\t\t\tif (sampleId in d.samples) Object.assign(s, d.samples[sampleId])\n\t\t\tconst _ref_ = this.hierClusterSamples.refs.bySampleId[sampleId] || {}\n\t\t\tif (!s._ref_) s._ref_ = _ref_\n\t\t\t// hierCluster refs.bySampleId will overwrite matrix reference properties with the same name\n\t\t\telse Object.assign(s._ref_, _ref_)\n\t\t}\n\n\t\t// combine this.hierClusterSamples.refs.byTermId into this.data.refs.byTermId\n\t\tconst t = this.hierClusterSamples.refs.byTermId\n\t\tfor (const $id of Object.keys(t)) {\n\t\t\td.refs.byTermId[$id] = Object.assign({}, d.refs.byTermId[$id] || {}, t[$id])\n\t\t}\n\t\tthis.data = { samples, lst, refs: d.refs, removedHierClusterTerms }\n\t}\n\n\tsetHierColorScale(c) {\n\t\tconst hc = this.settings.hierCluster\n\t\tconst scale = scaleLinear(colorScaleMap[hc.colorScale].domain, colorScaleMap[hc.colorScale].range).clamp(true)\n\t\tconst globalMinMaxes = []\n\t\tfor (const row of c.matrix) {\n\t\t\tglobalMinMaxes.push(...extent(row))\n\t\t}\n\t\tconst absMax = Math.min(hc.zScoreCap, Math.max(...extent(globalMinMaxes).map(Math.abs)))\n\t\t// if zScore transformation is not performed, should use min/max value from data\n\t\tconst [min, max] = hc.zScoreTransformation\n\t\t\t? [-absMax, absMax]\n\t\t\t: [Math.min(...globalMinMaxes), Math.max(...globalMinMaxes)]\n\t\t// what's purpose of assigning this.hierClusterValues{}, to signal something to matrix code?\n\t\tthis.hierClusterValues = { scale, min, max }\n\t}\n\n\tgetValueColor(value) {\n\t\tconst hc = this.settings.hierCluster\n\t\tif (hc.zScoreTransformation) {\n\t\t\tconst zScoreCap = this.settings.hierCluster.zScoreCap\n\t\t\treturn this.hierClusterValues.scale((value - -zScoreCap) / (zScoreCap * 2))\n\t\t} else {\n\t\t\treturn this.hierClusterValues.scale(value / this.hierClusterValues.max)\n\t\t}\n\t}\n\n\t/* returns list of gene terms as request parameter, e.g. {gene,chr,start,stop}\n\trequest parameter only need term but not tw, as it will simply fetch continuous sample values on terms without transform\n\n\tuse of this function is unfortunate because:\n\t\tthe incomplete migration of {name} to {gene} for gene-based term\n\t\tgeneset edit ui is hardcoded to return {name}\n\t\texisting plot states contain {name}\n\n\t!!! migration instruction !!!\n\t- term.name is for display only, if a term is gene-based, it has term.gene=str\n\t- 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\n\n\t*/\n\tgetClusterRowTermsAsParameter() {\n\t\tconst lst = this.hcTermGroup.lst.map(this.opts.app.vocabApi.getTwMinCopy)\n\t\t// this helps caching by having a more consistent URL string\n\t\tlst.sort((a, b) => (a.term.name < b.term.name ? -1 : 1))\n\t\treturn lst\n\t}\n}\n\nfor (const methods of [renderers, interactivity]) {\n\tfor (const methodName in methods) HierCluster.prototype[methodName] = methods[methodName]\n}\n\nexport const hierClusterInit = getCompInit(HierCluster)\nexport const componentInit = hierClusterInit\n"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAYO,IAAM,cAAN,MAAM,qBAAoB,OAAO;AAAA,EACvC,OAAO,OAAO;AAAA,EAEd,YAAY,MAAM;AACjB,UAAM,IAAI;AACV,SAAK,OAAO,aAAY;AACxB,SAAK,YAAY,aAAY;AAAA,EAC9B;AAAA,EAEA,MAAM,KAAK,UAAU;AACpB,UAAM,MAAM,KAAK,QAAQ;AAEzB,SAAK,oBAAoB,QAAQ;AAEjC,SAAK,WAAW,KAAK,eAAe;AACpC,SAAK,IAAI,aAAa,KAAK,IAAI,IAC7B,OAAO,MAAM,EACb,OAAO,UAAU,EACjB,KAAK,MAAM,KAAK,QAAQ,EAExB,KAAK,iBAAiB,gBAAgB,EACtC,OAAO,MAAM,EACb,KAAK,WAAW,OAAO;AAEzB,SAAK,IAAI,gBAAgB,KAAK,IAAI,IAChC,OAAO,KAAK,GAAG,EACf,KAAK,aAAa,QAAQ,KAAK,QAAQ,GAAG,EAC1C,OAAO,GAAG,EACV,KAAK,SAAS,wBAAwB,EACtC,KAAK,eAAe,4BAA4B,EAChD,GAAG,SAAS,WAAS;AACrB,YAAM,mBAAmB,KAAK,4BAA4B,KAAK;AAC/D,UAAI,kBAAkB;AACrB,aAAK,oBAAoB,KAAK,yBAAyB,gBAAgB;AACvE,aAAK,kBAAkB,KAAK,gBAAgB;AAE5C,cAAM,iBAAiB,KAAK,gBAAgB,WAAW,IAAI,eAAe,IAAI,gBAAgB;AAC9F,cAAM,4BAA4B,eAAe,SAAS,IAAI,OAAK,EAAE,IAAI;AACzE,aAAK,0BAA0B,2BAA2B,KAAK;AAAA,MAChE,OAAO;AAEN,eAAO,KAAK;AAAA,MACb;AAGA,UAAI,KAAK,uBAAuB;AAE/B,eAAO,KAAK;AACZ,aAAK,qBAAqB;AAAA,MAC3B,MAAO,MAAK,qBAAqB,KAAK;AAAA,IACvC,CAAC;AAEF,SAAK,IAAI,iBAAiB,KAAK,IAAI,IACjC,OAAO,KAAK,GAAG,EACf,KAAK,SAAS,wBAAwB,EACtC,KAAK,eAAe,6BAA6B,EACjD,GAAG,SAAS,WAAS;AACrB,YAAM,uBAAuB,KAAK,6BAA6B,KAAK;AACpE,UAAI,sBAAsB;AACzB,aAAK,wBAAwB,KAAK,yBAAyB,sBAAsB,IAAI;AACrF,aAAK,sBAAsB,KAAK,oBAAoB;AAEpD,cAAM,qBAAqB,KAAK,gBAAgB,WAAW,IAAI,eAAe,IAAI,oBAAoB;AAEtG,cAAM,8BAA8B,mBAAmB,SAAS,IAAI,OAAK,EAAE,IAAI;AAE/E,aAAK,uBAAuB,6BAA6B,KAAK;AAAA,MAC/D,OAAO;AAEN,eAAO,KAAK;AAAA,MACb;AAGA,UAAI,KAAK,mBAAmB;AAE3B,eAAO,KAAK;AACZ,aAAK,qBAAqB;AAAA,MAC3B,MAAO,MAAK,qBAAqB,MAAM;AAAA,IACxC,CAAC;AAAA,EAEH;AAAA,EAEA,MAAM,mBAAmB,QAAQ,CAAC,GAAG;AACpC,SAAK,iBAAiB,KAAK;AAC3B,UAAM,CAAC,GAAG,KAAK,IAAI,MAAM,KAAK,YAAY,CAAC,CAAC;AAC5C,QAAI,EAAE,MAAO,OAAM,EAAE;AACrB,SAAK,iBAAiB,gBAAgB,CAAC;AACvC,QAAI,CAAC,UAAU,KAAK,gBAAgB,KAAK,cAAc,GAAG;AAGzD,aAAO,KAAK;AACZ,aAAO,KAAK;AAAA,IACb;AACA,UAAM,IAAI,KAAK,SAAS;AAExB,QAAI,CAAC,EAAE,YAAY;AAElB,UAAI,EAAE,MAAM;AAEX,cAAM,4DAA4D,EAAE,IAAI;AAAA,MACzE;AAAA,IAED;AACA,SAAK,kBAAkB;AAEvB,UAAM,IAAI,KAAK,gBAAgB;AAC/B,SAAK,kBAAkB,CAAC;AAExB,UAAM,UAAU,CAAC;AAKjB,eAAW,CAAC,GAAG,MAAM,KAAK,EAAE,IAAI,MAAM,QAAQ,GAAG;AAChD,cAAQ,OAAO,IAAI,IAAI,EAAE,QAAQ,OAAO,KAAK;AAC7C,iBAAW,CAAC,GAAG,GAAG,KAAK,EAAE,IAAI,MAAM,QAAQ,GAAG;AAC7C,cAAM,KAAK,MAAM,KAAK,CAAAA,QAAMA,IAAG,QAAQ,IAAI,QAAQA,IAAG,OAAO,IAAI,IAAI;AAKrE,cAAM,QAAQ,EAAE,OAAO,CAAC,EAAE,CAAC;AAC3B,gBAAQ,OAAO,IAAI,EAAE,GAAG,GAAG,IAAI;AAAA,UAC9B,KAAK,GAAG,KAAK;AAAA,UACb,QAAQ;AAAA,YACP;AAAA,cACC,QAAQ,OAAO;AAAA,cACf,IAAI,aAAa,KAAK,MAAM,OAAO,QAAQ;AAAA,cAC3C,OAAO,EAAE;AAAA;AAAA;AAAA;AAAA,cAIT;AAAA;AAAA;AAAA,YAGD;AAAA,UACD;AAAA,QACD;AAAA,MACD;AAAA,IACD;AAEA,SAAK,kBACJ,KAAK,SAAS,YAAY,mBAAmB,aAC1C,MAAM,IAAI,OAAK,EAAE,KAAK,IAAI,IAC1B,KAAK,SAAS,YAAY,mBAAmB,WAC7C,MAAM,IAAI,OAAK,EAAE,KAAK,IAAI,EAAE,KAAK,IACjC,KAAK,OAAO,aAAa,0BACzB,EAAE,IAAI,MAAM,IAAI,SAAO,MAAM,KAAK,OAAK,EAAE,OAAO,IAAI,QAAQ,EAAE,KAAK,MAAM,IAAI,IAAI,GAAG,KAAK,IAAI,IAC7F,EAAE,IAAI,MAAM,IAAI,SAAO,MAAM,KAAK,OAAK,EAAE,OAAO,IAAI,IAAI,GAAG,KAAK,IAAI;AAExE,QAAI,KAAK,gBAAgB,SAAS,MAAS,EAAG,OAAM;AAEpD,SAAK,eAAe,CAAC,GAAG,MAAM;AAC7B,YAAM,IAAI,KAAK,gBAAgB,QAAQ,EAAE,GAAG,KAAK,IAAI;AACrD,YAAM,IAAI,KAAK,gBAAgB,QAAQ,EAAE,GAAG,KAAK,IAAI;AACrD,UAAI,KAAK,MAAM,KAAK,GAAI,QAAO;AAC/B,UAAI,KAAK,GAAI,QAAO;AACpB,UAAI,KAAK,GAAI,QAAO;AACpB,aAAO,IAAI;AAAA,IACZ;AAEA,SAAK,oBAAoB,EAAE,IAAI,MAAM,IAAI,SAAO,IAAI,IAAI;AACxD,SAAK,iBAAiB,CAAC,GAAG,MAAM;AAC/B,YAAM,IAAI,KAAK,kBAAkB,QAAQ,EAAE,MAAM;AACjD,YAAM,IAAI,KAAK,kBAAkB,QAAQ,EAAE,MAAM;AACjD,UAAI,KAAK,MAAM,KAAK,GAAI,QAAO;AAC/B,UAAI,KAAK,GAAI,QAAO;AACpB,UAAI,KAAK,GAAI,QAAO;AACpB,aAAO,IAAI;AAAA,IACZ;AAGA,UAAM,WAAW,CAAC;AAClB,eAAW,MAAM,OAAO;AACvB,UAAI,EAAE,SAAS,GAAG,KAAK,IAAI,EAAG,UAAS,GAAG,GAAG,IAAI,EAAE,SAAS,GAAG,KAAK,IAAI;AAAA,IACzE;AACA,SAAK,qBAAqB;AAAA,MACzB,MAAM,EAAE,UAAU,YAAY,EAAE,WAAW;AAAA,MAC3C,KAAK,EAAE,IAAI,MAAM,IAAI,CAAAC,OAAK,QAAQA,GAAE,IAAI,CAAC;AAAA,MACzC;AAAA,MACA,yBAAyB,EAAE;AAAA,IAC5B;AAAA,EACD;AAAA,EAEA,MAAM,cAAc;AAEnB,UAAM,OAAO,KAAK,iBAAiB,eAAe,KAAK,iBAAiB,KAAK,KAAK;AAClF,UAAM,QAAQ,KAAK,YAAY;AAC/B,UAAM,OAAO,MAAM,SAAS,kBAAkB,EAAE,MAAM,QAAQ,KAAK,IAAI,iBAAiB,EAAE,CAAC;AAE3F,WAAO,CAAC,MAAM,KAAK;AAAA,EACpB;AAAA,EAEA,iBAAiB,OAAO;AACvB,SAAK,cACJ,KAAK,OAAO,WAAW,KAAK,SAAO,IAAI,QAAQ,aAAa,KAC5D,KAAK,WAAW,KAAK,OAAK,EAAE,IAAI,QAAQ,aAAa,GAAG;AAEzD,UAAM,IAAI,MAAM,OAAO,SAAS;AAIhC,UAAM,yBAAyB;AAAA,MAC9B,MAAM;AAAA,MACN,IAAI;AAAA,MACJ,MAAM;AAAA,MACN,KAAK,MAAM,OAAO,kBAAkB,IAAI,OAAO,OAAK,CAAC,EAAE,IAAI,gBAAgB;AAAA,IAC5E;AACA,UAAM,QAAQ,KAAK,8BAA8B;AACjD,QAAI,CAAC,MAAM,OAAQ,OAAM;AAOzB,QAAI,CAAC,iBAAiB,KAAK,OAAK,EAAE,SAAS,EAAE,aAAa,EAAG,OAAM;AACnE,QAAI,CAAC,kBAAkB,KAAK,OAAK,EAAE,SAAS,EAAE,cAAc,EAAG,OAAM;AACrE,UAAM,OAAO;AAAA,MACZ,QAAQ,MAAM,MAAM;AAAA,MACpB,SAAS,MAAM,MAAM;AAAA,MACrB,UAAU,MAAM,OAAO;AAAA,MACvB,eAAe,EAAE;AAAA,MACjB,gBAAgB,EAAE;AAAA,MAClB,sBAAsB,EAAE;AAAA,MACxB;AAAA,MACA,QAAQ,cAAc,WAAW,CAAC,MAAM,QAAQ,sBAAsB,CAAC,CAAC;AAAA,MACxE,SAAS,MAAM;AAAA,IAChB;AACA,WAAO;AAAA,EACR;AAAA,EAEA,cAAc;AACb,QAAI,CAAC,KAAK,mBAAoB;AAC9B,UAAM,IAAI,KAAK;AACf,UAAM,0BAA0B,KAAK,mBAAmB;AACxD,UAAM,UAAU,CAAC;AACjB,UAAM,MAAM,CAAC;AAEb,eAAW,YAAY,KAAK,mBAAmB,SAAS;AACvD,YAAM,IAAI,KAAK,mBAAmB,QAAQ,QAAQ;AAClD,cAAQ,QAAQ,IAAI;AACpB,UAAI,KAAK,CAAC;AACV,UAAI,YAAY,EAAE,QAAS,QAAO,OAAO,GAAG,EAAE,QAAQ,QAAQ,CAAC;AAC/D,YAAM,QAAQ,KAAK,mBAAmB,KAAK,WAAW,QAAQ,KAAK,CAAC;AACpE,UAAI,CAAC,EAAE,MAAO,GAAE,QAAQ;AAAA,UAEnB,QAAO,OAAO,EAAE,OAAO,KAAK;AAAA,IAClC;AAGA,UAAM,IAAI,KAAK,mBAAmB,KAAK;AACvC,eAAW,OAAO,OAAO,KAAK,CAAC,GAAG;AACjC,QAAE,KAAK,SAAS,GAAG,IAAI,OAAO,OAAO,CAAC,GAAG,EAAE,KAAK,SAAS,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC;AAAA,IAC5E;AACA,SAAK,OAAO,EAAE,SAAS,KAAK,MAAM,EAAE,MAAM,wBAAwB;AAAA,EACnE;AAAA,EAEA,kBAAkB,GAAG;AACpB,UAAM,KAAK,KAAK,SAAS;AACzB,UAAM,QAAQ,OAAY,cAAc,GAAG,UAAU,EAAE,QAAQ,cAAc,GAAG,UAAU,EAAE,KAAK,EAAE,MAAM,IAAI;AAC7G,UAAM,iBAAiB,CAAC;AACxB,eAAW,OAAO,EAAE,QAAQ;AAC3B,qBAAe,KAAK,GAAG,OAAO,GAAG,CAAC;AAAA,IACnC;AACA,UAAM,SAAS,KAAK,IAAI,GAAG,WAAW,KAAK,IAAI,GAAG,OAAO,cAAc,EAAE,IAAI,KAAK,GAAG,CAAC,CAAC;AAEvF,UAAM,CAAC,KAAK,GAAG,IAAI,GAAG,uBACnB,CAAC,CAAC,QAAQ,MAAM,IAChB,CAAC,KAAK,IAAI,GAAG,cAAc,GAAG,KAAK,IAAI,GAAG,cAAc,CAAC;AAE5D,SAAK,oBAAoB,EAAE,OAAO,KAAK,IAAI;AAAA,EAC5C;AAAA,EAEA,cAAc,OAAO;AACpB,UAAM,KAAK,KAAK,SAAS;AACzB,QAAI,GAAG,sBAAsB;AAC5B,YAAM,YAAY,KAAK,SAAS,YAAY;AAC5C,aAAO,KAAK,kBAAkB,OAAO,QAAQ,CAAC,cAAc,YAAY,EAAE;AAAA,IAC3E,OAAO;AACN,aAAO,KAAK,kBAAkB,MAAM,QAAQ,KAAK,kBAAkB,GAAG;AAAA,IACvE;AAAA,EACD;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA,EAeA,gCAAgC;AAC/B,UAAM,MAAM,KAAK,YAAY,IAAI,IAAI,KAAK,KAAK,IAAI,SAAS,YAAY;AAExE,QAAI,KAAK,CAAC,GAAG,MAAO,EAAE,KAAK,OAAO,EAAE,KAAK,OAAO,KAAK,CAAE;AACvD,WAAO;AAAA,EACR;AACD;AAEA,WAAW,WAAW,CAAC,+BAAW,iCAAa,GAAG;AACjD,aAAW,cAAc,QAAS,aAAY,UAAU,UAAU,IAAI,QAAQ,UAAU;AACzF;AAEO,IAAM,kBAAkB,YAAY,WAAW;AAC/C,IAAM,gBAAgB;",
+  "names": ["tw", "c"]
+}

package/dist/chunk-5HXU7P5T.js

@@ -0,0 +1,276 @@
+import {
+  parsesample
+} from "./chunk-6ADTREVX.js";
+import {
+  dtfusionrna,
+  dtsv,
+  mclassfusionrna,
+  mclasssv
+} from "./chunk-6VQ4C735.js";
+
+// ../shared/utils/src/bulk.sv.js
+function parseheader(line, flag, issv) {
+  const header = line.toLowerCase().split("	");
+  if (header.length <= 1) return "invalid file header for fusions";
+  const htry = (...lst) => {
+    for (const a of lst) {
+      const j = header.indexOf(a);
+      if (j != -1) return j;
+    }
+    return -1;
+  };
+  let i = htry("gene_a", "gene1", "genea");
+  if (i == -1) return "gene_a missing from header";
+  header[i] = "gene1";
+  i = htry("gene_b", "gene2", "geneb");
+  if (i == -1) return "gene_b missing from header";
+  header[i] = "gene2";
+  i = htry("chr_a", "chr1", "chra");
+  if (i == -1) return "chr_a missing from header";
+  header[i] = "chr1";
+  i = htry("chr_b", "chr2", "chrb");
+  if (i == -1) return "chr_b missing from header";
+  header[i] = "chr2";
+  i = htry("pos_a", "position_a", "position1", "posa");
+  if (i == -1) return "pos_a missing from header";
+  header[i] = "position1";
+  i = htry("pos_b", "position_b", "position2", "posb");
+  if (i == -1) return "pos_b missing from header";
+  header[i] = "position2";
+  i = htry("isoform_a", "refseq_a", "refseq1", "isoform1", "sv_refseqa");
+  if (i == -1) return "isoform_a missing from header";
+  header[i] = "isoform1";
+  i = htry("isoform_b", "refseq_b", "refseq2", "isoform2", "sv_refseqb");
+  if (i == -1) return "isoform_b missing from header";
+  header[i] = "isoform2";
+  i = htry("strand_a", "orta");
+  if (i == -1) return "strand_a missing from header";
+  header[i] = "strand1";
+  i = htry("strand_b", "ortb");
+  if (i == -1) return "strand_b missing from header";
+  header[i] = "strand2";
+  i = htry("sample", "sample_name", "tumor_sample_barcode");
+  if (i != -1) header[i] = "sample";
+  i = htry("patient", "donor", "target_case_id");
+  if (i != -1) header[i] = "patient";
+  i = htry("sampletype", "sample type", "sample_type");
+  if (i != -1) header[i] = "sampletype";
+  i = htry("disease");
+  if (i != -1) header[i] = "disease";
+  i = htry("origin");
+  if (i != -1) header[i] = "origin";
+  if (issv) {
+    flag.sv.loaded = true;
+    flag.sv.header = header;
+  } else {
+    flag.fusion.loaded = true;
+    flag.fusion.header = header;
+  }
+  return false;
+}
+function parseline(i, line, flag, issv) {
+  if (line == "" || line[0] == "#") return;
+  const lst = line.split("	");
+  const m = {};
+  const header = issv ? flag.sv.header : flag.fusion.header;
+  const badlines = issv ? flag.sv.badlines : flag.fusion.badlines;
+  for (let j = 0; j < header.length; j++) {
+    m[header[j]] = lst[j];
+  }
+  if (!m.chr1) {
+    badlines.push([i, "missing chr1", lst]);
+    return;
+  }
+  if (m.chr1.toLowerCase().indexOf("chr") != 0) {
+    m.chr1 = "chr" + m.chr1;
+  }
+  if (!m.chr2) {
+    badlines.push([i, "missing chr2", lst]);
+    return;
+  }
+  if (m.chr2.toLowerCase().indexOf("chr") != 0) {
+    m.chr2 = "chr" + m.chr2;
+  }
+  let v = m.position1;
+  if (!v) {
+    badlines.push([i, "missing position1", lst]);
+    return;
+  }
+  let v2 = Number.parseInt(v);
+  if (Number.isNaN(v2) || v2 <= 0) {
+    badlines.push([i, "invalid value for position1", lst]);
+    return;
+  }
+  m.position1 = v2;
+  v = m.position2;
+  if (!v) {
+    badlines.push([i, "missing position2", lst]);
+    return;
+  }
+  v2 = Number.parseInt(v);
+  if (Number.isNaN(v2) || v2 <= 0) {
+    badlines.push([i, "invalid value for position2", lst]);
+    return;
+  }
+  m.position2 = v2;
+  if (parsesample(m, flag, i, lst, badlines)) {
+    return;
+  }
+  if (m.isoform1 && m.isoform1.indexOf(",") != -1) {
+    const lst2 = m.isoform1.split(",");
+    m.isoform1 = void 0;
+    for (const t of lst2) {
+      if (t != "") m.isoform1 = t;
+    }
+  }
+  if (m.isoform2 && m.isoform2.indexOf(",") != -1) {
+    const lst2 = m.isoform2.split(",");
+    m.isoform2 = void 0;
+    for (const t of lst2) {
+      if (t != "") m.isoform2 = t;
+    }
+  }
+  if (!m.gene1) {
+    m.isoform1 = void 0;
+  }
+  if (!m.gene2) {
+    m.isoform2 = void 0;
+  }
+  if (m.gene1) {
+    flag.good++;
+    const m2 = {
+      dt: issv ? dtsv : dtfusionrna,
+      class: issv ? mclasssv : mclassfusionrna,
+      isoform: m.isoform1,
+      mname: m.gene2 || m.chr2,
+      sample: m.sample,
+      patient: m.patient,
+      sampletype: m.sampletype,
+      origin: m.origin,
+      disease: m.disease,
+      pairlst: [
+        {
+          a: {
+            name: m.gene1,
+            isoform: m.isoform1,
+            strand: m.strand1,
+            chr: m.chr1,
+            position: m.position1
+          },
+          b: {
+            name: m.gene2,
+            isoform: m.isoform2,
+            strand: m.strand2,
+            chr: m.chr2,
+            position: m.position2
+          }
+        }
+      ]
+    };
+    const n = flag.geneToUpper ? m.gene1.toUpperCase() : m.gene1;
+    if (!flag.data[n]) {
+      flag.data[n] = [];
+    }
+    flag.data[n].push(m2);
+  }
+  if (m.gene2 && m.gene2 != m.gene1) {
+    flag.good++;
+    const m2 = {
+      dt: issv ? dtsv : dtfusionrna,
+      class: issv ? mclasssv : mclassfusionrna,
+      isoform: m.isoform2,
+      mname: m.gene1 || m.chr1,
+      sample: m.sample,
+      patient: m.patient,
+      sampletype: m.sampletype,
+      origin: m.origin,
+      disease: m.disease,
+      pairlst: [
+        {
+          a: {
+            name: m.gene1,
+            isoform: m.isoform1,
+            strand: m.strand1,
+            chr: m.chr1,
+            position: m.position1
+          },
+          b: {
+            name: m.gene2,
+            isoform: m.isoform2,
+            strand: m.strand2,
+            chr: m.chr2,
+            position: m.position2
+          }
+        }
+      ]
+    };
+    const n = flag.geneToUpper ? m.gene2.toUpperCase() : m.gene2;
+    if (!flag.data[n]) {
+      flag.data[n] = [];
+    }
+    flag.data[n].push(m2);
+  }
+}
+function duplicate(m) {
+  const n = {};
+  for (const k in m) {
+    if (k == "pairlst") continue;
+    const v = m[k];
+    const type = typeof v;
+    if (type == "object") {
+      continue;
+    }
+    n[k] = v;
+  }
+  if (m.pairlst) {
+    n.pairlst = [];
+    for (const pair of m.pairlst) {
+      const p = {};
+      for (const k in pair) {
+        if (k == "a" || k == "b" || k == "interstitial") {
+          continue;
+        }
+        p[k] = pair[k];
+      }
+      if (pair.a) {
+        p.a = {};
+        for (const k in pair.a) {
+          const v = pair.a[k];
+          if (typeof v == "object") {
+            continue;
+          }
+          p.a[k] = v;
+        }
+      }
+      if (pair.b) {
+        p.b = {};
+        for (const k in pair.b) {
+          const v = pair.b[k];
+          if (typeof v == "object") {
+            continue;
+          }
+          p.b[k] = v;
+        }
+      }
+      if (pair.interstitial) {
+        p.interstitial = {};
+        for (const k in pair.interstitial) {
+          const v = pair.interstitial[k];
+          if (typeof v == "object") {
+            continue;
+          }
+          p.interstitial[k] = v;
+        }
+      }
+      n.pairlst.push(p);
+    }
+  }
+  return n;
+}
+
+export {
+  parseheader,
+  parseline,
+  duplicate
+};
+//# sourceMappingURL=chunk-5HXU7P5T.js.map