heor-agent-mcp 0.2.0 → 0.7.0

package/README.md

@@ -126,6 +126,50 @@ Pass `project: "project-id"` to any tool and results are saved automatically.
 
 ---
 
+## Examples
+
+Copy-paste prompts to try in Claude Code, Claude Desktop, or the [web UI](https://web-michael-ns-projects.vercel.app).
+
+### Single-tool examples
+
+**Literature search**
+> Search the literature for tirzepatide cardiovascular outcomes in type 2 diabetes. Use PubMed, ClinicalTrials.gov, and NICE TAs.
+
+**Survival curve fitting**
+> Fit survival curves to this OS data from KEYNOTE-189: time 0 survival 1.0, time 6 survival 0.88, time 12 survival 0.72, time 18 survival 0.60, time 24 survival 0.51, time 36 survival 0.38. Use months.
+
+**Budget impact**
+> Estimate the 5-year NHS budget impact of semaglutide for obesity. 200,000 eligible patients, drug cost £1,200/year, comparator (orlistat) £250/year, uptake 15% year 1 to 40% year 5.
+
+**Cost-effectiveness model**
+> Build a CE model for semaglutide vs sitagliptin in T2D, NHS perspective, lifetime horizon, with PSA.
+
+**Indirect comparison (Bucher)**
+> I have two trials: SUSTAIN-1 showed semaglutide vs placebo HR 0.74 (0.58-0.95) for HbA1c, and AWARD-5 showed dulaglutide vs placebo HR 0.78 (0.65-0.93). Run a Bucher indirect comparison between semaglutide and dulaglutide.
+
+**MAIC (population-adjusted comparison)**
+> Run a MAIC between SUSTAIN-7 (N=300, semaglutide vs placebo, HR 0.74, CI 0.58-0.95, age 56±10, BMI 33±5) and AWARD-11 (N=600, dulaglutide vs placebo, HR 0.78, CI 0.65-0.93, age 58±9, BMI 35±6). Adjust for age and BMI.
+
+### Multi-tool workflows
+
+**Abstract screening workflow**
+> Search PubMed for pembrolizumab in NSCLC, then screen the results with population adults with NSCLC, intervention pembrolizumab, comparator chemotherapy, outcomes overall survival and PFS.
+
+**Evidence network + NMA feasibility**
+> Search for GLP-1 receptor agonists in T2D using PubMed, build an evidence network from the results, and assess NMA feasibility.
+
+**CE model with scenarios**
+> Build a CE model for dapagliflozin vs placebo in heart failure, NHS perspective, lifetime horizon, with PSA. Add scenarios: "20% price reduction" with drug cost 400, "10-year horizon" with time_horizon 10yr.
+
+### End-to-end HTA workflow
+
+**Full dossier preparation**
+> Create a project for semaglutide in obesity targeting NICE and ICER. Search literature for evidence, screen the results for adults with obesity comparing semaglutide to placebo for weight loss outcomes, then draft a NICE STA dossier using the screened results.
+
+This single prompt exercises: `project_create` → `literature_search` → `screen_abstracts` → `hta_dossier_prep` (with auto-GRADE).
+
+---
+
 ## Data Sources
 
 **41 sources across 9 categories.** Every `literature_search` call includes a source selection table showing used/not-used status and reason for each.

package/dist/formatters/comparisonMarkdown.d.ts

@@ -0,0 +1,3 @@
+import type { IndirectComparisonResult } from "../network/types.js";
+export declare function comparisonToMarkdown(result: IndirectComparisonResult): string;
+//# sourceMappingURL=comparisonMarkdown.d.ts.map

package/dist/formatters/comparisonMarkdown.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"comparisonMarkdown.d.ts","sourceRoot":"","sources":["../../src/formatters/comparisonMarkdown.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,wBAAwB,EAAE,MAAM,qBAAqB,CAAC;AAYpE,wBAAgB,oBAAoB,CAClC,MAAM,EAAE,wBAAwB,GAC/B,MAAM,CAgFR"}

package/dist/formatters/comparisonMarkdown.js

@@ -0,0 +1,68 @@
+function formatNum(n, decimals = 2) {
+    return n.toFixed(decimals);
+}
+function formatP(p) {
+    if (p < 0.001)
+        return "<0.001";
+    if (p < 0.01)
+        return formatNum(p, 3);
+    return formatNum(p, 3);
+}
+export function comparisonToMarkdown(result) {
+    const lines = [];
+    lines.push("## Indirect Treatment Comparisons");
+    lines.push("");
+    const methodLabel = result.method === "bucher"
+        ? "Bucher Method"
+        : result.method === "frequentist_nma"
+            ? "Frequentist Network Meta-Analysis"
+            : "Mixed Methods";
+    lines.push(`**Method:** ${methodLabel}`);
+    lines.push(`**Comparisons:** ${result.estimates.length}`);
+    lines.push("");
+    if (result.estimates.length === 0) {
+        lines.push("No indirect comparisons could be computed. Check that comparisons share a common comparator and measure the same outcome.");
+        return lines.join("\n");
+    }
+    // Group by outcome
+    const outcomes = new Map();
+    for (const e of result.estimates) {
+        if (!outcomes.has(e.outcome))
+            outcomes.set(e.outcome, []);
+        outcomes.get(e.outcome).push(e);
+    }
+    for (const [outcome, estimates] of outcomes) {
+        lines.push(`### ${outcome}`);
+        lines.push("");
+        const measure = estimates[0].measure;
+        const label = measure === "MD" ? "MD" : measure;
+        lines.push(`| Comparison | ${label} | 95% CI | p-value | Via | Method |`);
+        lines.push("|------------|---------|--------|---------|-----|--------|");
+        for (const e of estimates) {
+            const est = measure === "MD"
+                ? formatNum(e.estimate)
+                : formatNum(e.estimate);
+            const ci = `[${formatNum(e.ci_lower)}, ${formatNum(e.ci_upper)}]`;
+            const p = formatP(e.p_value);
+            const via = e.commonComparator;
+            const method = e.method === "bucher" ? "Bucher" : "Freq. NMA";
+            lines.push(`| ${e.intervention} vs ${e.comparator} | ${est} | ${ci} | ${p} | ${via} | ${method} |`);
+        }
+        lines.push("");
+    }
+    // Warnings
+    if (result.warnings.length > 0) {
+        lines.push("### Warnings");
+        for (const w of result.warnings) {
+            lines.push(`- ⚠️ ${w}`);
+        }
+        lines.push("");
+    }
+    // Limitations
+    lines.push("### Limitations");
+    for (const l of result.limitations) {
+        lines.push(`- ${l}`);
+    }
+    return lines.join("\n");
+}
+//# sourceMappingURL=comparisonMarkdown.js.map

package/dist/formatters/comparisonMarkdown.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"comparisonMarkdown.js","sourceRoot":"","sources":["../../src/formatters/comparisonMarkdown.ts"],"names":[],"mappings":"AAEA,SAAS,SAAS,CAAC,CAAS,EAAE,WAAmB,CAAC;IAChD,OAAO,CAAC,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;AAC7B,CAAC;AAED,SAAS,OAAO,CAAC,CAAS;IACxB,IAAI,CAAC,GAAG,KAAK;QAAE,OAAO,QAAQ,CAAC;IAC/B,IAAI,CAAC,GAAG,IAAI;QAAE,OAAO,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,OAAO,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACzB,CAAC;AAED,MAAM,UAAU,oBAAoB,CAClC,MAAgC;IAEhC,MAAM,KAAK,GAAa,EAAE,CAAC;IAE3B,KAAK,CAAC,IAAI,CAAC,mCAAmC,CAAC,CAAC;IAChD,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAEf,MAAM,WAAW,GACf,MAAM,CAAC,MAAM,KAAK,QAAQ;QACxB,CAAC,CAAC,eAAe;QACjB,CAAC,CAAC,MAAM,CAAC,MAAM,KAAK,iBAAiB;YACnC,CAAC,CAAC,mCAAmC;YACrC,CAAC,CAAC,eAAe,CAAC;IACxB,KAAK,CAAC,IAAI,CAAC,eAAe,WAAW,EAAE,CAAC,CAAC;IACzC,KAAK,CAAC,IAAI,CACR,oBAAoB,MAAM,CAAC,SAAS,CAAC,MAAM,EAAE,CAC9C,CAAC;IACF,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAEf,IAAI,MAAM,CAAC,SAAS,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAClC,KAAK,CAAC,IAAI,CACR,2HAA2H,CAC5H,CAAC;QACF,OAAO,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC1B,CAAC;IAED,mBAAmB;IACnB,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAmC,CAAC;IAC5D,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,SAAS,EAAE,CAAC;QACjC,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC;YAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,EAAE,EAAE,CAAC,CAAC;QAC1D,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,CAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACnC,CAAC;IAED,KAAK,MAAM,CAAC,OAAO,EAAE,SAAS,CAAC,IAAI,QAAQ,EAAE,CAAC;QAC5C,KAAK,CAAC,IAAI,CAAC,OAAO,OAAO,EAAE,CAAC,CAAC;QAC7B,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAEf,MAAM,OAAO,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC;QACrC,MAAM,KAAK,GAAG,OAAO,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC;QAEhD,KAAK,CAAC,IAAI,CACR,kBAAkB,KAAK,sCAAsC,CAC9D,CAAC;QACF,KAAK,CAAC,IAAI,CACR,4DAA4D,CAC7D,CAAC;QAEF,KAAK,MAAM,CAAC,IAAI,SAAS,EAAE,CAAC;YAC1B,MAAM,GAAG,GACP,OAAO,KAAK,IAAI;gBACd,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC;gBACvB,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;YAC5B,MAAM,EAAE,GAAG,IAAI,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC;YAClE,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;YAC7B,MAAM,GAAG,GAAG,CAAC,CAAC,gBAAgB,CAAC;YAC/B,MAAM,MAAM,GACV,CAAC,CAAC,MAAM,KAAK,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,WAAW,CAAC;YAEjD,KAAK,CAAC,IAAI,CACR,KAAK,CAAC,CAAC,YAAY,OAAO,CAAC,CAAC,UAAU,MAAM,GAAG,MAAM,EAAE,MAAM,CAAC,MAAM,GAAG,MAAM,MAAM,IAAI,CACxF,CAAC;QACJ,CAAC;QACD,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACjB,CAAC;IAED,WAAW;IACX,IAAI,MAAM,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC/B,KAAK,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC;QAC3B,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,QAAQ,EAAE,CAAC;YAChC,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;QAC1B,CAAC;QACD,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACjB,CAAC;IAED,cAAc;IACd,KAAK,CAAC,IAAI,CAAC,iBAAiB,CAAC,CAAC;IAC9B,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,WAAW,EAAE,CAAC;QACnC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;IACvB,CAAC;IAED,OAAO,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AAC1B,CAAC"}

package/dist/models/evppi.d.ts

@@ -0,0 +1,33 @@
+/**
+ * Expected Value of Partial Perfect Information (EVPPI)
+ *
+ * Per-parameter VOI analysis: shows which specific parameters are
+ * worth further research by computing how much value would be gained
+ * from resolving uncertainty in each parameter individually.
+ *
+ * Method: non-parametric regression (GAM approximation via binning)
+ * following Strong et al. (2014) "Estimating multiparameter partial
+ * expected value of perfect information from a probabilistic
+ * sensitivity analysis sample"
+ */
+export interface EVPPIResult {
+    parameter: string;
+    evppi: number;
+    evppi_proportion: number;
+}
+export interface PSAIterationData {
+    delta_cost: number;
+    delta_qaly: number;
+    params: Record<string, number>;
+}
+/**
+ * Compute EVPPI for each parameter using the binning method.
+ *
+ * For each parameter theta_j:
+ * 1. Sort iterations by theta_j
+ * 2. Bin into K groups
+ * 3. Within each bin, compute E[max_arm(NMB)]
+ * 4. EVPPI(theta_j) = E_bins[max(NMB_bin)] - max(E[NMB])
+ */
+export declare function computeEVPPI(iterations: PSAIterationData[], lambda: number, parameterNames: string[]): EVPPIResult[];
+//# sourceMappingURL=evppi.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"evppi.d.ts","sourceRoot":"","sources":["../../src/models/evppi.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,MAAM,WAAW,WAAW;IAC1B,SAAS,EAAE,MAAM,CAAC;IAClB,KAAK,EAAE,MAAM,CAAC;IACd,gBAAgB,EAAE,MAAM,CAAC;CAC1B;AAED,MAAM,WAAW,gBAAgB;IAC/B,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;CAChC;AAED;;;;;;;;GAQG;AACH,wBAAgB,YAAY,CAC1B,UAAU,EAAE,gBAAgB,EAAE,EAC9B,MAAM,EAAE,MAAM,EACd,cAAc,EAAE,MAAM,EAAE,GACvB,WAAW,EAAE,CAsEf"}

package/dist/models/evppi.js

@@ -0,0 +1,75 @@
+/**
+ * Expected Value of Partial Perfect Information (EVPPI)
+ *
+ * Per-parameter VOI analysis: shows which specific parameters are
+ * worth further research by computing how much value would be gained
+ * from resolving uncertainty in each parameter individually.
+ *
+ * Method: non-parametric regression (GAM approximation via binning)
+ * following Strong et al. (2014) "Estimating multiparameter partial
+ * expected value of perfect information from a probabilistic
+ * sensitivity analysis sample"
+ */
+/**
+ * Compute EVPPI for each parameter using the binning method.
+ *
+ * For each parameter theta_j:
+ * 1. Sort iterations by theta_j
+ * 2. Bin into K groups
+ * 3. Within each bin, compute E[max_arm(NMB)]
+ * 4. EVPPI(theta_j) = E_bins[max(NMB_bin)] - max(E[NMB])
+ */
+export function computeEVPPI(iterations, lambda, parameterNames) {
+    const N = iterations.length;
+    if (N < 20)
+        return [];
+    // Overall expected NMB for each arm
+    const nmb_intervention_mean = iterations.reduce((sum, it) => sum + (lambda * it.delta_qaly - it.delta_cost), 0) / N;
+    const nmb_comparator_mean = 0;
+    const max_e_nmb = Math.max(nmb_intervention_mean, nmb_comparator_mean);
+    // Total EVPI for reference
+    const e_max_nmb = iterations.reduce((sum, it) => {
+        const nmb_int = lambda * it.delta_qaly - it.delta_cost;
+        return sum + Math.max(nmb_int, 0);
+    }, 0) / N;
+    const totalEVPI = Math.max(0, e_max_nmb - max_e_nmb);
+    // Number of bins (Sturges' rule)
+    const K = Math.max(5, Math.min(30, Math.ceil(1 + 3.322 * Math.log10(N))));
+    const binSize = Math.ceil(N / K);
+    const results = [];
+    for (const paramName of parameterNames) {
+        // Check if parameter exists in iterations
+        if (iterations[0]?.params[paramName] === undefined)
+            continue;
+        // Sort by this parameter
+        const sorted = [...iterations].sort((a, b) => (a.params[paramName] ?? 0) - (b.params[paramName] ?? 0));
+        // Bin and compute conditional expectations
+        let e_max_conditional = 0;
+        for (let bin = 0; bin < K; bin++) {
+            const start = bin * binSize;
+            const end = Math.min(start + binSize, N);
+            const binN = end - start;
+            if (binN === 0)
+                continue;
+            // Expected NMB within this bin for intervention
+            let sumNMBInt = 0;
+            for (let i = start; i < end; i++) {
+                const it = sorted[i];
+                sumNMBInt += lambda * it.delta_qaly - it.delta_cost;
+            }
+            const meanNMBInt = sumNMBInt / binN;
+            // Max of conditional expectations
+            e_max_conditional += Math.max(meanNMBInt, 0) * (binN / N);
+        }
+        const evppi = Math.max(0, e_max_conditional - max_e_nmb);
+        results.push({
+            parameter: paramName,
+            evppi,
+            evppi_proportion: totalEVPI > 0 ? evppi / totalEVPI : 0,
+        });
+    }
+    // Sort by EVPPI descending
+    results.sort((a, b) => b.evppi - a.evppi);
+    return results;
+}
+//# sourceMappingURL=evppi.js.map

package/dist/models/evppi.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"evppi.js","sourceRoot":"","sources":["../../src/models/evppi.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAcH;;;;;;;;GAQG;AACH,MAAM,UAAU,YAAY,CAC1B,UAA8B,EAC9B,MAAc,EACd,cAAwB;IAExB,MAAM,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC;IAC5B,IAAI,CAAC,GAAG,EAAE;QAAE,OAAO,EAAE,CAAC;IAEtB,oCAAoC;IACpC,MAAM,qBAAqB,GACzB,UAAU,CAAC,MAAM,CACf,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,CAAC,EAC3D,CAAC,CACF,GAAG,CAAC,CAAC;IACR,MAAM,mBAAmB,GAAG,CAAC,CAAC;IAC9B,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,qBAAqB,EAAE,mBAAmB,CAAC,CAAC;IAEvE,2BAA2B;IAC3B,MAAM,SAAS,GACb,UAAU,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE;QAC5B,MAAM,OAAO,GAAG,MAAM,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,CAAC;QACvD,OAAO,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IACpC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;IACZ,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,GAAG,SAAS,CAAC,CAAC;IAErD,iCAAiC;IACjC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1E,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAEjC,MAAM,OAAO,GAAkB,EAAE,CAAC;IAElC,KAAK,MAAM,SAAS,IAAI,cAAc,EAAE,CAAC;QACvC,0CAA0C;QAC1C,IAAI,UAAU,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,SAAS,CAAC,KAAK,SAAS;YAAE,SAAS;QAE7D,yBAAyB;QACzB,MAAM,MAAM,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,IAAI,CACjC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,SAAS,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,SAAS,CAAC,IAAI,CAAC,CAAC,CAClE,CAAC;QAEF,2CAA2C;QAC3C,IAAI,iBAAiB,GAAG,CAAC,CAAC;QAE1B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;YACjC,MAAM,KAAK,GAAG,GAAG,GAAG,OAAO,CAAC;YAC5B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,OAAO,EAAE,CAAC,CAAC,CAAC;YACzC,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;YACzB,IAAI,IAAI,KAAK,CAAC;gBAAE,SAAS;YAEzB,gDAAgD;YAChD,IAAI,SAAS,GAAG,CAAC,CAAC;YAClB,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBACjC,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAE,CAAC;gBACtB,SAAS,IAAI,MAAM,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,CAAC;YACtD,CAAC;YACD,MAAM,UAAU,GAAG,SAAS,GAAG,IAAI,CAAC;YAEpC,kCAAkC;YAClC,iBAAiB,IAAI,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QAC5D,CAAC;QAED,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,iBAAiB,GAAG,SAAS,CAAC,CAAC;QAEzD,OAAO,CAAC,IAAI,CAAC;YACX,SAAS,EAAE,SAAS;YACpB,KAAK;YACL,gBAAgB,EAAE,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;SACxD,CAAC,CAAC;IACL,CAAC;IAED,2BAA2B;IAC3B,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC;IAE1C,OAAO,OAAO,CAAC;AACjB,CAAC"}

package/dist/models/modelUtils.d.ts

@@ -4,9 +4,14 @@
  */
 import type { CEModelParams } from "../providers/types.js";
 import type { MarkovParams } from "./markov.js";
+export declare function getTimeHorizonYears(horizon: CEModelParams["time_horizon"]): number;
 /**
- * Build MarkovParams from CEModelParams using a 2-state model
- * (On-Treatment, Off-Treatment).
+ * Build MarkovParams from CEModelParams using a 3-state model
+ * (On-Treatment, Off-Treatment, Dead).
+ *
+ * The Dead state is an absorbing state (utility=0, cost=0) — the cohort
+ * fraction that enters Dead never leaves, which correctly bounds life-year
+ * and QALY accumulation over the time horizon.
  */
 export declare function buildMarkovParamsFromCE(params: CEModelParams): MarkovParams;
 /**

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

@@ -1 +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"}
+{"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;AAK/E,wBAAgB,mBAAmB,CACjC,OAAO,EAAE,aAAa,CAAC,cAAc,CAAC,GACrC,MAAM,CAKR;AAED;;;;;;;GAOG;AACH,wBAAgB,uBAAuB,CAAC,MAAM,EAAE,aAAa,GAAG,YAAY,CAgG3E;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,CAgEA"}

package/dist/models/modelUtils.js

@@ -4,7 +4,7 @@
  */
 import { runMarkovModel } from "./markov.js";
 const DISCOUNT_RATE = 0.035;
-function getTimeHorizonYears(horizon) {
+export function getTimeHorizonYears(horizon) {
     if (horizon === "lifetime")
         return 40;
     if (horizon === "5yr")
@@ -14,8 +14,12 @@ function getTimeHorizonYears(horizon) {
     return Number(horizon);
 }
 /**
- * Build MarkovParams from CEModelParams using a 2-state model
- * (On-Treatment, Off-Treatment).
+ * Build MarkovParams from CEModelParams using a 3-state model
+ * (On-Treatment, Off-Treatment, Dead).
+ *
+ * The Dead state is an absorbing state (utility=0, cost=0) — the cohort
+ * fraction that enters Dead never leaves, which correctly bounds life-year
+ * and QALY accumulation over the time horizon.
  */
 export function buildMarkovParamsFromCE(params) {
     const years = getTimeHorizonYears(params.time_horizon);
@@ -29,50 +33,62 @@ export function buildMarkovParamsFromCE(params) {
     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));
+    // Background annual mortality rate (~2% baseline, reduced by mortalityReduction)
+    const baseMortality = 0.02;
+    const interventionMortality = Math.max(0.005, baseMortality * (1 - mortalityReduction));
+    const comparatorMortality = baseMortality;
+    // Transition: prob of staying in On-Treatment state
+    const probStayOnIntervention = Math.max(0.05, Math.min(0.93, 0.5 + efficacyDelta * 0.5));
+    const baselineProbStayOn = Math.max(0.05, Math.min(0.88, probStayOnIntervention * 0.7));
     const states = [
         { name: "On-Treatment", utility: utilityOn, cost_annual: costIntervention },
         { name: "Off-Treatment", utility: utilityOff, cost_annual: 0 },
+        { name: "Dead", utility: 0, cost_annual: 0 },
     ];
-    const statesComparator = [
-        { name: "On-Treatment", utility: utilityOn, cost_annual: costComparator },
-        { name: "Off-Treatment", utility: utilityOff, cost_annual: 0 },
-    ];
+    // Helper: ensure row sums to 1.0
+    function normalizeRow(row) {
+        const sum = Object.values(row).reduce((a, b) => a + b, 0);
+        if (Math.abs(sum - 1.0) < 1e-10)
+            return row;
+        const result = {};
+        for (const [k, v] of Object.entries(row)) {
+            result[k] = v / sum;
+        }
+        return result;
+    }
     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": normalizeRow({
+            "On-Treatment": probStayOnIntervention,
+            "Off-Treatment": 1 - probStayOnIntervention - interventionMortality,
+            Dead: interventionMortality,
+        }),
+        "Off-Treatment": normalizeRow({
             "On-Treatment": 0.05,
-            "Off-Treatment": 0.95,
-        },
+            "Off-Treatment": 0.95 - interventionMortality,
+            Dead: interventionMortality,
+        }),
+        Dead: { "On-Treatment": 0, "Off-Treatment": 0, Dead: 1 },
     };
     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": normalizeRow({
+            "On-Treatment": baselineProbStayOn,
+            "Off-Treatment": 1 - baselineProbStayOn - comparatorMortality,
+            Dead: comparatorMortality,
+        }),
+        "Off-Treatment": normalizeRow({
             "On-Treatment": 0.05,
-            "Off-Treatment": 0.95,
-        },
+            "Off-Treatment": 0.95 - comparatorMortality,
+            Dead: comparatorMortality,
+        }),
+        Dead: { "On-Treatment": 0, "Off-Treatment": 0, Dead: 1 },
     };
     return {
         states,
         transition_matrix_intervention,
         transition_matrix_comparator,
-        initial_cohort: { "On-Treatment": 1.0, "Off-Treatment": 0.0 },
+        initial_cohort: { "On-Treatment": 1.0, "Off-Treatment": 0.0, Dead: 0.0 },
         cycle_length_years,
         n_cycles,
         discount_rate_costs: DISCOUNT_RATE,
@@ -97,10 +113,12 @@ export function runMarkovAndComputeICER(params) {
     const statesIntervention = [
         { name: "On-Treatment", utility: utilityOn, cost_annual: costIntervention },
         { name: "Off-Treatment", utility: utilityOff, cost_annual: 0 },
+        { name: "Dead", utility: 0, cost_annual: 0 },
     ];
     const statesComparator = [
         { name: "On-Treatment", utility: utilityOn, cost_annual: costComparator },
         { name: "Off-Treatment", utility: utilityOff, cost_annual: 0 },
+        { name: "Dead", utility: 0, 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

package/dist/models/modelUtils.js.map

@@ -1 +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"}
+{"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,MAAM,UAAU,mBAAmB,CACjC,OAAsC;IAEtC,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;;;;;;;GAOG;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,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,iFAAiF;IACjF,MAAM,aAAa,GAAG,IAAI,CAAC;IAC3B,MAAM,qBAAqB,GAAG,IAAI,CAAC,GAAG,CACpC,KAAK,EACL,aAAa,GAAG,CAAC,CAAC,GAAG,kBAAkB,CAAC,CACzC,CAAC;IACF,MAAM,mBAAmB,GAAG,aAAa,CAAC;IAE1C,oDAAoD;IACpD,MAAM,sBAAsB,GAAG,IAAI,CAAC,GAAG,CACrC,IAAI,EACJ,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,GAAG,aAAa,GAAG,GAAG,CAAC,CAC1C,CAAC;IACF,MAAM,kBAAkB,GAAG,IAAI,CAAC,GAAG,CACjC,IAAI,EACJ,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,sBAAsB,GAAG,GAAG,CAAC,CAC7C,CAAC;IAEF,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;QAC9D,EAAE,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE;KAC7C,CAAC;IAEF,iCAAiC;IACjC,SAAS,YAAY,CAAC,GAA2B;QAC/C,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1D,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,KAAK;YAAE,OAAO,GAAG,CAAC;QAC5C,MAAM,MAAM,GAA2B,EAAE,CAAC;QAC1C,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,EAAE,CAAC;YACzC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;QACtB,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,MAAM,8BAA8B,GAAqB;QACvD,cAAc,EAAE,YAAY,CAAC;YAC3B,cAAc,EAAE,sBAAsB;YACtC,eAAe,EAAE,CAAC,GAAG,sBAAsB,GAAG,qBAAqB;YACnE,IAAI,EAAE,qBAAqB;SAC5B,CAAC;QACF,eAAe,EAAE,YAAY,CAAC;YAC5B,cAAc,EAAE,IAAI;YACpB,eAAe,EAAE,IAAI,GAAG,qBAAqB;YAC7C,IAAI,EAAE,qBAAqB;SAC5B,CAAC;QACF,IAAI,EAAE,EAAE,cAAc,EAAE,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE;KACzD,CAAC;IAEF,MAAM,4BAA4B,GAAqB;QACrD,cAAc,EAAE,YAAY,CAAC;YAC3B,cAAc,EAAE,kBAAkB;YAClC,eAAe,EAAE,CAAC,GAAG,kBAAkB,GAAG,mBAAmB;YAC7D,IAAI,EAAE,mBAAmB;SAC1B,CAAC;QACF,eAAe,EAAE,YAAY,CAAC;YAC5B,cAAc,EAAE,IAAI;YACpB,eAAe,EAAE,IAAI,GAAG,mBAAmB;YAC3C,IAAI,EAAE,mBAAmB;SAC1B,CAAC;QACF,IAAI,EAAE,EAAE,cAAc,EAAE,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE;KACzD,CAAC;IAEF,OAAO;QACL,MAAM;QACN,8BAA8B;QAC9B,4BAA4B;QAC5B,cAAc,EAAE,EAAE,cAAc,EAAE,GAAG,EAAE,eAAe,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;QACxE,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;QAC9D,EAAE,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE;KAC7C,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;QAC9D,EAAE,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE;KAC7C,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"}

package/dist/models/psa.d.ts

@@ -2,11 +2,14 @@
  * Probabilistic Sensitivity Analysis (PSA) using Monte Carlo simulation.
  */
 import type { CEModelParams } from "../providers/types.js";
+import { type EVPPIResult } from "./evppi.js";
 export type CEModelInputs = CEModelParams;
 export interface PSAParams {
     base_params: CEModelInputs;
     n_iterations: number;
     seed?: number;
+    /** WTP threshold for EVPI calculation. Defaults to perspective-appropriate value. */
+    evpi_lambda?: number;
 }
 export interface PSAIteration {
     delta_cost: number;
@@ -24,6 +27,7 @@ export interface PSAResult {
         prob_ce: number;
     }>;
     evpi: number;
+    evppi: EVPPIResult[];
     scatter_sample: PSAIteration[];
 }
 /**

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

@@ -1 +1 @@
-{"version":3,"file":"psa.d.ts","sourceRoot":"","sources":["../../src/models/psa.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AAO3D,MAAM,MAAM,aAAa,GAAG,aAAa,CAAC;AAE1C,MAAM,WAAW,SAAS;IACxB,WAAW,EAAE,aAAa,CAAC;IAC3B,YAAY,EAAE,MAAM,CAAC;IACrB,IAAI,CAAC,EAAE,MAAM,CAAC;CACf;AAED,MAAM,WAAW,YAAY;IAC3B,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,IAAI,EAAE,MAAM,CAAC;CACd;AAED,MAAM,WAAW,SAAS;IACxB,UAAU,EAAE,YAAY,EAAE,CAAC;IAC3B,SAAS,EAAE,MAAM,CAAC;IAClB,aAAa,EAAE,MAAM,CAAC;IACtB,aAAa,EAAE,MAAM,CAAC;IACtB,mBAAmB,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC5C,IAAI,EAAE,KAAK,CAAC;QAAE,GAAG,EAAE,MAAM,CAAC;QAAC,OAAO,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;IAC9C,IAAI,EAAE,MAAM,CAAC;IACb,cAAc,EAAE,YAAY,EAAE,CAAC;CAChC;AAiCD;;GAEG;AACH,wBAAgB,MAAM,CAAC,MAAM,EAAE,SAAS,GAAG,SAAS,CA8InD"}
+{"version":3,"file":"psa.d.ts","sourceRoot":"","sources":["../../src/models/psa.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AAM3D,OAAO,EAEL,KAAK,WAAW,EAEjB,MAAM,YAAY,CAAC;AAEpB,MAAM,MAAM,aAAa,GAAG,aAAa,CAAC;AAE1C,MAAM,WAAW,SAAS;IACxB,WAAW,EAAE,aAAa,CAAC;IAC3B,YAAY,EAAE,MAAM,CAAC;IACrB,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,qFAAqF;IACrF,WAAW,CAAC,EAAE,MAAM,CAAC;CACtB;AAED,MAAM,WAAW,YAAY;IAC3B,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,MAAM,CAAC;IACnB,IAAI,EAAE,MAAM,CAAC;CACd;AAED,MAAM,WAAW,SAAS;IACxB,UAAU,EAAE,YAAY,EAAE,CAAC;IAC3B,SAAS,EAAE,MAAM,CAAC;IAClB,aAAa,EAAE,MAAM,CAAC;IACtB,aAAa,EAAE,MAAM,CAAC;IACtB,mBAAmB,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC5C,IAAI,EAAE,KAAK,CAAC;QAAE,GAAG,EAAE,MAAM,CAAC;QAAC,OAAO,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;IAC9C,IAAI,EAAE,MAAM,CAAC;IACb,KAAK,EAAE,WAAW,EAAE,CAAC;IACrB,cAAc,EAAE,YAAY,EAAE,CAAC;CAChC;AAiCD;;GAEG;AACH,wBAAgB,MAAM,CAAC,MAAM,EAAE,SAAS,GAAG,SAAS,CA8KnD"}

package/dist/models/psa.js

@@ -3,6 +3,7 @@
  */
 import { betaSample, gammaSample, createSeededRng } from "./distributions.js";
 import { runMarkovAndComputeICER, } from "./modelUtils.js";
+import { computeEVPPI, } from "./evppi.js";
 const WTP_THRESHOLDS_PSA = {
     nhs_low: 25000,
     nhs_high: 35000,
@@ -37,6 +38,7 @@ function perturbGamma(value, cv, rng) {
 export function runPSA(params) {
     const rng = createSeededRng(params.seed ?? 42);
     const iterations = [];
+    const evppiData = [];
     for (let i = 0; i < params.n_iterations; i++) {
         // Sample perturbed parameters
         const perturbedParams = {
@@ -67,6 +69,27 @@ export function runPSA(params) {
         };
         const { delta_cost, delta_qaly, icer } = runMarkovAndComputeICER(perturbedParams);
         iterations.push({ delta_cost, delta_qaly, icer });
+        // Track parameter values for EVPPI
+        evppiData.push({
+            delta_cost,
+            delta_qaly,
+            params: {
+                efficacy_delta: perturbedParams.clinical_inputs.efficacy_delta,
+                drug_cost_annual: perturbedParams.cost_inputs.drug_cost_annual,
+                comparator_cost_annual: perturbedParams.cost_inputs.comparator_cost_annual,
+                ...(perturbedParams.clinical_inputs.mortality_reduction !== undefined
+                    ? {
+                        mortality_reduction: perturbedParams.clinical_inputs.mortality_reduction,
+                    }
+                    : {}),
+                ...(perturbedParams.utility_inputs
+                    ? {
+                        qaly_on_treatment: perturbedParams.utility_inputs.qaly_on_treatment,
+                        qaly_comparator: perturbedParams.utility_inputs.qaly_comparator,
+                    }
+                    : {}),
+            },
+        });
     }
     // Compute statistics
     const finiteICERs = iterations
@@ -96,8 +119,8 @@ export function runPSA(params) {
         ceac.push({ wtp, prob_ce: n_ce / iterations.length });
     }
     // EVPI: E[max_arm(NMB)] - max_arm(E[NMB])
-    // Using a reference WTP (midpoint of common thresholds: $50,000)
-    const lambda = 50000;
+    // Use perspective-appropriate WTP threshold (passed via evpi_lambda)
+    const lambda = params.evpi_lambda ?? 50000;
     const nmb_intervention_mean = iterations.reduce((sum, it) => sum + (lambda * it.delta_qaly - it.delta_cost), 0) / iterations.length;
     const nmb_comparator_mean = 0; // comparator is reference (delta = 0)
     const e_max_nmb = iterations.reduce((sum, it) => {
@@ -106,6 +129,9 @@ export function runPSA(params) {
     }, 0) / iterations.length;
     const max_e_nmb = Math.max(nmb_intervention_mean, nmb_comparator_mean);
     const evpi = Math.max(0, e_max_nmb - max_e_nmb);
+    // EVPPI: per-parameter partial value of information
+    const paramNames = evppiData.length > 0 ? Object.keys(evppiData[0].params) : [];
+    const evppiResults = computeEVPPI(evppiData, lambda, paramNames);
     const scatter_sample = iterations.slice(0, 500);
     return {
         iterations,
@@ -115,6 +141,7 @@ export function runPSA(params) {
         prob_cost_effective,
         ceac,
         evpi,
+        evppi: evppiResults,
         scatter_sample,
     };
 }

package/dist/models/psa.js.map

@@ -1 +1 @@
-{"version":3,"file":"psa.js","sourceRoot":"","sources":["../../src/models/psa.ts"],"names":[],"mappings":"AAAA;;GAEG;AAGH,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAC9E,OAAO,EAEL,uBAAuB,GACxB,MAAM,iBAAiB,CAAC;AA2BzB,MAAM,kBAAkB,GAA2B;IACjD,OAAO,EAAE,KAAK;IACd,QAAQ,EAAE,KAAK;IACf,YAAY,EAAE,MAAM;IACpB,YAAY,EAAE,MAAM;IACpB,QAAQ,EAAE,KAAK;CAChB,CAAC;AAEF;;;GAGG;AACH,SAAS,WAAW,CAAC,KAAa,EAAE,EAAU,EAAE,GAAiB;IAC/D,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC,CAAC;IACrD,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IACxC,4CAA4C;IAC5C,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC;IAC7C,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,WAAW,CAAC,CAAC;IACrD,OAAO,UAAU,CAAC,IAAI,EAAE,YAAY,EAAE,GAAG,CAAC,CAAC;AAC7C,CAAC;AAED;;;GAGG;AACH,SAAS,YAAY,CAAC,KAAa,EAAE,EAAU,EAAE,GAAiB;IAChE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IAChC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IACxC,OAAO,WAAW,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,MAAM,CAAC,MAAiB;IACtC,MAAM,GAAG,GAAG,eAAe,CAAC,MAAM,CAAC,IAAI,IAAI,EAAE,CAAC,CAAC;IAC/C,MAAM,UAAU,GAAmB,EAAE,CAAC;IAEtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,YAAY,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7C,8BAA8B;QAC9B,MAAM,eAAe,GAAkB;YACrC,GAAG,MAAM,CAAC,WAAW;YACrB,eAAe,EAAE;gBACf,GAAG,MAAM,CAAC,WAAW,CAAC,eAAe;gBACrC,cAAc,EAAE,WAAW,CACzB,IAAI,CAAC,GAAG,CACN,KAAK,EACL,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,WAAW,CAAC,eAAe,CAAC,cAAc,CAAC,CACnE,EACD,GAAG,EACH,GAAG,CACJ;gBACD,mBAAmB,EACjB,MAAM,CAAC,WAAW,CAAC,eAAe,CAAC,mBAAmB,KAAK,SAAS;oBAClE,CAAC,CAAC,WAAW,CACT,IAAI,CAAC,GAAG,CACN,KAAK,EACL,IAAI,CAAC,GAAG,CACN,KAAK,EACL,MAAM,CAAC,WAAW,CAAC,eAAe,CAAC,mBAAmB,CACvD,CACF,EACD,GAAG,EACH,GAAG,CACJ;oBACH,CAAC,CAAC,SAAS;aAChB;YACD,WAAW,EAAE;gBACX,gBAAgB,EAAE,YAAY,CAC5B,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,gBAAgB,EAC/C,GAAG,EACH,GAAG,CACJ;gBACD,sBAAsB,EAAE,YAAY,CAClC,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,sBAAsB,EACrD,GAAG,EACH,GAAG,CACJ;gBACD,UAAU,EACR,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,UAAU,KAAK,SAAS;oBACrD,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,UAAU,EAAE,GAAG,EAAE,GAAG,CAAC;oBACnE,CAAC,CAAC,SAAS;gBACf,OAAO,EACL,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,OAAO,KAAK,SAAS;oBAClD,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,OAAO,EAAE,GAAG,EAAE,GAAG,CAAC;oBAChE,CAAC,CAAC,SAAS;aAChB;YACD,cAAc,EAAE,MAAM,CAAC,WAAW,CAAC,cAAc;gBAC/C,CAAC,CAAC;oBACE,iBAAiB,EAAE,WAAW,CAC5B,MAAM,CAAC,WAAW,CAAC,cAAc,CAAC,iBAAiB,EACnD,IAAI,EACJ,GAAG,CACJ;oBACD,eAAe,EAAE,WAAW,CAC1B,MAAM,CAAC,WAAW,CAAC,cAAc,CAAC,eAAe,EACjD,IAAI,EACJ,GAAG,CACJ;iBACF;gBACH,CAAC,CAAC,SAAS;SACd,CAAC;QAEF,MAAM,EAAE,UAAU,EAAE,UAAU,EAAE,IAAI,EAAE,GACpC,uBAAuB,CAAC,eAAe,CAAC,CAAC;QAC3C,UAAU,CAAC,IAAI,CAAC,EAAE,UAAU,EAAE,UAAU,EAAE,IAAI,EAAE,CAAC,CAAC;IACpD,CAAC;IAED,qBAAqB;IACrB,MAAM,WAAW,GAAG,UAAU;SAC3B,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC;SACpB,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAChC,MAAM,WAAW,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAE3D,MAAM,SAAS,GACb,WAAW,CAAC,MAAM,GAAG,CAAC;QACpB,CAAC,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,WAAW,CAAC,MAAM;QAC7D,CAAC,CAAC,QAAQ,CAAC;IAEf,MAAM,aAAa,GACjB,WAAW,CAAC,MAAM,GAAG,CAAC;QACpB,CAAC,CAAC,CAAC,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,WAAW,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,IAAI,WAAW,CAAC,CAAC,CAAE,CAAC;QAC1E,CAAC,CAAC,CAAC,CAAC;IACR,MAAM,aAAa,GACjB,WAAW,CAAC,MAAM,GAAG,CAAC;QACpB,CAAC,CAAC,CAAC,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,WAAW,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC;YACpD,WAAW,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,CAAE,CAAC;QACvC,CAAC,CAAC,CAAC,CAAC;IAER,iDAAiD;IACjD,MAAM,mBAAmB,GAA2B,EAAE,CAAC;IACvD,KAAK,MAAM,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,MAAM,CAAC,OAAO,CAAC,kBAAkB,CAAC,EAAE,CAAC;QAC7D,MAAM,IAAI,GAAG,UAAU,CAAC,MAAM,CAC5B,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,GAAG,CAAC,CAChD,CAAC,MAAM,CAAC;QACT,mBAAmB,CAAC,IAAI,CAAC,GAAG,IAAI,GAAG,UAAU,CAAC,MAAM,CAAC;IACvD,CAAC;IAED,gDAAgD;IAChD,MAAM,IAAI,GAA4C,EAAE,CAAC;IACzD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,IAAI,MAAM,EAAE,GAAG,IAAI,IAAI,EAAE,CAAC;QAC7C,MAAM,IAAI,GAAG,UAAU,CAAC,MAAM,CAC5B,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,GAAG,CAAC,CAChD,CAAC,MAAM,CAAC;QACT,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,IAAI,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,CAAC;IACxD,CAAC;IAED,0CAA0C;IAC1C,iEAAiE;IACjE,MAAM,MAAM,GAAG,KAAK,CAAC;IACrB,MAAM,qBAAqB,GACzB,UAAU,CAAC,MAAM,CACf,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,CAAC,EAC3D,CAAC,CACF,GAAG,UAAU,CAAC,MAAM,CAAC;IACxB,MAAM,mBAAmB,GAAG,CAAC,CAAC,CAAC,sCAAsC;IACrE,MAAM,SAAS,GACb,UAAU,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE;QAC5B,MAAM,OAAO,GAAG,MAAM,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,CAAC;QACvD,OAAO,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,mBAAmB,CAAC,CAAC;IACtD,CAAC,EAAE,CAAC,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC;IAC5B,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,qBAAqB,EAAE,mBAAmB,CAAC,CAAC;IACvE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,GAAG,SAAS,CAAC,CAAC;IAEhD,MAAM,cAAc,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAEhD,OAAO;QACL,UAAU;QACV,SAAS;QACT,aAAa;QACb,aAAa;QACb,mBAAmB;QACnB,IAAI;QACJ,IAAI;QACJ,cAAc;KACf,CAAC;AACJ,CAAC"}
+{"version":3,"file":"psa.js","sourceRoot":"","sources":["../../src/models/psa.ts"],"names":[],"mappings":"AAAA;;GAEG;AAGH,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAC9E,OAAO,EAEL,uBAAuB,GACxB,MAAM,iBAAiB,CAAC;AACzB,OAAO,EACL,YAAY,GAGb,MAAM,YAAY,CAAC;AA8BpB,MAAM,kBAAkB,GAA2B;IACjD,OAAO,EAAE,KAAK;IACd,QAAQ,EAAE,KAAK;IACf,YAAY,EAAE,MAAM;IACpB,YAAY,EAAE,MAAM;IACpB,QAAQ,EAAE,KAAK;CAChB,CAAC;AAEF;;;GAGG;AACH,SAAS,WAAW,CAAC,KAAa,EAAE,EAAU,EAAE,GAAiB;IAC/D,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC,CAAC;IACrD,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IACxC,4CAA4C;IAC5C,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC;IAC7C,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,WAAW,CAAC,CAAC;IACrD,OAAO,UAAU,CAAC,IAAI,EAAE,YAAY,EAAE,GAAG,CAAC,CAAC;AAC7C,CAAC;AAED;;;GAGG;AACH,SAAS,YAAY,CAAC,KAAa,EAAE,EAAU,EAAE,GAAiB;IAChE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IAChC,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IACxC,OAAO,WAAW,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,MAAM,CAAC,MAAiB;IACtC,MAAM,GAAG,GAAG,eAAe,CAAC,MAAM,CAAC,IAAI,IAAI,EAAE,CAAC,CAAC;IAC/C,MAAM,UAAU,GAAmB,EAAE,CAAC;IACtC,MAAM,SAAS,GAAuB,EAAE,CAAC;IAEzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,YAAY,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7C,8BAA8B;QAC9B,MAAM,eAAe,GAAkB;YACrC,GAAG,MAAM,CAAC,WAAW;YACrB,eAAe,EAAE;gBACf,GAAG,MAAM,CAAC,WAAW,CAAC,eAAe;gBACrC,cAAc,EAAE,WAAW,CACzB,IAAI,CAAC,GAAG,CACN,KAAK,EACL,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,WAAW,CAAC,eAAe,CAAC,cAAc,CAAC,CACnE,EACD,GAAG,EACH,GAAG,CACJ;gBACD,mBAAmB,EACjB,MAAM,CAAC,WAAW,CAAC,eAAe,CAAC,mBAAmB,KAAK,SAAS;oBAClE,CAAC,CAAC,WAAW,CACT,IAAI,CAAC,GAAG,CACN,KAAK,EACL,IAAI,CAAC,GAAG,CACN,KAAK,EACL,MAAM,CAAC,WAAW,CAAC,eAAe,CAAC,mBAAmB,CACvD,CACF,EACD,GAAG,EACH,GAAG,CACJ;oBACH,CAAC,CAAC,SAAS;aAChB;YACD,WAAW,EAAE;gBACX,gBAAgB,EAAE,YAAY,CAC5B,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,gBAAgB,EAC/C,GAAG,EACH,GAAG,CACJ;gBACD,sBAAsB,EAAE,YAAY,CAClC,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,sBAAsB,EACrD,GAAG,EACH,GAAG,CACJ;gBACD,UAAU,EACR,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,UAAU,KAAK,SAAS;oBACrD,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,UAAU,EAAE,GAAG,EAAE,GAAG,CAAC;oBACnE,CAAC,CAAC,SAAS;gBACf,OAAO,EACL,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,OAAO,KAAK,SAAS;oBAClD,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,WAAW,CAAC,WAAW,CAAC,OAAO,EAAE,GAAG,EAAE,GAAG,CAAC;oBAChE,CAAC,CAAC,SAAS;aAChB;YACD,cAAc,EAAE,MAAM,CAAC,WAAW,CAAC,cAAc;gBAC/C,CAAC,CAAC;oBACE,iBAAiB,EAAE,WAAW,CAC5B,MAAM,CAAC,WAAW,CAAC,cAAc,CAAC,iBAAiB,EACnD,IAAI,EACJ,GAAG,CACJ;oBACD,eAAe,EAAE,WAAW,CAC1B,MAAM,CAAC,WAAW,CAAC,cAAc,CAAC,eAAe,EACjD,IAAI,EACJ,GAAG,CACJ;iBACF;gBACH,CAAC,CAAC,SAAS;SACd,CAAC;QAEF,MAAM,EAAE,UAAU,EAAE,UAAU,EAAE,IAAI,EAAE,GACpC,uBAAuB,CAAC,eAAe,CAAC,CAAC;QAC3C,UAAU,CAAC,IAAI,CAAC,EAAE,UAAU,EAAE,UAAU,EAAE,IAAI,EAAE,CAAC,CAAC;QAElD,mCAAmC;QACnC,SAAS,CAAC,IAAI,CAAC;YACb,UAAU;YACV,UAAU;YACV,MAAM,EAAE;gBACN,cAAc,EAAE,eAAe,CAAC,eAAe,CAAC,cAAc;gBAC9D,gBAAgB,EAAE,eAAe,CAAC,WAAW,CAAC,gBAAgB;gBAC9D,sBAAsB,EACpB,eAAe,CAAC,WAAW,CAAC,sBAAsB;gBACpD,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,mBAAmB,KAAK,SAAS;oBACnE,CAAC,CAAC;wBACE,mBAAmB,EACjB,eAAe,CAAC,eAAe,CAAC,mBAAmB;qBACtD;oBACH,CAAC,CAAC,EAAE,CAAC;gBACP,GAAG,CAAC,eAAe,CAAC,cAAc;oBAChC,CAAC,CAAC;wBACE,iBAAiB,EACf,eAAe,CAAC,cAAc,CAAC,iBAAiB;wBAClD,eAAe,EAAE,eAAe,CAAC,cAAc,CAAC,eAAe;qBAChE;oBACH,CAAC,CAAC,EAAE,CAAC;aACR;SACF,CAAC,CAAC;IACL,CAAC;IAED,qBAAqB;IACrB,MAAM,WAAW,GAAG,UAAU;SAC3B,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC;SACpB,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAChC,MAAM,WAAW,GAAG,CAAC,GAAG,WAAW,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAE3D,MAAM,SAAS,GACb,WAAW,CAAC,MAAM,GAAG,CAAC;QACpB,CAAC,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,WAAW,CAAC,MAAM;QAC7D,CAAC,CAAC,QAAQ,CAAC;IAEf,MAAM,aAAa,GACjB,WAAW,CAAC,MAAM,GAAG,CAAC;QACpB,CAAC,CAAC,CAAC,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,WAAW,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,IAAI,WAAW,CAAC,CAAC,CAAE,CAAC;QAC1E,CAAC,CAAC,CAAC,CAAC;IACR,MAAM,aAAa,GACjB,WAAW,CAAC,MAAM,GAAG,CAAC;QACpB,CAAC,CAAC,CAAC,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,WAAW,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC;YACpD,WAAW,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,CAAE,CAAC;QACvC,CAAC,CAAC,CAAC,CAAC;IAER,iDAAiD;IACjD,MAAM,mBAAmB,GAA2B,EAAE,CAAC;IACvD,KAAK,MAAM,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,MAAM,CAAC,OAAO,CAAC,kBAAkB,CAAC,EAAE,CAAC;QAC7D,MAAM,IAAI,GAAG,UAAU,CAAC,MAAM,CAC5B,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,GAAG,CAAC,CAChD,CAAC,MAAM,CAAC;QACT,mBAAmB,CAAC,IAAI,CAAC,GAAG,IAAI,GAAG,UAAU,CAAC,MAAM,CAAC;IACvD,CAAC;IAED,gDAAgD;IAChD,MAAM,IAAI,GAA4C,EAAE,CAAC;IACzD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,IAAI,MAAM,EAAE,GAAG,IAAI,IAAI,EAAE,CAAC;QAC7C,MAAM,IAAI,GAAG,UAAU,CAAC,MAAM,CAC5B,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,GAAG,CAAC,CAChD,CAAC,MAAM,CAAC;QACT,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,IAAI,GAAG,UAAU,CAAC,MAAM,EAAE,CAAC,CAAC;IACxD,CAAC;IAED,0CAA0C;IAC1C,qEAAqE;IACrE,MAAM,MAAM,GAAG,MAAM,CAAC,WAAW,IAAI,KAAK,CAAC;IAC3C,MAAM,qBAAqB,GACzB,UAAU,CAAC,MAAM,CACf,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,CAAC,EAC3D,CAAC,CACF,GAAG,UAAU,CAAC,MAAM,CAAC;IACxB,MAAM,mBAAmB,GAAG,CAAC,CAAC,CAAC,sCAAsC;IACrE,MAAM,SAAS,GACb,UAAU,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE;QAC5B,MAAM,OAAO,GAAG,MAAM,GAAG,EAAE,CAAC,UAAU,GAAG,EAAE,CAAC,UAAU,CAAC;QACvD,OAAO,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,mBAAmB,CAAC,CAAC;IACtD,CAAC,EAAE,CAAC,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC;IAC5B,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,qBAAqB,EAAE,mBAAmB,CAAC,CAAC;IACvE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,GAAG,SAAS,CAAC,CAAC;IAEhD,oDAAoD;IACpD,MAAM,UAAU,GACd,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;IAChE,MAAM,YAAY,GAAG,YAAY,CAAC,SAAS,EAAE,MAAM,EAAE,UAAU,CAAC,CAAC;IAEjE,MAAM,cAAc,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAEhD,OAAO;QACL,UAAU;QACV,SAAS;QACT,aAAa;QACb,aAAa;QACb,mBAAmB;QACnB,IAAI;QACJ,IAAI;QACJ,KAAK,EAAE,YAAY;QACnB,cAAc;KACf,CAAC;AACJ,CAAC"}

package/dist/models/survivalFitting.d.ts

@@ -0,0 +1,47 @@
+/**
+ * Survival Curve Fitting — fit parametric distributions to Kaplan-Meier data.
+ *
+ * Supports: Exponential, Weibull, Log-logistic, Log-normal, Gompertz.
+ * Uses maximum likelihood estimation via Newton-Raphson optimization.
+ * Model selection via AIC/BIC.
+ *
+ * References:
+ * - Latimer NR. NICE DSU TSD 14: Survival analysis (2013)
+ * - Collett D. Modelling Survival Data in Medical Research (2015)
+ */
+export interface KMDataPoint {
+    time: number;
+    survival: number;
+    n_at_risk?: number;
+    n_events?: number;
+}
+export type DistributionName = "exponential" | "weibull" | "log_logistic" | "log_normal" | "gompertz";
+export interface FittedDistribution {
+    name: DistributionName;
+    params: Record<string, number>;
+    aic: number;
+    bic: number;
+    log_likelihood: number;
+    survival_at: (t: number) => number;
+    hazard_at: (t: number) => number;
+    median_survival: number;
+    mean_survival_restricted: number;
+}
+export interface SurvivalFitResult {
+    fits: FittedDistribution[];
+    best_aic: FittedDistribution;
+    best_bic: FittedDistribution;
+    km_data: KMDataPoint[];
+    time_unit: string;
+    extrapolations: Array<{
+        time: number;
+        km_observed?: number;
+        exponential: number;
+        weibull: number;
+        log_logistic: number;
+        log_normal: number;
+        gompertz: number;
+    }>;
+}
+export declare function fitSurvivalCurves(data: KMDataPoint[], timeUnit?: string): SurvivalFitResult;
+//# sourceMappingURL=survivalFitting.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"survivalFitting.d.ts","sourceRoot":"","sources":["../../src/models/survivalFitting.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH,MAAM,WAAW,WAAW;IAC1B,IAAI,EAAE,MAAM,CAAC;IACb,QAAQ,EAAE,MAAM,CAAC;IACjB,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,EAAE,MAAM,CAAC;CACnB;AAED,MAAM,MAAM,gBAAgB,GACxB,aAAa,GACb,SAAS,GACT,cAAc,GACd,YAAY,GACZ,UAAU,CAAC;AAEf,MAAM,WAAW,kBAAkB;IACjC,IAAI,EAAE,gBAAgB,CAAC;IACvB,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/B,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;IACZ,cAAc,EAAE,MAAM,CAAC;IACvB,WAAW,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,CAAC;IACnC,SAAS,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,CAAC;IACjC,eAAe,EAAE,MAAM,CAAC;IACxB,wBAAwB,EAAE,MAAM,CAAC;CAClC;AAED,MAAM,WAAW,iBAAiB;IAChC,IAAI,EAAE,kBAAkB,EAAE,CAAC;IAC3B,QAAQ,EAAE,kBAAkB,CAAC;IAC7B,QAAQ,EAAE,kBAAkB,CAAC;IAC7B,OAAO,EAAE,WAAW,EAAE,CAAC;IACvB,SAAS,EAAE,MAAM,CAAC;IAClB,cAAc,EAAE,KAAK,CAAC;QACpB,IAAI,EAAE,MAAM,CAAC;QACb,WAAW,CAAC,EAAE,MAAM,CAAC;QACrB,WAAW,EAAE,MAAM,CAAC;QACpB,OAAO,EAAE,MAAM,CAAC;QAChB,YAAY,EAAE,MAAM,CAAC;QACrB,UAAU,EAAE,MAAM,CAAC;QACnB,QAAQ,EAAE,MAAM,CAAC;KAClB,CAAC,CAAC;CACJ;AAyXD,wBAAgB,iBAAiB,CAC/B,IAAI,EAAE,WAAW,EAAE,EACnB,QAAQ,GAAE,MAAiB,GAC1B,iBAAiB,CAiDnB"}