qesuite 1.0.53 → 1.0.55

package/dist/index.d.mts

@@ -175,8 +175,9 @@ declare const Capability: {
     * @customfunction
     * @param data The data to evaluate.
     * @param spec The specification to evaluate the capability against.
+    * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
     */
-    Analysis(data: number[], spec: Specification): {
+    Analysis(data: number[], spec: Specification, subgroupSize?: number): {
         Cp: number;
         Cpk: number;
         Pp: number;
@@ -187,15 +188,49 @@ declare const Capability: {
     * @customfunction
     * @param data The data to evaluate.
     * @param spec The specification to evaluate the capability against.
+    * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
     */
-    Cp(data: number[], spec: Specification): number;
+    Cp(data: number[], spec: Specification, subgroupSize?: number): number;
     /**
     * Evaluates the capability index Cpk(overall).
     * @customfunction
     * @param data The data to evaluate.
     * @param spec The specification to evaluate the capability against.
+    * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
     */
-    Cpk(data: number[], spec: Specification): number;
+    Cpk(data: number[], spec: Specification, subgroupSize?: number): number;
+    /**
+    * Evaluates the capability index Cpk(overall).
+    * @customfunction
+    * @param data The data to evaluate.
+    * @param spec The specification to evaluate the capability against.
+    * @param target The target value of the process
+    */
+    Cpm(data: number[], spec: Specification, target: number): number;
+    /**
+    * Calculates the defective parts per million.
+    * @customfunction
+    * @param data The data to evaluate.
+    * @param spec The specification to evaluate against.
+    * @param confidenceInterval [OPTIONAL] The desired confidence level of the analysis. DEFAULT = 0.95 (95%)
+    */
+    DPPM(data: number[], spec: Specification, confidenceInterval?: number): {
+        LSL: {
+            PPM: number;
+            LB: number;
+            UB: number;
+        };
+        USL: {
+            PPM: number;
+            LB: number;
+            UB: number;
+        };
+        Total: {
+            PPM: number;
+            LB: number;
+            UB: number;
+        };
+    };
     /**
     * Evaluates the capability index Pp(within).
     * @customfunction
@@ -621,7 +656,7 @@ declare const Distributions: {
         Chebyshev(data: number[], spec: Specification): number;
         JarqueBera(data: number[]): {
             JB: number;
-            p: number;
+            p: string | number;
         };
         Kurtosis(data: any): number;
         KolmogorovSmirnov(data: any): {
@@ -647,14 +682,14 @@ declare const Distributions: {
         ShapiroWilk: {
             Royston(data: number[]): {
                 W: number;
-                p: number;
+                p: string | number;
             };
             Test(data: number[], algorithmic?: boolean): {
                 W: number;
-                p: number;
+                p: string | number;
             };
             constants: number[][];
-            p(W: number, N: number): number;
+            p(W: number, N: number): number | "<0.005";
             pChart: {
                 wValues: number[][];
                 pValues: number[];

package/dist/index.d.ts

@@ -175,8 +175,9 @@ declare const Capability: {
     * @customfunction
     * @param data The data to evaluate.
     * @param spec The specification to evaluate the capability against.
+    * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
     */
-    Analysis(data: number[], spec: Specification): {
+    Analysis(data: number[], spec: Specification, subgroupSize?: number): {
         Cp: number;
         Cpk: number;
         Pp: number;
@@ -187,15 +188,49 @@ declare const Capability: {
     * @customfunction
     * @param data The data to evaluate.
     * @param spec The specification to evaluate the capability against.
+    * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
     */
-    Cp(data: number[], spec: Specification): number;
+    Cp(data: number[], spec: Specification, subgroupSize?: number): number;
     /**
     * Evaluates the capability index Cpk(overall).
     * @customfunction
     * @param data The data to evaluate.
     * @param spec The specification to evaluate the capability against.
+    * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
     */
-    Cpk(data: number[], spec: Specification): number;
+    Cpk(data: number[], spec: Specification, subgroupSize?: number): number;
+    /**
+    * Evaluates the capability index Cpk(overall).
+    * @customfunction
+    * @param data The data to evaluate.
+    * @param spec The specification to evaluate the capability against.
+    * @param target The target value of the process
+    */
+    Cpm(data: number[], spec: Specification, target: number): number;
+    /**
+    * Calculates the defective parts per million.
+    * @customfunction
+    * @param data The data to evaluate.
+    * @param spec The specification to evaluate against.
+    * @param confidenceInterval [OPTIONAL] The desired confidence level of the analysis. DEFAULT = 0.95 (95%)
+    */
+    DPPM(data: number[], spec: Specification, confidenceInterval?: number): {
+        LSL: {
+            PPM: number;
+            LB: number;
+            UB: number;
+        };
+        USL: {
+            PPM: number;
+            LB: number;
+            UB: number;
+        };
+        Total: {
+            PPM: number;
+            LB: number;
+            UB: number;
+        };
+    };
     /**
     * Evaluates the capability index Pp(within).
     * @customfunction
@@ -621,7 +656,7 @@ declare const Distributions: {
         Chebyshev(data: number[], spec: Specification): number;
         JarqueBera(data: number[]): {
             JB: number;
-            p: number;
+            p: string | number;
         };
         Kurtosis(data: any): number;
         KolmogorovSmirnov(data: any): {
@@ -647,14 +682,14 @@ declare const Distributions: {
         ShapiroWilk: {
             Royston(data: number[]): {
                 W: number;
-                p: number;
+                p: string | number;
             };
             Test(data: number[], algorithmic?: boolean): {
                 W: number;
-                p: number;
+                p: string | number;
             };
             constants: number[][];
-            p(W: number, N: number): number;
+            p(W: number, N: number): number | "<0.005";
             pChart: {
                 wValues: number[][];
                 pValues: number[];

package/dist/index.js

@@ -134,7 +134,7 @@ var Specification = class {
   full_tolerance;
   // methods
   ParseSpec(spec) {
-    let string = spec.replace(/\,/g, ".");
+    let string = spec.replace(/\,/g, ".").toUpperCase().replace(/MIN/g, ">").replace(/MAX/g, "<").replace(/\-\-/g, "\u2014");
     let nominal;
     let USL;
     let LSL;
@@ -142,6 +142,39 @@ var Specification = class {
     let tol_minus;
     let tol;
     let units;
+    if (string.indexOf("\u2014") !== -1) {
+      let bookends = string.split("\u2014").map((c) => {
+        return Number(c);
+      });
+      USL = Math.max(...bookends);
+      LSL = Math.min(...bookends);
+      tol = USL - LSL;
+      return {
+        nominal: NaN,
+        USL,
+        LSL,
+        units: "",
+        tolerance: tol
+      };
+    }
+    if (string.indexOf("<") !== -1) {
+      return {
+        nominal: NaN,
+        USL: Number(string.substring(string.indexOf("<") + 1)),
+        LSL: NaN,
+        units: "",
+        tolerance: NaN
+      };
+    }
+    if (string.indexOf(">") !== -1) {
+      return {
+        nominal: NaN,
+        USL: NaN,
+        LSL: Number(string.substring(string.indexOf(">") + 1)),
+        units: "",
+        tolerance: NaN
+      };
+    }
     if (string.indexOf("\xB1") !== -1) {
       nominal = Number(string.substring(0, string.indexOf("\xB1")).replace("\xD8", ""));
       let after = this.SplitUnits(string.substring(string.indexOf("\xB1") + 1));
@@ -154,16 +187,18 @@ var Specification = class {
       tol_plus = Number(this.SplitUnits(string.substring(string.indexOf("+") + 1, string.indexOf("/") > 0 ? string.indexOf("/") : string.indexOf("-") > 0 ? string.indexOf("-") : void 0))[0]);
       nominal = Number(this.SplitUnits(string)[0].substring(0, string.indexOf("+")));
       if (string.indexOf("-") === -1) {
-        tol_minus = 0;
+        tol_minus = NaN;
       }
     }
     if (string.indexOf("-") !== -1) {
       tol_minus = Number(this.SplitUnits(string.substring(string.indexOf("/") > 0 ? Math.max(string.indexOf("/"), string.indexOf("-")) + 1 : string.indexOf("-") + 1))[0]);
       if (string.indexOf("+") === -1) {
-        tol_plus = 0;
+        tol_plus = NaN;
         nominal = Number(this.SplitUnits(string.substring(0, string.indexOf("/") > 0 ? Math.min(string.indexOf("/"), string.indexOf("-")) : string.indexOf("-")))[0]);
       }
     }
+    USL = nominal + tol_plus;
+    LSL = nominal - tol_minus;
     if (string.indexOf("+-") !== -1 || string.indexOf("+/-") !== -1) {
       nominal = Number(string.substring(0, string.indexOf("+")));
       let after = this.SplitUnits(string.substring(string.indexOf("-") + 1));
@@ -519,13 +554,35 @@ function RoundTO(number, digits) {
   }
   let str = number.toString();
   let decimalIndex = str.indexOf(".");
+  let roundFinalDigit = Number(str.substring(decimalIndex + digits + 2, decimalIndex + digits + 3)) >= 5;
   let returner = Number(str.substring(0, decimalIndex + digits + 1));
+  if (roundFinalDigit) {
+    returner = returner + Math.pow(10, -1 * digits);
+  }
   let i = 0;
   while (returner === 0) {
     i++;
     returner = Number(str.substring(0, decimalIndex + digits + 1 + i));
   }
-  return returner;
+  if (Math.abs(number) < 1e-4) {
+    let f_i = 0;
+    let l_i = 0;
+    let d_i = 0;
+    let i2 = 0;
+    while (d_i <= digits) {
+      if (str[i2] != "0" && str[i2] != ".") {
+        d_i += 1;
+        if (f_i == 0) {
+          f_i = i2;
+        }
+        ;
+      }
+      i2 += 1;
+    }
+    l_i = i2;
+    return str[f_i] + "." + str.substring(f_i + 1, f_i + digits) + "e-" + (f_i - decimalIndex);
+  }
+  return Number(`${returner}`.substring(0, decimalIndex + digits + 1));
 }
 var Capability = {
   /**
@@ -533,11 +590,12 @@ var Capability = {
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Analysis(data, spec) {
+  Analysis(data, spec, subgroupSize = 2) {
     return {
-      Cp: Capability.Cp(data, spec),
-      Cpk: Capability.Cpk(data, spec),
+      Cp: Capability.Cp(data, spec, subgroupSize),
+      Cpk: Capability.Cpk(data, spec, subgroupSize),
       Pp: Capability.Pp(data, spec),
       Ppk: Capability.Ppk(data, spec)
     };
@@ -547,22 +605,28 @@ var Capability = {
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Cp(data, spec) {
-    let standarddeviation = StDev.W(data, 2);
-    return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+  Cp(data, spec, subgroupSize = 2) {
+    let standarddeviation = StDev.W(data, subgroupSize);
+    if (!Number.isNaN(spec.USL) && !Number.isNaN(spec.LSL)) {
+      return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+    } else {
+      return NaN;
+    }
   },
   /**
   * Evaluates the capability index Cpk(overall).
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Cpk(data, spec) {
+  Cpk(data, spec, subgroupSize = 2) {
     let mean = Mean(data);
-    let standarddeviation = StDev.W(data, 2);
-    let AbsUSL = spec.USL - mean;
-    let AbsLSL = mean - spec.LSL;
+    let standarddeviation = StDev.W(data, subgroupSize);
+    let AbsUSL = !Number.isNaN(spec.USL) ? spec.USL - mean : Infinity;
+    let AbsLSL = !Number.isNaN(spec.LSL) ? mean - spec.LSL : Infinity;
     let min;
     if (AbsUSL < AbsLSL) {
       min = AbsUSL;
@@ -572,6 +636,61 @@ var Capability = {
     return Number(min / (3 * standarddeviation));
   },
   /**
+  * Evaluates the capability index Cpk(overall).
+  * @customfunction
+  * @param data The data to evaluate.
+  * @param spec The specification to evaluate the capability against.
+  * @param target The target value of the process
+  */
+  Cpm(data, spec, target) {
+    return Math.min(spec.USL - target, target - spec.LSL) / (3 * Math.sqrt(SumSquaredDeviation(data, target) / data.length));
+  },
+  /**
+  * Calculates the defective parts per million.
+  * @customfunction
+  * @param data The data to evaluate.
+  * @param spec The specification to evaluate against.
+  * @param confidenceInterval [OPTIONAL] The desired confidence level of the analysis. DEFAULT = 0.95 (95%)
+  */
+  DPPM(data, spec, confidenceInterval = 0.95) {
+    let mean = Mean(data);
+    let std = StDev.S(data);
+    let a = 1 - confidenceInterval;
+    let Z_1_a2 = Distributions.Normal.inv(1 - a / 2);
+    let Z_lsl = (mean - spec.LSL) / std;
+    let Z_usl = (spec.USL - mean) / std;
+    let v = data.length - 1;
+    let UL = Z_lsl + Z_1_a2 * Math.sqrt(1 / data.length + Math.pow(Z_lsl, 2) / (2 * v));
+    let LL = Z_lsl - Z_1_a2 * Math.sqrt(1 / data.length + Math.pow(Z_lsl, 2) / (2 * v));
+    let UU = Z_usl + Z_1_a2 * Math.sqrt(1 / data.length + Math.pow(Z_usl, 2) / (2 * v));
+    let LU = Z_usl - Z_1_a2 * Math.sqrt(1 / data.length + Math.pow(Z_usl, 2) / (2 * v));
+    let PPM_L = (1 - Distributions.Normal.cdf((mean - spec.LSL) / std, 0, 1)) * 1e6;
+    let PPM_U = (1 - Distributions.Normal.cdf((spec.USL - mean) / std, 0, 1)) * 1e6;
+    let p_hat = 1 - Distributions.Normal.cdf((mean - spec.LSL) / std, 0, 1) + (1 - Distributions.Normal.cdf((mean - spec.USL) / std, 0, 1));
+    let ln_p = Math.log(p_hat / (1 - p_hat));
+    let part_p_u = 1 / (p_hat * (1 - p_hat)) * (1 / std) * (Distributions.Normal.pdf((spec.USL - mean) / std) - Distributions.Normal.pdf((spec.LSL - mean) / std));
+    let part_p_v = 0.5 * (1 / (p_hat * (1 - p_hat))) * ((spec.USL - mean) / Math.pow(std, 3) * Distributions.Normal.pdf((spec.USL - mean) / std) - (spec.LSL - mean) / Math.pow(std, 3) * Distributions.Normal.pdf((spec.LSL - mean) / std));
+    let V_hat = Math.pow(part_p_u, 2) * Math.pow(std, 2) / mean + Math.pow(part_p_v, 2) * 2 * Math.pow(std, 4) / v;
+    let returner = {
+      LSL: {
+        PPM: PPM_L,
+        LB: 1e6 * (1 - Distributions.Normal.cdf(UL, 0, 1)),
+        UB: 1e6 * (1 - Distributions.Normal.cdf(LL, 0, 1))
+      },
+      USL: {
+        PPM: PPM_U,
+        LB: 1e6 * (1 - Distributions.Normal.cdf(UU, 0, 1)),
+        UB: 1e6 * (1 - Distributions.Normal.cdf(LU, 0, 1))
+      },
+      Total: {
+        PPM: PPM_L + PPM_U,
+        LB: Math.exp(ln_p - Z_1_a2 * Math.sqrt(V_hat)) / (1 + Math.exp(ln_p - Z_1_a2 * Math.sqrt(V_hat))) * 1e6,
+        UB: Math.exp(ln_p + Z_1_a2 * Math.sqrt(V_hat)) / (1 + Math.exp(ln_p + Z_1_a2 * Math.sqrt(V_hat))) * 1e6
+      }
+    };
+    return returner;
+  },
+  /**
   * Evaluates the capability index Pp(within).
   * @customfunction
   * @param data The data to evaluate.
@@ -579,7 +698,11 @@ var Capability = {
   */
   Pp(data, spec) {
     let standarddeviation = StDev.S(data);
-    return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+    if (!Number.isNaN(spec.USL) && !Number.isNaN(spec.LSL)) {
+      return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+    } else {
+      return NaN;
+    }
   },
   /**
   * Evaluates the capability index Ppk(within).
@@ -590,8 +713,8 @@ var Capability = {
   Ppk(data, spec) {
     let mean = Mean(data);
     let standarddeviation = StDev.S(data);
-    let AbsUSL = spec.USL - mean;
-    let AbsLSL = mean - spec.LSL;
+    let AbsUSL = !Number.isNaN(spec.USL) ? spec.USL - mean : Infinity;
+    let AbsLSL = !Number.isNaN(spec.LSL) ? mean - spec.LSL : Infinity;
     let min;
     if (AbsUSL < AbsLSL) {
       min = AbsUSL;
@@ -784,9 +907,11 @@ function CreateCapabilityPlot(data, specification) {
   ];
   AddInfoTable(ctx, overallInfo, 3 * chartSettings.width / 4, chartSettings.margins.top, chartSettings.width / 4, chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom, "Overall");
   const mean = Mean(data);
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value], chartSettings.width / 4, chartSettings.margins.top, chartSettings.width / 2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Overall");
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value], chartSettings.width / 4, chartSettings.margins.top + 40 + (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Within");
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [specification.LSL, specification.USL], chartSettings.width / 4, chartSettings.margins.top + 80 + 2 * (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Specs");
+  let plotMin = !isNaN(specification.LSL) ? specification.LSL : Math.max(mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value, mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value);
+  let plotMax = !isNaN(specification.USL) ? specification.USL : Math.min(mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value, mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value);
+  AddLinePlot(ctx, plotMin, plotMax, mean, [mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value], chartSettings.width / 4 + 25, chartSettings.margins.top, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Overall");
+  AddLinePlot(ctx, plotMin, plotMax, mean, [mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value], chartSettings.width / 4 + 25, chartSettings.margins.top + 40 + (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Within");
+  AddLinePlot(ctx, plotMin, plotMax, mean, [specification.LSL, specification.USL], chartSettings.width / 4 + 25, chartSettings.margins.top + 80 + 2 * (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Specs");
   return canvas;
 }
 var Beta = {
@@ -1634,9 +1759,10 @@ var Distributions = {
       let S = Distributions.Normal.Skewness(data);
       let K = Distributions.Normal.Kurtosis(data);
       let JB = N / 6 * (Math.pow(S, 2) + 0.25 * Math.pow(K, 2));
+      let p = Distributions.ChiSq.RightTail(JB, 2);
       return {
         JB,
-        p: Distributions.ChiSq.RightTail(JB, 2)
+        p: p > 1 || p < 5e-3 ? "<0.005" : p
       };
     },
     Kurtosis(data) {
@@ -1669,7 +1795,7 @@ var Distributions = {
       let p = GoodnessOfFit.AndersonDarling.p(KS, N);
       return {
         KS,
-        p
+        p: p > 1 || p < 5e-3 ? "<0.005" : p
       };
     },
     Lilliefors: {
@@ -1686,6 +1812,9 @@ var Distributions = {
         if (p === 1) {
           p = Distributions.Normal.Lilliefors.pApprox(KS, N);
         }
+        if (!Number.isNaN(p)) {
+          p = Number(p) > 1 || Number(p) < 5e-3 ? "<0.005" : p;
+        }
         return p;
       },
       pApprox(KS, N) {
@@ -1827,9 +1956,10 @@ var Distributions = {
         let p_mu = 38915e-7 * Math.pow(Math.log(N), 3) - 0.083751 * Math.pow(Math.log(N), 2) - 0.31082 * Math.log(N) - 1.5861;
         let p_sigma = Math.exp(30302e-7 * Math.pow(Math.log(N), 2) - 0.082676 * Math.log(N) - 0.4803);
         let z = (Math.log(1 - W) - p_mu) / p_sigma;
+        let p = 1 - Distributions.Normal.cdf(z, 0, 1);
         return {
           W,
-          p: 1 - Distributions.Normal.cdf(z, 0, 1)
+          p: p < 5e-3 ? "<0.005" : p
         };
       },
       Test(data, algorithmic) {
@@ -1898,7 +2028,8 @@ var Distributions = {
         let w2 = wArray[inflection];
         let p1 = Distributions.Normal.ShapiroWilk.pChart.pValues[inflection - 1];
         let p2 = Distributions.Normal.ShapiroWilk.pChart.pValues[inflection];
-        return p1 + (p2 - p1) / (w2 - w1) * (W - w1);
+        let p = p1 + (p2 - p1) / (w2 - w1) * (W - w1);
+        return p < 5e-3 ? "<0.005" : p;
       },
       pChart: {
         wValues: [
@@ -2300,7 +2431,7 @@ var GoodnessOfFit = {
       if (ADStar <= 0.2) {
         p = 1 - Math.exp(-13.436 + 101.14 * ADStar - 223.73 * Math.pow(ADStar, 2));
       }
-      return p;
+      return p > 1 || p < 5e-3 ? "<0.005" : p;
     },
     /**
     * Evaluates the Anderson-Darling goodness-of-fit test for the specified data and cumulative distribution function.
@@ -2523,11 +2654,10 @@ function Median(data) {
 }
 function MovingRange(data, w = 2) {
   let movingRange = [];
-  for (let i = 1; i < data.length; i++) {
-    let h = i - w + 1;
+  for (let i = 0; i < data.length + 1 - w; i++) {
     let values = [];
     for (let j = 0; j < w; j++) {
-      values.push(Number(data[h + j]));
+      values.push(Number(data[i + j]));
     }
     movingRange.push(CalculateRange(values));
   }
@@ -3815,7 +3945,6 @@ var ANOVA = {
       chartSet.width += 1e3;
       chartSet.axis.y.tickOnly = true;
       let graph = CreateSplitGraph(individualCharts, `Main Effects Plot of ${Object.keys(response)[0]}`, chartSet);
-      console.log(graph);
       return graph;
     },
     Table(factors, response) {
@@ -3931,7 +4060,6 @@ var ANOVA = {
             });
           }
         });
-        console.log(k, SS);
       });
       return anovaTable;
     }

package/dist/index.mjs

@@ -42,7 +42,7 @@ var Specification = class {
   full_tolerance;
   // methods
   ParseSpec(spec) {
-    let string = spec.replace(/\,/g, ".");
+    let string = spec.replace(/\,/g, ".").toUpperCase().replace(/MIN/g, ">").replace(/MAX/g, "<").replace(/\-\-/g, "\u2014");
     let nominal;
     let USL;
     let LSL;
@@ -50,6 +50,39 @@ var Specification = class {
     let tol_minus;
     let tol;
     let units;
+    if (string.indexOf("\u2014") !== -1) {
+      let bookends = string.split("\u2014").map((c) => {
+        return Number(c);
+      });
+      USL = Math.max(...bookends);
+      LSL = Math.min(...bookends);
+      tol = USL - LSL;
+      return {
+        nominal: NaN,
+        USL,
+        LSL,
+        units: "",
+        tolerance: tol
+      };
+    }
+    if (string.indexOf("<") !== -1) {
+      return {
+        nominal: NaN,
+        USL: Number(string.substring(string.indexOf("<") + 1)),
+        LSL: NaN,
+        units: "",
+        tolerance: NaN
+      };
+    }
+    if (string.indexOf(">") !== -1) {
+      return {
+        nominal: NaN,
+        USL: NaN,
+        LSL: Number(string.substring(string.indexOf(">") + 1)),
+        units: "",
+        tolerance: NaN
+      };
+    }
     if (string.indexOf("\xB1") !== -1) {
       nominal = Number(string.substring(0, string.indexOf("\xB1")).replace("\xD8", ""));
       let after = this.SplitUnits(string.substring(string.indexOf("\xB1") + 1));
@@ -62,16 +95,18 @@ var Specification = class {
       tol_plus = Number(this.SplitUnits(string.substring(string.indexOf("+") + 1, string.indexOf("/") > 0 ? string.indexOf("/") : string.indexOf("-") > 0 ? string.indexOf("-") : void 0))[0]);
       nominal = Number(this.SplitUnits(string)[0].substring(0, string.indexOf("+")));
       if (string.indexOf("-") === -1) {
-        tol_minus = 0;
+        tol_minus = NaN;
       }
     }
     if (string.indexOf("-") !== -1) {
       tol_minus = Number(this.SplitUnits(string.substring(string.indexOf("/") > 0 ? Math.max(string.indexOf("/"), string.indexOf("-")) + 1 : string.indexOf("-") + 1))[0]);
       if (string.indexOf("+") === -1) {
-        tol_plus = 0;
+        tol_plus = NaN;
         nominal = Number(this.SplitUnits(string.substring(0, string.indexOf("/") > 0 ? Math.min(string.indexOf("/"), string.indexOf("-")) : string.indexOf("-")))[0]);
       }
     }
+    USL = nominal + tol_plus;
+    LSL = nominal - tol_minus;
     if (string.indexOf("+-") !== -1 || string.indexOf("+/-") !== -1) {
       nominal = Number(string.substring(0, string.indexOf("+")));
       let after = this.SplitUnits(string.substring(string.indexOf("-") + 1));
@@ -427,13 +462,35 @@ function RoundTO(number, digits) {
   }
   let str = number.toString();
   let decimalIndex = str.indexOf(".");
+  let roundFinalDigit = Number(str.substring(decimalIndex + digits + 2, decimalIndex + digits + 3)) >= 5;
   let returner = Number(str.substring(0, decimalIndex + digits + 1));
+  if (roundFinalDigit) {
+    returner = returner + Math.pow(10, -1 * digits);
+  }
   let i = 0;
   while (returner === 0) {
     i++;
     returner = Number(str.substring(0, decimalIndex + digits + 1 + i));
   }
-  return returner;
+  if (Math.abs(number) < 1e-4) {
+    let f_i = 0;
+    let l_i = 0;
+    let d_i = 0;
+    let i2 = 0;
+    while (d_i <= digits) {
+      if (str[i2] != "0" && str[i2] != ".") {
+        d_i += 1;
+        if (f_i == 0) {
+          f_i = i2;
+        }
+        ;
+      }
+      i2 += 1;
+    }
+    l_i = i2;
+    return str[f_i] + "." + str.substring(f_i + 1, f_i + digits) + "e-" + (f_i - decimalIndex);
+  }
+  return Number(`${returner}`.substring(0, decimalIndex + digits + 1));
 }
 var Capability = {
   /**
@@ -441,11 +498,12 @@ var Capability = {
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Analysis(data, spec) {
+  Analysis(data, spec, subgroupSize = 2) {
     return {
-      Cp: Capability.Cp(data, spec),
-      Cpk: Capability.Cpk(data, spec),
+      Cp: Capability.Cp(data, spec, subgroupSize),
+      Cpk: Capability.Cpk(data, spec, subgroupSize),
       Pp: Capability.Pp(data, spec),
       Ppk: Capability.Ppk(data, spec)
     };
@@ -455,22 +513,28 @@ var Capability = {
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Cp(data, spec) {
-    let standarddeviation = StDev.W(data, 2);
-    return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+  Cp(data, spec, subgroupSize = 2) {
+    let standarddeviation = StDev.W(data, subgroupSize);
+    if (!Number.isNaN(spec.USL) && !Number.isNaN(spec.LSL)) {
+      return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+    } else {
+      return NaN;
+    }
   },
   /**
   * Evaluates the capability index Cpk(overall).
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Cpk(data, spec) {
+  Cpk(data, spec, subgroupSize = 2) {
     let mean = Mean(data);
-    let standarddeviation = StDev.W(data, 2);
-    let AbsUSL = spec.USL - mean;
-    let AbsLSL = mean - spec.LSL;
+    let standarddeviation = StDev.W(data, subgroupSize);
+    let AbsUSL = !Number.isNaN(spec.USL) ? spec.USL - mean : Infinity;
+    let AbsLSL = !Number.isNaN(spec.LSL) ? mean - spec.LSL : Infinity;
     let min;
     if (AbsUSL < AbsLSL) {
       min = AbsUSL;
@@ -480,6 +544,61 @@ var Capability = {
     return Number(min / (3 * standarddeviation));
   },
   /**
+  * Evaluates the capability index Cpk(overall).
+  * @customfunction
+  * @param data The data to evaluate.
+  * @param spec The specification to evaluate the capability against.
+  * @param target The target value of the process
+  */
+  Cpm(data, spec, target) {
+    return Math.min(spec.USL - target, target - spec.LSL) / (3 * Math.sqrt(SumSquaredDeviation(data, target) / data.length));
+  },
+  /**
+  * Calculates the defective parts per million.
+  * @customfunction
+  * @param data The data to evaluate.
+  * @param spec The specification to evaluate against.
+  * @param confidenceInterval [OPTIONAL] The desired confidence level of the analysis. DEFAULT = 0.95 (95%)
+  */
+  DPPM(data, spec, confidenceInterval = 0.95) {
+    let mean = Mean(data);
+    let std = StDev.S(data);
+    let a = 1 - confidenceInterval;
+    let Z_1_a2 = Distributions.Normal.inv(1 - a / 2);
+    let Z_lsl = (mean - spec.LSL) / std;
+    let Z_usl = (spec.USL - mean) / std;
+    let v = data.length - 1;
+    let UL = Z_lsl + Z_1_a2 * Math.sqrt(1 / data.length + Math.pow(Z_lsl, 2) / (2 * v));
+    let LL = Z_lsl - Z_1_a2 * Math.sqrt(1 / data.length + Math.pow(Z_lsl, 2) / (2 * v));
+    let UU = Z_usl + Z_1_a2 * Math.sqrt(1 / data.length + Math.pow(Z_usl, 2) / (2 * v));
+    let LU = Z_usl - Z_1_a2 * Math.sqrt(1 / data.length + Math.pow(Z_usl, 2) / (2 * v));
+    let PPM_L = (1 - Distributions.Normal.cdf((mean - spec.LSL) / std, 0, 1)) * 1e6;
+    let PPM_U = (1 - Distributions.Normal.cdf((spec.USL - mean) / std, 0, 1)) * 1e6;
+    let p_hat = 1 - Distributions.Normal.cdf((mean - spec.LSL) / std, 0, 1) + (1 - Distributions.Normal.cdf((mean - spec.USL) / std, 0, 1));
+    let ln_p = Math.log(p_hat / (1 - p_hat));
+    let part_p_u = 1 / (p_hat * (1 - p_hat)) * (1 / std) * (Distributions.Normal.pdf((spec.USL - mean) / std) - Distributions.Normal.pdf((spec.LSL - mean) / std));
+    let part_p_v = 0.5 * (1 / (p_hat * (1 - p_hat))) * ((spec.USL - mean) / Math.pow(std, 3) * Distributions.Normal.pdf((spec.USL - mean) / std) - (spec.LSL - mean) / Math.pow(std, 3) * Distributions.Normal.pdf((spec.LSL - mean) / std));
+    let V_hat = Math.pow(part_p_u, 2) * Math.pow(std, 2) / mean + Math.pow(part_p_v, 2) * 2 * Math.pow(std, 4) / v;
+    let returner = {
+      LSL: {
+        PPM: PPM_L,
+        LB: 1e6 * (1 - Distributions.Normal.cdf(UL, 0, 1)),
+        UB: 1e6 * (1 - Distributions.Normal.cdf(LL, 0, 1))
+      },
+      USL: {
+        PPM: PPM_U,
+        LB: 1e6 * (1 - Distributions.Normal.cdf(UU, 0, 1)),
+        UB: 1e6 * (1 - Distributions.Normal.cdf(LU, 0, 1))
+      },
+      Total: {
+        PPM: PPM_L + PPM_U,
+        LB: Math.exp(ln_p - Z_1_a2 * Math.sqrt(V_hat)) / (1 + Math.exp(ln_p - Z_1_a2 * Math.sqrt(V_hat))) * 1e6,
+        UB: Math.exp(ln_p + Z_1_a2 * Math.sqrt(V_hat)) / (1 + Math.exp(ln_p + Z_1_a2 * Math.sqrt(V_hat))) * 1e6
+      }
+    };
+    return returner;
+  },
+  /**
   * Evaluates the capability index Pp(within).
   * @customfunction
   * @param data The data to evaluate.
@@ -487,7 +606,11 @@ var Capability = {
   */
   Pp(data, spec) {
     let standarddeviation = StDev.S(data);
-    return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+    if (!Number.isNaN(spec.USL) && !Number.isNaN(spec.LSL)) {
+      return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+    } else {
+      return NaN;
+    }
   },
   /**
   * Evaluates the capability index Ppk(within).
@@ -498,8 +621,8 @@ var Capability = {
   Ppk(data, spec) {
     let mean = Mean(data);
     let standarddeviation = StDev.S(data);
-    let AbsUSL = spec.USL - mean;
-    let AbsLSL = mean - spec.LSL;
+    let AbsUSL = !Number.isNaN(spec.USL) ? spec.USL - mean : Infinity;
+    let AbsLSL = !Number.isNaN(spec.LSL) ? mean - spec.LSL : Infinity;
     let min;
     if (AbsUSL < AbsLSL) {
       min = AbsUSL;
@@ -692,9 +815,11 @@ function CreateCapabilityPlot(data, specification) {
   ];
   AddInfoTable(ctx, overallInfo, 3 * chartSettings.width / 4, chartSettings.margins.top, chartSettings.width / 4, chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom, "Overall");
   const mean = Mean(data);
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value], chartSettings.width / 4, chartSettings.margins.top, chartSettings.width / 2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Overall");
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value], chartSettings.width / 4, chartSettings.margins.top + 40 + (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Within");
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [specification.LSL, specification.USL], chartSettings.width / 4, chartSettings.margins.top + 80 + 2 * (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Specs");
+  let plotMin = !isNaN(specification.LSL) ? specification.LSL : Math.max(mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value, mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value);
+  let plotMax = !isNaN(specification.USL) ? specification.USL : Math.min(mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value, mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value);
+  AddLinePlot(ctx, plotMin, plotMax, mean, [mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value], chartSettings.width / 4 + 25, chartSettings.margins.top, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Overall");
+  AddLinePlot(ctx, plotMin, plotMax, mean, [mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value], chartSettings.width / 4 + 25, chartSettings.margins.top + 40 + (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Within");
+  AddLinePlot(ctx, plotMin, plotMax, mean, [specification.LSL, specification.USL], chartSettings.width / 4 + 25, chartSettings.margins.top + 80 + 2 * (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Specs");
   return canvas;
 }
 var Beta = {
@@ -1542,9 +1667,10 @@ var Distributions = {
       let S = Distributions.Normal.Skewness(data);
       let K = Distributions.Normal.Kurtosis(data);
       let JB = N / 6 * (Math.pow(S, 2) + 0.25 * Math.pow(K, 2));
+      let p = Distributions.ChiSq.RightTail(JB, 2);
       return {
         JB,
-        p: Distributions.ChiSq.RightTail(JB, 2)
+        p: p > 1 || p < 5e-3 ? "<0.005" : p
       };
     },
     Kurtosis(data) {
@@ -1577,7 +1703,7 @@ var Distributions = {
       let p = GoodnessOfFit.AndersonDarling.p(KS, N);
       return {
         KS,
-        p
+        p: p > 1 || p < 5e-3 ? "<0.005" : p
       };
     },
     Lilliefors: {
@@ -1594,6 +1720,9 @@ var Distributions = {
         if (p === 1) {
           p = Distributions.Normal.Lilliefors.pApprox(KS, N);
         }
+        if (!Number.isNaN(p)) {
+          p = Number(p) > 1 || Number(p) < 5e-3 ? "<0.005" : p;
+        }
         return p;
       },
       pApprox(KS, N) {
@@ -1735,9 +1864,10 @@ var Distributions = {
         let p_mu = 38915e-7 * Math.pow(Math.log(N), 3) - 0.083751 * Math.pow(Math.log(N), 2) - 0.31082 * Math.log(N) - 1.5861;
         let p_sigma = Math.exp(30302e-7 * Math.pow(Math.log(N), 2) - 0.082676 * Math.log(N) - 0.4803);
         let z = (Math.log(1 - W) - p_mu) / p_sigma;
+        let p = 1 - Distributions.Normal.cdf(z, 0, 1);
         return {
           W,
-          p: 1 - Distributions.Normal.cdf(z, 0, 1)
+          p: p < 5e-3 ? "<0.005" : p
         };
       },
       Test(data, algorithmic) {
@@ -1806,7 +1936,8 @@ var Distributions = {
         let w2 = wArray[inflection];
         let p1 = Distributions.Normal.ShapiroWilk.pChart.pValues[inflection - 1];
         let p2 = Distributions.Normal.ShapiroWilk.pChart.pValues[inflection];
-        return p1 + (p2 - p1) / (w2 - w1) * (W - w1);
+        let p = p1 + (p2 - p1) / (w2 - w1) * (W - w1);
+        return p < 5e-3 ? "<0.005" : p;
       },
       pChart: {
         wValues: [
@@ -2208,7 +2339,7 @@ var GoodnessOfFit = {
       if (ADStar <= 0.2) {
         p = 1 - Math.exp(-13.436 + 101.14 * ADStar - 223.73 * Math.pow(ADStar, 2));
       }
-      return p;
+      return p > 1 || p < 5e-3 ? "<0.005" : p;
     },
     /**
     * Evaluates the Anderson-Darling goodness-of-fit test for the specified data and cumulative distribution function.
@@ -2431,11 +2562,10 @@ function Median(data) {
 }
 function MovingRange(data, w = 2) {
   let movingRange = [];
-  for (let i = 1; i < data.length; i++) {
-    let h = i - w + 1;
+  for (let i = 0; i < data.length + 1 - w; i++) {
     let values = [];
     for (let j = 0; j < w; j++) {
-      values.push(Number(data[h + j]));
+      values.push(Number(data[i + j]));
     }
     movingRange.push(CalculateRange(values));
   }
@@ -3723,7 +3853,6 @@ var ANOVA = {
       chartSet.width += 1e3;
       chartSet.axis.y.tickOnly = true;
       let graph = CreateSplitGraph(individualCharts, `Main Effects Plot of ${Object.keys(response)[0]}`, chartSet);
-      console.log(graph);
       return graph;
     },
     Table(factors, response) {
@@ -3839,7 +3968,6 @@ var ANOVA = {
             });
           }
         });
-        console.log(k, SS);
       });
       return anovaTable;
     }

package/index.ts

@@ -50,7 +50,7 @@ export class Specification{
   ParseSpec(spec: string) {
       // This function parses specs and returns an array of values, in this order:
       // Nominal, USL, LSL, Units, Tolerance(if UT and LT differ then {UT, LT})
-      let string = spec.replace(/\,/g, ".");
+      let string = spec.replace(/\,/g, ".").toUpperCase().replace(/MIN/g, ">").replace(/MAX/g, "<").replace(/\-\-/g,"—");
       let nominal;
       let USL;
       let LSL;
@@ -58,7 +58,39 @@ export class Specification{
       let tol_minus;
       let tol;
       let units;
-    
+
+      if (string.indexOf("—") !== -1){
+        let bookends = string.split("—").map(c => {return Number(c)});
+        USL = Math.max(...bookends);
+        LSL = Math.min(...bookends);
+        tol = USL - LSL;
+        
+        return {
+          nominal: NaN,
+          USL: USL,
+          LSL: LSL,
+          units: '',
+          tolerance: tol
+        }
+      }
+      if (string.indexOf("<") !== -1){    
+        return {
+          nominal: NaN,
+          USL: Number(string.substring(string.indexOf("<")+1)),
+          LSL: NaN,
+          units: '',
+          tolerance: NaN
+        }
+      }
+      if (string.indexOf(">") !== -1){        
+        return {
+          nominal: NaN,
+          USL: NaN,
+          LSL: Number(string.substring(string.indexOf(">")+1)),
+          units: '',
+          tolerance: NaN
+        }
+      }
       if (string.indexOf("±") !== -1) {
         nominal = Number(string.substring(0, string.indexOf("±")).replace("Ø", ""));
         let after = this.SplitUnits(string.substring(string.indexOf("±") + 1));
@@ -71,16 +103,18 @@ export class Specification{
         tol_plus = Number(this.SplitUnits(string.substring(string.indexOf("+") + 1, string.indexOf("/") > 0 ? string.indexOf("/") : string.indexOf("-") > 0 ? string.indexOf("-") : undefined))[0]);
         nominal = Number(this.SplitUnits(string)[0].substring(0, string.indexOf("+")));
         if(string.indexOf("-") === -1){
-          tol_minus = 0;
+          tol_minus = NaN;
         }
       }
       if (string.indexOf("-") !== -1) {
         tol_minus = Number(this.SplitUnits(string.substring(string.indexOf("/") > 0 ? Math.max(string.indexOf("/"), string.indexOf("-")) + 1 : string.indexOf("-") + 1))[0]);
         if(string.indexOf("+") === -1){
-          tol_plus = 0;
+          tol_plus = NaN;
           nominal = Number(this.SplitUnits(string.substring(0, string.indexOf("/") > 0 ? Math.min(string.indexOf("/"), string.indexOf("-")) : string.indexOf("-")))[0])
         }
       }
+      USL = nominal + tol_plus
+      LSL = nominal - tol_minus
       if (string.indexOf("+-") !== -1 || string.indexOf("+/-") !== -1) {
         nominal = Number(string.substring(0, string.indexOf("+")));
         let after = this.SplitUnits(string.substring(string.indexOf("-") + 1));
@@ -488,14 +522,35 @@ export function RoundTO(number: number, digits: number){
   }
   let str = number.toString();
   let decimalIndex = str.indexOf('.');
+  let roundFinalDigit = Number(str.substring(decimalIndex + digits + 2,decimalIndex + digits + 3)) >= 5
+  
   let returner = Number(str.substring(0, decimalIndex + digits + 1));
-
+  if(roundFinalDigit){
+    returner = returner + Math.pow(10,-1*digits)
+  }
+  
   let i = 0;
   while(returner === 0){
       i++;
       returner =  Number(str.substring(0, decimalIndex + digits + 1 + i));
   }
-  return returner
+  if(Math.abs(number) < .0001){
+    let f_i = 0;
+    let l_i = 0;
+    let d_i = 0;
+    let i = 0;
+    while(d_i <= digits){
+      if(str[i] != "0" && str[i] != "."){
+        d_i += 1;
+        if(f_i == 0){f_i = i};
+      }
+      i += 1
+    }
+    l_i = i;
+    
+    return str[f_i] + "." + str.substring(f_i + 1, f_i + digits) + "e-" + (f_i - decimalIndex); 
+  }
+  return Number(`${returner}`.substring(0, decimalIndex + digits + 1))
 }
 
 // End Classes
@@ -509,11 +564,12 @@ export const Capability = {
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Analysis(data: number[], spec: Specification){
+  Analysis(data: number[], spec: Specification, subgroupSize: number = 2){
       return {
-        Cp: Capability.Cp(data, spec),
-        Cpk: Capability.Cpk(data, spec),
+        Cp: Capability.Cp(data, spec, subgroupSize),
+        Cpk: Capability.Cpk(data, spec, subgroupSize),
         Pp: Capability.Pp(data, spec),
         Ppk: Capability.Ppk(data, spec),
       }
@@ -523,23 +579,28 @@ export const Capability = {
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Cp(data: number[], spec: Specification){      
-    let standarddeviation = StDev.W(data, 2);
-    
-    return Number((spec.USL - spec.LSL) / (6 * standarddeviation));      
+  Cp(data: number[], spec: Specification, subgroupSize: number = 2){      
+    let standarddeviation = StDev.W(data, subgroupSize);
+    if(!Number.isNaN(spec.USL) && !Number.isNaN(spec.LSL)){
+      return Number((spec.USL - spec.LSL) / (6 * standarddeviation));
+    }else{
+      return NaN
+    }
   },
   /**
   * Evaluates the capability index Cpk(overall).
   * @customfunction
   * @param data The data to evaluate.
   * @param spec The specification to evaluate the capability against.
+  * @param subgroupSize The size of the subgroup to use for calculating within standard deviation
   */
-  Cpk(data: number[], spec: Specification){      
+  Cpk(data: number[], spec: Specification, subgroupSize: number = 2){      
     let mean = Mean(data);
-    let standarddeviation = StDev.W(data,2);
-    let AbsUSL = spec.USL - mean;
-    let AbsLSL = mean - spec.LSL;
+    let standarddeviation = StDev.W(data, subgroupSize);
+    let AbsUSL = !Number.isNaN(spec.USL) ? spec.USL - mean : Infinity;
+    let AbsLSL = !Number.isNaN(spec.LSL) ? mean - spec.LSL : Infinity;
     let min;
     if (AbsUSL < AbsLSL) {
       min = AbsUSL;
@@ -551,6 +612,64 @@ export const Capability = {
   
   },
   /**
+  * Evaluates the capability index Cpk(overall).
+  * @customfunction
+  * @param data The data to evaluate.
+  * @param spec The specification to evaluate the capability against.
+  * @param target The target value of the process
+  */
+  Cpm(data: number[], spec: Specification, target: number){
+    return Math.min(spec.USL - target, target - spec.LSL)/(3*Math.sqrt(SumSquaredDeviation(data, target)/data.length))  
+  },
+  /**
+  * Calculates the defective parts per million.
+  * @customfunction
+  * @param data The data to evaluate.
+  * @param spec The specification to evaluate against.
+  * @param confidenceInterval [OPTIONAL] The desired confidence level of the analysis. DEFAULT = 0.95 (95%)
+  */
+  DPPM(data: number[], spec: Specification, confidenceInterval: number = 0.95){
+    let mean = Mean(data)
+    let std = StDev.S(data)
+    let a = 1 - confidenceInterval
+    let Z_1_a2 = Distributions.Normal.inv(1 - a/2);
+    let Z_lsl = (mean - spec.LSL)/std
+    let Z_usl = (spec.USL - mean)/std
+    let v = data.length - 1
+    let UL = Z_lsl + Z_1_a2 * Math.sqrt(1/data.length + Math.pow(Z_lsl, 2)/(2*v))
+    let LL = Z_lsl - Z_1_a2 * Math.sqrt(1/data.length + Math.pow(Z_lsl, 2)/(2*v))
+    let UU = Z_usl + Z_1_a2 * Math.sqrt(1/data.length + Math.pow(Z_usl, 2)/(2*v))
+    let LU = Z_usl - Z_1_a2 * Math.sqrt(1/data.length + Math.pow(Z_usl, 2)/(2*v))
+    let PPM_L = (1 - Distributions.Normal.cdf((mean - spec.LSL)/std, 0, 1)) * 1000000
+    let PPM_U = (1 - Distributions.Normal.cdf((spec.USL - mean)/std, 0, 1)) * 1000000
+
+    let p_hat = (1 - Distributions.Normal.cdf((mean - spec.LSL)/std, 0, 1)) + (1 - Distributions.Normal.cdf((mean - spec.USL)/std, 0, 1))
+    let ln_p = Math.log(p_hat/(1 - p_hat));
+    let part_p_u = (1/(p_hat * (1 - p_hat))) * (1/std) * (Distributions.Normal.pdf((spec.USL - mean)/std) - Distributions.Normal.pdf((spec.LSL - mean)/std))
+    let part_p_v = .5 * (1/(p_hat * (1 - p_hat))) * ( (spec.USL - mean)/Math.pow(std, 3) * Distributions.Normal.pdf((spec.USL - mean)/std) - (spec.LSL - mean)/Math.pow(std, 3) * Distributions.Normal.pdf((spec.LSL - mean)/std))
+    let V_hat = Math.pow(part_p_u,2) * Math.pow(std, 2)/mean + Math.pow(part_p_v, 2) * 2 * Math.pow(std, 4)/v
+
+    let returner = {
+      LSL: {
+        PPM: PPM_L,
+        LB: 1000000 * (1 - Distributions.Normal.cdf(UL, 0, 1)),
+        UB: 1000000 * (1 - Distributions.Normal.cdf(LL, 0, 1)),
+      },
+      USL: {
+        PPM: PPM_U,
+        LB: 1000000 * (1 - Distributions.Normal.cdf(UU, 0, 1)),
+        UB: 1000000 * (1 - Distributions.Normal.cdf(LU, 0, 1)),
+      },
+      Total: {
+        PPM: PPM_L + PPM_U,
+        LB: Math.exp(ln_p - Z_1_a2*Math.sqrt(V_hat))/(1 + Math.exp(ln_p - Z_1_a2*Math.sqrt(V_hat))) * 1000000,
+        UB: Math.exp(ln_p + Z_1_a2*Math.sqrt(V_hat))/(1 + Math.exp(ln_p + Z_1_a2*Math.sqrt(V_hat))) * 1000000
+      }
+    }
+
+    return returner
+  },
+  /**
   * Evaluates the capability index Pp(within).
   * @customfunction
   * @param data The data to evaluate.
@@ -559,7 +678,12 @@ export const Capability = {
   Pp(data: number[], spec: Specification) {
     let standarddeviation = StDev.S(data);
   
-    return Number((spec.USL-spec.LSL) / (6 * standarddeviation));
+    if(!Number.isNaN(spec.USL) && !Number.isNaN(spec.LSL)){
+      return Number((spec.USL-spec.LSL) / (6 * standarddeviation));
+    }else{
+      return NaN
+    }
+    
   
   },
   /**
@@ -571,8 +695,8 @@ export const Capability = {
   Ppk(data: number[], spec: Specification) {
     let mean = Mean(data);
     let standarddeviation = StDev.S(data);
-    let AbsUSL = spec.USL - mean;
-    let AbsLSL = mean - spec.LSL;
+    let AbsUSL = !Number.isNaN(spec.USL) ? spec.USL - mean : Infinity;
+    let AbsLSL = !Number.isNaN(spec.LSL) ? mean - spec.LSL : Infinity;
     let min;
     if (AbsUSL < AbsLSL) {
       min = AbsUSL;
@@ -809,13 +933,17 @@ export function CreateCapabilityPlot(data: number[], specification: Specificatio
   AddInfoTable(ctx, overallInfo, 3*chartSettings.width/4, chartSettings.margins.top, chartSettings.width/4, chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom, "Overall");
   
   const mean = Mean(data);
+  let plotMin = !isNaN(specification.LSL) ? specification.LSL : Math.max(mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value, mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value);
+  let plotMax = !isNaN(specification.USL) ? specification.USL : Math.min(mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value, mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value);
+
   // Add Overall Line Plot
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value], chartSettings.width/4, chartSettings.margins.top, chartSettings.width/2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom)/3 - 10, "Overall")
-  
+  AddLinePlot(ctx, plotMin, plotMax, mean, [mean - 3 * overallInfo[0].value, mean + 3 * overallInfo[0].value], chartSettings.width / 4 + 25, chartSettings.margins.top, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Overall");
+
   // Add Within Line Plot
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value], chartSettings.width/4, chartSettings.margins.top + 40 + (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom)/3, chartSettings.width/2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom)/3 - 10, "Within")
+  AddLinePlot(ctx, plotMin, plotMax, mean, [mean - 3 * withinInfo[0].value, mean + 3 * withinInfo[0].value], chartSettings.width / 4 + 25, chartSettings.margins.top + 40 + (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Within");
+
   // Add Specs Line Plot
-  AddLinePlot(ctx, specification.LSL, specification.USL, mean, [specification.LSL, specification.USL], chartSettings.width/4, chartSettings.margins.top + 80 + 2*(chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom)/3, chartSettings.width/2, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom)/3 - 10, "Specs")
+  AddLinePlot(ctx, plotMin, plotMax, mean, [specification.LSL, specification.USL], chartSettings.width / 4 + 25, chartSettings.margins.top + 80 + 2 * (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3, chartSettings.width / 2 - 50, (chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom) / 3 - 10, "Specs");
   return canvas;
 }
 
@@ -1713,13 +1841,13 @@ export const Distributions = {
       let K = Distributions.Normal.Kurtosis(data);
       
       let JB = (N/6)*(Math.pow(S, 2) + .25 * Math.pow(K, 2));
+      let p = Distributions.ChiSq.RightTail(JB, 2)
       return {
         JB: JB,
-        p: Distributions.ChiSq.RightTail(JB, 2)
+        p: p > 1 || p < 0.005 ? "<0.005" : p
       }
     },
     Kurtosis(data){
-      // SRC: https://www.omnicalculator.com/statistics/skewness#what-are-skewness-and-kurtosis
       let N = data.length;
       let mean = Mean(data);
       let numsum = 0;
@@ -1748,7 +1876,7 @@ export const Distributions = {
       let p = GoodnessOfFit.AndersonDarling.p(KS, N);      
       return {
         KS: KS,
-        p: p
+        p: p > 1 || p < 0.005 ? "<0.005" : p
       };
     },
     Lilliefors:{
@@ -1765,6 +1893,9 @@ export const Distributions = {
         if(p === 1){
           p = Distributions.Normal.Lilliefors.pApprox(KS, N)
         }
+        if(!Number.isNaN(p)){
+          p = Number(p) > 1 || Number(p) < 0.005 ? "<0.005" : p
+        }
         return p;
       },
       pApprox(KS: number, N: number){
@@ -1910,10 +2041,10 @@ export const Distributions = {
         let p_mu = 0.0038915*Math.pow(Math.log(N),3)-0.083751*Math.pow(Math.log(N),2)-0.31082*Math.log(N)-1.5861;
         let p_sigma = Math.exp(0.0030302*Math.pow(Math.log(N),2)-0.082676*Math.log(N)-0.4803);
         let z = (Math.log(1 - W) - p_mu)/p_sigma;
-        
+        let p = 1 - Distributions.Normal.cdf(z, 0, 1)
         return {
           W: W,
-          p: 1 - Distributions.Normal.cdf(z, 0, 1)
+          p: p < .005 ? "<0.005" : p
         }
       },
       Test(data: number[], algorithmic?: boolean){
@@ -1986,7 +2117,8 @@ export const Distributions = {
         let w2 = wArray[inflection];
         let p1 = Distributions.Normal.ShapiroWilk.pChart.pValues[inflection - 1]
         let p2 = Distributions.Normal.ShapiroWilk.pChart.pValues[inflection]
-        return p1 + ((p2-p1)/(w2-w1))*(W-w1);
+        let p = p1 + ((p2-p1)/(w2-w1))*(W-w1); 
+        return p < .005 ? "<0.005" : p
       },
       pChart:{
         wValues: [[0.753,0.756,0.767,0.789,0.959,0.998,0.999,1,1,],
@@ -2377,7 +2509,7 @@ export const GoodnessOfFit = {
       if (ADStar <= .2) {
         p = 1 - Math.exp(-13.436 + 101.14 * ADStar - 223.73 * Math.pow(ADStar, 2))
       }
-      return p;
+      return p > 1 || p < 0.005 ? "<0.005" : p;
     },   
     /**
     * Evaluates the Anderson-Darling goodness-of-fit test for the specified data and cumulative distribution function.
@@ -2749,12 +2881,10 @@ export function Median(data: any[]){
 export function MovingRange(data: any[], w: number = 2){
   let movingRange: number[] = [];
   
-  for(let i = 1; i < data.length; i++){
-    let h = i - w + 1;
-
+  for(let i = 0; i < data.length + 1 - w; i++){  
     let values: number[] = [];
     for(let j = 0; j < w; j++){
-      values.push(Number(data[h + j]));
+      values.push(Number(data[i + j]));
     }
 
     movingRange.push(CalculateRange(values))
@@ -4126,7 +4256,6 @@ export const ANOVA = {
       chartSet.width += 1000
       chartSet.axis.y.tickOnly = true
       let graph = CreateSplitGraph(individualCharts, `Main Effects Plot of ${Object.keys(response)[0]}`, chartSet)
-      console.log(graph)
       return graph
     },
     Table(factors: any, response: any){
@@ -4257,7 +4386,6 @@ export const ANOVA = {
             })
           }
         })
-        console.log(k, SS);
       })
 
 

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "qesuite",
-  "version": "1.0.53",
+  "version": "1.0.55",
   "description": "Performs advanced statistical analysis of data. Specifically designed for engineering statistical analysis",
   "main": "dist/index.js",
   "module": "dist/index.mjs",

package/versions/1_0_54.md

@@ -0,0 +1,19 @@
+ParseSpec() - added support for single-sided and range specifications
+    New Specification Types
+        - nom-tol_minus : singled-sided spec.; all values are to be greater than nom - tol_minus
+        - MIN# : singled-sided spec.; all values are to be greater than the number(#); not case sensitive
+        - <# : singled-sided spec.; all values are to be less than the number(#)
+        - nom+tol_plus : singled-sided spec.; all values are to be less than nom + tol_plus
+        - MAX# : singled-sided spec.; all values are to be less than the number(#); not case sensitive
+        - ># : singled-sided spec.; all values are to be greater than the number(#)
+        - a--b: range spec.; all values are to be between numbers a & b 
+        - a—b: range spec.; all values are to be between numbers a & b 
+    Affected QE.FN()'s
+        - QE.CAPABILITYANALYSIS()
+        - QE.CP()
+        - QE.CPK()
+        - QE.PP()
+        - QE.PPK()
+        - QE.CHEBYSHEV()
+
+p-values: fixed bug that allowed p-values > 1

package/versions/1_0_55.md

@@ -0,0 +1,2 @@
+Bug Fix: Capability plots did not have missing bounds when using single sided tolerances
+Added rounding on final decimal of RoundTO() and added support for Scientific notation