ml-time-graph 1.0.0 → 1.0.5

package/dist/analyze/index.js

@@ -1,574 +1,15 @@
-// src/analyze/processor.ts
-var SeriesProcessor = class {
-  /**
-   * Standard interpolation for scalar DataPoints.
-   */
-  static interpolateDataPoint(p1, p2, t) {
-    return {
-      time: p1.time + t * (p2.time - p1.time),
-      value: (p1.value ?? 0) + t * ((p2.value ?? 0) - (p1.value ?? 0))
-    };
-  }
-  /**
-   * Standard interpolation for AggregatedPoints (interpolates min, max, avg and count).
-   */
-  static interpolateAggregatedPoint(p1, p2, t) {
-    const lerp = (v1, v2) => v1 !== null && v2 !== null ? v1 + t * (v2 - v1) : null;
-    return {
-      time: p1.time + t * (p2.time - p1.time),
-      min: lerp(p1.min, p2.min),
-      max: lerp(p1.max, p2.max),
-      avg: lerp(p1.avg, p2.avg),
-      count: Math.round(p1.count + t * (p2.count - p1.count))
-    };
-  }
-  /**
-   * Splits a data array into contiguous runs based on null values or time jumps.
-   *
-   * @param data The raw data points.
-   * @param isNull A predicate to identify "gap" points (e.g. value === null).
-   * @param gapThreshold Max time distance between points before a new run starts.
-   */
-  static getRuns(data, isNull, gapThreshold = 0) {
-    const sorted = [...data].sort((a, b) => a.time - b.time);
-    const runs = [];
-    let current = [];
-    let prev = null;
-    for (const p of sorted) {
-      const isPointNull = isNull(p);
-      const isJump = gapThreshold > 0 && prev && p.time - prev.time > gapThreshold;
-      if (isPointNull || isJump) {
-        if (current.length) {
-          runs.push(current);
-          current = [];
-        }
-      }
-      if (!isPointNull) {
-        current.push(p);
-      }
-      prev = p;
-    }
-    if (current.length) {
-      runs.push(current);
-    }
-    return runs;
-  }
-  /**
-   * Splits a contiguous run into sub-segments at the given boundary values.
-   * Inserts interpolated points at every boundary crossing so segments
-   * meet exactly at the boundary.
-   *
-   * @param run A gap-free array of points.
-   * @param boundaries Values at which to split the run.
-   * @param getValue Function to extract the numeric value used for splitting.
-   * @param interpolate Function to create an interpolated point between p1 and p2 at factor t [0..1].
-   */
-  static splitByBoundaries(run, boundaries, getValue, interpolate) {
-    if (run.length === 0) return [];
-    if (boundaries.length === 0) {
-      return [{ data: run, zoneIndex: 0 }];
-    }
-    const bs = [...boundaries].sort((a, b) => a - b);
-    const out = [];
-    const getZone = (v) => {
-      let idx = 0;
-      for (let i = 0; i < bs.length; i++) {
-        if (v >= bs[i]) idx = i + 1;
-        else break;
-      }
-      return idx;
-    };
-    let currentSeg = [run[0]];
-    for (let i = 1; i < run.length; i++) {
-      const p1 = run[i - 1];
-      const p2 = run[i];
-      const v1 = getValue(p1);
-      const v2 = getValue(p2);
-      let crossed;
-      if (v2 > v1) {
-        crossed = bs.filter((b) => b > v1 && b <= v2);
-      } else if (v2 < v1) {
-        crossed = bs.filter((b) => b >= v2 && b < v1).reverse();
-      } else {
-        crossed = [];
-      }
-      for (const b of crossed) {
-        const t = (b - v1) / (v2 - v1);
-        const pInt = interpolate(p1, p2, t);
-        currentSeg.push(pInt);
-        out.push({ data: currentSeg, zoneIndex: getZone((v1 + b) / 2) });
-        currentSeg = [pInt];
-      }
-      currentSeg.push(p2);
-    }
-    if (currentSeg.length > 0) {
-      const vStart = getValue(currentSeg[0]);
-      const vEnd = getValue(currentSeg[currentSeg.length - 1]);
-      out.push({ data: currentSeg, zoneIndex: getZone((vStart + vEnd) / 2) });
-    }
-    return out;
-  }
-  /**
-   * Splits a contiguous run into two groups: those below and those at/above a threshold.
-   * Internally uses splitByBoundaries to ensure exact intersection points.
-   */
-  static splitByThreshold(run, threshold, getValue, interpolate) {
-    const segments = this.splitByBoundaries(run, [threshold], getValue, interpolate);
-    const result = { above: [], below: [] };
-    for (const seg of segments) {
-      if (seg.zoneIndex === 0) result.below.push(seg.data);
-      else result.above.push(seg.data);
-    }
-    return result;
-  }
-};
-
-// src/analyze/aggregator.ts
-var MS_PER_HOUR = 36e5;
-var MS_PER_DAY = 864e5;
-var PRODUCT_PROFILES = {
-  pharma_cold: { limitLow: 2, limitHigh: 8, activationEnergy: 83144 },
-  pharma_ambient: {
-    limitLow: 15,
-    limitHigh: 25,
-    activationEnergy: 83144
-  },
-  blood: { limitLow: 2, limitHigh: 6, activationEnergy: 83144 },
-  food_chilled: { limitLow: 0, limitHigh: 4 },
-  // Lebensmittel nutzen MKT seltener, dH entfällt meist
-  frozen: { limitLow: -40, limitHigh: -18 }
-};
-var R = 8.314462618;
-var CELSIUS_TO_KELVIN = 273.15;
-function createAggrWithMkt(time, rawPointsInSlot, thresholds, minCountForMkt = 1) {
-  const validPoints = rawPointsInSlot.filter(
-    (pt) => pt.value !== null && pt.value !== void 0 && !isNaN(pt.value)
-  );
-  const values = validPoints.map((pt) => pt.value);
-  if (!values.length) {
-    return {
-      time,
-      min: null,
-      max: null,
-      avg: null,
-      mkt: null,
-      deltaMkt: null,
-      stdDev: null,
-      minutesAboveHigh: thresholds ? 0 : null,
-      minutesBelowLow: thresholds ? 0 : null,
-      count: 0
-    };
-  }
-  const mn = Math.min(...values);
-  const mx = Math.max(...values);
-  const avg = values.reduce((a, b) => a + b, 0) / values.length;
-  let minutesAboveHigh = thresholds ? 0 : null;
-  let minutesBelowLow = thresholds ? 0 : null;
-  if (thresholds) {
-    for (let i = 0; i < rawPointsInSlot.length; i++) {
-      const current = rawPointsInSlot[i];
-      if (current.value !== null && current.value !== void 0) {
-        const next = rawPointsInSlot[i + 1];
-        const durationMs = next ? next.time - current.time : 0;
-        const durationMinutes = durationMs / (1e3 * 60);
-        if (current.value > thresholds.limitHigh)
-          minutesAboveHigh += durationMinutes;
-        if (current.value < thresholds.limitLow)
-          minutesBelowLow += durationMinutes;
-      }
-    }
-  }
-  let sumSquareDiff = 0;
-  let sumExp = 0;
-  const deltaH = thresholds?.activationEnergy ?? 83144;
-  const dh_r = deltaH / R;
-  for (const cels of values) {
-    sumSquareDiff += Math.pow(cels - avg, 2);
-    const kelvin = cels + CELSIUS_TO_KELVIN;
-    sumExp += Math.exp(-dh_r / kelvin);
-  }
-  const variance = sumSquareDiff / values.length;
-  const stdDev2 = Math.sqrt(variance);
-  let mktCelsius = null;
-  let deltaMkt = null;
-  if (values.length >= minCountForMkt) {
-    const meanExp = sumExp / values.length;
-    const mktKelvin = dh_r / -Math.log(meanExp);
-    mktCelsius = Number((mktKelvin - CELSIUS_TO_KELVIN).toFixed(2));
-    deltaMkt = Number((mktCelsius - avg).toFixed(2));
-  }
-  return {
-    time,
-    min: Number(mn.toFixed(2)),
-    max: Number(mx.toFixed(2)),
-    avg: Number(avg.toFixed(2)),
-    mkt: mktCelsius,
-    deltaMkt,
-    stdDev: Number(stdDev2.toFixed(2)),
-    minutesAboveHigh: minutesAboveHigh !== null ? Number(minutesAboveHigh.toFixed(1)) : null,
-    minutesBelowLow: minutesBelowLow !== null ? Number(minutesBelowLow.toFixed(1)) : null,
-    count: values.length
-  };
-}
-function aggregateBySlot(data, mode = "hourly", customInterval, thresholds, minCountForMkt = 1) {
-  if (!data.length) return [];
-  const interval = mode === "hourly" ? MS_PER_HOUR : mode === "daily" ? MS_PER_DAY : customInterval ?? 60 * 60 * 1e3;
-  const sorted = [...data].sort((a, b) => a.time - b.time);
-  const slots = [];
-  let slotStart = sorted[0].time;
-  let pointsInSlot = [];
-  for (const pt of sorted) {
-    while (pt.time - slotStart >= interval) {
-      slots.push(
-        createAggrWithMkt(slotStart, pointsInSlot, thresholds, minCountForMkt)
-      );
-      slotStart += interval;
-      pointsInSlot = [];
-    }
-    pointsInSlot.push(pt);
-  }
-  slots.push(
-    createAggrWithMkt(slotStart, pointsInSlot, thresholds, minCountForMkt)
-  );
-  return slots;
-}
-function createAggr(time, values) {
-  const mn = Math.min(...values);
-  const mx = Math.max(...values);
-  const avg = values.reduce((a, b) => a + b, 0) / values.length;
-  return { time, min: mn, max: mx, avg, count: values.length };
-}
-function detectGaps(data, minGapMs = 6e4) {
-  if (data.length < 2) return [];
-  const sorted = [...data].sort((a, b) => a.time - b.time);
-  const gaps = [];
-  for (let i = 1; i < sorted.length; i++) {
-    const gap = sorted[i].time - sorted[i - 1].time;
-    if (gap > minGapMs) {
-      gaps.push({
-        startTime: sorted[i - 1].time,
-        endTime: sorted[i].time
-      });
-    }
-  }
-  return gaps;
-}
-function analyzeLongTermTrends(aggregatedData) {
-  const insights = [];
-  let totalMinutesTooHot = 0;
-  let maxStdDev = 0;
-  let slotsWithData = 0;
-  let highDeltaMktCount = 0;
-  for (const slot of aggregatedData) {
-    if (slot.count > 0) {
-      slotsWithData++;
-      if (slot.minutesAboveHigh) totalMinutesTooHot += slot.minutesAboveHigh;
-      if (slot.stdDev && slot.stdDev > maxStdDev) maxStdDev = slot.stdDev;
-      if (slot.deltaMkt && slot.deltaMkt > 2) {
-        highDeltaMktCount++;
-      }
-    }
-  }
-  if (slotsWithData === 0) return [];
-  if (totalMinutesTooHot > 120) {
-    insights.push({
-      type: "critical",
-      message: `Kritische Gesamtbelastung: Das Produkt war im gesamten Zeitraum insgesamt ${Math.round(totalMinutesTooHot)} Minuten zu warm.`,
-      metric: "Total Excursion Time"
-    });
-  }
-  if (highDeltaMktCount > slotsWithData * 0.15) {
-    insights.push({
-      type: "warning",
-      message: "H\xE4ufige thermische Schocks erkannt. Die kinetische Temperatur weicht oft stark vom Schnitt ab (m\xF6gliche regelm\xE4\xDFige T\xFCr\xF6ffnungen).",
-      metric: "Delta MKT Instability"
-    });
-  }
-  if (maxStdDev > 4) {
-    insights.push({
-      type: "warning",
-      message: `Hohe Instabilit\xE4t gemessen. Die maximale Standardabweichung lag bei ${maxStdDev}\xB0C. Das K\xFChlsystem regelt unsauber.`,
-      metric: "Standard Deviation Peak"
-    });
-  }
-  return insights;
-}
-function downsample(data, target) {
-  if (data.length <= target) return data;
-  const result = [];
-  const step = Math.ceil(data.length / target);
-  result.push(data[0]);
-  for (let i = 1; i < data.length; i += step) {
-    result.push(data[Math.min(i, data.length - 1)]);
-  }
-  return result;
-}
-function aggregateWithStats(data, mode = "hourly", opts) {
-  const thresholds = {
-    limitLow: opts.low,
-    limitHigh: opts.high,
-    activationEnergy: opts.activationEnergy
-  };
-  return aggregateBySlot(
-    data,
-    mode,
-    opts.interval,
-    thresholds,
-    opts.minCountForMkt ?? 1
-  );
-}
-
-// src/analyze/stats.ts
-var StatsAggregator = class {
-  /** Compute statistics from data points */
-  static compute(data) {
-    const values = data.map((d) => d.value).filter((v) => v !== null).sort((a, b) => a - b);
-    if (values.length === 0) {
-      return { min: NaN, max: NaN, avg: NaN, mean: NaN, median: NaN, stdDev: NaN, count: 0 };
-    }
-    const count = values.length;
-    const sum = values.reduce((a, b) => a + b, 0);
-    const mean = sum / count;
-    const median = count % 2 === 1 ? values[Math.floor(count / 2)] : (values[count / 2 - 1] + values[count / 2]) / 2;
-    const variance = values.reduce((a, b) => a + Math.pow(b - mean, 2), 0) / count;
-    const stdDev2 = Math.sqrt(variance);
-    return {
-      min: values[0],
-      max: values[count - 1],
-      avg: mean,
-      mean,
-      median,
-      stdDev: stdDev2,
-      count
-    };
-  }
-  /** Compute stats for a specific time range (viewport-scoped) */
-  static computeInRange(data, startTime, endTime) {
-    const filtered = data.filter((d) => d.time >= startTime && d.time <= endTime);
-    return this.compute(filtered);
-  }
-};
-
-// src/analyze/moving_avg.ts
-var MovingAvg = class _MovingAvg {
-  #data;
-  #windowSize;
-  #type;
-  #timeScale;
-  #valueScale;
-  #stroke;
-  #strokeWidth;
-  constructor(config) {
-    this.#data = [...config.data].sort((a, b) => a.time - b.time);
-    this.#windowSize = config.windowSize;
-    this.#type = config.type ?? "simple";
-    this.#timeScale = config.timeScale;
-    this.#valueScale = config.valueScale;
-    this.#stroke = config.stroke ?? "#e53e3e";
-    this.#strokeWidth = config.strokeWidth ?? 2;
-  }
-  /** Compute simple moving average points */
-  static simple(data, windowSize) {
-    const sorted = [...data].filter((d) => d.value !== null).sort((a, b) => a.time - b.time);
-    const result = [];
-    for (let i = 0; i < sorted.length; i++) {
-      const start = Math.max(0, i - windowSize + 1);
-      const window = sorted.slice(start, i + 1);
-      const avg = window.reduce((s, d) => s + d.value, 0) / window.length;
-      result.push({ time: sorted[i].time, value: avg });
-    }
-    return result;
-  }
-  /** Compute exponential moving average (EMA) */
-  static exponential(data, windowSize) {
-    const sorted = [...data].filter((d) => d.value !== null).sort((a, b) => a.time - b.time);
-    if (sorted.length === 0) return [];
-    const multiplier = 2 / (windowSize + 1);
-    const result = [];
-    let ema = sorted[0].value;
-    result.push({ time: sorted[0].time, value: ema });
-    for (let i = 1; i < sorted.length; i++) {
-      ema = sorted[i].value * multiplier + ema * (1 - multiplier);
-      result.push({ time: sorted[i].time, value: ema });
-    }
-    return result;
-  }
-  /** Get computed moving average as pixel points */
-  points() {
-    const raw = this.#type === "exponential" ? _MovingAvg.exponential(this.#data, this.#windowSize) : _MovingAvg.simple(this.#data, this.#windowSize);
-    return raw.map((p) => ({
-      x: this.#timeScale.map(p.time),
-      y: this.#valueScale.map(p.value)
-    }));
-  }
-  /** Render MA as a path command */
-  render() {
-    const pts = this.points();
-    if (pts.length < 2)
-      return pts.map((p) => ({
-        type: "circle",
-        cx: p.x,
-        cy: p.y,
-        r: this.#strokeWidth,
-        fill: this.#stroke
-      }));
-    return [
-      {
-        type: "path",
-        points: pts,
-        smoothing: true,
-        stroke: this.#stroke,
-        strokeWidth: this.#strokeWidth
-      }
-    ];
-  }
-};
-
-// src/analyze/mkt.ts
-var R2 = 8.314;
-var C_TO_K = 273.15;
-var DEFAULT_ACTIVATION_ENERGY = 83144;
-function mkt(samples, activationEnergy = DEFAULT_ACTIVATION_ENERGY) {
-  const dhR = activationEnergy / R2;
-  let sumExp = 0;
-  let count = 0;
-  for (const v of samples) {
-    if (v === null) continue;
-    sumExp += Math.exp(-dhR / (v + C_TO_K));
-    count++;
-  }
-  if (count === 0) return null;
-  const meanExp = sumExp / count;
-  const mktKelvin = dhR / -Math.log(meanExp);
-  return mktKelvin - C_TO_K;
-}
-function rollingMkt(data, windowMs, activationEnergy = DEFAULT_ACTIVATION_ENERGY) {
-  const out = new Array(data.length);
-  let left = 0;
-  for (let i = 0; i < data.length; i++) {
-    const t = data[i].time;
-    const wStart = t - windowMs;
-    while (left < i && data[left].time < wStart) left++;
-    const slice = data.slice(left, i + 1).map((p) => p.value);
-    const value = mkt(slice, activationEnergy);
-    out[i] = { time: t, value, synthetic: true };
-  }
-  return out;
-}
-
-// src/analyze/std_dev.ts
-function stdDev(samples) {
-  let sum = 0;
-  let count = 0;
-  for (const v of samples) {
-    if (v === null) continue;
-    sum += v;
-    count++;
-  }
-  if (count === 0) return null;
-  const mean = sum / count;
-  let sqSum = 0;
-  for (const v of samples) {
-    if (v === null) continue;
-    const d = v - mean;
-    sqSum += d * d;
-  }
-  return Math.sqrt(sqSum / count);
-}
-function sampleStdDev(samples) {
-  let sum = 0;
-  let count = 0;
-  for (const v of samples) {
-    if (v === null) continue;
-    sum += v;
-    count++;
-  }
-  if (count < 2) return null;
-  const mean = sum / count;
-  let sqSum = 0;
-  for (const v of samples) {
-    if (v === null) continue;
-    const d = v - mean;
-    sqSum += d * d;
-  }
-  return Math.sqrt(sqSum / (count - 1));
-}
-function rollingStdDev(data, windowMs, sample = false) {
-  const out = new Array(data.length);
-  let left = 0;
-  const compute = sample ? sampleStdDev : stdDev;
-  for (let i = 0; i < data.length; i++) {
-    const t = data[i].time;
-    const wStart = t - windowMs;
-    while (left < i && data[left].time < wStart) left++;
-    const slice = data.slice(left, i + 1).map((p) => p.value);
-    const value = compute(slice);
-    out[i] = { time: t, value, synthetic: true };
-  }
-  return out;
-}
-
-// src/analyze/limits.ts
-function computeLimitExcursions(data, opts) {
-  const { high, low } = opts;
-  const result = {
-    msAboveHigh: 0,
-    msBelowLow: 0,
-    globalMax: null,
-    globalMin: null,
-    excursions: []
-  };
-  let openEpisode = null;
-  for (let i = 0; i < data.length; i++) {
-    const p = data[i];
-    if (p.value === null) {
-      if (openEpisode) {
-        result.excursions.push(openEpisode);
-        openEpisode = null;
-      }
-      continue;
-    }
-    if (result.globalMax === null || p.value > result.globalMax) result.globalMax = p.value;
-    if (result.globalMin === null || p.value < result.globalMin) result.globalMin = p.value;
-    const aboveHigh = high !== void 0 && p.value > high;
-    const belowLow = low !== void 0 && p.value < low;
-    const side = aboveHigh ? "above" : belowLow ? "below" : null;
-    const prev = data[i - 1];
-    if (prev && prev.value !== null && side !== null) {
-      const dt = p.time - prev.time;
-      if (side === "above") result.msAboveHigh += dt;
-      else result.msBelowLow += dt;
-    }
-    if (side === null) {
-      if (openEpisode) {
-        result.excursions.push(openEpisode);
-        openEpisode = null;
-      }
-      continue;
-    }
-    if (openEpisode && openEpisode.side === side) {
-      openEpisode.endTime = p.time;
-      openEpisode.durationMs = openEpisode.endTime - openEpisode.startTime;
-      if (side === "above" && p.value > openEpisode.extremum) openEpisode.extremum = p.value;
-      else if (side === "below" && p.value < openEpisode.extremum) openEpisode.extremum = p.value;
-    } else {
-      if (openEpisode) result.excursions.push(openEpisode);
-      openEpisode = {
-        startTime: p.time,
-        endTime: p.time,
-        side,
-        extremum: p.value,
-        durationMs: 0
-      };
-    }
-  }
-  if (openEpisode) result.excursions.push(openEpisode);
-  return result;
-}
+var H=class{static interpolateDataPoint(n,i,t){return {time:n.time+t*(i.time-n.time),value:(n.value??0)+t*((i.value??0)-(n.value??0))}}static interpolateAggregatedPoint(n,i,t){let u=(e,l)=>e!==null&&l!==null?e+t*(l-e):null;return {time:n.time+t*(i.time-n.time),min:u(n.min,i.min),max:u(n.max,i.max),avg:u(n.avg,i.avg),count:Math.round(n.count+t*(i.count-n.count))}}static getRuns(n,i,t=0){let u=[...n].sort((r,o)=>r.time-o.time),e=[],l=[],a=null;for(let r of u){let o=i(r),m=t>0&&a&&r.time-a.time>t;(o||m)&&l.length&&(e.push(l),l=[]),o||l.push(r),a=r;}return l.length&&e.push(l),e}static splitByBoundaries(n,i,t,u){if(n.length===0)return [];if(i.length===0)return [{data:n,zoneIndex:0}];let e=[...i].sort((o,m)=>o-m),l=[],a=o=>{let m=0;for(let s=0;s<e.length&&o>=e[s];s++)m=s+1;return m},r=[n[0]];for(let o=1;o<n.length;o++){let m=n[o-1],s=n[o],h=t(m),g=t(s),d;g>h?d=e.filter(c=>c>h&&c<=g):g<h?d=e.filter(c=>c>=g&&c<h).reverse():d=[];for(let c of d){let p=(c-h)/(g-h),x=u(m,s,p);r.push(x),l.push({data:r,zoneIndex:a((h+c)/2)}),r=[x];}r.push(s);}if(r.length>0){let o=t(r[0]),m=t(r[r.length-1]);l.push({data:r,zoneIndex:a((o+m)/2)});}return l}static splitByThreshold(n,i,t,u){let e=this.splitByBoundaries(n,[i],t,u),l={above:[],below:[]};for(let a of e)a.zoneIndex===0?l.below.push(a.data):l.above.push(a.data);return l}},L={pharma_cold:{limitLow:2,limitHigh:8,activationEnergy:83144},pharma_ambient:{limitLow:15,limitHigh:25,activationEnergy:83144},blood:{limitLow:2,limitHigh:6,activationEnergy:83144},food_chilled:{limitLow:0,limitHigh:4},frozen:{limitLow:-40,limitHigh:-18}},T=8.314462618,w=273.15;function N(n,i,t,u=1){let e=i.filter(f=>f.value!==null&&f.value!==void 0&&!isNaN(f.value)).map(f=>f.value);if(!e.length)return {time:n,min:null,max:null,avg:null,mkt:null,deltaMkt:null,stdDev:null,minutesAboveHigh:t?0:null,minutesBelowLow:t?0:null,count:0};let l=Math.min(...e),a=Math.max(...e),r=e.reduce((f,v)=>f+v,0)/e.length,o=t?0:null,m=t?0:null;if(t)for(let f=0;f<i.length;f++){let v=i[f];if(v.value!==null&&v.value!==void 0){let b=i[f+1],M=(b?b.time-v.time:0)/(1e3*60);v.value>t.limitHigh&&(o+=M),v.value<t.limitLow&&(m+=M);}}let s=0,h=0,g=(t?.activationEnergy??83144)/T;for(let f of e){s+=Math.pow(f-r,2);let v=f+w;h+=Math.exp(-g/v);}let d=s/e.length,c=Math.sqrt(d),p=null,x=null;if(e.length>=u){let f=h/e.length,v=g/-Math.log(f);p=Number((v-w).toFixed(2)),x=Number((p-r).toFixed(2));}return {time:n,min:Number(l.toFixed(2)),max:Number(a.toFixed(2)),avg:Number(r.toFixed(2)),mkt:p,deltaMkt:x,stdDev:Number(c.toFixed(2)),minutesAboveHigh:o!==null?Number(o.toFixed(1)):null,minutesBelowLow:m!==null?Number(m.toFixed(1)):null,count:e.length}}function D(n,i="hourly",t,u,e=1){if(!n.length)return [];let l=i==="hourly"?36e5:i==="daily"?864e5:t??3600*1e3,a=[...n].sort((s,h)=>s.time-h.time),r=[],o=a[0].time,m=[];for(let s of a){for(;s.time-o>=l;)r.push(N(o,m,u,e)),o+=l,m=[];m.push(s);}return r.push(N(o,m,u,e)),r}function A(n,i){let t=Math.min(...i),u=Math.max(...i),e=i.reduce((l,a)=>l+a,0)/i.length;return {time:n,min:t,max:u,avg:e,count:i.length}}function B(n,i=6e4){if(n.length<2)return [];let t=[...n].sort((e,l)=>e.time-l.time),u=[];for(let e=1;e<t.length;e++)t[e].time-t[e-1].time>i&&u.push({startTime:t[e-1].time,endTime:t[e].time});return u}function I(n){let i=[],t=0,u=0,e=0,l=0;for(let a of n)a.count>0&&(e++,a.minutesAboveHigh&&(t+=a.minutesAboveHigh),a.stdDev&&a.stdDev>u&&(u=a.stdDev),a.deltaMkt&&a.deltaMkt>2&&l++);return e===0?[]:(t>120&&i.push({type:"critical",message:`Kritische Gesamtbelastung: Das Produkt war im gesamten Zeitraum insgesamt ${Math.round(t)} Minuten zu warm.`,metric:"Total Excursion Time"}),l>e*.15&&i.push({type:"warning",message:"H\xE4ufige thermische Schocks erkannt. Die kinetische Temperatur weicht oft stark vom Schnitt ab (m\xF6gliche regelm\xE4\xDFige T\xFCr\xF6ffnungen).",metric:"Delta MKT Instability"}),u>4&&i.push({type:"warning",message:`Hohe Instabilit\xE4t gemessen. Die maximale Standardabweichung lag bei ${u}\xB0C. Das K\xFChlsystem regelt unsauber.`,metric:"Standard Deviation Peak"}),i)}function z(n,i){if(n.length<=i)return n;let t=[],u=Math.ceil(n.length/i);t.push(n[0]);for(let e=1;e<n.length;e+=u)t.push(n[Math.min(e,n.length-1)]);return t}function C(n,i="hourly",t){let u={limitLow:t.low,limitHigh:t.high,activationEnergy:t.activationEnergy};return D(n,i,t.interval,u,t.minCountForMkt??1)}var P=class{static compute(n){let i=n.map(r=>r.value).filter(r=>r!==null).sort((r,o)=>r-o);if(i.length===0)return {min:NaN,max:NaN,avg:NaN,mean:NaN,median:NaN,stdDev:NaN,count:0};let t=i.length,u=i.reduce((r,o)=>r+o,0)/t,e=t%2===1?i[Math.floor(t/2)]:(i[t/2-1]+i[t/2])/2,l=i.reduce((r,o)=>r+Math.pow(o-u,2),0)/t,a=Math.sqrt(l);return {min:i[0],max:i[t-1],avg:u,mean:u,median:e,stdDev:a,count:t}}static computeInRange(n,i,t){let u=n.filter(e=>e.time>=i&&e.time<=t);return this.compute(u)}};function _(n,i){let t=typeof n=="number"?n:n.mean,u=typeof n=="number"?i??NaN:n.stdDev;return !Number.isFinite(t)||!Number.isFinite(u)||t===0?null:Math.abs(u/t)*100}var R=83144;function F(n,i=83144){let t=i/8.314,u=0,e=0;for(let a of n)a!==null&&(u+=Math.exp(-t/(a+273.15)),e++);if(e===0)return null;let l=u/e;return t/-Math.log(l)-273.15}function q(n,i,t=83144){let u=new Array(n.length),e=0;for(let l=0;l<n.length;l++){let a=n[l].time,r=a-i;for(;e<l&&n[e].time<r;)e++;let o=n.slice(e,l+1).map(s=>s.value),m=F(o,t);u[l]={time:a,value:m,synthetic:true};}return u}function k(n){let i=0,t=0;for(let l of n)l!==null&&(i+=l,t++);if(t===0)return null;let u=i/t,e=0;for(let l of n){if(l===null)continue;let a=l-u;e+=a*a;}return Math.sqrt(e/t)}function S(n){let i=0,t=0;for(let l of n)l!==null&&(i+=l,t++);if(t<2)return null;let u=i/t,e=0;for(let l of n){if(l===null)continue;let a=l-u;e+=a*a;}return Math.sqrt(e/(t-1))}function K(n,i,t=false){let u=new Array(n.length),e=0,l=t?S:k;for(let a=0;a<n.length;a++){let r=n[a].time,o=r-i;for(;e<a&&n[e].time<o;)e++;let m=n.slice(e,a+1).map(h=>h.value),s=l(m);u[a]={time:r,value:s,synthetic:true};}return u}function G(n,i){let{high:t,low:u}=i,e={msAboveHigh:0,msBelowLow:0,globalMax:null,globalMin:null,excursions:[]},l=null;for(let a=0;a<n.length;a++){let r=n[a];if(r.value===null){l&&(e.excursions.push(l),l=null);continue}(e.globalMax===null||r.value>e.globalMax)&&(e.globalMax=r.value),(e.globalMin===null||r.value<e.globalMin)&&(e.globalMin=r.value);let o=t!==void 0&&r.value>t,m=u!==void 0&&r.value<u,s=o?"above":m?"below":null,h=n[a-1];if(h&&h.value!==null&&s!==null){let g=r.time-h.time;s==="above"?e.msAboveHigh+=g:e.msBelowLow+=g;}if(s===null){l&&(e.excursions.push(l),l=null);continue}l&&l.side===s?(l.endTime=r.time,l.durationMs=l.endTime-l.startTime,(s==="above"&&r.value>l.extremum||s==="below"&&r.value<l.extremum)&&(l.extremum=r.value)):(l&&e.excursions.push(l),l={startTime:r.time,endTime:r.time,side:s,extremum:r.value,durationMs:0});}return l&&e.excursions.push(l),e}function O(n){if(n.length===0)return [];let i=new Map;for(let u of n)for(let e of u.data){if(e.value===null)continue;let l=i.get(e.time);l||(l=[],i.set(e.time,l)),l.push({name:u.name,value:e.value});}let t=[];for(let u of [...i.keys()].sort((e,l)=>e-l)){let e=i.get(u);if(e.length===0){t.push({time:u,min:null,max:null,delta:null,minSensor:null,maxSensor:null});continue}let l=e[0],a=e[0];for(let r=1;r<e.length;r++)e[r].value<l.value&&(l=e[r]),e[r].value>a.value&&(a=e[r]);t.push({time:u,min:l.value,max:a.value,delta:e.length<2?null:a.value-l.value,minSensor:l.name,maxSensor:a.name});}return t}function U(n){let i=[];for(let e of n){let l=0,a=0,r=1/0,o=-1/0;for(let m of e.data)m.value!==null&&(l+=m.value,a++,m.value<r&&(r=m.value),m.value>o&&(o=m.value));a!==0&&i.push({name:e.name,mean:l/a,min:r,max:o,count:a});}if(i.length===0)return {sensors:[],hottest:null,coldest:null,meanDelta:null};let t=i[0],u=i[0];for(let e=1;e<i.length;e++)i[e].mean>t.mean&&(t=i[e]),i[e].mean<u.mean&&(u=i[e]);return {sensors:i,hottest:t,coldest:u,meanDelta:t.mean-u.mean}}function y(n){let i=0,t=0,u=0,e=0,l=0;for(let m of n)m.value!==null&&(i++,t+=m.time,u+=m.value,e+=m.time*m.value,l+=m.time*m.time);if(i<2)return null;let a=i*l-t*t;if(a===0)return null;let r=(i*e-t*u)/a,o=(u-r*t)/i;return {slope:r,intercept:o,n:i}}function E(n,i){if(n.length===0)return null;let t=n[n.length-1].time,u=t-i,e=n.filter(l=>l.time>=u&&l.time<=t);return y(e)}function V(n,i){let t=[],u=0;for(let e=0;e<n.length;e++){let l=n[e].time,a=l-i;for(;u<n.length&&n[u].time<a;)u++;let r=y(n.slice(u,e+1));t.push({time:l,fit:r});}return t}function W(n,i,t={}){let u=t.lookbackMs??9e5,e=t.side??"either";if(n.length===0)return null;let l=n[n.length-1];if(l.value===null)return null;let a=E(n,u);if(!a||a.slope===0)return null;let r=a.slope>0&&l.value<i,o=a.slope<0&&l.value>i;if(e==="above"&&!r||e==="below"&&!o||e==="either"&&!r&&!o)return null;let m=(i-a.intercept)/a.slope,s=m-l.time;return s<=0?null:{msUntil:s,eta:m,fit:a,currentValue:l.value}}function Y(n,i={}){let t=i.refTempCelsius??121.11,u=i.zValueKelvin??10,e=0;for(let l=0;l<n.length-1;l++){let a=n[l],r=n[l+1];if(a.value===null||r.value===null)continue;let o=(r.time-a.time)/6e4;if(o<=0)continue;let m=Math.pow(10,(a.value-t)/u),s=Math.pow(10,(r.value-t)/u);e+=(m+s)/2*o;}return e}/*!
+* MLTimeGraph — Copyright (c) 2026 Michael Lechner
+* MIT with Attribution: free use incl. commercial requires visible credit to
+* "Michael Lechner". Commercial license (no attribution) on request. See LICENSE.
+*/
 /*!
- * MLTimeGraph — Copyright (c) 2026 Michael Lechner
+ * ml-time-analyze — Copyright (c) 2026 Michael Lechner
  * MIT with Attribution: free use incl. commercial requires visible credit to
  * "Michael Lechner". Commercial license (no attribution) on request. See LICENSE.
  */
-
-export { DEFAULT_ACTIVATION_ENERGY, MovingAvg, PRODUCT_PROFILES, SeriesProcessor, StatsAggregator, aggregateBySlot, aggregateWithStats, analyzeLongTermTrends, computeLimitExcursions, createAggr, detectGaps, downsample, mkt, rollingMkt, rollingStdDev, sampleStdDev, stdDev };
+/*!
+ * MLTimeGraph — Copyright (c) 2026 Michael Lechner
+ * MIT with Attribution: free use incl. commercial requires visible credit to
+ * "Michael Lechner". Commercial license (no attribution) on request. See LICENSE.
+ */export{R as DEFAULT_ACTIVATION_ENERGY,L as PRODUCT_PROFILES,H as SeriesProcessor,P as StatsAggregator,D as aggregateBySlot,C as aggregateWithStats,I as analyzeLongTermTrends,G as computeLimitExcursions,A as createAggr,E as currentTrend,B as detectGaps,z as downsample,Y as f0Sterilization,U as hotColdSpots,y as linearFit,F as mkt,W as predictTimeToThreshold,q as rollingMkt,K as rollingStdDev,V as rollingTrend,S as sampleStdDev,O as spatialDelta,k as stdDev,_ as varianceCoefficient};