@sjcrh/proteinpaint-client 2.204.0 → 2.205.0

package/dist/chunk-6JBQLOJW.js

@@ -0,0 +1,1339 @@
+import {
+  ColorScale,
+  axisstyle,
+  drawBoxplot,
+  first_genetrack_tolist,
+  makeSsmLink,
+  sayerror
+} from "./chunk-J7JDCNLU.js";
+import {
+  Menu
+} from "./chunk-ELJX3QIQ.js";
+import {
+  axisBottom,
+  axisTop
+} from "./chunk-YLJOZP4P.js";
+import {
+  linear,
+  log
+} from "./chunk-HDTFYTEL.js";
+import {
+  roundValue
+} from "./chunk-TLT4YIG3.js";
+import {
+  rgb_default
+} from "./chunk-Q5RDQNIT.js";
+
+// plots/regression/estimateMsg.ts
+var refGrp_NA = "NA";
+function getEstimateMsg(arg) {
+  const {
+    est,
+    tw,
+    tw2,
+    categoryKey,
+    categoryKey2,
+    isIntercept,
+    isUnivariate,
+    outcomeTw,
+    independentTws,
+    termdbConfig,
+    getIndependentInput
+  } = arg;
+  const regtype = arg.regressionType;
+  const category = tw?.term?.values && tw.term.values[categoryKey] ? tw.term.values[categoryKey].label : categoryKey;
+  const category2 = tw2?.term?.values && tw2.term.values[categoryKey2] ? tw2.term.values[categoryKey2].label : categoryKey2;
+  const refGrp = tw?.term?.values && tw.term.values[tw.refGrp] ? tw.term.values[tw.refGrp].label : tw?.refGrp;
+  const refGrp2 = tw2?.term?.values && tw2.term.values[tw2.refGrp] ? tw2.term.values[tw2.refGrp].label : tw2?.refGrp;
+  if (!Number.isFinite(est) || regtype != "linear" && est <= 0) {
+    return "The estimate of this variable is not available.";
+  }
+  let msg;
+  if (regtype == "linear") {
+    msg = tw2 ? getInteractionMsg() : `Mean ${styleVariable(outcomeTw)} is`;
+    if (isIntercept) {
+      const baselines = getBaselines(independentTws);
+      return `${msg} ${est} ${unitsOf(outcomeTw)} when ${joinVariables(baselines)}.`;
+    }
+    msg += est == 0 ? " no different " : ` ${Math.abs(est)} ${unitsOf(outcomeTw)} ${est < 0 ? "lower" : "higher"} `;
+  } else if (regtype == "logistic") {
+    msg = tw2 ? getInteractionMsg() : `Odds of ${styleVariable(outcomeTw, outcomeTw.nonRefGrp)} is`;
+    if (isIntercept) {
+      const baselines = getBaselines(independentTws);
+      return `${msg} ${est} when ${joinVariables(baselines)}.`;
+    }
+    msg += getRatioMsg();
+  } else if (regtype == "cox") {
+    msg = tw2 ? getInteractionMsg() : `Hazard (instantaneous rate) of ${styleVariable(outcomeTw, outcomeTw.eventLabel)} is`;
+    msg += getRatioMsg();
+  } else {
+    throw "regression type not recognized";
+  }
+  const interactions = [];
+  const interactionsBaselines = [];
+  if (tw.interactions?.length && !tw2) {
+    for (const tid of tw.interactions) {
+      const t = getIndependentInput(tid).term;
+      if (t.term.snps) {
+        for (const snp of t.term.snps) interactions.push(snp.snpid);
+      } else {
+        interactions.push(tid);
+      }
+    }
+    if (!interactions.length) throw "interactions[] is empty";
+    const interactingTws = independentTws.filter((t) => interactions.includes(t.$id || t.id));
+    interactionsBaselines.push(...getBaselines(interactingTws));
+  }
+  if (category) {
+    msg += `in ${joinVariables([styleVariable(tw, category), ...interactionsBaselines])} compared to ${joinVariables([
+      styleVariable(tw, refGrp),
+      ...interactionsBaselines
+    ])}`;
+  } else if (tw.q.mode == "continuous") {
+    msg += `for every ${oneUnitOf(tw)} increase of ${styleVariable(tw)}`;
+    if (interactionsBaselines.length) msg += ` when ${joinVariables(interactionsBaselines)}`;
+  } else if (tw.q.geneticModel === 0) {
+    msg += `for every additional ${tw.effectAllele} allele of ${styleVariable(tw)}`;
+    if (interactionsBaselines.length) msg += ` when ${joinVariables(interactionsBaselines)}`;
+  } else if (tw.q.geneticModel == 1 || tw.q.geneticModel == 2) {
+    const gts = Object.keys(tw.gt2count);
+    const testGts = gts.filter((gt) => {
+      if (tw.q.geneticModel == 1) {
+        return gt.includes(tw.effectAllele);
+      } else {
+        return gt.replace(/[^a-zA-Z]/g, "").split("").every((c) => c == tw.effectAllele);
+      }
+    });
+    const refGts = gts.filter((gt) => !testGts.includes(gt));
+    msg += `in ${joinVariables([
+      styleVariable(tw, testGts.join(", ")),
+      ...interactionsBaselines
+    ])} compared to ${joinVariables([styleVariable(tw, refGts.join(", ")), ...interactionsBaselines])}`;
+  }
+  const tids = [tw.$id || tw.id];
+  if (tw.interactions?.length) {
+    if (tw2) tids.push(tw2.$id || tw2.id);
+    else tids.push(...interactions);
+  }
+  const covariates = independentTws.filter((t) => !tids.includes(t.$id || t.id)).map((t) => styleVariable(t));
+  if (covariates.length && !isUnivariate) msg += `, adjusting for ${joinVariables(covariates)}`;
+  msg += ".";
+  if (regtype == "cox") {
+    if (outcomeTw.q.timeScale == "age") {
+      msg += " Time is measured as attained age during follow-up.";
+    } else {
+      const unit = termdbConfig?.timeUnit ? ` in ${termdbConfig.timeUnit}` : "";
+      const start = termdbConfig?.cohortStartTimeMsg ? ` from ${termdbConfig.cohortStartTimeMsg}` : "";
+      if (unit || start) msg += ` Time is measured${unit}${start}.`;
+    }
+  }
+  return msg;
+  function oneUnitOf(tw3) {
+    const u = tw3?.term?.valueConversion?.toUnit;
+    return u ? `1 ${u}` : "one unit";
+  }
+  function unitsOf(tw3) {
+    const u = tw3?.term?.valueConversion?.toUnit;
+    return u ? `${u}s` : "units";
+  }
+  function getRatioMsg() {
+    if (est == 1) return " no different ";
+    return est > 1 ? ` ${est} times higher ` : ` ${roundValue(1 / est, 3)} times lower `;
+  }
+  function styleVariable(tw3, category3) {
+    const spans = [
+      `<span class="term_name_btn sja_filter_tag_btn" style="padding: 3px 6px; margin: 2.5px 0px; border-radius: ${category3 ? "6px 0px 0px 6px" : "6px"};">${tw3.term.name.length < 40 ? tw3.term.name : tw3.term.name.substring(0, 35) + " ..."}</span>`
+    ];
+    if (category3) {
+      spans.push(
+        `<span class="ts_summary_btn sja_filter_tag_btn" style="padding: 3px 6px; margin: 2.5px 0px; border-radius: 0px 6px 6px 0px; font-style: italic;">${category3}</span>`
+      );
+    }
+    return `<div style="display: inline; white-space: nowrap; font-size: 0.9em">${spans.join("")}</div>`;
+  }
+  function getInteractionMsg() {
+    let msg2 = regtype == "linear" ? `The difference in mean ${styleVariable(outcomeTw)}` : regtype == "logistic" ? `The difference in odds of ${styleVariable(outcomeTw, outcomeTw.nonRefGrp)}` : `The difference in hazard (instantaneous rate) of ${styleVariable(outcomeTw, outcomeTw.eventLabel)}`;
+    if (category2) {
+      msg2 += ` between ${styleVariable(tw2, category2)} and ${styleVariable(tw2, refGrp2)} is`;
+    } else if (tw2.q.mode == "continuous") {
+      msg2 += ` for every ${oneUnitOf(tw2)} increase of ${styleVariable(tw2)} is`;
+    } else if (tw2.q.geneticModel === 0) {
+      msg2 += ` for every additional ${tw2.effectAllele} allele of ${styleVariable(tw2)} is`;
+    } else if (tw2.q.geneticModel == 1 || tw2.q.geneticModel == 2) {
+      const gts = Object.keys(tw2.gt2count);
+      const testGts = gts.filter((gt) => {
+        if (tw2.q.geneticModel == 1) {
+          return gt.includes(tw2.effectAllele);
+        } else {
+          return gt.replace(/[^a-zA-Z]/g, "").split("").every((c) => c == tw2.effectAllele);
+        }
+      });
+      const refGts = gts.filter((gt) => !testGts.includes(gt));
+      msg2 += ` between ${styleVariable(tw2, testGts.join(", "))} and ${styleVariable(tw2, refGts.join(", "))} is`;
+    }
+    return msg2;
+  }
+  function getBaselines(tws) {
+    const baselines = tws.map((tw3) => {
+      if (tw3.q.mode != "spline" && "refGrp" in tw3 && tw3.refGrp != refGrp_NA) {
+        const refGrp3 = tw3?.term?.values && tw3.term.values[tw3.refGrp] ? tw3.term.values[tw3.refGrp].label : tw3?.refGrp;
+        return styleVariable(tw3, refGrp3);
+      } else if (tw3.q.mode == "continuous" || tw3.q.mode == "spline") {
+        return styleVariable(tw3, "0");
+      } else if (tw3.q.geneticModel === 0) {
+        return styleVariable(tw3, `No ${tw3.effectAllele} alleles`);
+      } else if (tw3.q.geneticModel == 1 || tw3.q.geneticModel == 2) {
+        const gts = Object.keys(tw3.gt2count);
+        const refGts = gts.filter((gt) => {
+          if (tw3.q.geneticModel == 1) {
+            return !gt.includes(tw3.effectAllele);
+          } else {
+            return !gt.replace(/[^a-zA-Z]/g, "").split("").every((c) => c == tw3.effectAllele);
+          }
+        });
+        return styleVariable(tw3, refGts.join(", "));
+      }
+    });
+    return baselines.filter(Boolean);
+  }
+  function joinVariables(variables) {
+    if (!variables.length) return "";
+    else if (variables.length == 1) return variables[0];
+    else if (variables.length == 2) return variables.join(" and ");
+    else return `${variables.slice(0, -1).join(", ")}, and ${variables.slice(-1)}`;
+  }
+}
+
+// plots/regression/regression.results.ts
+var forestcolor = "#126e08";
+var boxplotcolor = forestcolor;
+var graytextcolor = "#555";
+var graytextopacity = 0.75;
+var RegressionResults = class {
+  constructor(opts) {
+    this.opts = opts;
+    this.app = opts.app;
+    this.parent = opts.parent;
+    this.vocabApi = this.parent.vocabApi;
+    this.type = "regression";
+    setRenderers(this);
+    const holder = this.opts.holder;
+    holder.append("div").style("margin-top", "10px").style("padding-top", "20px").style("font-size", "1.2em").style("opacity", graytextopacity).html("Results");
+    this.dom = {
+      holder,
+      err_div: holder.append("div"),
+      snplocusBlockDiv: holder.append("div"),
+      // is where newDiv() and displayResult_oneset() writes to
+      oneSetResultDiv: holder.append("div").style("margin", "10px"),
+      tip: new Menu({ padding: "9px" })
+    };
+  }
+  async main() {
+    try {
+      this.config = this.parent.config;
+      this.state = this.parent.state;
+      if (!this.state.formIsComplete || this.parent.inputs.hasError || this.config.hasUnsubmittedEdits && !this.hasUnsubmittedEdits_nullify_singleuse) {
+        delete this.snplocusBlock;
+        this.dom.snplocusBlockDiv.selectAll("*").remove();
+        this.dom.holder.style("display", "none");
+        return;
+      }
+      delete this.hasUnsubmittedEdits_nullify_singleuse;
+      if (this.snplocusBlock) {
+        this.snplocusBlock.cloakOn();
+      }
+      this.parent.inputs.dom.submitBtn.text("Running...");
+      this.parent.dom.inputs.style("opacity", 0.5).style("pointer-events", "none");
+      const data = await this.vocabApi.getRegressionData(this.getDataRequestOpts());
+      if (data.error) throw data.error;
+      this.dom.err_div.style("display", "none");
+      this.dom.oneSetResultDiv.selectAll("*").remove();
+      this.dom.holder.style("display", "block");
+      await this.displayResult(data);
+      const results_y = this.dom.holder.node().getBoundingClientRect().top + window.scrollY;
+      const nav_height = document.querySelector(".sjpp-nav")?.getBoundingClientRect().height;
+      if (!nav_height || isNaN(nav_height)) throw new Error("cannot get nav height");
+      window.scroll({ behavior: "smooth", top: results_y - nav_height });
+    } catch (e) {
+      this.hasError = true;
+      this.dom.holder.style("display", "block");
+      this.dom.err_div.style("display", "block");
+      sayerror(this.dom.err_div, "Error: " + (e.error || e));
+      this.parent.inputs.dom.submitBtn.property("disabled", true);
+      if (e.stack) console.log(e.stack);
+      console.error(e);
+    } finally {
+      this.parent.dom.inputs.style("opacity", 1).style("pointer-events", "auto");
+    }
+  }
+  // creates an opts object for the vocabApi.getRegressionData()
+  getDataRequestOpts() {
+    const c = this.config;
+    const opts = {
+      regressionType: c.regressionType,
+      outcome: c.outcome,
+      independent: c.independent,
+      includeUnivariate: c.includeUnivariate
+    };
+    opts.filter = this.parent.filter;
+    opts.filter0 = this.state.termfilter.filter0;
+    return opts;
+  }
+  getIndependentInput(tid) {
+    for (const i of this.parent.inputs.independent.inputLst) {
+      if (!i.term) continue;
+      if (i.term.term && i.term.term.snps) {
+        for (const snp of i.term.term.snps) {
+          if (snp.snpid == tid) {
+            const tw = {
+              id: tid,
+              q: {
+                geneticModel: i.term.q.geneticModel
+              },
+              term: {
+                id: tid,
+                name: tid
+              },
+              interactions: i.term.interactions,
+              effectAllele: i.term.q.snp2effAle[tid],
+              gt2count: snp.gt2count
+            };
+            if (i.term.q.snp2refGrp) {
+              tw.refGrp = i.term.q.snp2refGrp[tid];
+            }
+            if (snp.mlst) {
+              const m = snp.mlst.find((j) => j.alt == i.term.q.snp2effAle[tid]);
+              if (m) {
+                tw.term.name = m.mname;
+              } else {
+                tw.term.name = snp.mlst[0].mname;
+              }
+            }
+            return { term: tw };
+          }
+        }
+      }
+      if (i.term.$id == tid) return i;
+    }
+    return {
+      term: {
+        id: tid,
+        q: { mode: "continuous" },
+        term: { name: tid }
+      }
+    };
+  }
+};
+function setRenderers(self) {
+  self.displayResult = async (result) => {
+    const snplocusInput = self.parent.inputs.independent.inputLst.find((i) => i.term && i.term.term.type == "snplocus");
+    if (snplocusInput) {
+      if (!self.snplocusBlock) {
+        self.dom.snplocusBlockDiv.append("div").style("margin-top", "30px").style("opacity", graytextopacity).text("Click on a variant within the browser to view its regression results");
+        self.snplocusBlock = await createGenomebrowser(self, snplocusInput, result.resultLst);
+      } else {
+        await updateMds3Tk(self, snplocusInput, result.resultLst);
+      }
+      return;
+    }
+    delete self.snplocusBlock;
+    self.dom.snplocusBlockDiv.selectAll("*").remove();
+    if (!result.resultLst[0] || !result.resultLst[0].data) throw "result is not [ {data:{}} ]";
+    self.displayResult_oneset(result.resultLst[0].data);
+  };
+  self.displayResult_oneset = (result) => {
+    self.dom.oneSetResultDiv.selectAll("*").remove();
+    self.dom.LDresultDiv = self.dom.oneSetResultDiv.append("div");
+    self.mayshow_warn(result);
+    if (result.sampleSize) self.newDiv("Sample size:", result.sampleSize);
+    if (result.eventCnt) self.newDiv("Number of events:", result.eventCnt);
+    self.mayshow_headerRow(result);
+    self.mayshow_splinePlots(result);
+    self.mayshow_residuals(result);
+    self.mayshow_coefficients(result);
+    self.mayshow_coxDisclaimer();
+    self.mayshow_totalSnpEffect(result);
+    self.mayshow_type3(result);
+    self.mayshow_nonlinearity(result);
+    self.mayshow_tests(result);
+    self.mayshow_other(result);
+    self.mayshow_fisher(result);
+    self.mayshow_wilcoxon(result);
+    self.mayshow_cuminc(result);
+  };
+  self.newDiv = (label, label2, getrow) => {
+    const div = self.dom.oneSetResultDiv.append("div").style("margin", "20px 0px 10px 0px").attr("name", label);
+    const row = div.append("div");
+    row.append("span").style("text-decoration", "underline").text(label);
+    if (label2) {
+      row.append("span").html(label2).style("margin-left", "5px");
+    }
+    return getrow ? row : div.append("div").style("margin-left", "20px");
+  };
+  self.mayshow_warn = (result) => {
+    if (!result.warnings) return;
+    const div = self.newDiv("Warnings");
+    const warnings = new Set(result.warnings);
+    for (const line of warnings) {
+      div.append("p").style("margin", "5px").text(line);
+    }
+  };
+  self.mayshow_headerRow = (result) => {
+    if (!result.headerRow) return;
+    const k = result.headerRow.k;
+    const v = result.headerRow.v;
+    const snplocusInput = self.parent.inputs.independent.inputLst.find((i) => i.term && i.term.term.type == "snplocus");
+    if (snplocusInput) {
+      const snp = snplocusInput.term.term.snps.find((snp2) => snp2.snpid == v.snpid);
+      const m = snp.mlst[0];
+      m.chr = snp.chr;
+      const row = self.newDiv(k, null, true);
+      const snpLabelDom = row.append("span").text(`${m.chr}:${m.pos + 1} ${m.ref && m.alt ? m.ref + ">" + m.alt : ""}`).style("margin-left", "5px");
+      const urlConfig = self.app.vocabApi.termdbConfig.urlTemplates?.ssm || self.app.vocabApi.termdbConfig.queries?.snvindel?.ssmUrl;
+      if (urlConfig) {
+        const separateUrls = makeSsmLink(urlConfig, m, snpLabelDom, self.parent.genomeObj.name);
+        if (separateUrls?.length) {
+          row.append("span").style("margin-left", "10px").html(separateUrls.join(" "));
+        }
+      }
+      let labels;
+      const gt_label = `Genotypes: ${v.gtcounts.join(", ")}`;
+      if (v.monomorphic) {
+        labels = [gt_label];
+      } else {
+        const effale_label = `Effect allele: ${v.effAle}`;
+        const af_label = `Allele frequency: ${v.af}`;
+        labels = [effale_label, af_label, gt_label];
+      }
+      row.append("span").html(`&nbsp;&#65372;&nbsp;${labels.join("&nbsp;&#65372;&nbsp;")}`);
+    } else {
+      self.newDiv(k, v);
+    }
+  };
+  self.mayshow_splinePlots = (result) => {
+    if (!result.splinePlots) return;
+    const div = self.newDiv("Cubic spline plots");
+    div.style("display", "flex").style("align-items", "center");
+    result.splinePlots.sort((a, b) => {
+      if (a.type == "univariate" && b.type == "multivariate") return -1;
+      if (a.type == "multivariate" && b.type == "univariate") return 1;
+      return 0;
+    });
+    for (const plot of result.splinePlots) {
+      const plotDiv = div.append("div").style("margin", "0px 50px 5px 0px");
+      plotDiv.append("img").attr("src", plot.src).attr("width", 670);
+    }
+  };
+  self.mayshow_residuals = (result) => {
+    if (!result.residuals) return;
+    const div = self.newDiv(result.residuals.label);
+    const table = div.append("table").style("border-spacing", "8px").attr("name", "sjpp-residuals-table");
+    const tr1 = table.append("tr").style("opacity", graytextopacity);
+    const tr2 = table.append("tr");
+    for (let i = 0; i < result.residuals.header.length; i++) {
+      tr1.append("td").text(result.residuals.header[i]);
+      tr2.append("td").text(result.residuals.rows[i]);
+    }
+  };
+  self.mayshow_cuminc = async (result) => {
+    if (!result.cuminc) return;
+    const holder = self.newDiv(
+      "Cumulative incidence test:"
+      /*, 'p-value = ' + result.cuminc.pvalue*/
+    );
+    const _ = await import("./Cuminc-6F2C5C4E.js");
+    const plotter = new _.Cuminc({
+      holder,
+      config: {
+        term: self.config.outcome,
+        term2: {
+          term: {
+            name: "Variant",
+            values: {
+              1: { key: 1, label: "Has minor allele" },
+              2: { key: 2, label: "No minor allele" }
+            }
+          }
+        }
+      }
+    });
+    if (result.cuminc.ci_data) {
+      plotter.main(result.cuminc.ci_data);
+    } else {
+      holder.append("div").style("margin", "20px").text(result.cuminc.msg);
+    }
+  };
+  self.mayshow_wilcoxon = (result) => {
+    if (!result.wilcoxon) return;
+    const div = self.newDiv("Wilcoxon rank sum test:", "p-value = " + result.wilcoxon.pvalue);
+    if (result.wilcoxon.boxplots) {
+      const bs = result.wilcoxon.boxplots;
+      const boxplotHeight = 20, boxplotWidth = 400, leftLabelWidth = 160, axisheight = 40, labpad = 20, vpad = 10;
+      const scale = linear().domain([bs.minv, bs.maxv]).range([0, boxplotWidth]);
+      const svg = div.append("svg").style("margin-top", "10px").attr("width", leftLabelWidth + labpad + boxplotWidth + 10).attr("height", vpad * 3 + boxplotHeight * 2 + axisheight);
+      const g = svg.append("g").attr("transform", `translate(${leftLabelWidth + labpad},${vpad})`);
+      drawBoxplot({
+        g: g.append("g"),
+        bp: bs.hasEff,
+        scale,
+        rowheight: boxplotHeight,
+        color: boxplotcolor,
+        labpad
+      });
+      drawBoxplot({
+        g: g.append("g").attr("transform", `translate(0,${boxplotHeight + vpad})`),
+        bp: bs.noEff,
+        scale,
+        rowheight: boxplotHeight,
+        color: boxplotcolor,
+        labpad
+      });
+      {
+        const axisg = g.append("g").attr("transform", `translate(0,${boxplotHeight * 2 + vpad * 2})`);
+        const axis = axisBottom(scale);
+        axisstyle({
+          axis: axisg.call(axis),
+          color: boxplotcolor,
+          showline: true
+        });
+        axisg.append("text").text(self.config.outcome.term.name).attr("font-size", 15).attr("x", boxplotWidth / 2).attr("y", axisheight - 5).attr("text-anchor", "middle").attr("fill", boxplotcolor);
+      }
+    }
+  };
+  self.mayshow_fisher = (result) => {
+    if (!result.fisher) return;
+    const div = self.newDiv(
+      result.fisher.isChi ? "Chi-square test:" : "Fisher's exact test:",
+      "p-value = " + result.fisher.pvalue
+    );
+    const table = div.append("table").style("margin", "20px").style("border-spacing", "5px").style("border-collapse", "separate");
+    for (const r of result.fisher.rows) {
+      const tr = table.append("tr");
+      for (const c of r) {
+        tr.append("td").text(c);
+      }
+    }
+  };
+  self.mayshow_coefficients = (result) => {
+    if (!result.coefficients) {
+      if (result.coefficients_uni && result.coefficients_multi) {
+        self.mayshow_coefficients_uniMulti(result);
+      }
+      return;
+    }
+    const div = self.newDiv(result.coefficients.label);
+    const table = div.append("table").style("border-spacing", "0px").attr("data-testid", "sjpp-regression-resultCoefficientTable");
+    const forestPlotter = self.getForestPlotter(result.coefficients.terms, result.coefficients.interactions);
+    {
+      const tr = table.append("tr");
+      const header = result.coefficients.header;
+      tr.append("td").text(header.shift()).style("padding", "8px");
+      tr.append("td").text(header.shift()).style("padding", "8px");
+      if (self.config.regressionType == "cox") {
+        header.shift();
+        header.shift();
+      }
+      header.splice(1, 2, "95% CI");
+      self.fillDataHeaders(header, tr, void 0, void 0, forestPlotter);
+    }
+    self.independentTws = Object.keys(result.coefficients.terms).map((tid) => self.getIndependentInput(tid).term);
+    let varcount = 0;
+    const intercept = result.coefficients.intercept;
+    if (intercept) {
+      const tr = table.append("tr").style("background", ++varcount % 2 ? "#eee" : "none");
+      tr.append("td").text(intercept.shift()).style("padding", "8px");
+      tr.append("td").text(intercept.shift()).style("padding", "8px");
+      tr.append("td");
+      self.fillCoefDataCols({ tr, cols: intercept, isIntercept: true });
+    }
+    let rowcolor;
+    for (const tid in result.coefficients.terms) {
+      const termdata = result.coefficients.terms[tid];
+      const tw = self.getIndependentInput(tid).term;
+      rowcolor = ++varcount % 2 ? "#eee" : "none";
+      let tr = table.append("tr").style("background", rowcolor);
+      const termNameTd = tr.append("td").style("padding", "8px");
+      fillCoefficientTermname(tw, termNameTd);
+      if (termdata.fields) {
+        const cols = termdata.fields;
+        {
+          const td = tr.append("td").style("padding", "8px");
+          fillColumn2coefficientsTable(td, tw);
+        }
+        if (self.config.regressionType == "cox") {
+          cols.shift();
+          cols.shift();
+        }
+        forestPlotter(tr.append("td"), cols);
+        self.fillCoefDataCols({ tr, cols, tw });
+      } else if (termdata.categories) {
+        const orderedCategories = [];
+        const input = self.getIndependentInput(tid);
+        if (input.orderedLabels) {
+          for (const k of input.orderedLabels) {
+            if (termdata.categories[k]) orderedCategories.push(k);
+          }
+        }
+        for (const k in termdata.categories) {
+          if (!orderedCategories.includes(k)) orderedCategories.push(k);
+        }
+        termNameTd.attr("rowspan", orderedCategories.length).style("vertical-align", "top");
+        let isfirst = true;
+        for (const k of orderedCategories) {
+          if (!isfirst) {
+            tr = table.append("tr").style("background", rowcolor);
+          }
+          const cols = termdata.categories[k];
+          const td = tr.append("td").style("padding", "8px");
+          fillColumn2coefficientsTable(td, tw, k);
+          if (self.config.regressionType == "cox") {
+            if (tw.q.mode == "spline") {
+              cols.shift();
+              cols.shift();
+            } else {
+              const [samplesize_ref, samplesize_c] = cols.shift().split("/");
+              const [eventcnt_ref, eventcnt_c] = cols.shift().split("/");
+              if (isfirst) {
+                const variableBottomDiv = termNameTd.select(".sjpcb-coef-variable-bottom");
+                variableBottomDiv.style("align-items", "baseline");
+                const refGrpDiv = variableBottomDiv.selectAll("div").filter((d, i) => i === 1);
+                refGrpDiv.append("div").html(`n=${samplesize_ref}<br>events=${eventcnt_ref}`);
+              }
+              td.append("div").style("font-size", ".8em").html(`n=${samplesize_c}<br>events=${eventcnt_c}`);
+            }
+          }
+          forestPlotter(tr.append("td"), cols);
+          self.fillCoefDataCols({ tr, cols, tw, categoryKey: k });
+          isfirst = false;
+        }
+      } else {
+        tr.append("td").text("ERROR: no .fields[] or .categories{}");
+      }
+    }
+    for (const i of result.coefficients.interactions) {
+      rowcolor = ++varcount % 2 ? "#eee" : "none";
+      let tr = table.append("tr").style("background", rowcolor);
+      const term1 = self.getIndependentInput(i.term1).term;
+      const term2 = self.getIndependentInput(i.term2).term;
+      {
+        const td = tr.append("td").style("padding", "8px");
+        fillTdName(td.append("div"), term1 ? term1.term.name + " : " : i.term1 + " : ");
+        fillTdName(td.append("div"), term2 ? term2.term.name : i.term2);
+        td.attr("rowspan", i.categories.length).style("vertical-align", "top");
+      }
+      let isfirst = true;
+      for (const c of i.categories) {
+        if (!isfirst) tr = table.append("tr").style("background", rowcolor);
+        const td = tr.append("td").style("padding", "8px");
+        fillColumn2coefficientsTable(td.append("div"), term1, c.category1);
+        fillColumn2coefficientsTable(td.append("div"), term2, c.category2);
+        const cols = c.lst;
+        if (self.config.regressionType == "cox") {
+          cols.shift();
+          cols.shift();
+        }
+        forestPlotter(tr.append("td"), cols);
+        self.fillCoefDataCols({ tr, cols, tw: term1, tw2: term2, categoryKey: c.category1, categoryKey2: c.category2 });
+        isfirst = false;
+      }
+    }
+  };
+  self.mayshow_coefficients_uniMulti = (result) => {
+    if (!result.coefficients_uni || !result.coefficients_multi) return;
+    const div = self.newDiv(result.coefficients_uni.label);
+    div.style("margin-bottom", "200px");
+    const table = div.append("table").style("border-spacing", "0px").attr("data-testid", "sjpp-regression-resultCoefficientTable");
+    const forestPlotter_uni = self.getForestPlotter(result.coefficients_uni.terms, result.coefficients_uni.interactions);
+    const forestPlotter_multi = self.getForestPlotter(
+      result.coefficients_multi.terms,
+      result.coefficients_multi.interactions
+    );
+    {
+      const tr_label = table.append("tr").style("opacity", graytextopacity);
+      const tr = table.append("tr");
+      const header_uni = result.coefficients_uni.header;
+      const header_multi = result.coefficients_multi.header;
+      tr.append("td").text(header_uni.shift()).style("padding", "8px");
+      tr_label.append("td").style("padding", "8px");
+      header_multi.shift();
+      tr.append("td").text(header_uni.shift()).style("padding", "8px");
+      tr_label.append("td").style("padding", "8px");
+      header_multi.shift();
+      if (self.config.regressionType == "cox") {
+        header_uni.shift();
+        header_uni.shift();
+        header_multi.shift();
+        header_multi.shift();
+      }
+      header_uni.splice(1, 2, "95% CI");
+      header_multi.splice(1, 2, "95% CI");
+      self.fillDataHeaders(header_uni, tr, tr_label, "Univariate", forestPlotter_uni);
+      tr.append("td").style("width", "2px");
+      tr_label.append("td").style("width", "2px");
+      self.fillDataHeaders(header_multi, tr, tr_label, "Multivariable-adjusted", forestPlotter_multi);
+    }
+    self.independentTws = Object.keys(result.coefficients_uni.terms).map((tid) => self.getIndependentInput(tid).term);
+    let varcount = 0, rowcolor;
+    for (const tid in result.coefficients_uni.terms) {
+      const termdata = result.coefficients_uni.terms[tid];
+      const termdata_multi = result.coefficients_multi.terms[tid];
+      const tw = self.getIndependentInput(tid).term;
+      rowcolor = ++varcount % 2 ? "#eee" : "none";
+      let tr = table.append("tr").style("background", rowcolor);
+      const termNameTd = tr.append("td").style("padding", "8px");
+      fillCoefficientTermname(tw, termNameTd);
+      if (termdata.fields) {
+        const cols = termdata.fields;
+        const cols_multi = termdata_multi.fields;
+        {
+          const td = tr.append("td").style("padding", "8px");
+          fillColumn2coefficientsTable(td, tw);
+        }
+        if (self.config.regressionType == "cox") {
+          cols.shift();
+          cols.shift();
+          cols_multi.shift();
+          cols_multi.shift();
+        }
+        forestPlotter_uni(tr.append("td"), cols);
+        self.fillCoefDataCols({ tr, cols, tw, isUnivariate: true });
+        tr.append("td").style("width", "2px");
+        forestPlotter_multi(tr.append("td"), cols_multi);
+        self.fillCoefDataCols({ tr, cols: cols_multi, tw });
+      } else if (termdata.categories) {
+        const orderedCategories = [];
+        const input = self.getIndependentInput(tid);
+        if (input.orderedLabels) {
+          for (const k of input.orderedLabels) {
+            if (termdata.categories[k]) orderedCategories.push(k);
+          }
+        }
+        for (const k in termdata.categories) {
+          if (!orderedCategories.includes(k)) orderedCategories.push(k);
+        }
+        termNameTd.attr("rowspan", orderedCategories.length).style("vertical-align", "top");
+        let isfirst = true;
+        for (const k of orderedCategories) {
+          if (!isfirst) {
+            tr = table.append("tr").style("background", rowcolor);
+          }
+          const cols = termdata.categories[k];
+          const cols_multi = termdata_multi.categories[k];
+          const td = tr.append("td").style("padding", "8px");
+          fillColumn2coefficientsTable(td, tw, k);
+          if (self.config.regressionType == "cox") {
+            if (tw.q.mode == "spline") {
+              cols.shift();
+              cols.shift();
+              cols_multi.shift();
+              cols_multi.shift();
+            } else {
+              const [samplesize_ref, samplesize_c] = cols.shift().split("/");
+              const [eventcnt_ref, eventcnt_c] = cols.shift().split("/");
+              if (isfirst) {
+                const variableBottomDiv = termNameTd.select(".sjpcb-coef-variable-bottom");
+                variableBottomDiv.style("align-items", "baseline");
+                const refGrpDiv = variableBottomDiv.selectAll("div").filter((d, i) => i === 1);
+                refGrpDiv.append("div").html(`n=${samplesize_ref}<br>events=${eventcnt_ref}`);
+              }
+              td.append("div").style("font-size", ".8em").html(`n=${samplesize_c}<br>events=${eventcnt_c}`);
+              cols_multi.shift();
+              cols_multi.shift();
+            }
+          }
+          forestPlotter_uni(tr.append("td"), cols);
+          self.fillCoefDataCols({ tr, cols, tw, categoryKey: k, isUnivariate: true });
+          tr.append("td").style("width", "2px");
+          forestPlotter_multi(tr.append("td"), cols_multi);
+          self.fillCoefDataCols({ tr, cols: cols_multi, tw, categoryKey: k });
+          isfirst = false;
+        }
+      } else {
+        tr.append("td").text("ERROR: no .fields[] or .categories{}");
+      }
+    }
+  };
+  self.fillDataHeaders = (header, tr, tr_label, label, forestPlotter) => {
+    const startColN = tr.selectAll("td").size();
+    const forestTd = tr.append("td");
+    if (forestPlotter) {
+      forestTd.style("vertical-align", "bottom");
+      forestPlotter(forestTd);
+    }
+    header.forEach((h, i, arr) => {
+      if (i === 0) {
+        const est = h;
+        const estTd = tr.append("td").style("padding", "8px").text(est);
+        const estInfo = estTd.append("sup").style("cursor", "default").html("&nbsp;&#9432;");
+        estInfo.on("mouseover", (event) => {
+          const tip = self.dom.tip.clear();
+          tip.d.append("div").text("Hover over each value to view explanation of the result");
+          tip.showunder(event.target);
+        });
+        estInfo.on("mouseout", () => self.dom.tip.hide());
+      } else {
+        const td = tr.append("td").text(h).style("padding", "8px");
+        if (i === arr.length - 1) td.style("font-style", "italic");
+      }
+    });
+    if (tr_label) {
+      const endColN = tr.selectAll("td").size();
+      tr_label.append("td").attr("colspan", endColN - startColN).style("padding", "0px 8px").style("text-align", "center").append("div").text(label).style("border-bottom", "1px solid").style("padding", "5px");
+    }
+  };
+  self.fillCoefDataCols = (arg) => {
+    const { tr, cols, tw } = arg;
+    const est = cols.shift();
+    const estSpan = tr.append("td").style("padding", "8px").style("cursor", "default").append("span").text(est);
+    estSpan.on("mouseover", (event) => {
+      if (tw && tw.q.mode == "spline") return;
+      const tip = self.dom.tip.clear();
+      let estimateMsg = self.getEstimateMsg(Object.assign({ est: Number(est) }, arg));
+      if (tw) {
+        const pvalue = Number(cols[cols.length - 1]);
+        estimateMsg += `<br><br><span style="font-style: italic">This association is ${pvalue < 0.05 ? "statistically significant (P < 0.05)" : "not statistically significant (P \u2265 0.05)</span>"}.`;
+      }
+      tip.d.append("div").style("max-width", "500px").html(estimateMsg);
+      tip.showunder(event.target);
+    });
+    estSpan.on("mouseout", () => self.dom.tip.hide());
+    tr.append("td").html(`${cols.shift()} &ndash; ${cols.shift()}`).style("padding", "8px");
+    for (const v of cols) tr.append("td").text(v).style("padding", "8px");
+  };
+  self.getEstimateMsg = (arg) => getEstimateMsg(
+    Object.assign(
+      {
+        regressionType: self.config.regressionType,
+        outcomeTw: self.config.outcome,
+        independentTws: self.independentTws,
+        termdbConfig: self.app.vocabApi.termdbConfig,
+        getIndependentInput: (tid) => self.getIndependentInput(tid)
+      },
+      arg
+    )
+  );
+  self.mayshow_coxDisclaimer = () => {
+    const disclaimer = self.app.vocabApi.termdbConfig.regression?.settings?.coxDisclaimer;
+    if (disclaimer && self.config.regressionType == "cox") {
+      self.dom.oneSetResultDiv.append("div").style("white-space", "wrap").attr("data-testid", "sjpp-regression-result-coxDisclaimer").style("margin", "20px 0px 20px 10px").style("font-size", ".8em").style("text-align", "left").text(disclaimer);
+    }
+  };
+  self.mayshow_totalSnpEffect = (result) => {
+    if (!result.totalSnpEffect) return;
+    const div = self.newDiv(result.totalSnpEffect.label);
+    const table = div.append("table").style("border-spacing", "0px");
+    {
+      const tr2 = table.append("tr").style("opacity", graytextopacity);
+      for (const v of result.totalSnpEffect.header) {
+        tr2.append("td").text(v).style("padding", "8px");
+      }
+    }
+    const tr = table.append("tr").style("background", "#eee");
+    for (const v of result.totalSnpEffect.lst) {
+      tr.append("td").text(v).style("padding", "8px");
+    }
+    const snp = self.getIndependentInput(result.totalSnpEffect.snp).term;
+    const interactions = result.totalSnpEffect.interactions.map((interaction) => {
+      return {
+        t1: self.getIndependentInput(interaction.term1).term,
+        t2: self.getIndependentInput(interaction.term2).term
+      };
+    });
+    const bottomInfo = `Total: total effect of removing the snp (${snp.term.name}) and its interactions (${interactions.map((interaction) => interaction.t1.term.name + " : " + interaction.t2.term.name).join(" ; ")}) from the model`;
+    div.append("div").style("margin", "20px 0px 20px 10px").style("font-size", ".8em").style("text-align", "left").style("color", graytextcolor).text(bottomInfo);
+  };
+  self.mayshow_type3 = (result) => {
+    if (!result.type3 || self.app.vocabApi.termdbConfig.regression?.settings?.hideType3) return;
+    const div = self.newDiv(result.type3.label);
+    const table = div.append("table").style("border-spacing", "0px");
+    {
+      const tr = table.append("tr").style("opacity", graytextopacity);
+      for (const v of result.type3.header) {
+        tr.append("td").text(v).style("padding", "8px");
+      }
+    }
+    if (self.config.regressionType != "cox") {
+      const tr = table.append("tr").style("background", "#eee");
+      for (const v of result.type3.intercept) {
+        tr.append("td").text(v).style("padding", "8px");
+      }
+    }
+    let rowcount = self.config.regressionType == "cox" ? 1 : 0;
+    for (const tid in result.type3.terms) {
+      const termdata = result.type3.terms[tid];
+      const tw = self.getIndependentInput(tid).term;
+      const tr = table.append("tr").style("background", rowcount++ % 2 ? "#eee" : "none");
+      const termNameTd = tr.append("td").style("padding", "8px");
+      fillTdName(termNameTd, tw.term.name);
+      for (const v of termdata) {
+        tr.append("td").text(v).style("padding", "8px");
+      }
+    }
+    for (const row of result.type3.interactions) {
+      const tr = table.append("tr").style("background", rowcount++ % 2 ? "#eee" : "none");
+      const t1 = self.getIndependentInput(row.term1).term;
+      const t2 = self.getIndependentInput(row.term2).term;
+      const td = tr.append("td").style("padding", "8px");
+      fillTdName(td.append("div"), t1.term.name + " : ");
+      fillTdName(td.append("div"), t2.term.name);
+      for (const v of row.lst) {
+        tr.append("td").text(v).style("padding", "8px");
+      }
+    }
+  };
+  self.mayshow_nonlinearity = (result) => {
+    if (!result.nonlinearity) return;
+    const div = self.newDiv(result.nonlinearity.label);
+    const table = div.append("table").style("border-spacing", "0px");
+    {
+      const tr = table.append("tr").style("opacity", graytextopacity);
+      for (const v of result.nonlinearity.header) {
+        tr.append("td").text(v).style("padding", "8px");
+      }
+    }
+    let rowcount = 0;
+    for (const tid in result.nonlinearity.terms) {
+      const tw = self.getIndependentInput(tid).term;
+      const tr = table.append("tr").style("background", rowcount++ % 2 ? "none" : "#eee");
+      fillTdName(tr.append("td").style("padding", "8px"), tw.term.name);
+      for (const v of result.nonlinearity.terms[tid]) {
+        tr.append("td").text(v).style("padding", "8px");
+      }
+    }
+    div.append("div").style("margin", "10px 0px 0px 0px").style("font-size", ".8em").style("color", graytextcolor).text(
+      "Comparison of the cubic spline fit against a linear fit of the same variable. A small p-value indicates the effect of the variable departs from linearity."
+    );
+  };
+  self.mayshow_tests = (result) => {
+    if (!result.tests || self.app.vocabApi.termdbConfig.regression?.settings?.hideTests) return;
+    const div = self.newDiv(result.tests.label);
+    const table = div.append("table").style("border-spacing", "0px");
+    const header = table.append("tr").style("opacity", graytextopacity);
+    for (const cell of result.tests.header) {
+      header.append("td").text(cell).style("padding", "8px");
+    }
+    let rowcount = 0;
+    for (const row of result.tests.rows) {
+      const tr = table.append("tr").style("background", rowcount++ % 2 ? "none" : "#eee");
+      for (const cell of row) {
+        tr.append("td").text(cell).style("padding", "8px");
+      }
+    }
+  };
+  self.mayshow_other = (result) => {
+    if (!result.other) return;
+    const div = self.newDiv(result.other.label);
+    const table = div.append("table").style("border-spacing", "8px");
+    for (let i = 0; i < result.other.header.length; i++) {
+      const tr = table.append("tr");
+      tr.append("td").style("opacity", graytextopacity).text(result.other.header[i]);
+      tr.append("td").text(result.other.rows[i]);
+    }
+  };
+  self.getForestPlotter = (terms, interactions) => {
+    let midIdx, CIlow, CIhigh, axislab, baselineValue, capMin, capMax;
+    if (self.config.regressionType == "linear") {
+      midIdx = 0;
+      CIlow = 1;
+      CIhigh = 2;
+      axislab = "Beta value";
+      baselineValue = 0;
+      capMin = null;
+      capMax = null;
+    } else if (self.config.regressionType == "logistic") {
+      midIdx = 0;
+      CIlow = 1;
+      CIhigh = 2;
+      axislab = "Odds ratio";
+      baselineValue = 1;
+      capMin = 0.1;
+      capMax = 10;
+    } else if (self.config.regressionType == "cox") {
+      midIdx = 0;
+      CIlow = 1;
+      CIhigh = 2;
+      axislab = "Hazard ratio";
+      baselineValue = 1;
+      capMin = 0.1;
+      capMax = 10;
+    } else {
+      throw "unknown regressionType";
+    }
+    const values = [];
+    for (const tid in terms) {
+      const d = terms[tid];
+      if (d.fields) {
+        numbers2array(d.fields);
+      } else {
+        for (const k in d.categories) {
+          numbers2array(d.categories[k]);
+        }
+      }
+    }
+    for (const i of interactions) {
+      for (const k of i.categories) {
+        numbers2array(k.lst);
+      }
+    }
+    if (values.length == 0) {
+      return () => {
+      };
+    }
+    values.sort((a, b) => a - b);
+    if (capMin == null) {
+      capMin = values[0];
+      capMax = values[values.length - 1];
+    }
+    const width = 180;
+    const height = 20;
+    const xleftpad = 10, xrightpad = 10;
+    const scale = get_scale(values);
+    return (td, lst) => {
+      if (!scale) {
+        return;
+      }
+      const fontsize = 12;
+      const svg = td.append("svg").attr("width", width + xleftpad + xrightpad).attr("height", lst ? height : height + fontsize);
+      const g = svg.append("g").attr("transform", `translate(${xleftpad},${lst ? 0 : height + fontsize - 1})`);
+      if (!lst) {
+        const tickFormat = self.config.regressionType == "logistic" ? ".1r" : void 0;
+        const axis = axisTop(scale).ticks(4, tickFormat);
+        axisstyle({
+          axis: g.call(axis),
+          color: forestcolor,
+          showline: true
+        });
+        g.append("text").attr("fill", forestcolor).text(axislab).attr("x", width / 2).attr("y", -height);
+        return;
+      }
+      {
+        const x = scale(baselineValue);
+        g.append("line").attr("x1", x).attr("y1", 0).attr("x2", x).attr("y2", height).attr("stroke", "#ccc");
+      }
+      const mid = Number(lst[midIdx]), cilow = Number(lst[CIlow]), cihigh = Number(lst[CIhigh]);
+      if (Number.isNaN(mid)) {
+        return;
+      }
+      g.append("circle").attr("cx", scale(Math.min(Math.max(mid, capMin), capMax))).attr("cy", height / 2).attr("r", 3).attr("fill", forestcolor);
+      if (Number.isNaN(cilow) || Number.isNaN(cihigh)) {
+        return;
+      }
+      g.append("line").attr("x1", scale(Math.min(Math.max(cilow, capMin), capMax))).attr("y1", height / 2).attr("x2", scale(Math.min(Math.max(cihigh, capMin), capMax))).attr("y2", height / 2).attr("stroke", forestcolor);
+    };
+    function numbers2array(_lst) {
+      const lst = self.config.regressionType == "cox" ? _lst.slice(2) : _lst;
+      const m = Number(lst[midIdx]);
+      if (!Number.isNaN(m)) values.push(m);
+      const l = Number(lst[CIlow]), h = Number(lst[CIhigh]);
+      if (!Number.isNaN(l) && !Number.isNaN(h)) {
+        values.push(l);
+        values.push(h);
+      }
+    }
+    function get_scale(values2) {
+      if (self.config.regressionType == "logistic") {
+        let i = 0;
+        while (values2[i] <= 0) {
+          i++;
+        }
+        if (i >= values2.length || values2[i] <= 0) {
+          return;
+        }
+        const min = values2[i];
+        const max = values2[values2.length - 1];
+        return log().domain([Math.max(min, capMin), Math.min(max, capMax)]).range([0, width]).nice();
+      }
+      if (self.config.regressionType == "linear" || self.config.regressionType == "cox") {
+        return linear().domain([Math.max(values2[0], capMin), Math.min(values2[values2.length - 1], capMax)]).range([0, width]);
+      }
+      throw "unknown type";
+    }
+  };
+}
+function fillTdName(td, name) {
+  if (name.length < 40) {
+    td.text(name);
+  } else {
+    td.text(name.substring(0, 35) + " ...").attr("aria-label", name);
+  }
+}
+function fillCoefficientTermname(tw, td) {
+  fillTdName(td, tw.term.name || tw.term.id || "");
+  const hasRefGrp = "refGrp" in tw && tw.refGrp != refGrp_NA && tw.q.mode != "spline";
+  if (hasRefGrp || tw.effectAllele) {
+    const bottomDiv = td.append("div").attr("class", "sjpcb-coef-variable-bottom").style("display", "flex").style("align-items", "center").style("margin-top", "2px").style("font-size", ".8em");
+    let label;
+    if (hasRefGrp) {
+      label = tw.term.values && tw.term.values[tw.refGrp] ? tw.term.values[tw.refGrp].label : tw.refGrp;
+    } else {
+      label = tw.effectAllele;
+    }
+    bottomDiv.append("div").style("padding", "1px 5px").style("border", "1px solid #aaa").style("border-radius", "10px").style("font-size", ".7em").text(hasRefGrp ? "REF" : "EFFECT ALLELE");
+    bottomDiv.append("div").style("padding", "1px 3px").text(label);
+  }
+  const toUnit = tw.term.valueConversion?.toUnit;
+  if (toUnit && (tw.q.mode == "continuous" || tw.q.mode == "spline")) {
+    td.append("div").style("margin-top", "2px").style("font-size", ".7em").style("opacity", 0.6).text(`per ${toUnit}`);
+  }
+}
+function make_mds3_variants(tw, resultLst, regressionType) {
+  const mlst = [];
+  for (const snp of tw.term.snps) {
+    const m = {
+      chr: snp.chr,
+      pos: snp.pos,
+      ssm_id: snp.snpid
+      // needed for highlighting dot
+    };
+    mlst.push(m);
+    const effAle = tw.q.snp2effAle[snp.snpid];
+    const m2 = snp.mlst.find((i) => i.alt == effAle);
+    if (m2) {
+      Object.assign(m, m2);
+    } else {
+      Object.assign(m, snp.mlst[0]);
+    }
+    m.regressionPvalue = "NA";
+    m.mlpv = 0;
+    const thisresult = resultLst.find((i) => i.id == snp.snpid);
+    if (!thisresult) {
+      m.regressionResult = {
+        data: {
+          err: ["No result for this variant at " + snp.snpid]
+        }
+      };
+      continue;
+    }
+    m.regressionResult = thisresult;
+    const d = thisresult.data;
+    if (!d) throw ".data{} missing";
+    if (d.type3) {
+      const v = getSnpPvalueFromRegressionResults(d, snp.snpid);
+      if (v == void 0) {
+      } else {
+        m.regressionPvalue = v;
+        m.mlpv = -Math.log10(v);
+      }
+      if (!d.coefficients || !d.coefficients.terms) throw ".data.coefficients.terms{} missing";
+      const r = d.coefficients.terms[snp.snpid];
+      if (!r) throw "snp missing from data.coefficients.terms{}";
+      if (Array.isArray(r.fields)) {
+        m.regressionEstimate = regressionType == "cox" ? r.fields[2] : r.fields[0];
+      } else if (r.categories) {
+        const lst = [];
+        for (const gt in r.categories) {
+          lst.push(`${gt}:${regressionType == "cox" ? r.categories[gt][2] : r.categories[gt][0]}`);
+        }
+        m.regressionEstimate = " " + lst.join(" ");
+      } else {
+        throw "unknown way to get snp estimates from coefficients table";
+      }
+    } else if (d.fisher) {
+      m.regressionPvalue = d.fisher.pvalue;
+      m.mlpv = -Math.log10(d.fisher.pvalue);
+      m.shape = "filledTriangle";
+    } else if (d.wilcoxon) {
+      m.regressionPvalue = d.wilcoxon.pvalue;
+      m.mlpv = -Math.log10(d.wilcoxon.pvalue);
+      m.shape = "filledTriangle";
+    } else if (d.cuminc) {
+      m.regressionPvalue = d.cuminc.pvalue;
+      m.mlpv = -Math.log10(d.cuminc.pvalue);
+      m.shape = "filledTriangle";
+    } else {
+      m.shape = "emptyCircle";
+    }
+  }
+  return mlst;
+}
+async function createGenomebrowser(self, input, resultLst) {
+  const arg = {
+    holder: self.dom.snplocusBlockDiv,
+    genome: self.parent.genomeObj,
+    chr: input.term.q.chr,
+    start: input.term.q.start,
+    stop: input.term.q.stop,
+    nobox: true,
+    tklst: [],
+    onCoordinateChange: async (rglst) => {
+      for (const t of self.snplocusBlock.tklst) {
+        if (t.type == "mds3") delete t.skewer.hlssmid;
+      }
+      const { chr, start, stop } = rglst[0];
+      const overrideTw = {
+        term: {
+          id: input.term.term.id,
+          type: "snplocus"
+        },
+        q: JSON.parse(JSON.stringify(input.term.q))
+      };
+      overrideTw.q.chr = chr;
+      overrideTw.q.start = start;
+      overrideTw.q.stop = stop;
+      const _2 = await import("./snplocus-TL25OOPE.js");
+      await _2.fillTW(overrideTw, self.vocabApi);
+      self.hasUnsubmittedEdits_nullify_singleuse = true;
+      input.pill.runCallback(overrideTw);
+    }
+  };
+  arg.tklst.push({
+    type: "mds3",
+    // tkt.mds3
+    name: "Variants",
+    skewerModes: [
+      {
+        type: "numeric",
+        byAttribute: "mlpv",
+        // corresponds to the "mlpv" attribute in m{}, can be anything
+        label: "-log10 p-value",
+        inuse: true,
+        tooltipPrintValue: (m) => getMtooltipValues(m, self.config.regressionType)
+      }
+    ],
+    custom_variants: make_mds3_variants(input.term, resultLst, self.config.regressionType),
+    legend: {
+      customShapeLabels: {
+        filledCircle: "common variants analyzed by model-fitting",
+        filledTriangle: "rare variants analyzed by " + (self.config.regressionType == "linear" ? "Wilcoxon rank sum test" : self.config.regressionType == "logistic" ? "Fisher's exact test" : "Cumulative incidence test"),
+        emptyCircle: "monomorphic variants skipped"
+      }
+    },
+    click_snvindel: async (m) => {
+      self.displayResult_oneset(structuredClone(m.regressionResult.data));
+      await mayCheckLD(m, input, self);
+      const result_y = self.dom.oneSetResultDiv.node().getBoundingClientRect().top + window.scrollY;
+      const nav_height = document.querySelector(".sjpp-nav")?.getBoundingClientRect().height || 0;
+      window.scroll({ behavior: "smooth", top: result_y - nav_height });
+    }
+  });
+  first_genetrack_tolist(self.parent.genomeObj, arg.tklst);
+  const _ = await import("./block-CR75JHV3.js");
+  return new _.Block(arg);
+}
+async function updateMds3Tk(self, input, resultLst) {
+  const tk = self.snplocusBlock.tklst.find((i) => i.type == "mds3");
+  tk.custom_variants = make_mds3_variants(input.term, resultLst, self.config.regressionType);
+  const r = self.snplocusBlock.rglst[0];
+  if (r.chr == input.term.q.chr && r.start == input.term.q.start && r.stop == input.term.q.stop) {
+    tk.load();
+  } else {
+    await self.snplocusBlock.jump_1basedcoordinate(input.term.q);
+  }
+  self.snplocusBlock.cloakOff();
+}
+var LDcolor0 = "#2E6594";
+var LDcolor1 = "#ff0000";
+var LDcolorScale = rgb_default(LDcolor0, LDcolor1);
+async function mayCheckLD(m, input, self) {
+  if (!input.term.q.restrictAncestry) {
+    return;
+  }
+  const tk = self.snplocusBlock.tklst.find((i) => i.type == "mds3");
+  if (!tk || !tk.skewer || !tk.skewer.nmg) return;
+  for (const m2 of tk.custom_variants) delete m2.regressionR2;
+  const wait = self.dom.LDresultDiv.append("span").text("Loading LD data...");
+  try {
+    const data = await self.vocabApi.getLDdata(input.term.q.restrictAncestry.name, m);
+    if (data.error) throw data.error;
+    if (data.nodata || !data.lst || data.lst.length == 0) {
+      wait.text("No LD data");
+      tk.skewer.nmg.selectAll(".sja_aa_disk_fill").attr("fill", (m2) => m2.shapeCircle ? "none" : tk.color4disc(m2));
+      return;
+    }
+    tk.skewer.nmg.selectAll(".sja_aa_disk_fill").attr("fill", (m2) => {
+      if (m2.pos == m.pos && m2.ref == m.ref && m2.alt == m.alt) {
+        return LDcolor1;
+      }
+      for (const i of data.lst) {
+        if (i.pos == m2.pos && i.alleles == m2.ref + "." + m2.alt) {
+          m2.regressionR2 = i.r2;
+          return LDcolorScale(i.r2);
+        }
+      }
+      return LDcolorScale(0);
+    });
+    wait.html(input.term.q.restrictAncestry.name + " LD r<sup>2</sup>");
+    showLDlegend(self.dom.LDresultDiv, LDcolorScale);
+  } catch (e) {
+    wait.text("Error: " + (e.message || e));
+  }
+}
+function showLDlegend(div, colorScale) {
+  const colorbardiv = div.append("span").style("margin-left", "10px");
+  const colorlst = [];
+  for (let i = 0; i <= 1; i += 0.1) {
+    colorlst.push(colorScale(i));
+  }
+  const axisheight = 20;
+  const barheight = 15;
+  const xpad = 10;
+  const axiswidth = 150;
+  const domain = colorlst.map((_, i) => i / (colorlst.length - 1));
+  new ColorScale({
+    holder: colorbardiv,
+    domain,
+    topTicks: true,
+    width: xpad * 2 + axiswidth,
+    height: axisheight + barheight,
+    barheight,
+    barwidth: axiswidth,
+    fontSize: 12,
+    colors: colorlst,
+    position: `${xpad},${axisheight}`,
+    tickSize: 6
+  });
+}
+function getMtooltipValues(m, regressionType) {
+  const lst = [{ k: "p-value", v: m.regressionPvalue }];
+  if (m.regressionResult.AFstr) {
+    lst.push({ k: "AF", v: m.regressionResult.AFstr });
+  }
+  if (m.regressionEstimate) {
+    if (regressionType == "linear") lst.push({ k: "beta", v: m.regressionEstimate });
+    else if (regressionType == "logistic") lst.push({ k: "odds ratio", v: m.regressionEstimate });
+    else if (regressionType == "cox") lst.push({ k: "hazard ratio", v: m.regressionEstimate });
+    else throw "unknown regression type";
+  }
+  if (m.regressionR2) {
+    lst.push({ k: "LD r2", v: m.regressionR2 });
+  }
+  return lst;
+}
+function getSnpPvalueFromRegressionResults(d, snpid) {
+  let str;
+  if (d.totalSnpEffect) {
+    str = d.totalSnpEffect.lst[d.totalSnpEffect.lst.length - 1];
+  } else {
+    if (!d.type3.terms) throw ".data{type3:{terms}} missing";
+    if (!d.type3.terms[snpid]) throw snpid + " missing in type3.terms{}";
+    if (!Array.isArray(d.type3.terms[snpid])) throw `type3.terms[${snpid}] not array`;
+    str = d.type3.terms[snpid][d.type3.terms[snpid].length - 1];
+  }
+  const v = Number(str);
+  if (Number.isFinite(v)) {
+    return v;
+  }
+  return void 0;
+}
+function fillColumn2coefficientsTable(div, tw, categoryKey) {
+  if (categoryKey) {
+    div.text(tw && tw.term.values && tw.term.values[categoryKey] ? tw.term.values[categoryKey].label : categoryKey);
+    return;
+  }
+  div.style("opacity", graytextopacity);
+  if ("geneticModel" in tw.q) {
+    const v = tw.q.geneticModel;
+    div.text(v == 0 ? "(additive)" : v == 1 ? "(dominant)" : "(recessive)");
+    return;
+  }
+  if (tw.q.mode) {
+    div.text("(" + tw.q.mode + ")");
+    return;
+  }
+}
+
+export {
+  RegressionResults,
+  showLDlegend
+};
+//# sourceMappingURL=chunk-6JBQLOJW.js.map