qesuite 1.0.36 → 1.0.38

package/dist/index.d.mts

@@ -618,7 +618,11 @@ declare const Distributions: {
             Critical(data: number[]): number;
         };
         ShapiroWilk: {
-            Test(data: number[], expanded?: boolean): {
+            Royston(data: number[]): {
+                W: number;
+                p: number;
+            };
+            Test(data: number[], algorithmic?: boolean): {
                 W: number;
                 p: number;
             };
@@ -963,6 +967,7 @@ declare function CreateSummaryChart(charts: HTMLCanvasElement[], title: string):
 declare function CreateStackedChart(charts: Chart[], title?: string, chartSettings?: ChartSettings): HTMLCanvasElement;
 declare function GetFunctionValues(callback: Function, range: number[], count: number): DataSet;
 declare function QQPlotChart(data: number[], distribution?: string): HTMLCanvasElement;
+declare function CreateControlChart(data: number[], chartTitle?: string, UCL?: number, LCL?: number): HTMLCanvasElement;
 /** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 */
 /**
@@ -1145,4 +1150,4 @@ declare function GRRByPartChart(data: GRRData): HTMLCanvasElement;
 declare function GRRPartxOperatorChart(data: GRRData): HTMLCanvasElement;
 declare function CompileGRRObject(operators: any, parts: any, measurements: any): GRRData;
 
-export { AnovaTableTwoWay, AverageMovingRange, AxisSettings, Beta, CanvasDrawSettings, Capability, Chart, ChartSettings, CompileGRRObject, CreateBoxandWhiskerGraph, CreateCapabilityHistogram, CreateCapabilityPlot, CreateCategoricalBarGraph, CreateContinuousBarGraph, CreateScatterPlot, CreateSplitGraph, CreateStackedChart, CreateSummaryChart, CustomGridline, DataCollection, DataSet, Distributions, ERF, EngFont, FixedPointIteration, G1, G1Graphs, GRR, GRRByPartChart, GRRData, GRRGraphs, GRROperator, GRRPartxOperatorChart, GRRReplication, GageEvaluation, Gamma, GetFunctionValues, GetHistogramDataset, GoodnessOfFit, IndividualDistributionIdentification, IndividualDistributionPlots, IndividualValuePlot, LastObservationsChart, Line, Margin, MaximumLikelihoodParameters, Mean, Median, MovingRange, MovingRangeChart, NewtonRaphson, ObjectToArray, ParameterizedCDF, Point, QQPlot, QQPlotChart, Rang, RoundTO, SerializeData, Specification, StDev, Sum, UnbiasingConstant, VarianceComponents };
+export { AnovaTableTwoWay, AverageMovingRange, AxisSettings, Beta, CanvasDrawSettings, Capability, Chart, ChartSettings, CompileGRRObject, CreateBoxandWhiskerGraph, CreateCapabilityHistogram, CreateCapabilityPlot, CreateCategoricalBarGraph, CreateContinuousBarGraph, CreateControlChart, CreateScatterPlot, CreateSplitGraph, CreateStackedChart, CreateSummaryChart, CustomGridline, DataCollection, DataSet, Distributions, ERF, EngFont, FixedPointIteration, G1, G1Graphs, GRR, GRRByPartChart, GRRData, GRRGraphs, GRROperator, GRRPartxOperatorChart, GRRReplication, GageEvaluation, Gamma, GetFunctionValues, GetHistogramDataset, GoodnessOfFit, IndividualDistributionIdentification, IndividualDistributionPlots, IndividualValuePlot, LastObservationsChart, Line, Margin, MaximumLikelihoodParameters, Mean, Median, MovingRange, MovingRangeChart, NewtonRaphson, ObjectToArray, ParameterizedCDF, Point, QQPlot, QQPlotChart, Rang, RoundTO, SerializeData, Specification, StDev, Sum, UnbiasingConstant, VarianceComponents };

package/dist/index.d.ts

@@ -618,7 +618,11 @@ declare const Distributions: {
             Critical(data: number[]): number;
         };
         ShapiroWilk: {
-            Test(data: number[], expanded?: boolean): {
+            Royston(data: number[]): {
+                W: number;
+                p: number;
+            };
+            Test(data: number[], algorithmic?: boolean): {
                 W: number;
                 p: number;
             };
@@ -963,6 +967,7 @@ declare function CreateSummaryChart(charts: HTMLCanvasElement[], title: string):
 declare function CreateStackedChart(charts: Chart[], title?: string, chartSettings?: ChartSettings): HTMLCanvasElement;
 declare function GetFunctionValues(callback: Function, range: number[], count: number): DataSet;
 declare function QQPlotChart(data: number[], distribution?: string): HTMLCanvasElement;
+declare function CreateControlChart(data: number[], chartTitle?: string, UCL?: number, LCL?: number): HTMLCanvasElement;
 /** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 */
 /**
@@ -1145,4 +1150,4 @@ declare function GRRByPartChart(data: GRRData): HTMLCanvasElement;
 declare function GRRPartxOperatorChart(data: GRRData): HTMLCanvasElement;
 declare function CompileGRRObject(operators: any, parts: any, measurements: any): GRRData;
 
-export { AnovaTableTwoWay, AverageMovingRange, AxisSettings, Beta, CanvasDrawSettings, Capability, Chart, ChartSettings, CompileGRRObject, CreateBoxandWhiskerGraph, CreateCapabilityHistogram, CreateCapabilityPlot, CreateCategoricalBarGraph, CreateContinuousBarGraph, CreateScatterPlot, CreateSplitGraph, CreateStackedChart, CreateSummaryChart, CustomGridline, DataCollection, DataSet, Distributions, ERF, EngFont, FixedPointIteration, G1, G1Graphs, GRR, GRRByPartChart, GRRData, GRRGraphs, GRROperator, GRRPartxOperatorChart, GRRReplication, GageEvaluation, Gamma, GetFunctionValues, GetHistogramDataset, GoodnessOfFit, IndividualDistributionIdentification, IndividualDistributionPlots, IndividualValuePlot, LastObservationsChart, Line, Margin, MaximumLikelihoodParameters, Mean, Median, MovingRange, MovingRangeChart, NewtonRaphson, ObjectToArray, ParameterizedCDF, Point, QQPlot, QQPlotChart, Rang, RoundTO, SerializeData, Specification, StDev, Sum, UnbiasingConstant, VarianceComponents };
+export { AnovaTableTwoWay, AverageMovingRange, AxisSettings, Beta, CanvasDrawSettings, Capability, Chart, ChartSettings, CompileGRRObject, CreateBoxandWhiskerGraph, CreateCapabilityHistogram, CreateCapabilityPlot, CreateCategoricalBarGraph, CreateContinuousBarGraph, CreateControlChart, CreateScatterPlot, CreateSplitGraph, CreateStackedChart, CreateSummaryChart, CustomGridline, DataCollection, DataSet, Distributions, ERF, EngFont, FixedPointIteration, G1, G1Graphs, GRR, GRRByPartChart, GRRData, GRRGraphs, GRROperator, GRRPartxOperatorChart, GRRReplication, GageEvaluation, Gamma, GetFunctionValues, GetHistogramDataset, GoodnessOfFit, IndividualDistributionIdentification, IndividualDistributionPlots, IndividualValuePlot, LastObservationsChart, Line, Margin, MaximumLikelihoodParameters, Mean, Median, MovingRange, MovingRangeChart, NewtonRaphson, ObjectToArray, ParameterizedCDF, Point, QQPlot, QQPlotChart, Rang, RoundTO, SerializeData, Specification, StDev, Sum, UnbiasingConstant, VarianceComponents };

package/dist/index.js

@@ -33,6 +33,7 @@ __export(QESuite_exports, {
   CreateCapabilityPlot: () => CreateCapabilityPlot,
   CreateCategoricalBarGraph: () => CreateCategoricalBarGraph,
   CreateContinuousBarGraph: () => CreateContinuousBarGraph,
+  CreateControlChart: () => CreateControlChart,
   CreateScatterPlot: () => CreateScatterPlot,
   CreateSplitGraph: () => CreateSplitGraph,
   CreateStackedChart: () => CreateStackedChart,
@@ -1670,7 +1671,37 @@ var Distributions = {
       }
     },
     ShapiroWilk: {
-      Test(data, expanded) {
+      Royston(data) {
+        let N = data.length;
+        let mean = Mean(data);
+        let sorted = [...data].sort();
+        let u = 1 / Math.sqrt(N);
+        let m_i = [];
+        let M = 0;
+        for (let i = 0; i < N; i++) {
+          let m = Distributions.Normal.inv((i + 1 - 0.375) / (N + 0.25));
+          m_i[i] = m;
+          M += Math.pow(m, 2);
+        }
+        let a_i = [];
+        a_i[N - 1] = -2.706056 * Math.pow(u, 5) + 4.434685 * Math.pow(u, 4) - 2.07119 * Math.pow(u, 3) - 0.147981 * Math.pow(u, 2) + 0.221157 * u + m_i[N - 1] / Math.sqrt(M);
+        a_i[N - 2] = -3.582633 * Math.pow(u, 5) + 5.682633 * Math.pow(u, 4) - 1.752461 * Math.pow(u, 3) - 0.293762 * Math.pow(u, 2) + 0.042981 * u + m_i[N - 2] / Math.sqrt(M);
+        a_i[0] = -1 * a_i[N - 1];
+        a_i[1] = -1 * a_i[N - 2];
+        let eps = (M - 2 * Math.pow(m_i[N - 1], 2) - 2 * Math.pow(m_i[N - 2], 2)) / (1 - 2 * Math.pow(a_i[N - 1], 2) - 2 * Math.pow(a_i[N - 2], 2));
+        for (let i = 2; i < N - 2; i++) {
+          a_i[i] = m_i[i] / Math.sqrt(eps);
+        }
+        let W = Math.pow(Correlate(a_i, sorted), 2);
+        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;
+        return {
+          W,
+          p: 1 - Distributions.Normal.cdf(z, 0, 1)
+        };
+      },
+      Test(data, algorithmic) {
         let N = data.length;
         let mean = Mean(data);
         let sorted = [...data].sort();
@@ -1679,39 +1710,21 @@ var Distributions = {
         sorted.forEach((s) => {
           SS += Math.pow(s - mean, 2);
         });
-        if (expanded || N > 50) {
-          let u = 1 / Math.sqrt(N);
-          let m_i = [];
-          let M = 0;
-          for (let i = 0; i < N; i++) {
-            let m = Distributions.Normal.inv((i + 1 - 0.375) / (N + 0.25));
-            m_i[i] = m;
-            M += Math.pow(m, 2);
-          }
-          let a_i = [];
-          a_i[N - 1] = -2.706056 * Math.pow(u, 5) + 4.434685 * Math.pow(u, 4) - 2.07119 * Math.pow(u, 3) - 0.147981 * Math.pow(u, 2) + 0.221157 * u + m_i[N - 1] / Math.sqrt(M);
-          a_i[N - 2] = -3.582633 * Math.pow(u, 5) + 5.682633 * Math.pow(u, 4) - 1.752461 * Math.pow(u, 3) - 0.293762 * Math.pow(u, 2) + 0.042981 * u + m_i[N - 2] / Math.sqrt(M);
-          a_i[0] = -1 * a_i[N - 1];
-          a_i[1] = -1 * a_i[N - 2];
-          let eps = (M - 2 * Math.pow(m_i[N - 1], 2) - 2 * Math.pow(m_i[N - 2], 2)) / (1 - 2 * Math.pow(a_i[N - 1], 2) - 2 * Math.pow(a_i[N - 2], 2));
-          for (let i = 2; i < N - 2; i++) {
-            a_i[i] = m_i[i] / Math.sqrt(eps);
-          }
-          let W2 = Math.pow(Correlate(a_i, sorted), 2);
-          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 - W2) - p_mu) / p_sigma;
-          return {
-            W: W2,
-            p: 1 - Distributions.Normal.cdf(z, 0, 1)
-          };
-        }
         let b = 0;
         for (let i = 0; i < Math.floor(N / 2); i++) {
           let a = Distributions.Normal.ShapiroWilk.constants[i][N - 2];
           b += a * (invSort[i] - sorted[i]);
         }
         let W = Math.pow(b, 2) / SS;
+        if (algorithmic) {
+          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;
+          return {
+            W,
+            p: 1 - Distributions.Normal.cdf(z, 0, 1)
+          };
+        }
         let p = Distributions.Normal.ShapiroWilk.p(W, N);
         return {
           W,
@@ -3179,7 +3192,7 @@ function CreateScatterPlot(data, title, chartSettings) {
   AddPlotBorder(ctx, chartSettings);
   return canvas;
 }
-function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSettings("red", "black", 1, 1), controlTests = [1, 2, 3, 4, 5, 6, 7, 8]) {
+function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSettings("red", "red", 1, 1), controlTests = [1]) {
   if (CL === void 0) {
     CL = Mean([LCL, UCL]);
   }
@@ -3200,7 +3213,7 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
     if (controlTests.includes(2) && (v.y > CL && Test2 >= 0 || v.y < CL && Test2 <= 0)) {
       Test2 += v.y < CL ? -1 : 1;
       if (Test2 >= 8 || Test2 <= -8) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("orange", "red", 1, 1);
       }
     } else {
       Test2 = 0;
@@ -3209,7 +3222,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
       if (controlTests.includes(3) && (v.y > data.values[i - 1].y && Test3 >= 0 || v.y < data.values[i - 1].y && Test3 <= 0)) {
         Test3 += v.y < data.values[i - 1].y ? -1 : 1;
         if (Test3 >= 5 || Test3 <= -6) {
-          v.display = display;
+          v.display = new CanvasDrawSettings("yellow", "red", 1, 1);
+          ;
         }
       } else {
         Test3 = 0;
@@ -3221,7 +3235,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
           Test4 = 0;
         }
         if (Math.abs(Test4) >= 13) {
-          v.display = display;
+          v.display = new CanvasDrawSettings("purple", "red", 1, 1);
+          ;
         }
       }
       if (controlTests.includes(4) && v.y > data.values[i - 1].y) {
@@ -3231,7 +3246,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
           Test4 = 0;
         }
         if (Math.abs(Test4) >= 13) {
-          v.display = display;
+          v.display = new CanvasDrawSettings("purple", "red", 1, 1);
+          ;
         }
       }
     }
@@ -3240,7 +3256,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
     if (controlTests.includes(5) && v.y > 2 * std + CL) {
       Test5[2] = 1;
       if (Test5[0] === 1 || Test5[1] === 1) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("white", "red", 1, 1);
+        ;
       }
     }
     if (controlTests.includes(5) && v.y < CL - 2 * std) {
@@ -3258,7 +3275,7 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
       Test6[5] = 1;
       let plus = Test6.filter((t) => t === 1);
       if (plus.length >= 4) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("grey", "red", 1, 1);
       }
     } else if (controlTests.includes(6) && v.y < CL - std) {
       Test6[5] = -1;
@@ -3272,7 +3289,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
     if (controlTests.includes(7) && v.y >= CL - std && v.y <= CL + std) {
       Test7 += 1;
       if (Test7 >= 15) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("pink", "red", 1, 1);
+        ;
       }
     } else {
       Test7 = 0;
@@ -3280,7 +3298,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
     if (controlTests.includes(8) && (v.y < CL - std || v.y > CL + std)) {
       Test8 += 1;
       if (Test8 >= 8) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("brown", "red", 1, 1);
+        ;
       }
     } else {
       Test8 = 0;
@@ -3634,12 +3653,33 @@ function QQPlotChart(data, distribution = "Normal") {
   let chart = CreateStackedChart([middleLineChart, new Chart([dataset], "scatter")], `${distribution[0].toUpperCase()}${distribution.substring(1)} Probability Plot`, chartSettings);
   if (distribution.toLowerCase() === "normal" || distribution.toLowerCase().includes("boxcox")) {
     let AD = GoodnessOfFit.AndersonDarling.Test(data, GoodnessOfFit.Distributions.Normal);
-    let RJ = Distributions.Normal.RyanJoiner.Test(data);
-    let KS = Distributions.Normal.KolmogorovSmirnov(data);
     let JB = Distributions.Normal.JarqueBera(data);
+    let KS = Distributions.Normal.KolmogorovSmirnov(data);
+    let LF = Distributions.Normal.Lilliefors.Test(data);
+    let RJ = Distributions.Normal.RyanJoiner.Test(data);
+    let SWR = Distributions.Normal.ShapiroWilk.Royston(data);
     let ctx = chart.getContext("2d");
     if (ctx) {
-      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }, { key: "", value: "" }, { key: "RJ:", value: RJ.RJ }, { key: "crit.:", value: RJ.critical }, { key: "", value: "" }, { key: "KS:", value: KS.KS }, { key: "p-value:", value: KS.p }, { key: "", value: "" }, { key: "JB:", value: JB.JB }, { key: "p-value:", value: JB.p }], 1075, 150, 175, 100, "", 30);
+      AddInfoTable(ctx, [
+        { key: "AD:", value: AD.AD },
+        { key: "p-value:", value: AD.p },
+        { key: "", value: "" },
+        { key: "JB:", value: JB.JB },
+        { key: "p-value:", value: JB.p },
+        { key: "", value: "" },
+        { key: "KS:", value: KS.KS },
+        { key: "p-value:", value: KS.p },
+        { key: "", value: "" },
+        { key: "Lf:", value: LF.KS },
+        { key: "p-value:", value: LF.p },
+        { key: "", value: "" },
+        { key: "RJ:", value: RJ.RJ },
+        { key: "crit.:", value: RJ.critical },
+        { key: "", value: "" },
+        { key: "SW(R):", value: SWR.W },
+        { key: "p-value:", value: SWR.p },
+        { key: "", value: "" }
+      ], 1075, 150, 175, 100, "", 30);
     }
   } else {
     let dist = distribution.replace("-", "").replace("arameter", "").replace(/\s/g, "");
@@ -3685,6 +3725,30 @@ function QQDistributionModificationX(p, distribution, data) {
   }
   return p;
 }
+function CreateControlChart(data, chartTitle, UCL, LCL) {
+  let centerline = Mean(data);
+  let std = StDev.S(data);
+  if (UCL === void 0) {
+    UCL = centerline + 3 * std;
+  }
+  if (LCL === void 0) {
+    LCL = centerline - 3 * std;
+  }
+  let dataSet = new DataSet(data.map((d, i) => {
+    return new Point(i, d);
+  }));
+  dataSet.line.type = "straight";
+  let customGridlines = [
+    new CustomGridline(centerline, `x\u0304 = ${Math.round(centerline * 1e3) / 1e3}`, new CanvasDrawSettings("", "green"), void 0, true),
+    new CustomGridline(LCL, `LCL = ${Math.round(LCL * 1e3) / 1e3}`, new CanvasDrawSettings("", "maroon"), void 0, true),
+    new CustomGridline(UCL, `UCL = ${Math.round(UCL * 1e3) / 1e3}`, new CanvasDrawSettings("", "maroon"), void 0, true)
+  ];
+  let chartSettings = new ChartSettings();
+  chartSettings.axis.y.customGridlines = customGridlines;
+  MarkOutOfControl(dataSet, UCL, LCL, centerline, std);
+  let graph = CreateScatterPlot([dataSet], chartTitle, chartSettings);
+  return graph;
+}
 function AnovaTableTwoWay(data, interaction) {
   let ANOVA;
   if (!interaction) {
@@ -4364,6 +4428,7 @@ function SplitObjectArrayByProp(array, prop) {
   CreateCapabilityPlot,
   CreateCategoricalBarGraph,
   CreateContinuousBarGraph,
+  CreateControlChart,
   CreateScatterPlot,
   CreateSplitGraph,
   CreateStackedChart,

package/dist/index.mjs

@@ -1584,7 +1584,37 @@ var Distributions = {
       }
     },
     ShapiroWilk: {
-      Test(data, expanded) {
+      Royston(data) {
+        let N = data.length;
+        let mean = Mean(data);
+        let sorted = [...data].sort();
+        let u = 1 / Math.sqrt(N);
+        let m_i = [];
+        let M = 0;
+        for (let i = 0; i < N; i++) {
+          let m = Distributions.Normal.inv((i + 1 - 0.375) / (N + 0.25));
+          m_i[i] = m;
+          M += Math.pow(m, 2);
+        }
+        let a_i = [];
+        a_i[N - 1] = -2.706056 * Math.pow(u, 5) + 4.434685 * Math.pow(u, 4) - 2.07119 * Math.pow(u, 3) - 0.147981 * Math.pow(u, 2) + 0.221157 * u + m_i[N - 1] / Math.sqrt(M);
+        a_i[N - 2] = -3.582633 * Math.pow(u, 5) + 5.682633 * Math.pow(u, 4) - 1.752461 * Math.pow(u, 3) - 0.293762 * Math.pow(u, 2) + 0.042981 * u + m_i[N - 2] / Math.sqrt(M);
+        a_i[0] = -1 * a_i[N - 1];
+        a_i[1] = -1 * a_i[N - 2];
+        let eps = (M - 2 * Math.pow(m_i[N - 1], 2) - 2 * Math.pow(m_i[N - 2], 2)) / (1 - 2 * Math.pow(a_i[N - 1], 2) - 2 * Math.pow(a_i[N - 2], 2));
+        for (let i = 2; i < N - 2; i++) {
+          a_i[i] = m_i[i] / Math.sqrt(eps);
+        }
+        let W = Math.pow(Correlate(a_i, sorted), 2);
+        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;
+        return {
+          W,
+          p: 1 - Distributions.Normal.cdf(z, 0, 1)
+        };
+      },
+      Test(data, algorithmic) {
         let N = data.length;
         let mean = Mean(data);
         let sorted = [...data].sort();
@@ -1593,39 +1623,21 @@ var Distributions = {
         sorted.forEach((s) => {
           SS += Math.pow(s - mean, 2);
         });
-        if (expanded || N > 50) {
-          let u = 1 / Math.sqrt(N);
-          let m_i = [];
-          let M = 0;
-          for (let i = 0; i < N; i++) {
-            let m = Distributions.Normal.inv((i + 1 - 0.375) / (N + 0.25));
-            m_i[i] = m;
-            M += Math.pow(m, 2);
-          }
-          let a_i = [];
-          a_i[N - 1] = -2.706056 * Math.pow(u, 5) + 4.434685 * Math.pow(u, 4) - 2.07119 * Math.pow(u, 3) - 0.147981 * Math.pow(u, 2) + 0.221157 * u + m_i[N - 1] / Math.sqrt(M);
-          a_i[N - 2] = -3.582633 * Math.pow(u, 5) + 5.682633 * Math.pow(u, 4) - 1.752461 * Math.pow(u, 3) - 0.293762 * Math.pow(u, 2) + 0.042981 * u + m_i[N - 2] / Math.sqrt(M);
-          a_i[0] = -1 * a_i[N - 1];
-          a_i[1] = -1 * a_i[N - 2];
-          let eps = (M - 2 * Math.pow(m_i[N - 1], 2) - 2 * Math.pow(m_i[N - 2], 2)) / (1 - 2 * Math.pow(a_i[N - 1], 2) - 2 * Math.pow(a_i[N - 2], 2));
-          for (let i = 2; i < N - 2; i++) {
-            a_i[i] = m_i[i] / Math.sqrt(eps);
-          }
-          let W2 = Math.pow(Correlate(a_i, sorted), 2);
-          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 - W2) - p_mu) / p_sigma;
-          return {
-            W: W2,
-            p: 1 - Distributions.Normal.cdf(z, 0, 1)
-          };
-        }
         let b = 0;
         for (let i = 0; i < Math.floor(N / 2); i++) {
           let a = Distributions.Normal.ShapiroWilk.constants[i][N - 2];
           b += a * (invSort[i] - sorted[i]);
         }
         let W = Math.pow(b, 2) / SS;
+        if (algorithmic) {
+          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;
+          return {
+            W,
+            p: 1 - Distributions.Normal.cdf(z, 0, 1)
+          };
+        }
         let p = Distributions.Normal.ShapiroWilk.p(W, N);
         return {
           W,
@@ -3093,7 +3105,7 @@ function CreateScatterPlot(data, title, chartSettings) {
   AddPlotBorder(ctx, chartSettings);
   return canvas;
 }
-function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSettings("red", "black", 1, 1), controlTests = [1, 2, 3, 4, 5, 6, 7, 8]) {
+function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSettings("red", "red", 1, 1), controlTests = [1]) {
   if (CL === void 0) {
     CL = Mean([LCL, UCL]);
   }
@@ -3114,7 +3126,7 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
     if (controlTests.includes(2) && (v.y > CL && Test2 >= 0 || v.y < CL && Test2 <= 0)) {
       Test2 += v.y < CL ? -1 : 1;
       if (Test2 >= 8 || Test2 <= -8) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("orange", "red", 1, 1);
       }
     } else {
       Test2 = 0;
@@ -3123,7 +3135,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
       if (controlTests.includes(3) && (v.y > data.values[i - 1].y && Test3 >= 0 || v.y < data.values[i - 1].y && Test3 <= 0)) {
         Test3 += v.y < data.values[i - 1].y ? -1 : 1;
         if (Test3 >= 5 || Test3 <= -6) {
-          v.display = display;
+          v.display = new CanvasDrawSettings("yellow", "red", 1, 1);
+          ;
         }
       } else {
         Test3 = 0;
@@ -3135,7 +3148,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
           Test4 = 0;
         }
         if (Math.abs(Test4) >= 13) {
-          v.display = display;
+          v.display = new CanvasDrawSettings("purple", "red", 1, 1);
+          ;
         }
       }
       if (controlTests.includes(4) && v.y > data.values[i - 1].y) {
@@ -3145,7 +3159,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
           Test4 = 0;
         }
         if (Math.abs(Test4) >= 13) {
-          v.display = display;
+          v.display = new CanvasDrawSettings("purple", "red", 1, 1);
+          ;
         }
       }
     }
@@ -3154,7 +3169,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
     if (controlTests.includes(5) && v.y > 2 * std + CL) {
       Test5[2] = 1;
       if (Test5[0] === 1 || Test5[1] === 1) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("white", "red", 1, 1);
+        ;
       }
     }
     if (controlTests.includes(5) && v.y < CL - 2 * std) {
@@ -3172,7 +3188,7 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
       Test6[5] = 1;
       let plus = Test6.filter((t) => t === 1);
       if (plus.length >= 4) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("grey", "red", 1, 1);
       }
     } else if (controlTests.includes(6) && v.y < CL - std) {
       Test6[5] = -1;
@@ -3186,7 +3202,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
     if (controlTests.includes(7) && v.y >= CL - std && v.y <= CL + std) {
       Test7 += 1;
       if (Test7 >= 15) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("pink", "red", 1, 1);
+        ;
       }
     } else {
       Test7 = 0;
@@ -3194,7 +3211,8 @@ function MarkOutOfControl(data, UCL, LCL, CL, std, display = new CanvasDrawSetti
     if (controlTests.includes(8) && (v.y < CL - std || v.y > CL + std)) {
       Test8 += 1;
       if (Test8 >= 8) {
-        v.display = display;
+        v.display = new CanvasDrawSettings("brown", "red", 1, 1);
+        ;
       }
     } else {
       Test8 = 0;
@@ -3548,12 +3566,33 @@ function QQPlotChart(data, distribution = "Normal") {
   let chart = CreateStackedChart([middleLineChart, new Chart([dataset], "scatter")], `${distribution[0].toUpperCase()}${distribution.substring(1)} Probability Plot`, chartSettings);
   if (distribution.toLowerCase() === "normal" || distribution.toLowerCase().includes("boxcox")) {
     let AD = GoodnessOfFit.AndersonDarling.Test(data, GoodnessOfFit.Distributions.Normal);
-    let RJ = Distributions.Normal.RyanJoiner.Test(data);
-    let KS = Distributions.Normal.KolmogorovSmirnov(data);
     let JB = Distributions.Normal.JarqueBera(data);
+    let KS = Distributions.Normal.KolmogorovSmirnov(data);
+    let LF = Distributions.Normal.Lilliefors.Test(data);
+    let RJ = Distributions.Normal.RyanJoiner.Test(data);
+    let SWR = Distributions.Normal.ShapiroWilk.Royston(data);
     let ctx = chart.getContext("2d");
     if (ctx) {
-      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }, { key: "", value: "" }, { key: "RJ:", value: RJ.RJ }, { key: "crit.:", value: RJ.critical }, { key: "", value: "" }, { key: "KS:", value: KS.KS }, { key: "p-value:", value: KS.p }, { key: "", value: "" }, { key: "JB:", value: JB.JB }, { key: "p-value:", value: JB.p }], 1075, 150, 175, 100, "", 30);
+      AddInfoTable(ctx, [
+        { key: "AD:", value: AD.AD },
+        { key: "p-value:", value: AD.p },
+        { key: "", value: "" },
+        { key: "JB:", value: JB.JB },
+        { key: "p-value:", value: JB.p },
+        { key: "", value: "" },
+        { key: "KS:", value: KS.KS },
+        { key: "p-value:", value: KS.p },
+        { key: "", value: "" },
+        { key: "Lf:", value: LF.KS },
+        { key: "p-value:", value: LF.p },
+        { key: "", value: "" },
+        { key: "RJ:", value: RJ.RJ },
+        { key: "crit.:", value: RJ.critical },
+        { key: "", value: "" },
+        { key: "SW(R):", value: SWR.W },
+        { key: "p-value:", value: SWR.p },
+        { key: "", value: "" }
+      ], 1075, 150, 175, 100, "", 30);
     }
   } else {
     let dist = distribution.replace("-", "").replace("arameter", "").replace(/\s/g, "");
@@ -3599,6 +3638,30 @@ function QQDistributionModificationX(p, distribution, data) {
   }
   return p;
 }
+function CreateControlChart(data, chartTitle, UCL, LCL) {
+  let centerline = Mean(data);
+  let std = StDev.S(data);
+  if (UCL === void 0) {
+    UCL = centerline + 3 * std;
+  }
+  if (LCL === void 0) {
+    LCL = centerline - 3 * std;
+  }
+  let dataSet = new DataSet(data.map((d, i) => {
+    return new Point(i, d);
+  }));
+  dataSet.line.type = "straight";
+  let customGridlines = [
+    new CustomGridline(centerline, `x\u0304 = ${Math.round(centerline * 1e3) / 1e3}`, new CanvasDrawSettings("", "green"), void 0, true),
+    new CustomGridline(LCL, `LCL = ${Math.round(LCL * 1e3) / 1e3}`, new CanvasDrawSettings("", "maroon"), void 0, true),
+    new CustomGridline(UCL, `UCL = ${Math.round(UCL * 1e3) / 1e3}`, new CanvasDrawSettings("", "maroon"), void 0, true)
+  ];
+  let chartSettings = new ChartSettings();
+  chartSettings.axis.y.customGridlines = customGridlines;
+  MarkOutOfControl(dataSet, UCL, LCL, centerline, std);
+  let graph = CreateScatterPlot([dataSet], chartTitle, chartSettings);
+  return graph;
+}
 function AnovaTableTwoWay(data, interaction) {
   let ANOVA;
   if (!interaction) {
@@ -4277,6 +4340,7 @@ export {
   CreateCapabilityPlot,
   CreateCategoricalBarGraph,
   CreateContinuousBarGraph,
+  CreateControlChart,
   CreateScatterPlot,
   CreateSplitGraph,
   CreateStackedChart,

package/index.ts

@@ -1748,7 +1748,44 @@ export const Distributions = {
       }
     },
     ShapiroWilk:{
-      Test(data: number[], expanded?: boolean){
+      Royston(data: number[]){
+        let N = data.length;
+        let mean = Mean(data);
+        
+        let sorted = [...data].sort();
+
+        let u = 1/Math.sqrt(N);
+        let m_i: number[] = [];
+        let M = 0;
+        for(let i = 0; i < N; i++){
+          let m = Distributions.Normal.inv((i+1-0.375)/(N+.25));
+          m_i[i] = m;
+          M += Math.pow(m,2);
+        }
+        
+        let a_i: number[] = [];
+        a_i[N - 1] = -2.706056*Math.pow(u,5)+4.434685*Math.pow(u,4)-2.07119*Math.pow(u,3)-0.147981*Math.pow(u,2)+0.221157*u+m_i[N - 1]/Math.sqrt(M);
+        a_i[N - 2] = -3.582633*Math.pow(u,5)+5.682633*Math.pow(u,4)-1.752461*Math.pow(u,3)-0.293762*Math.pow(u,2)+0.042981*u+m_i[N - 2]/Math.sqrt(M);
+        a_i[0] = -1 * a_i[N-1];
+        a_i[1] = -1 * a_i[N-2];
+        let eps = (M-2*Math.pow(m_i[N-1],2)-2*Math.pow(m_i[N-2],2))/(1-2*Math.pow(a_i[N-1],2)-2*Math.pow(a_i[N-2],2))
+        
+        for(let i = 2; i < N-2; i++){
+          a_i[i] = m_i[i]/Math.sqrt(eps);
+        }
+        
+        let W = Math.pow(Correlate(a_i, sorted), 2);
+
+        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;
+        
+        return {
+          W: W,
+          p: 1 - Distributions.Normal.cdf(z, 0, 1)
+        }
+      },
+      Test(data: number[], algorithmic?: boolean){
         let N = data.length;
         let mean = Mean(data);
         
@@ -1759,46 +1796,25 @@ export const Distributions = {
         sorted.forEach(s => {
           SS += Math.pow(s - mean, 2);
         })
-        // ShapiroWilk Royston Expanded Approximation
-        if(expanded || N > 50){
-          let u = 1/Math.sqrt(N);
-          let m_i: number[] = [];
-          let M = 0;
-          for(let i = 0; i < N; i++){
-            let m = Distributions.Normal.inv((i+1-0.375)/(N+.25));
-            m_i[i] = m;
-            M += Math.pow(m,2);
-          }
-          
-          let a_i: number[] = [];
-          a_i[N - 1] = -2.706056*Math.pow(u,5)+4.434685*Math.pow(u,4)-2.07119*Math.pow(u,3)-0.147981*Math.pow(u,2)+0.221157*u+m_i[N - 1]/Math.sqrt(M);
-          a_i[N - 2] = -3.582633*Math.pow(u,5)+5.682633*Math.pow(u,4)-1.752461*Math.pow(u,3)-0.293762*Math.pow(u,2)+0.042981*u+m_i[N - 2]/Math.sqrt(M);
-          a_i[0] = -1 * a_i[N-1];
-          a_i[1] = -1 * a_i[N-2];
-          let eps = (M-2*Math.pow(m_i[N-1],2)-2*Math.pow(m_i[N-2],2))/(1-2*Math.pow(a_i[N-1],2)-2*Math.pow(a_i[N-2],2))
-          
-          for(let i = 2; i < N-2; i++){
-            a_i[i] = m_i[i]/Math.sqrt(eps);
-          }
-          
-          let W = Math.pow(Correlate(a_i, sorted), 2);
+        // ShapiroWilk
+        let b = 0;
+        for(let i = 0; i < Math.floor(N/2); i++){
+          let a = Distributions.Normal.ShapiroWilk.constants[i][N-2];
+          b += a*(invSort[i]-sorted[i]);
+        }
+        let W = Math.pow(b, 2)/SS;
 
+        if(algorithmic){
           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;
-          
+
           return {
             W: W,
             p: 1 - Distributions.Normal.cdf(z, 0, 1)
           }
         }
-        // ShapiroWilk
-        let b = 0;
-        for(let i = 0; i < Math.floor(N/2); i++){
-          let a = Distributions.Normal.ShapiroWilk.constants[i][N-2];
-          b += a*(invSort[i]-sorted[i]);
-        }
-        let W = Math.pow(b, 2)/SS;
+
         let p = Distributions.Normal.ShapiroWilk.p(W, N);
         return {
           W: W,
@@ -3558,7 +3574,7 @@ export function CreateScatterPlot(data: DataSet[], title?: string, chartSettings
 * @param std Standard Deviation of Data
 * @param display Display settings for out-of-control points. DEFAULT: RED SQUARE
 */
-function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std?: any, display = new CanvasDrawSettings("red", "black", 1, 1), controlTests: number[] = [1,2,3,4,5,6,7,8]) {
+function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std?: any, display = new CanvasDrawSettings("red", "red", 1, 1), controlTests: number[] = [1]) {
   if(CL === undefined){CL = Mean([LCL, UCL])}
   if(std === undefined){std = (CL - LCL)/3}
   let Test2 = 0;
@@ -3577,7 +3593,7 @@ function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std
     if(controlTests.includes(2) && ((v.y > CL && Test2 >= 0) || (v.y < CL && Test2 <= 0))){
       Test2 += v.y < CL ? -1 : 1;
       if(Test2 >= 8 || Test2 <= -8){
-        v.display = display;
+        v.display = new CanvasDrawSettings("orange", "red", 1, 1);
       }
     }else{
       Test2 = 0;
@@ -3588,7 +3604,7 @@ function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std
       if(controlTests.includes(3) && ((v.y > data.values[i-1].y && Test3 >= 0) || (v.y < data.values[i-1].y && Test3 <= 0))){
         Test3 += v.y < data.values[i-1].y ? -1 : 1;
         if(Test3 >= 5 || Test3 <= -6){
-          v.display = display;
+          v.display = new CanvasDrawSettings("yellow", "red", 1, 1);;
         }
       }else{
         Test3 = 0;
@@ -3600,7 +3616,7 @@ function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std
           Test4 = -1*(Test4 + 1);
         }else{Test4 = 0}
         if(Math.abs(Test4) >= 13){
-          v.display = display;
+          v.display = new CanvasDrawSettings("purple", "red", 1, 1);;
         }
       }      
       if(controlTests.includes(4) && v.y > data.values[i-1].y){
@@ -3608,7 +3624,7 @@ function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std
           Test4 = -1*(Test4 - 1);
         }else{Test4 = 0}
         if(Math.abs(Test4) >= 13){
-          v.display = display;
+          v.display = new CanvasDrawSettings("purple", "red", 1, 1);;
         }
       }
     }
@@ -3619,7 +3635,7 @@ function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std
     if (controlTests.includes(5) && v.y > 2*std+CL) {
       Test5[2] =  1;
       if(Test5[0] === 1 || Test5[1] === 1){
-        v.display = display;
+        v.display = new CanvasDrawSettings("white", "red", 1, 1);;
       }
     }
     if (controlTests.includes(5) && v.y < CL-2*std) {
@@ -3639,7 +3655,7 @@ function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std
       Test6[5] =  1;
       let plus = Test6.filter(t => t === 1);
       if(plus.length >= 4){
-        v.display = display
+        v.display = new CanvasDrawSettings("grey", "red", 1, 1);
       }
     }else if (controlTests.includes(6) && v.y < CL-std) {
       Test6[5] =  -1;
@@ -3653,17 +3669,17 @@ function MarkOutOfControl(data: DataSet, UCL: number, LCL: number, CL?: any, std
     if(controlTests.includes(7) && v.y >= CL - std && v.y <= CL + std){
       Test7 += 1;
       if(Test7 >= 15){
-        v.display = display;
+        v.display = new CanvasDrawSettings("pink", "red", 1, 1);;
       }
     }else{
       Test7 = 0
     }
 
-    // Test7: Fifteen points in a row within 1σ of center line (either side)
+    // Test8: Fifteen points in a row within 1σ of center line (either side)
     if(controlTests.includes(8) && (v.y < CL - std || v.y > CL + std)){
       Test8 += 1;
       if(Test8 >= 8){
-        v.display = display;
+        v.display = new CanvasDrawSettings("brown", "red", 1, 1);;
       }
     }else{
       Test8 = 0
@@ -4149,12 +4165,20 @@ export function QQPlotChart(data: number[], distribution: string = 'Normal'){
     
   if (distribution.toLowerCase() === "normal" || distribution.toLowerCase().includes("boxcox")) {
     let AD = GoodnessOfFit.AndersonDarling.Test(data, GoodnessOfFit.Distributions.Normal);
-    let RJ = Distributions.Normal.RyanJoiner.Test(data);
-    let KS = Distributions.Normal.KolmogorovSmirnov(data);
     let JB = Distributions.Normal.JarqueBera(data);
+    let KS = Distributions.Normal.KolmogorovSmirnov(data);
+    let LF = Distributions.Normal.Lilliefors.Test(data);
+    let RJ = Distributions.Normal.RyanJoiner.Test(data);
+    let SWR = Distributions.Normal.ShapiroWilk.Royston(data);
     let ctx = chart.getContext("2d");
     if(ctx){
-      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }, { key: "", value: '' }, { key: "RJ:", value: RJ.RJ }, { key: "crit.:", value: RJ.critical }, { key: "", value: "" }, { key: "KS:", value: KS.KS }, { key: "p-value:", value: KS.p }, { key: "", value: "" }, { key: "JB:", value: JB.JB }, { key: "p-value:", value: JB.p }], 1075, 150, 175, 100, "", 30);
+      AddInfoTable(ctx, [
+        { key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }, { key: "", value: '' }, 
+        { key: "JB:", value: JB.JB }, { key: "p-value:", value: JB.p }, { key: "", value: "" }, 
+        { key: "KS:", value: KS.KS }, { key: "p-value:", value: KS.p }, { key: "", value: "" },  
+        { key: "Lf:", value: LF.KS }, { key: "p-value:", value: LF.p }, { key: "", value: "" },  
+        { key: "RJ:", value: RJ.RJ }, { key: "crit.:", value: RJ.critical }, { key: "", value: "" }, 
+        { key: "SW(R):", value: SWR.W }, { key: "p-value:", value: SWR.p }, { key: "", value: "" },], 1075, 150, 175, 100, "", 30);
     }
   }else{
     let dist = distribution.replace("-", "").replace("arameter", "").replace(/\s/g, "");
@@ -4222,7 +4246,7 @@ function QQDistributionModificationX(p: number, distribution: string, data: numb
   return p;
 }
 
-function CreateControlChart(data: number[], chartTitle?: string, UCL?: number, LCL?: number){
+export function CreateControlChart(data: number[], chartTitle?: string, UCL?: number, LCL?: number){
   let centerline = Mean(data);
   let std = StDev.S(data);
 

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "qesuite",
-  "version": "1.0.36",
+  "version": "1.0.38",
   "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_37.md

@@ -0,0 +1,14 @@
+- Separated RoystonShapiroWilk from ShapiroWilk.Test()
+- - ShapiroWilk.Test() can now evaluate p-value using the approximately normal distribution of the p-value using ln(1-W)
+- - ShapiroWilk.Royston() performs Royston Expanded test
+
+- Control Chart Tests are now identified by the color of the point fill
+    - Test 1: red
+    - Test 2: orange
+    - Test 3: yellow
+    - Test 4: purple
+    - Test 5: black
+    - Test 6: grey
+    - Test 7: pink
+    - Test 8: brown
+- Control chart defaults to only perform test 1

package/versions/1_0_38.md

@@ -0,0 +1,3 @@
+- include export for CreateControlChart();
+- include Shapiro-Wilk (Royston) and Lilliefors in QQPlot for normality
+- alphabetized the order of normality results table on QQPlot by test name