formlab-mcp 0.6.33 → 0.6.35

package/README.md

@@ -51,7 +51,7 @@ Ask Claude (or another MCP client) things like:
 | `compare_batches` | Reproducibility of 2+ runs (ideally of one formula): each run's **yield %**, actual produced mass, cost/kg; per-ingredient **drift** across the runs vs the formula's proposed wt-%; and the biggest **outlier** run. `basis`: `wt_percent` (default) or `amount` |
 | `get_batch_pivot` | Production analytics — aggregate batches by `group_by` (formula / status / month / prepared_by / project) × `metric` (count, avg_yield_pct, total_produced_kg, avg_cost_per_batch, total_samples, total_tests), with share + total for additive metrics. Optional `batch_uids` scope + `status` filter |
 | `get_project_pivot` | Portfolio analytics — aggregate projects by `group_by` (status / phase / priority / business_unit / site / customer / lead, or `cf:<custom field>`) × `metric` (count, total_formulas / batches / samples, total_cost, avg_formulas / batches / cost_per_batch), rolling child activity + cost up by any project attribute. Optional `project_uids` scope |
-| `get_stability` | Grounded stability / shelf-life analysis for a formulation, sample or batch × parameter (auto-detected): time-series, drift (per-month slope + R²), an I-chart (mean ±3σ + out-of-control points), spec status, and an ICH-Q1E-flavored projected shelf life (point + 95%-CI crossing), split by storage condition |
+| `get_stability` | Grounded stability / shelf-life analysis for a formulation, sample or batch × parameter (auto-detected): time-series, drift (per-month slope + R²), an I-chart (mean ±3σ + out-of-control points), spec status, and an ICH-Q1E-flavored projected shelf life (point + 95%-CI crossing), all on **one point per test date per run condition** (replicates averaged, never separate time points), split by storage condition |
 | `get_design_space` | COMPOSITION design space for one project: which ingredients vary and over what observed wt-% ranges (ranked), an occupied-region summary ("where you've been"), the biggest untested **interior gap** as ready-to-seed DOE ranges (a combination you could have made but skipped), and an optional standardised **2-component PCA** of all varying ingredients — loadings, variance-explained, scree, and a full-dimensional gap. Matches the in-app Design Space Viewer exactly. Answers *"where haven't I explored?"* / *"what should I formulate next?"* |
 | `list_doe_designs` | Saved DOE designs — the *recipe* behind each batch of runs: design type, factors + ranges, constraints, run count, D-efficiency. Filters: `project`, `design_type`, `constrained_only`, `name_contains` |
 | `get_doe_design` | One design in full: constraints in plain language, the model it was optimised for, generation settings (run budget, replicates, centre points, seed) and the complete run matrix with the formulation each run became |

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "formlab-mcp",
-  "version": "0.6.33",
+  "version": "0.6.35",
   "mcpName": "io.github.juliu1980/formlab-mcp",
   "description": "Read-only Model Context Protocol server for FormLab \u2014 lets Claude (and other MCP clients) read and analyze your FormLab data, from a local export file OR your live cloud workspace.",
   "type": "module",

package/server.json

@@ -2,7 +2,7 @@
   "$schema": "https://static.modelcontextprotocol.io/schemas/2025-12-11/server.schema.json",
   "name": "io.github.juliu1980/formlab-mcp",
   "description": "Read-only MCP for FormLab — let Claude query your formulation lab. Free reads a local export; Pro connects to your live cloud workspace with a dedicated read-only token.",
-  "version": "0.6.33",
+  "version": "0.6.35",
   "websiteUrl": "https://formvix.com/mcp",
   "repository": {
     "url": "https://github.com/juliu1980/FormLab",
@@ -13,7 +13,7 @@
     {
       "registryType": "npm",
       "identifier": "formlab-mcp",
-      "version": "0.6.33",
+      "version": "0.6.35",
       "runtimeHint": "npx",
       "transport": {
         "type": "stdio"

package/tools/analytics.js

@@ -226,7 +226,7 @@ const get_test_result = {
 const get_doe_matrix = {
   definition: {
     name: 'get_doe_matrix',
-    description: 'Build a Design-of-Experiments style pivot matrix: rows = formulations (or samples / batches / tests), columns = composition wt-% per ingredient + aggregated test-parameter values. Same shape as the FormLab DOE Matrix view. A parameter measured under two or more RUN CONDITIONS (e.g. storage 25 °C vs 40 °C) is one column per condition — "pH · 40 °C" — and the bare parameter column is the base arm (readings with no condition, else the condition with the most readings); they are different responses, never averaged together. Returns CSV when format=csv (default) for compact LLM consumption, or structured JSON when format=json.',
+    description: 'Build a Design-of-Experiments style pivot matrix: rows = formulations (or samples / batches / tests), columns = composition wt-% per ingredient + aggregated test-parameter values. Same shape as the FormLab DOE Matrix view. A parameter measured under two or more RUN CONDITIONS (e.g. storage 25 °C vs 40 °C) is one column per condition — "pH · 40 °C" — and the bare parameter column is the base arm (readings with no condition, else the condition with the most readings, on a tie the one declared as the default on the Test Method / Test Panel); they are different responses, never averaged together. Returns CSV when format=csv (default) for compact LLM consumption, or structured JSON when format=json.',
     inputSchema: {
       type: 'object',
       properties: {
@@ -453,7 +453,34 @@ function _mcpCondTags(conds) {
   const lead = tags.map(t => t.split(/\s+[\/·]\s+/)[0]);
   return (new Set(lead).size === lead.length) ? lead : tags;
 }
-// items: [{ parameter, m }] → Map(parameter → Map(condSig → column name)),
+// Is this reading taken under the DEFAULT run condition its Test Method or
+// Test Panel declares? Port of the app's flCondDefaultScore (js/utils.js):
+// the tiebreak for the primary arm when two arms have equally long records,
+// instead of the condition's name ("25 °C" sorting before "40 °C").
+function _mcpCondDefaultScore(m, t, db) {
+  const c = m && m.conditions;
+  if (!c || typeof c !== 'object') return 0;
+  const decls = [];
+  const ps = db.parameters || [];
+  const nm = String(m.parameter || '').trim().toLowerCase();
+  const method = (m.parameterId && ps.find(p => p && p.id === m.parameterId)) || (nm ? ps.find(p => p && String(p.name || '').trim().toLowerCase() === nm) : null);
+  if (method && Array.isArray(method.conditions)) decls.push(...method.conditions);
+  const tmpl = t && t.templateId ? (db.templates || []).find(x => x && x.id === t.templateId) : null;
+  if (tmpl && Array.isArray(tmpl.conditions)) decls.push(...tmpl.conditions);
+  const lc = {};
+  Object.keys(c).forEach(k => { lc[k.toLowerCase()] = c[k]; });
+  let matched = 0;
+  for (const d of decls) {
+    if (!d || !d.key || d.default == null || d.default === '') continue;
+    const v = lc[String(d.key).toLowerCase()];
+    if (v == null || v === '') continue;
+    if (String(v).trim() !== String(d.default).trim()) return 0;
+    matched++;
+  }
+  return matched ? 1 : 0;
+}
+
+// items: [{ parameter, m, test? }] → Map(parameter → Map(condSig → column name)),
 // only for parameters with two or more arms.
 function _mcpArmNames(items) {
   const arms = new Map();
@@ -461,14 +488,26 @@ function _mcpArmNames(items) {
     const sig = _mcpCondSig(it.m);
     if (!arms.has(it.parameter)) arms.set(it.parameter, new Map());
     const a = arms.get(it.parameter);
-    const cur = a.get(sig) || { n: 0, m: it.m };
+    const cur = a.get(sig) || { n: 0, m: it.m, t: it.test || null };
     cur.n++; a.set(sig, cur);
   }
   const out = new Map();
+  let testOf = null;   // measurement → its report, built only if a tie needs a panel default
+  const score = (arm) => {
+    const { db } = getStore();
+    let t = arm.t;
+    if (!t) {
+      if (!testOf) { testOf = new Map(); (db.testResults || []).forEach(r => (r.measurements || []).forEach(x => testOf.set(x, r))); }
+      t = testOf.get(arm.m) || null;
+    }
+    return _mcpCondDefaultScore(arm.m, t, db);
+  };
   arms.forEach((a, param) => {
     if (a.size < 2) return;
     const keys = [...a.keys()];
-    const base = a.has('') ? '' : keys.slice().sort((x, y) => (a.get(y).n - a.get(x).n) || x.localeCompare(y))[0];
+    // No condition wins outright; among conditioned arms, most readings, then
+    // the declared default condition, then by name — as in the app.
+    const base = a.has('') ? '' : keys.slice().sort((x, y) => (a.get(y).n - a.get(x).n) || (score(a.get(y)) - score(a.get(x))) || x.localeCompare(y))[0];
     const ordered = [base, ...keys.filter(k => k !== base).sort((x, y) => x.localeCompare(y))];
     const tags = _mcpCondTags(ordered.map(k => k ? _mcpCondLabel(a.get(k).m) : ''));
     out.set(param, new Map(ordered.map((k, i) => [k, i === 0 ? param : `${param} · ${tags[i]}`])));
@@ -1088,7 +1127,7 @@ function _shelfLifeAssess(sl, targetMonths, now) {
 const get_stability = {
   definition: {
     name: 'get_stability',
-    description: 'Grounded stability / shelf-life analysis for one entity × parameter — the deterministic numbers behind the Stability & Trends view. Pass a formulation (pools all its samples), a sample (that sample over time), or a batch (pools its samples) by id or UID, plus a parameter name. Returns the time series, drift (per-month slope + R²), an individuals control chart (mean ± 3σ + out-of-control points), spec status, and an ICH-flavored projected shelf life (linear extrapolation to the spec crossing: point estimate + the conservative 95%-confidence-bound crossing). When measurements carry run conditions (e.g. Storage 25 °C vs 40 °C) it splits into one series per condition — a true accelerated-vs-ambient comparison. Use these figures directly; do not re-derive the regression.',
+    description: 'Grounded stability / shelf-life analysis for one entity × parameter — the deterministic numbers behind the Stability & Trends view. Pass a formulation (pools all its samples), a sample (that sample over time), or a batch (pools its samples) by id or UID, plus a parameter name. Returns the time series, drift (per-month slope + R²), an individuals control chart (mean ± 3σ + out-of-control points), spec status, and an ICH-flavored projected shelf life (linear extrapolation to the spec crossing: point estimate + the conservative 95%-confidence-bound crossing). When measurements carry run conditions (e.g. Storage 25 °C vs 40 °C) it splits into one series per condition — a true accelerated-vs-ambient comparison. A point is one TEST DATE: every reading taken that day under that condition is averaged (points = test dates, readings = the raw count), so replicates never count as separate time points in the drift, control chart or shelf-life fit — the same points as the app. Use these figures directly; do not re-derive the regression.',
     inputSchema: {
       type: 'object',
       properties: {
@@ -1143,7 +1182,7 @@ const get_stability = {
         }
         const cond = (m.conditions && Object.keys(m.conditions).length)
           ? Object.entries(m.conditions).map(([k, v]) => `${k} ${v}`).join(' · ') : '';
-        points.push({ date, value, cond, spec: sp });
+        points.push({ date, day: dStr, value, cond, spec: sp });
       });
     });
     const specList = [...specByKey.values()].map(e => ({ min: e.min, max: e.max, readings: e.readings, sources: [...e.sources] }));
@@ -1154,7 +1193,21 @@ const get_stability = {
       const pdef = (db.parameters || []).find(p => (p.name || '').trim().toLowerCase() === paramLc);
       if (pdef && pdef.acceptanceCriteria) { specMin = pdef.acceptanceCriteria.min ?? null; specMax = pdef.acceptanceCriteria.max ?? null; }
     }
-    points.sort((a, b) => a.date - b.date);
+    // ONE POINT PER TEST DATE PER RUN CONDITION — the mean of every reading
+    // taken that day under that condition — exactly as the app's Stability view
+    // (FormLab build 24k). Replicates are not separate time points: counting a
+    // triplicate as three consecutive observations took the control chart's
+    // sigma from within-day differences (Clarifying Shampoo's ambient pH: 10
+    // out-of-control points on 53 readings, 1 on its 18 test dates) and gave
+    // the drift and shelf-life fits degrees of freedom the study does not have.
+    const byDay = new Map();
+    points.forEach(p => { const k = p.day + '\u0000' + p.cond; if (!byDay.has(k)) byDay.set(k, []); byDay.get(k).push(p); });
+    points.length = 0;
+    byDay.forEach(g => {
+      const value = g.reduce((a, r) => a + r.value, 0) / g.length;
+      points.push({ date: new Date(g[0].day), value, cond: g[0].cond, spec: (g.find(r => r.spec) || {}).spec || null, readings: g.length });
+    });
+    points.sort((a, b) => (a.date - b.date) || a.cond.localeCompare(b.cond));
     if (!points.length) return { error: `No dated "${args.parameter}" measurements on ${entity.uid || entity.id}.` };
 
     const spec = { min: specMin, max: specMax };
@@ -1169,7 +1222,7 @@ const get_stability = {
     const analyze = (pts) => {
       const ys = pts.map(p => p.value);
       const n = ys.length;
-      const out = { points: n, from: pts[0].date.toISOString().slice(0, 10), to: pts[n - 1].date.toISOString().slice(0, 10), values: ys.map(v => _round(v, 6)), latest: _round(ys[n - 1], 6) };
+      const out = { points: n, readings: pts.reduce((a, p) => a + (p.readings || 1), 0), from: pts[0].date.toISOString().slice(0, 10), to: pts[n - 1].date.toISOString().slice(0, 10), values: ys.map(v => _round(v, 6)), latest: _round(ys[n - 1], 6) };
       if (n < 2) { out.note = 'Fewer than 2 points — no trend/shelf-life.'; return out; }
       // Drift on real elapsed months (not point order), like the shelf-life
       // projection and the app's Stability view: uneven ICH schedules
@@ -1229,7 +1282,7 @@ const get_stability = {
   },
 };
 
-export { _evalPassFail, _resolveSpec, _shelfLife, _shelfLifeAssess, _mcpCondTags, _mcpArmNames };
+export { _evalPassFail, _resolveSpec, _shelfLife, _shelfLifeAssess, _mcpCondTags, _mcpArmNames, _mcpCondDefaultScore };
 export const tools = {
   list_test_results,
   get_test_result,