heor-agent-mcp 0.1.0

package/dist/knowledge/yaml.js

@@ -0,0 +1,131 @@
+/**
+ * Parse YAML frontmatter from a markdown string.
+ * Returns { frontmatter, body } where body is the content after the frontmatter.
+ * If no frontmatter present, returns { frontmatter: {}, body: content }.
+ */
+export function parseFrontmatter(content) {
+    if (!content.startsWith("---\n") && !content.startsWith("---\r\n")) {
+        return { frontmatter: {}, body: content };
+    }
+    const lines = content.split(/\r?\n/);
+    // Find closing ---
+    let endIdx = -1;
+    for (let i = 1; i < lines.length; i++) {
+        if (lines[i] === "---") {
+            endIdx = i;
+            break;
+        }
+    }
+    if (endIdx === -1)
+        return { frontmatter: {}, body: content };
+    const yamlLines = lines.slice(1, endIdx);
+    const rawBody = lines.slice(endIdx + 1).join("\n");
+    // Strip a single leading newline that stringifyFrontmatter inserts
+    const body = rawBody.startsWith("\n") ? rawBody.slice(1) : rawBody;
+    const frontmatter = parseYamlLines(yamlLines);
+    return { frontmatter, body };
+}
+function parseYamlLines(lines) {
+    const result = {};
+    let i = 0;
+    while (i < lines.length) {
+        const line = lines[i];
+        if (!line.trim() || line.trim().startsWith("#")) {
+            i++;
+            continue;
+        }
+        // Key: value pattern
+        const match = line.match(/^([a-zA-Z_][a-zA-Z0-9_]*):\s*(.*)$/);
+        if (!match) {
+            i++;
+            continue;
+        }
+        const key = match[1];
+        const rawValue = match[2].trim();
+        if (rawValue === "" || rawValue === "|" || rawValue === ">") {
+            // Multi-line or array follows. Check next line indentation.
+            const items = [];
+            let j = i + 1;
+            while (j < lines.length && lines[j].startsWith("  ")) {
+                const itemLine = lines[j].trim();
+                if (itemLine.startsWith("- ")) {
+                    items.push(parseScalar(itemLine.slice(2).trim()));
+                }
+                j++;
+            }
+            if (items.length > 0)
+                result[key] = items;
+            else
+                result[key] = "";
+            i = j;
+        }
+        else {
+            result[key] = parseScalar(rawValue);
+            i++;
+        }
+    }
+    return result;
+}
+function parseScalar(value) {
+    const trimmed = value.trim();
+    if (trimmed === "null" || trimmed === "~" || trimmed === "")
+        return null;
+    if (trimmed === "true")
+        return true;
+    if (trimmed === "false")
+        return false;
+    // Number
+    if (/^-?\d+$/.test(trimmed))
+        return parseInt(trimmed, 10);
+    if (/^-?\d*\.\d+$/.test(trimmed))
+        return parseFloat(trimmed);
+    // Quoted string
+    if ((trimmed.startsWith('"') && trimmed.endsWith('"')) ||
+        (trimmed.startsWith("'") && trimmed.endsWith("'"))) {
+        return trimmed.slice(1, -1);
+    }
+    return trimmed;
+}
+/**
+ * Serialize a Frontmatter object + body into a markdown string with YAML frontmatter.
+ */
+export function stringifyFrontmatter(frontmatter, body) {
+    const lines = ["---"];
+    for (const [key, value] of Object.entries(frontmatter)) {
+        if (value === undefined || value === null)
+            continue;
+        if (Array.isArray(value)) {
+            if (value.length === 0)
+                continue;
+            lines.push(`${key}:`);
+            value.forEach((item) => lines.push(`  - ${stringifyScalar(item)}`));
+        }
+        else if (value instanceof Date) {
+            lines.push(`${key}: ${value.toISOString()}`);
+        }
+        else if (typeof value === "object") {
+            // Skip nested objects for now (v1 limitation)
+            continue;
+        }
+        else {
+            lines.push(`${key}: ${stringifyScalar(value)}`);
+        }
+    }
+    lines.push("---");
+    lines.push("");
+    lines.push(body);
+    return lines.join("\n");
+}
+function stringifyScalar(value) {
+    if (typeof value === "boolean")
+        return value ? "true" : "false";
+    if (typeof value === "number")
+        return String(value);
+    // String — quote if contains special chars
+    const str = String(value);
+    if (/[:#\[\]{}&*!|>'"%@`,\n]/.test(str) || str.includes("  ")) {
+        return `"${str.replace(/"/g, '\\"')}"`;
+    }
+    return str;
+}
+//# sourceMappingURL=yaml.js.map

package/dist/knowledge/yaml.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"yaml.js","sourceRoot":"","sources":["../../src/knowledge/yaml.ts"],"names":[],"mappings":"AAYA;;;;GAIG;AACH,MAAM,UAAU,gBAAgB,CAAC,OAAe;IAI9C,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,OAAO,CAAC,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,SAAS,CAAC,EAAE,CAAC;QACnE,OAAO,EAAE,WAAW,EAAE,EAAE,EAAE,IAAI,EAAE,OAAO,EAAE,CAAC;IAC5C,CAAC;IACD,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IACrC,mBAAmB;IACnB,IAAI,MAAM,GAAG,CAAC,CAAC,CAAC;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,KAAK,EAAE,CAAC;YACvB,MAAM,GAAG,CAAC,CAAC;YACX,MAAM;QACR,CAAC;IACH,CAAC;IACD,IAAI,MAAM,KAAK,CAAC,CAAC;QAAE,OAAO,EAAE,WAAW,EAAE,EAAE,EAAE,IAAI,EAAE,OAAO,EAAE,CAAC;IAE7D,MAAM,SAAS,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACzC,MAAM,OAAO,GAAG,KAAK,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACnD,mEAAmE;IACnE,MAAM,IAAI,GAAG,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC;IACnE,MAAM,WAAW,GAAG,cAAc,CAAC,SAAS,CAAC,CAAC;IAC9C,OAAO,EAAE,WAAW,EAAE,IAAI,EAAE,CAAC;AAC/B,CAAC;AAED,SAAS,cAAc,CAAC,KAAe;IACrC,MAAM,MAAM,GAAgB,EAAE,CAAC;IAC/B,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,OAAO,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;QACxB,MAAM,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACtB,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,IAAI,CAAC,IAAI,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,EAAE,CAAC;YAChD,CAAC,EAAE,CAAC;YACJ,SAAS;QACX,CAAC;QAED,qBAAqB;QACrB,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,oCAAoC,CAAC,CAAC;QAC/D,IAAI,CAAC,KAAK,EAAE,CAAC;YACX,CAAC,EAAE,CAAC;YACJ,SAAS;QACX,CAAC;QACD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACrB,MAAM,QAAQ,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;QAEjC,IAAI,QAAQ,KAAK,EAAE,IAAI,QAAQ,KAAK,GAAG,IAAI,QAAQ,KAAK,GAAG,EAAE,CAAC;YAC5D,4DAA4D;YAC5D,MAAM,KAAK,GAAa,EAAE,CAAC;YAC3B,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACd,OAAO,CAAC,GAAG,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,IAAI,CAAC,EAAE,CAAC;gBACrD,MAAM,QAAQ,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;gBACjC,IAAI,QAAQ,CAAC,UAAU,CAAC,IAAI,CAAC,EAAE,CAAC;oBAC9B,KAAK,CAAC,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAW,CAAC,CAAC;gBAC9D,CAAC;gBACD,CAAC,EAAE,CAAC;YACN,CAAC;YACD,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC;gBAAE,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;;gBACrC,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;YACtB,CAAC,GAAG,CAAC,CAAC;QACR,CAAC;aAAM,CAAC;YACN,MAAM,CAAC,GAAG,CAAC,GAAG,WAAW,CAAC,QAAQ,CAAC,CAAC;YACpC,CAAC,EAAE,CAAC;QACN,CAAC;IACH,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,WAAW,CAAC,KAAa;IAChC,MAAM,OAAO,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC;IAC7B,IAAI,OAAO,KAAK,MAAM,IAAI,OAAO,KAAK,GAAG,IAAI,OAAO,KAAK,EAAE;QAAE,OAAO,IAAI,CAAC;IACzE,IAAI,OAAO,KAAK,MAAM;QAAE,OAAO,IAAI,CAAC;IACpC,IAAI,OAAO,KAAK,OAAO;QAAE,OAAO,KAAK,CAAC;IACtC,SAAS;IACT,IAAI,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC;QAAE,OAAO,QAAQ,CAAC,OAAO,EAAE,EAAE,CAAC,CAAC;IAC1D,IAAI,cAAc,CAAC,IAAI,CAAC,OAAO,CAAC;QAAE,OAAO,UAAU,CAAC,OAAO,CAAC,CAAC;IAC7D,gBAAgB;IAChB,IACE,CAAC,OAAO,CAAC,UAAU,CAAC,GAAG,CAAC,IAAI,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;QAClD,CAAC,OAAO,CAAC,UAAU,CAAC,GAAG,CAAC,IAAI,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,EAClD,CAAC;QACD,OAAO,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC9B,CAAC;IACD,OAAO,OAAO,CAAC;AACjB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,oBAAoB,CAClC,WAAwB,EACxB,IAAY;IAEZ,MAAM,KAAK,GAAa,CAAC,KAAK,CAAC,CAAC;IAChC,KAAK,MAAM,CAAC,GAAG,EAAE,KAAK,CAAC,IAAI,MAAM,CAAC,OAAO,CAAC,WAAW,CAAC,EAAE,CAAC;QACvD,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,IAAI;YAAE,SAAS;QACpD,IAAI,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC;YACzB,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;gBAAE,SAAS;YACjC,KAAK,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YACtB,KAAK,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,eAAe,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QACtE,CAAC;aAAM,IAAI,KAAK,YAAY,IAAI,EAAE,CAAC;YACjC,KAAK,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,KAAK,CAAC,WAAW,EAAE,EAAE,CAAC,CAAC;QAC/C,CAAC;aAAM,IAAI,OAAO,KAAK,KAAK,QAAQ,EAAE,CAAC;YACrC,8CAA8C;YAC9C,SAAS;QACX,CAAC;aAAM,CAAC;YACN,KAAK,CAAC,IAAI,CACR,GAAG,GAAG,KAAK,eAAe,CAAC,KAAkC,CAAC,EAAE,CACjE,CAAC;QACJ,CAAC;IACH,CAAC;IACD,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IAClB,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACf,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACjB,OAAO,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AAC1B,CAAC;AAED,SAAS,eAAe,CAAC,KAAgC;IACvD,IAAI,OAAO,KAAK,KAAK,SAAS;QAAE,OAAO,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC;IAChE,IAAI,OAAO,KAAK,KAAK,QAAQ;QAAE,OAAO,MAAM,CAAC,KAAK,CAAC,CAAC;IACpD,2CAA2C;IAC3C,MAAM,GAAG,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;IAC1B,IAAI,yBAAyB,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;QAC9D,OAAO,IAAI,GAAG,CAAC,OAAO,CAAC,IAAI,EAAE,KAAK,CAAC,GAAG,CAAC;IACzC,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC"}

package/dist/models/distributions.d.ts

@@ -0,0 +1,28 @@
+/**
+ * Pure functions for sampling from statistical distributions (used in PSA Monte Carlo).
+ * All functions use method of moments parameterization: convert mean/variance to distribution params.
+ */
+/**
+ * Simple LCG (Linear Congruential Generator) seeded RNG.
+ * Returns values in [0, 1).
+ */
+export declare function createSeededRng(seed: number): () => number;
+/**
+ * Beta distribution sample using method of moments.
+ * Mean and variance parameterize Beta(alpha, beta).
+ * Output clamped to [0.001, 0.999].
+ */
+export declare function betaSample(mean: number, variance: number, rng: () => number): number;
+/**
+ * Gamma distribution sample using method of moments.
+ * shape = mean^2/variance, scale = variance/mean
+ * Returns positive value.
+ */
+export declare function gammaSample(mean: number, variance: number, rng: () => number): number;
+/**
+ * Log-Normal distribution sample using method of moments.
+ * params on log scale derived from mean/variance.
+ * Returns positive value.
+ */
+export declare function logNormalSample(mean: number, variance: number, rng: () => number): number;
+//# sourceMappingURL=distributions.d.ts.map

package/dist/models/distributions.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"distributions.d.ts","sourceRoot":"","sources":["../../src/models/distributions.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH;;;GAGG;AACH,wBAAgB,eAAe,CAAC,IAAI,EAAE,MAAM,GAAG,MAAM,MAAM,CAO1D;AAcD;;;;GAIG;AACH,wBAAgB,UAAU,CAAC,IAAI,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,MAAM,GAAG,MAAM,CAYpF;AAsCD;;;;GAIG;AACH,wBAAgB,WAAW,CAAC,IAAI,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,MAAM,GAAG,MAAM,CAIrF;AAED;;;;GAIG;AACH,wBAAgB,eAAe,CAAC,IAAI,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,MAAM,GAAG,MAAM,CAQzF"}

package/dist/models/distributions.js

@@ -0,0 +1,99 @@
+/**
+ * Pure functions for sampling from statistical distributions (used in PSA Monte Carlo).
+ * All functions use method of moments parameterization: convert mean/variance to distribution params.
+ */
+/**
+ * Simple LCG (Linear Congruential Generator) seeded RNG.
+ * Returns values in [0, 1).
+ */
+export function createSeededRng(seed) {
+    let state = seed >>> 0; // ensure unsigned 32-bit int
+    return function () {
+        // LCG parameters from Numerical Recipes
+        state = (Math.imul(1664525, state) + 1013904223) >>> 0;
+        return state / 0x100000000;
+    };
+}
+/**
+ * Box-Muller transform: produces a standard normal sample from two uniform samples.
+ */
+function boxMullerNormal(rng) {
+    let u1 = rng();
+    let u2 = rng();
+    // Guard against log(0)
+    while (u1 === 0)
+        u1 = rng();
+    while (u2 === 0)
+        u2 = rng();
+    return Math.sqrt(-2 * Math.log(u1)) * Math.cos(2 * Math.PI * u2);
+}
+/**
+ * Beta distribution sample using method of moments.
+ * Mean and variance parameterize Beta(alpha, beta).
+ * Output clamped to [0.001, 0.999].
+ */
+export function betaSample(mean, variance, rng) {
+    // Method of moments
+    const factor = mean * (1 - mean) / variance - 1;
+    const alpha = mean * factor;
+    const beta = (1 - mean) * factor;
+    // Sample Beta via Gamma ratio: Beta = Gamma(alpha) / (Gamma(alpha) + Gamma(beta))
+    const x = gammaSampleShape(alpha, 1, rng);
+    const y = gammaSampleShape(beta, 1, rng);
+    const result = x / (x + y);
+    return Math.max(0.001, Math.min(0.999, result));
+}
+/**
+ * Gamma distribution sample using Marsaglia-Tsang method.
+ * shape >= 1 case — returns Gamma(shape, scale=1) sample.
+ */
+function gammaSampleShape(shape, _scale, rng) {
+    if (shape < 1) {
+        // Boost: Gamma(shape) = Gamma(shape+1) * U^(1/shape)
+        const u = rng();
+        return gammaSampleShape(shape + 1, 1, rng) * Math.pow(u, 1 / shape);
+    }
+    // Marsaglia-Tsang
+    const d = shape - 1 / 3;
+    const c = 1 / Math.sqrt(9 * d);
+    while (true) {
+        let x;
+        let v;
+        do {
+            x = boxMullerNormal(rng);
+            v = 1 + c * x;
+        } while (v <= 0);
+        v = v * v * v;
+        const u = rng();
+        if (u < 1 - 0.0331 * (x * x) * (x * x)) {
+            return d * v;
+        }
+        if (Math.log(u) < 0.5 * x * x + d * (1 - v + Math.log(v))) {
+            return d * v;
+        }
+    }
+}
+/**
+ * Gamma distribution sample using method of moments.
+ * shape = mean^2/variance, scale = variance/mean
+ * Returns positive value.
+ */
+export function gammaSample(mean, variance, rng) {
+    const shape = (mean * mean) / variance;
+    const scale = variance / mean;
+    return gammaSampleShape(shape, 1, rng) * scale;
+}
+/**
+ * Log-Normal distribution sample using method of moments.
+ * params on log scale derived from mean/variance.
+ * Returns positive value.
+ */
+export function logNormalSample(mean, variance, rng) {
+    // Log-scale parameters
+    const sigmaSquared = Math.log(1 + variance / (mean * mean));
+    const mu = Math.log(mean) - sigmaSquared / 2;
+    const sigma = Math.sqrt(sigmaSquared);
+    const z = boxMullerNormal(rng);
+    return Math.exp(mu + sigma * z);
+}
+//# sourceMappingURL=distributions.js.map

package/dist/models/distributions.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"distributions.js","sourceRoot":"","sources":["../../src/models/distributions.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH;;;GAGG;AACH,MAAM,UAAU,eAAe,CAAC,IAAY;IAC1C,IAAI,KAAK,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,6BAA6B;IACrD,OAAO;QACL,wCAAwC;QACxC,KAAK,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC;QACvD,OAAO,KAAK,GAAG,WAAW,CAAC;IAC7B,CAAC,CAAC;AACJ,CAAC;AAED;;GAEG;AACH,SAAS,eAAe,CAAC,GAAiB;IACxC,IAAI,EAAE,GAAG,GAAG,EAAE,CAAC;IACf,IAAI,EAAE,GAAG,GAAG,EAAE,CAAC;IACf,uBAAuB;IACvB,OAAO,EAAE,KAAK,CAAC;QAAE,EAAE,GAAG,GAAG,EAAE,CAAC;IAC5B,OAAO,EAAE,KAAK,CAAC;QAAE,EAAE,GAAG,GAAG,EAAE,CAAC;IAC5B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;AACnE,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,UAAU,CAAC,IAAY,EAAE,QAAgB,EAAE,GAAiB;IAC1E,oBAAoB;IACpB,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC;IAChD,MAAM,KAAK,GAAG,IAAI,GAAG,MAAM,CAAC;IAC5B,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,MAAM,CAAC;IAEjC,kFAAkF;IAClF,MAAM,CAAC,GAAG,gBAAgB,CAAC,KAAK,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAC1C,MAAM,CAAC,GAAG,gBAAgB,CAAC,IAAI,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IACzC,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAE3B,OAAO,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAC;AAClD,CAAC;AAED;;;GAGG;AACH,SAAS,gBAAgB,CAAC,KAAa,EAAE,MAAc,EAAE,GAAiB;IACxE,IAAI,KAAK,GAAG,CAAC,EAAE,CAAC;QACd,qDAAqD;QACrD,MAAM,CAAC,GAAG,GAAG,EAAE,CAAC;QAChB,OAAO,gBAAgB,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC;IACtE,CAAC;IAED,kBAAkB;IAClB,MAAM,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;IACxB,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAE/B,OAAO,IAAI,EAAE,CAAC;QACZ,IAAI,CAAS,CAAC;QACd,IAAI,CAAS,CAAC;QACd,GAAG,CAAC;YACF,CAAC,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;YACzB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE;QAEjB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACd,MAAM,CAAC,GAAG,GAAG,EAAE,CAAC;QAEhB,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YACvC,OAAO,CAAC,GAAG,CAAC,CAAC;QACf,CAAC;QAED,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;YAC1D,OAAO,CAAC,GAAG,CAAC,CAAC;QACf,CAAC;IACH,CAAC;AACH,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,WAAW,CAAC,IAAY,EAAE,QAAgB,EAAE,GAAiB;IAC3E,MAAM,KAAK,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,QAAQ,CAAC;IACvC,MAAM,KAAK,GAAG,QAAQ,GAAG,IAAI,CAAC;IAC9B,OAAO,gBAAgB,CAAC,KAAK,EAAE,CAAC,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC;AACjD,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,eAAe,CAAC,IAAY,EAAE,QAAgB,EAAE,GAAiB;IAC/E,uBAAuB;IACvB,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC;IAC5D,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,YAAY,GAAG,CAAC,CAAC;IAC7C,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;IAEtC,MAAM,CAAC,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;IAC/B,OAAO,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC;AAClC,CAAC"}

package/dist/models/index.d.ts

@@ -0,0 +1,11 @@
+export { runMarkovModel } from "./markov.js";
+export type { MarkovParams, MarkovState, TransitionMatrix, MarkovRunResult } from "./markov.js";
+export { runPartSA } from "./partsa.js";
+export type { PartSAParams, PartSAResult, SurvivalParams } from "./partsa.js";
+export { runPSA } from "./psa.js";
+export type { PSAParams, PSAResult, PSAIteration } from "./psa.js";
+export { runOWSA, buildDefaultOWSAParameters } from "./owsa.js";
+export type { OWSAParameter, OWSAResult } from "./owsa.js";
+export { betaSample, gammaSample, logNormalSample, createSeededRng } from "./distributions.js";
+export { buildMarkovParamsFromCE, runMarkovAndComputeICER } from "./modelUtils.js";
+//# sourceMappingURL=index.d.ts.map

package/dist/models/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/models/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAC7C,YAAY,EAAE,YAAY,EAAE,WAAW,EAAE,gBAAgB,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAEhG,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,YAAY,EAAE,YAAY,EAAE,YAAY,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAE9E,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAClC,YAAY,EAAE,SAAS,EAAE,SAAS,EAAE,YAAY,EAAE,MAAM,UAAU,CAAC;AAEnE,OAAO,EAAE,OAAO,EAAE,0BAA0B,EAAE,MAAM,WAAW,CAAC;AAChE,YAAY,EAAE,aAAa,EAAE,UAAU,EAAE,MAAM,WAAW,CAAC;AAE3D,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAE/F,OAAO,EAAE,uBAAuB,EAAE,uBAAuB,EAAE,MAAM,iBAAiB,CAAC"}

package/dist/models/index.js

@@ -0,0 +1,7 @@
+export { runMarkovModel } from "./markov.js";
+export { runPartSA } from "./partsa.js";
+export { runPSA } from "./psa.js";
+export { runOWSA, buildDefaultOWSAParameters } from "./owsa.js";
+export { betaSample, gammaSample, logNormalSample, createSeededRng } from "./distributions.js";
+export { buildMarkovParamsFromCE, runMarkovAndComputeICER } from "./modelUtils.js";
+//# sourceMappingURL=index.js.map

package/dist/models/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/models/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAG7C,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAGxC,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAGlC,OAAO,EAAE,OAAO,EAAE,0BAA0B,EAAE,MAAM,WAAW,CAAC;AAGhE,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAE/F,OAAO,EAAE,uBAAuB,EAAE,uBAAuB,EAAE,MAAM,iBAAiB,CAAC"}

package/dist/models/markov.d.ts

@@ -0,0 +1,42 @@
+/**
+ * Multi-state Markov cohort model with half-cycle correction and discounting.
+ */
+export interface MarkovState {
+    name: string;
+    utility: number;
+    cost_annual: number;
+}
+export interface TransitionMatrix {
+    [fromState: string]: {
+        [toState: string]: number;
+    };
+}
+export interface MarkovParams {
+    states: MarkovState[];
+    transition_matrix_intervention: TransitionMatrix;
+    transition_matrix_comparator: TransitionMatrix;
+    initial_cohort: {
+        [stateName: string]: number;
+    };
+    cycle_length_years: number;
+    n_cycles: number;
+    discount_rate_costs: number;
+    discount_rate_outcomes: number;
+}
+export interface MarkovRunResult {
+    total_cost: number;
+    total_qaly: number;
+    total_lys: number;
+    state_trace: Array<{
+        cycle: number;
+        distribution: Record<string, number>;
+    }>;
+}
+/**
+ * Run Markov model for both intervention and comparator arms.
+ */
+export declare function runMarkovModel(params: MarkovParams): {
+    intervention: MarkovRunResult;
+    comparator: MarkovRunResult;
+};
+//# sourceMappingURL=markov.d.ts.map

package/dist/models/markov.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"markov.d.ts","sourceRoot":"","sources":["../../src/models/markov.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,MAAM,WAAW,WAAW;IAC1B,IAAI,EAAE,MAAM,CAAC;IACb,OAAO,EAAE,MAAM,CAAC;IAChB,WAAW,EAAE,MAAM,CAAC;CACrB;AAED,MAAM,WAAW,gBAAgB;IAE/B,CAAC,SAAS,EAAE,MAAM,GAAG;QAAE,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM,CAAA;KAAE,CAAC;CACpD;AAED,MAAM,WAAW,YAAY;IAC3B,MAAM,EAAE,WAAW,EAAE,CAAC;IACtB,8BAA8B,EAAE,gBAAgB,CAAC;IACjD,4BAA4B,EAAE,gBAAgB,CAAC;IAC/C,cAAc,EAAE;QAAE,CAAC,SAAS,EAAE,MAAM,GAAG,MAAM,CAAA;KAAE,CAAC;IAChD,kBAAkB,EAAE,MAAM,CAAC;IAC3B,QAAQ,EAAE,MAAM,CAAC;IACjB,mBAAmB,EAAE,MAAM,CAAC;IAC5B,sBAAsB,EAAE,MAAM,CAAC;CAChC;AAED,MAAM,WAAW,eAAe;IAC9B,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,SAAS,EAAE,MAAM,CAAC;IAClB,WAAW,EAAE,KAAK,CAAC;QAAE,KAAK,EAAE,MAAM,CAAC;QAAC,YAAY,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAA;KAAE,CAAC,CAAC;CAC7E;AA2FD;;GAEG;AACH,wBAAgB,cAAc,CAC5B,MAAM,EAAE,YAAY,GACnB;IAAE,YAAY,EAAE,eAAe,CAAC;IAAC,UAAU,EAAE,eAAe,CAAA;CAAE,CAIhE"}

package/dist/models/markov.js

@@ -0,0 +1,82 @@
+/**
+ * Multi-state Markov cohort model with half-cycle correction and discounting.
+ */
+/**
+ * Run one arm of a Markov cohort simulation.
+ * Uses half-cycle correction: add half a cycle's value at cycle 0 and final cycle.
+ */
+function runArm(params, transitionMatrix) {
+    const { states, initial_cohort, cycle_length_years, n_cycles, discount_rate_costs, discount_rate_outcomes } = params;
+    const stateNames = states.map(s => s.name);
+    // Initialize cohort distribution
+    let cohort = {};
+    for (const name of stateNames) {
+        cohort[name] = initial_cohort[name] ?? 0;
+    }
+    let total_cost = 0;
+    let total_qaly = 0;
+    let total_lys = 0;
+    const state_trace = [];
+    // Half-cycle correction at cycle 0 (t=0)
+    const discountCost0 = 1 / Math.pow(1 + discount_rate_costs, 0);
+    const discountOutcome0 = 1 / Math.pow(1 + discount_rate_outcomes, 0);
+    for (const state of states) {
+        const prop = cohort[state.name] ?? 0;
+        // Half cycle at start
+        total_cost += prop * state.cost_annual * cycle_length_years * 0.5 * discountCost0;
+        total_qaly += prop * state.utility * cycle_length_years * 0.5 * discountOutcome0;
+        // Life years: any non-dead state contributes (assume "Dead" states have utility ~0 — but we track LYs by cohort alive)
+        if (state.utility > 0) {
+            total_lys += prop * cycle_length_years * 0.5 * discountOutcome0;
+        }
+    }
+    state_trace.push({ cycle: 0, distribution: { ...cohort } });
+    for (let cycle = 1; cycle <= n_cycles; cycle++) {
+        // Apply transition matrix: new cohort
+        const newCohort = {};
+        for (const name of stateNames) {
+            newCohort[name] = 0;
+        }
+        for (const fromName of stateNames) {
+            const fromProp = cohort[fromName] ?? 0;
+            if (fromProp <= 0)
+                continue;
+            const row = transitionMatrix[fromName];
+            if (!row) {
+                // No transitions defined: stay in same state
+                newCohort[fromName] = (newCohort[fromName] ?? 0) + fromProp;
+                continue;
+            }
+            for (const toName of stateNames) {
+                const prob = row[toName] ?? 0;
+                newCohort[toName] = (newCohort[toName] ?? 0) + fromProp * prob;
+            }
+        }
+        cohort = newCohort;
+        // Discount factors (time t = cycle in years from start)
+        const t = cycle * cycle_length_years;
+        const discountCost = 1 / Math.pow(1 + discount_rate_costs, t);
+        const discountOutcome = 1 / Math.pow(1 + discount_rate_outcomes, t);
+        const isLastCycle = cycle === n_cycles;
+        const cycleFraction = isLastCycle ? 0.5 : 1.0; // half-cycle at last cycle
+        for (const state of states) {
+            const prop = cohort[state.name] ?? 0;
+            total_cost += prop * state.cost_annual * cycle_length_years * cycleFraction * discountCost;
+            total_qaly += prop * state.utility * cycle_length_years * cycleFraction * discountOutcome;
+            if (state.utility > 0) {
+                total_lys += prop * cycle_length_years * cycleFraction * discountOutcome;
+            }
+        }
+        state_trace.push({ cycle, distribution: { ...cohort } });
+    }
+    return { total_cost, total_qaly, total_lys, state_trace };
+}
+/**
+ * Run Markov model for both intervention and comparator arms.
+ */
+export function runMarkovModel(params) {
+    const intervention = runArm(params, params.transition_matrix_intervention);
+    const comparator = runArm(params, params.transition_matrix_comparator);
+    return { intervention, comparator };
+}
+//# sourceMappingURL=markov.js.map

package/dist/models/markov.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"markov.js","sourceRoot":"","sources":["../../src/models/markov.ts"],"names":[],"mappings":"AAAA;;GAEG;AA+BH;;;GAGG;AACH,SAAS,MAAM,CACb,MAAoB,EACpB,gBAAkC;IAElC,MAAM,EAAE,MAAM,EAAE,cAAc,EAAE,kBAAkB,EAAE,QAAQ,EAAE,mBAAmB,EAAE,sBAAsB,EAAE,GAAG,MAAM,CAAC;IACrH,MAAM,UAAU,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAE3C,iCAAiC;IACjC,IAAI,MAAM,GAA2B,EAAE,CAAC;IACxC,KAAK,MAAM,IAAI,IAAI,UAAU,EAAE,CAAC;QAC9B,MAAM,CAAC,IAAI,CAAC,GAAG,cAAc,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC3C,CAAC;IAED,IAAI,UAAU,GAAG,CAAC,CAAC;IACnB,IAAI,UAAU,GAAG,CAAC,CAAC;IACnB,IAAI,SAAS,GAAG,CAAC,CAAC;IAElB,MAAM,WAAW,GAAmE,EAAE,CAAC;IAEvF,yCAAyC;IACzC,MAAM,aAAa,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,mBAAmB,EAAE,CAAC,CAAC,CAAC;IAC/D,MAAM,gBAAgB,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,sBAAsB,EAAE,CAAC,CAAC,CAAC;IAErE,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACrC,sBAAsB;QACtB,UAAU,IAAI,IAAI,GAAG,KAAK,CAAC,WAAW,GAAG,kBAAkB,GAAG,GAAG,GAAG,aAAa,CAAC;QAClF,UAAU,IAAI,IAAI,GAAG,KAAK,CAAC,OAAO,GAAG,kBAAkB,GAAG,GAAG,GAAG,gBAAgB,CAAC;QACjF,uHAAuH;QACvH,IAAI,KAAK,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC;YACtB,SAAS,IAAI,IAAI,GAAG,kBAAkB,GAAG,GAAG,GAAG,gBAAgB,CAAC;QAClE,CAAC;IACH,CAAC;IAED,WAAW,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,YAAY,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,CAAC;IAE5D,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,IAAI,QAAQ,EAAE,KAAK,EAAE,EAAE,CAAC;QAC/C,sCAAsC;QACtC,MAAM,SAAS,GAA2B,EAAE,CAAC;QAC7C,KAAK,MAAM,IAAI,IAAI,UAAU,EAAE,CAAC;YAC9B,SAAS,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACtB,CAAC;QAED,KAAK,MAAM,QAAQ,IAAI,UAAU,EAAE,CAAC;YAClC,MAAM,QAAQ,GAAG,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;YACvC,IAAI,QAAQ,IAAI,CAAC;gBAAE,SAAS;YAE5B,MAAM,GAAG,GAAG,gBAAgB,CAAC,QAAQ,CAAC,CAAC;YACvC,IAAI,CAAC,GAAG,EAAE,CAAC;gBACT,6CAA6C;gBAC7C,SAAS,CAAC,QAAQ,CAAC,GAAG,CAAC,SAAS,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,GAAG,QAAQ,CAAC;gBAC5D,SAAS;YACX,CAAC;YAED,KAAK,MAAM,MAAM,IAAI,UAAU,EAAE,CAAC;gBAChC,MAAM,IAAI,GAAG,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;gBAC9B,SAAS,CAAC,MAAM,CAAC,GAAG,CAAC,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC;YACjE,CAAC;QACH,CAAC;QAED,MAAM,GAAG,SAAS,CAAC;QAEnB,wDAAwD;QACxD,MAAM,CAAC,GAAG,KAAK,GAAG,kBAAkB,CAAC;QACrC,MAAM,YAAY,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,mBAAmB,EAAE,CAAC,CAAC,CAAC;QAC9D,MAAM,eAAe,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,sBAAsB,EAAE,CAAC,CAAC,CAAC;QAEpE,MAAM,WAAW,GAAG,KAAK,KAAK,QAAQ,CAAC;QACvC,MAAM,aAAa,GAAG,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,2BAA2B;QAE1E,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;YAC3B,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YACrC,UAAU,IAAI,IAAI,GAAG,KAAK,CAAC,WAAW,GAAG,kBAAkB,GAAG,aAAa,GAAG,YAAY,CAAC;YAC3F,UAAU,IAAI,IAAI,GAAG,KAAK,CAAC,OAAO,GAAG,kBAAkB,GAAG,aAAa,GAAG,eAAe,CAAC;YAC1F,IAAI,KAAK,CAAC,OAAO,GAAG,CAAC,EAAE,CAAC;gBACtB,SAAS,IAAI,IAAI,GAAG,kBAAkB,GAAG,aAAa,GAAG,eAAe,CAAC;YAC3E,CAAC;QACH,CAAC;QAED,WAAW,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,YAAY,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,CAAC;IAC3D,CAAC;IAED,OAAO,EAAE,UAAU,EAAE,UAAU,EAAE,SAAS,EAAE,WAAW,EAAE,CAAC;AAC5D,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,cAAc,CAC5B,MAAoB;IAEpB,MAAM,YAAY,GAAG,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,8BAA8B,CAAC,CAAC;IAC3E,MAAM,UAAU,GAAG,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,4BAA4B,CAAC,CAAC;IACvE,OAAO,EAAE,YAAY,EAAE,UAAU,EAAE,CAAC;AACtC,CAAC"}

package/dist/models/modelUtils.d.ts

@@ -0,0 +1,29 @@
+/**
+ * Shared utilities for building Markov params from CEModelParams
+ * and computing ICERs.
+ */
+import type { CEModelParams } from "../providers/types.js";
+import type { MarkovParams } from "./markov.js";
+/**
+ * Build MarkovParams from CEModelParams using a 2-state model
+ * (On-Treatment, Off-Treatment).
+ */
+export declare function buildMarkovParamsFromCE(params: CEModelParams): MarkovParams;
+/**
+ * Run Markov model from CEModelParams and return delta_cost, delta_qaly, ICER.
+ *
+ * Runs the model twice — once with intervention state costs and once with
+ * comparator state costs — so each arm uses the correct drug cost.
+ */
+export declare function runMarkovAndComputeICER(params: CEModelParams): {
+    delta_cost: number;
+    delta_qaly: number;
+    icer: number;
+    intervention_cost: number;
+    comparator_cost: number;
+    intervention_qaly: number;
+    comparator_qaly: number;
+    intervention_lys: number;
+    comparator_lys: number;
+};
+//# sourceMappingURL=modelUtils.d.ts.map

package/dist/models/modelUtils.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"modelUtils.d.ts","sourceRoot":"","sources":["../../src/models/modelUtils.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAEH,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AAC3D,OAAO,KAAK,EAAE,YAAY,EAAiC,MAAM,aAAa,CAAC;AAY/E;;;GAGG;AACH,wBAAgB,uBAAuB,CAAC,MAAM,EAAE,aAAa,GAAG,YAAY,CAqF3E;AAED;;;;;GAKG;AACH,wBAAgB,uBAAuB,CAAC,MAAM,EAAE,aAAa,GAAG;IAC9D,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,IAAI,EAAE,MAAM,CAAC;IACb,iBAAiB,EAAE,MAAM,CAAC;IAC1B,eAAe,EAAE,MAAM,CAAC;IACxB,iBAAiB,EAAE,MAAM,CAAC;IAC1B,eAAe,EAAE,MAAM,CAAC;IACxB,gBAAgB,EAAE,MAAM,CAAC;IACzB,cAAc,EAAE,MAAM,CAAC;CACxB,CA8DA"}

package/dist/models/modelUtils.js

@@ -0,0 +1,139 @@
+/**
+ * Shared utilities for building Markov params from CEModelParams
+ * and computing ICERs.
+ */
+import { runMarkovModel } from "./markov.js";
+const DISCOUNT_RATE = 0.035;
+function getTimeHorizonYears(horizon) {
+    if (horizon === "lifetime")
+        return 40;
+    if (horizon === "5yr")
+        return 5;
+    if (horizon === "10yr")
+        return 10;
+    return Number(horizon);
+}
+/**
+ * Build MarkovParams from CEModelParams using a 2-state model
+ * (On-Treatment, Off-Treatment).
+ */
+export function buildMarkovParamsFromCE(params) {
+    const years = getTimeHorizonYears(params.time_horizon);
+    const cycle_length_years = 1;
+    const n_cycles = years;
+    // Utility values
+    const utilityOn = params.utility_inputs?.qaly_on_treatment ?? 0.75;
+    const utilityOff = params.utility_inputs?.qaly_comparator ?? 0.7;
+    // Costs
+    const costIntervention = params.cost_inputs.drug_cost_annual + (params.cost_inputs.admin_cost ?? 0);
+    const costComparator = params.cost_inputs.comparator_cost_annual +
+        (params.cost_inputs.admin_cost ?? 0);
+    // Efficacy delta used to derive annual probability of staying on treatment
+    // Higher efficacy → more cycles on treatment → better health
+    const efficacyDelta = Math.max(0, Math.min(0.999, params.clinical_inputs.efficacy_delta));
+    const mortalityReduction = params.clinical_inputs.mortality_reduction ?? 0;
+    // Transition: prob of staying in On-Treatment state for intervention
+    // Base: driven by efficacy (0.4 efficacy → ~0.8 prob of staying on treatment, modulated)
+    const probStayOnIntervention = Math.max(0.05, Math.min(0.95, 0.5 + efficacyDelta * 0.5));
+    // For comparator: lower stay probability (baseline)
+    const baselineProbStayOn = Math.max(0.05, Math.min(0.9, probStayOnIntervention * 0.7));
+    // Mortality reduction affects off-treatment → death transition
+    // (simplified: mortality reduction lowers probability of going to off-treatment)
+    const mortalityEffect = Math.max(0, Math.min(0.3, mortalityReduction * 0.3));
+    const states = [
+        { name: "On-Treatment", utility: utilityOn, cost_annual: costIntervention },
+        { name: "Off-Treatment", utility: utilityOff, cost_annual: 0 },
+    ];
+    const statesComparator = [
+        { name: "On-Treatment", utility: utilityOn, cost_annual: costComparator },
+        { name: "Off-Treatment", utility: utilityOff, cost_annual: 0 },
+    ];
+    const transition_matrix_intervention = {
+        "On-Treatment": {
+            "On-Treatment": Math.min(0.95, probStayOnIntervention + mortalityEffect),
+            "Off-Treatment": Math.max(0.05, 1 - probStayOnIntervention - mortalityEffect),
+        },
+        "Off-Treatment": {
+            "On-Treatment": 0.05,
+            "Off-Treatment": 0.95,
+        },
+    };
+    const transition_matrix_comparator = {
+        "On-Treatment": {
+            "On-Treatment": Math.min(0.9, baselineProbStayOn),
+            "Off-Treatment": Math.max(0.1, 1 - baselineProbStayOn),
+        },
+        "Off-Treatment": {
+            "On-Treatment": 0.05,
+            "Off-Treatment": 0.95,
+        },
+    };
+    return {
+        states,
+        transition_matrix_intervention,
+        transition_matrix_comparator,
+        initial_cohort: { "On-Treatment": 1.0, "Off-Treatment": 0.0 },
+        cycle_length_years,
+        n_cycles,
+        discount_rate_costs: DISCOUNT_RATE,
+        discount_rate_outcomes: DISCOUNT_RATE,
+    };
+}
+/**
+ * Run Markov model from CEModelParams and return delta_cost, delta_qaly, ICER.
+ *
+ * Runs the model twice — once with intervention state costs and once with
+ * comparator state costs — so each arm uses the correct drug cost.
+ */
+export function runMarkovAndComputeICER(params) {
+    const baseParams = buildMarkovParamsFromCE(params);
+    // Utility values (shared across arms)
+    const utilityOn = params.utility_inputs?.qaly_on_treatment ?? 0.75;
+    const utilityOff = params.utility_inputs?.qaly_comparator ?? 0.7;
+    // Costs per arm
+    const costIntervention = params.cost_inputs.drug_cost_annual + (params.cost_inputs.admin_cost ?? 0);
+    const costComparator = params.cost_inputs.comparator_cost_annual +
+        (params.cost_inputs.admin_cost ?? 0);
+    const statesIntervention = [
+        { name: "On-Treatment", utility: utilityOn, cost_annual: costIntervention },
+        { name: "Off-Treatment", utility: utilityOff, cost_annual: 0 },
+    ];
+    const statesComparator = [
+        { name: "On-Treatment", utility: utilityOn, cost_annual: costComparator },
+        { name: "Off-Treatment", utility: utilityOff, cost_annual: 0 },
+    ];
+    // Run intervention arm: use intervention states + intervention transition matrix.
+    // We pass intervention states in both arms' params; we only use the intervention
+    // arm result from this call.
+    const interventionRun = runMarkovModel({
+        ...baseParams,
+        states: statesIntervention,
+    });
+    // Run comparator arm: use comparator states + comparator transition matrix.
+    // We only use the comparator arm result from this call.
+    const comparatorRun = runMarkovModel({
+        ...baseParams,
+        states: statesComparator,
+    });
+    const intervention = interventionRun.intervention;
+    const comparator = comparatorRun.comparator;
+    const delta_cost = intervention.total_cost - comparator.total_cost;
+    const delta_qaly = intervention.total_qaly - comparator.total_qaly;
+    const icer = delta_qaly > 0
+        ? delta_cost / delta_qaly
+        : delta_qaly < 0
+            ? -Infinity
+            : Infinity;
+    return {
+        delta_cost,
+        delta_qaly,
+        icer,
+        intervention_cost: intervention.total_cost,
+        comparator_cost: comparator.total_cost,
+        intervention_qaly: intervention.total_qaly,
+        comparator_qaly: comparator.total_qaly,
+        intervention_lys: intervention.total_lys,
+        comparator_lys: comparator.total_lys,
+    };
+}
+//# sourceMappingURL=modelUtils.js.map

package/dist/models/modelUtils.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"modelUtils.js","sourceRoot":"","sources":["../../src/models/modelUtils.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAIH,OAAO,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAE7C,MAAM,aAAa,GAAG,KAAK,CAAC;AAE5B,SAAS,mBAAmB,CAAC,OAAsC;IACjE,IAAI,OAAO,KAAK,UAAU;QAAE,OAAO,EAAE,CAAC;IACtC,IAAI,OAAO,KAAK,KAAK;QAAE,OAAO,CAAC,CAAC;IAChC,IAAI,OAAO,KAAK,MAAM;QAAE,OAAO,EAAE,CAAC;IAClC,OAAO,MAAM,CAAC,OAAO,CAAC,CAAC;AACzB,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,uBAAuB,CAAC,MAAqB;IAC3D,MAAM,KAAK,GAAG,mBAAmB,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC;IACvD,MAAM,kBAAkB,GAAG,CAAC,CAAC;IAC7B,MAAM,QAAQ,GAAG,KAAK,CAAC;IAEvB,iBAAiB;IACjB,MAAM,SAAS,GAAG,MAAM,CAAC,cAAc,EAAE,iBAAiB,IAAI,IAAI,CAAC;IACnE,MAAM,UAAU,GAAG,MAAM,CAAC,cAAc,EAAE,eAAe,IAAI,GAAG,CAAC;IAEjE,QAAQ;IACR,MAAM,gBAAgB,GACpB,MAAM,CAAC,WAAW,CAAC,gBAAgB,GAAG,CAAC,MAAM,CAAC,WAAW,CAAC,UAAU,IAAI,CAAC,CAAC,CAAC;IAC7E,MAAM,cAAc,GAClB,MAAM,CAAC,WAAW,CAAC,sBAAsB;QACzC,CAAC,MAAM,CAAC,WAAW,CAAC,UAAU,IAAI,CAAC,CAAC,CAAC;IAEvC,2EAA2E;IAC3E,6DAA6D;IAC7D,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAC5B,CAAC,EACD,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,eAAe,CAAC,cAAc,CAAC,CACvD,CAAC;IACF,MAAM,kBAAkB,GAAG,MAAM,CAAC,eAAe,CAAC,mBAAmB,IAAI,CAAC,CAAC;IAE3E,qEAAqE;IACrE,yFAAyF;IACzF,MAAM,sBAAsB,GAAG,IAAI,CAAC,GAAG,CACrC,IAAI,EACJ,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,GAAG,aAAa,GAAG,GAAG,CAAC,CAC1C,CAAC;IACF,oDAAoD;IACpD,MAAM,kBAAkB,GAAG,IAAI,CAAC,GAAG,CACjC,IAAI,EACJ,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,sBAAsB,GAAG,GAAG,CAAC,CAC5C,CAAC;IAEF,+DAA+D;IAC/D,iFAAiF;IACjF,MAAM,eAAe,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,kBAAkB,GAAG,GAAG,CAAC,CAAC,CAAC;IAE7E,MAAM,MAAM,GAAkB;QAC5B,EAAE,IAAI,EAAE,cAAc,EAAE,OAAO,EAAE,SAAS,EAAE,WAAW,EAAE,gBAAgB,EAAE;QAC3E,EAAE,IAAI,EAAE,eAAe,EAAE,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,CAAC,EAAE;KAC/D,CAAC;IAEF,MAAM,gBAAgB,GAAkB;QACtC,EAAE,IAAI,EAAE,cAAc,EAAE,OAAO,EAAE,SAAS,EAAE,WAAW,EAAE,cAAc,EAAE;QACzE,EAAE,IAAI,EAAE,eAAe,EAAE,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,CAAC,EAAE;KAC/D,CAAC;IAEF,MAAM,8BAA8B,GAAqB;QACvD,cAAc,EAAE;YACd,cAAc,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,sBAAsB,GAAG,eAAe,CAAC;YACxE,eAAe,EAAE,IAAI,CAAC,GAAG,CACvB,IAAI,EACJ,CAAC,GAAG,sBAAsB,GAAG,eAAe,CAC7C;SACF;QACD,eAAe,EAAE;YACf,cAAc,EAAE,IAAI;YACpB,eAAe,EAAE,IAAI;SACtB;KACF,CAAC;IAEF,MAAM,4BAA4B,GAAqB;QACrD,cAAc,EAAE;YACd,cAAc,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,kBAAkB,CAAC;YACjD,eAAe,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,kBAAkB,CAAC;SACvD;QACD,eAAe,EAAE;YACf,cAAc,EAAE,IAAI;YACpB,eAAe,EAAE,IAAI;SACtB;KACF,CAAC;IAEF,OAAO;QACL,MAAM;QACN,8BAA8B;QAC9B,4BAA4B;QAC5B,cAAc,EAAE,EAAE,cAAc,EAAE,GAAG,EAAE,eAAe,EAAE,GAAG,EAAE;QAC7D,kBAAkB;QAClB,QAAQ;QACR,mBAAmB,EAAE,aAAa;QAClC,sBAAsB,EAAE,aAAa;KACtC,CAAC;AACJ,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,uBAAuB,CAAC,MAAqB;IAW3D,MAAM,UAAU,GAAG,uBAAuB,CAAC,MAAM,CAAC,CAAC;IAEnD,sCAAsC;IACtC,MAAM,SAAS,GAAG,MAAM,CAAC,cAAc,EAAE,iBAAiB,IAAI,IAAI,CAAC;IACnE,MAAM,UAAU,GAAG,MAAM,CAAC,cAAc,EAAE,eAAe,IAAI,GAAG,CAAC;IAEjE,gBAAgB;IAChB,MAAM,gBAAgB,GACpB,MAAM,CAAC,WAAW,CAAC,gBAAgB,GAAG,CAAC,MAAM,CAAC,WAAW,CAAC,UAAU,IAAI,CAAC,CAAC,CAAC;IAC7E,MAAM,cAAc,GAClB,MAAM,CAAC,WAAW,CAAC,sBAAsB;QACzC,CAAC,MAAM,CAAC,WAAW,CAAC,UAAU,IAAI,CAAC,CAAC,CAAC;IAEvC,MAAM,kBAAkB,GAAkB;QACxC,EAAE,IAAI,EAAE,cAAc,EAAE,OAAO,EAAE,SAAS,EAAE,WAAW,EAAE,gBAAgB,EAAE;QAC3E,EAAE,IAAI,EAAE,eAAe,EAAE,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,CAAC,EAAE;KAC/D,CAAC;IAEF,MAAM,gBAAgB,GAAkB;QACtC,EAAE,IAAI,EAAE,cAAc,EAAE,OAAO,EAAE,SAAS,EAAE,WAAW,EAAE,cAAc,EAAE;QACzE,EAAE,IAAI,EAAE,eAAe,EAAE,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,CAAC,EAAE;KAC/D,CAAC;IAEF,kFAAkF;IAClF,iFAAiF;IACjF,6BAA6B;IAC7B,MAAM,eAAe,GAAG,cAAc,CAAC;QACrC,GAAG,UAAU;QACb,MAAM,EAAE,kBAAkB;KAC3B,CAAC,CAAC;IAEH,4EAA4E;IAC5E,wDAAwD;IACxD,MAAM,aAAa,GAAG,cAAc,CAAC;QACnC,GAAG,UAAU;QACb,MAAM,EAAE,gBAAgB;KACzB,CAAC,CAAC;IAEH,MAAM,YAAY,GAAG,eAAe,CAAC,YAAY,CAAC;IAClD,MAAM,UAAU,GAAG,aAAa,CAAC,UAAU,CAAC;IAE5C,MAAM,UAAU,GAAG,YAAY,CAAC,UAAU,GAAG,UAAU,CAAC,UAAU,CAAC;IACnE,MAAM,UAAU,GAAG,YAAY,CAAC,UAAU,GAAG,UAAU,CAAC,UAAU,CAAC;IACnE,MAAM,IAAI,GACR,UAAU,GAAG,CAAC;QACZ,CAAC,CAAC,UAAU,GAAG,UAAU;QACzB,CAAC,CAAC,UAAU,GAAG,CAAC;YACd,CAAC,CAAC,CAAC,QAAQ;YACX,CAAC,CAAC,QAAQ,CAAC;IAEjB,OAAO;QACL,UAAU;QACV,UAAU;QACV,IAAI;QACJ,iBAAiB,EAAE,YAAY,CAAC,UAAU;QAC1C,eAAe,EAAE,UAAU,CAAC,UAAU;QACtC,iBAAiB,EAAE,YAAY,CAAC,UAAU;QAC1C,eAAe,EAAE,UAAU,CAAC,UAAU;QACtC,gBAAgB,EAAE,YAAY,CAAC,SAAS;QACxC,cAAc,EAAE,UAAU,CAAC,SAAS;KACrC,CAAC;AACJ,CAAC"}