qesuite 1.0.58 → 1.0.59

package/dist/index.d.mts

@@ -1,5 +1,5 @@
 /** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- */
+*/
 declare class Specification {
     constructor(spec: any, USL?: number, LSL?: number, tolerance?: number, tol_plus?: number, tol_minus?: number, units?: string);
     nominal: number;
@@ -1211,8 +1211,10 @@ declare const ANOVA: {
     };
 };
 declare const MANOVA: {
+    WilksLambda(responses: any, factors: any, interactions?: any[]): any;
     Report(responses: any, factors: any, interactions?: any[]): any;
-    HypothesisMatrix(responses: any, factor: any): void;
+    ErrorMatrix(responses: any, factor: any): any;
+    HypothesisMatrix(responses: any, factor: any): any;
 };
 /**
 * Calculates the average moving range of the array supplied.
@@ -1717,6 +1719,18 @@ declare const OutlierTesting: {
         rejectNull: boolean;
         chart: HTMLCanvasElement;
     };
+    /**
+    * Performs Dixon's Test to test for and identify (if one exists) an outlier. H_0: data has no outliers | H_A: there is one outlier
+    * @param data The array of numbers to evaluate
+    * @param tails OPTIONAL: Whether to perform single or two-tailed Grubbs' Test. DEFAULT = 0: two-tailed | -1: lower tailed (checks if smallest value is an outlier) | 1: upper tailed (checks if largest value is an outlier)
+    * @param alpha OPTIONAL: The significance value to evaluate outliers. DEFAULT = 0.05
+   */
+    Dixon(data: any[], i?: number, j?: number, alpha?: number, name?: string): {
+        test: string;
+        r: number;
+        p: number;
+        chart: HTMLCanvasElement;
+    };
 };
 declare const Probability: {
     Combinations(n: number, r: number): number;
@@ -1890,7 +1904,7 @@ declare function CreateLayeredChart(charts: Chart[], title?: string, chartSettin
 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, controlTests?: any): HTMLCanvasElement;
-declare function GeneralizedLinearModel(dependent: number[], factors: any, covariates?: any, interactions?: any[][]): {
+declare function GeneralizedLinearModel(dependent: number[], factors: any, covariates?: any, interactions?: any[][], simple?: boolean): {
     ANOVA: any;
     Coefficients: {};
     ModelSummary: {
@@ -2276,5 +2290,6 @@ declare function GRRByPartChart(data: GRRData): HTMLCanvasElement;
 declare function NestedGRRByPartChart(data: GRRData): HTMLCanvasElement;
 declare function GRRPartxOperatorChart(data: GRRData): HTMLCanvasElement;
 declare function CompileGRRObject(operators: any, parts: any, measurements: any): GRRData;
+declare function Unique(values: any[]): any[];
 
-export { ANOMEChartEMP, ANOMRChartEMP, ANOVA, AnovaTableTwoWay, AttributeAgreementAnalysis, AverageMovingRange, AxisSettings, Beta, CanvasDrawSettings, Capability, Chart, ChartSettings, CohenKappa, CompileGRRObject, Correlation, Count, CreateBarChart, CreateBoxPlot, CreateBoxandWhiskerGraph, CreateCapabilityHistogram, CreateCapabilityPlot, CreateCategoricalBarGraph, CreateContinuousBarGraph, CreateControlChart, CreateLayeredChart, CreateMatrixChart, CreateScatterPlot, CreateSplitGraph, CreateStackedChart, CreateSummaryChart, CustomDraw, CustomGridline, CustomGridlines, DataCollection, DataSet, DescriptiveStatistics, Distributions, ERF, EngFont, EquivalenceTesting, EvaluateMeasurementProcess, FixedPointIteration, FleissKappa, G1, G1Graphs, GRR, GRRByPartChart, GRRData, GRRGraphs, GRROperator, GRRPartxOperatorChart, GRRReplication, GageEvaluation, GageLinearityAndBias, Gamma, GeneralizedLinearModel, GetFunctionValues, GetHistogramDataset, GoodnessOfFit, HypothesisTesting, IndividualDistributionIdentification, IndividualDistributionPlots, IndividualValuePlot, LastObservationsChart, Line, MANOVA, Margin, MaximumLikelihoodParameters, Mean, Median, MedianMovingRange, MovingRange, MovingRangeChart, MultipleLinearRegression, NestedGRR, NestedGRRByPartChart, NewtonRaphson, NonlinearRootFinding, NonparametricTesting, ObjectToArray, OutlierTesting, ParameterizedCDF, PlochhammerSymbol, Point, Probability, Product, QQPlot, QQPlotChart, Quartile, Rang, SerializeData, SimpleLinearRegression, Specification, StDev, StandardError, Sum, SumSq, SumSquaredDeviation, TrimmedMean, UnbiasingConstant, Variance, VarianceComponents };
+export { ANOMEChartEMP, ANOMRChartEMP, ANOVA, AnovaTableTwoWay, AttributeAgreementAnalysis, AverageMovingRange, AxisSettings, Beta, CanvasDrawSettings, Capability, Chart, ChartSettings, CohenKappa, CompileGRRObject, Correlation, Count, CreateBarChart, CreateBoxPlot, CreateBoxandWhiskerGraph, CreateCapabilityHistogram, CreateCapabilityPlot, CreateCategoricalBarGraph, CreateContinuousBarGraph, CreateControlChart, CreateLayeredChart, CreateMatrixChart, CreateScatterPlot, CreateSplitGraph, CreateStackedChart, CreateSummaryChart, CustomDraw, CustomGridline, CustomGridlines, DataCollection, DataSet, DescriptiveStatistics, Distributions, ERF, EngFont, EquivalenceTesting, EvaluateMeasurementProcess, FixedPointIteration, FleissKappa, G1, G1Graphs, GRR, GRRByPartChart, GRRData, GRRGraphs, GRROperator, GRRPartxOperatorChart, GRRReplication, GageEvaluation, GageLinearityAndBias, Gamma, GeneralizedLinearModel, GetFunctionValues, GetHistogramDataset, GoodnessOfFit, HypothesisTesting, IndividualDistributionIdentification, IndividualDistributionPlots, IndividualValuePlot, LastObservationsChart, Line, MANOVA, Margin, MaximumLikelihoodParameters, Mean, Median, MedianMovingRange, MovingRange, MovingRangeChart, MultipleLinearRegression, NestedGRR, NestedGRRByPartChart, NewtonRaphson, NonlinearRootFinding, NonparametricTesting, ObjectToArray, OutlierTesting, ParameterizedCDF, PlochhammerSymbol, Point, Probability, Product, QQPlot, QQPlotChart, Quartile, Rang, SerializeData, SimpleLinearRegression, Specification, StDev, StandardError, Sum, SumSq, SumSquaredDeviation, TrimmedMean, UnbiasingConstant, Unique, Variance, VarianceComponents };

package/dist/index.d.ts

@@ -1,5 +1,5 @@
 /** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- */
+*/
 declare class Specification {
     constructor(spec: any, USL?: number, LSL?: number, tolerance?: number, tol_plus?: number, tol_minus?: number, units?: string);
     nominal: number;
@@ -1211,8 +1211,10 @@ declare const ANOVA: {
     };
 };
 declare const MANOVA: {
+    WilksLambda(responses: any, factors: any, interactions?: any[]): any;
     Report(responses: any, factors: any, interactions?: any[]): any;
-    HypothesisMatrix(responses: any, factor: any): void;
+    ErrorMatrix(responses: any, factor: any): any;
+    HypothesisMatrix(responses: any, factor: any): any;
 };
 /**
 * Calculates the average moving range of the array supplied.
@@ -1717,6 +1719,18 @@ declare const OutlierTesting: {
         rejectNull: boolean;
         chart: HTMLCanvasElement;
     };
+    /**
+    * Performs Dixon's Test to test for and identify (if one exists) an outlier. H_0: data has no outliers | H_A: there is one outlier
+    * @param data The array of numbers to evaluate
+    * @param tails OPTIONAL: Whether to perform single or two-tailed Grubbs' Test. DEFAULT = 0: two-tailed | -1: lower tailed (checks if smallest value is an outlier) | 1: upper tailed (checks if largest value is an outlier)
+    * @param alpha OPTIONAL: The significance value to evaluate outliers. DEFAULT = 0.05
+   */
+    Dixon(data: any[], i?: number, j?: number, alpha?: number, name?: string): {
+        test: string;
+        r: number;
+        p: number;
+        chart: HTMLCanvasElement;
+    };
 };
 declare const Probability: {
     Combinations(n: number, r: number): number;
@@ -1890,7 +1904,7 @@ declare function CreateLayeredChart(charts: Chart[], title?: string, chartSettin
 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, controlTests?: any): HTMLCanvasElement;
-declare function GeneralizedLinearModel(dependent: number[], factors: any, covariates?: any, interactions?: any[][]): {
+declare function GeneralizedLinearModel(dependent: number[], factors: any, covariates?: any, interactions?: any[][], simple?: boolean): {
     ANOVA: any;
     Coefficients: {};
     ModelSummary: {
@@ -2276,5 +2290,6 @@ declare function GRRByPartChart(data: GRRData): HTMLCanvasElement;
 declare function NestedGRRByPartChart(data: GRRData): HTMLCanvasElement;
 declare function GRRPartxOperatorChart(data: GRRData): HTMLCanvasElement;
 declare function CompileGRRObject(operators: any, parts: any, measurements: any): GRRData;
+declare function Unique(values: any[]): any[];
 
-export { ANOMEChartEMP, ANOMRChartEMP, ANOVA, AnovaTableTwoWay, AttributeAgreementAnalysis, AverageMovingRange, AxisSettings, Beta, CanvasDrawSettings, Capability, Chart, ChartSettings, CohenKappa, CompileGRRObject, Correlation, Count, CreateBarChart, CreateBoxPlot, CreateBoxandWhiskerGraph, CreateCapabilityHistogram, CreateCapabilityPlot, CreateCategoricalBarGraph, CreateContinuousBarGraph, CreateControlChart, CreateLayeredChart, CreateMatrixChart, CreateScatterPlot, CreateSplitGraph, CreateStackedChart, CreateSummaryChart, CustomDraw, CustomGridline, CustomGridlines, DataCollection, DataSet, DescriptiveStatistics, Distributions, ERF, EngFont, EquivalenceTesting, EvaluateMeasurementProcess, FixedPointIteration, FleissKappa, G1, G1Graphs, GRR, GRRByPartChart, GRRData, GRRGraphs, GRROperator, GRRPartxOperatorChart, GRRReplication, GageEvaluation, GageLinearityAndBias, Gamma, GeneralizedLinearModel, GetFunctionValues, GetHistogramDataset, GoodnessOfFit, HypothesisTesting, IndividualDistributionIdentification, IndividualDistributionPlots, IndividualValuePlot, LastObservationsChart, Line, MANOVA, Margin, MaximumLikelihoodParameters, Mean, Median, MedianMovingRange, MovingRange, MovingRangeChart, MultipleLinearRegression, NestedGRR, NestedGRRByPartChart, NewtonRaphson, NonlinearRootFinding, NonparametricTesting, ObjectToArray, OutlierTesting, ParameterizedCDF, PlochhammerSymbol, Point, Probability, Product, QQPlot, QQPlotChart, Quartile, Rang, SerializeData, SimpleLinearRegression, Specification, StDev, StandardError, Sum, SumSq, SumSquaredDeviation, TrimmedMean, UnbiasingConstant, Variance, VarianceComponents };
+export { ANOMEChartEMP, ANOMRChartEMP, ANOVA, AnovaTableTwoWay, AttributeAgreementAnalysis, AverageMovingRange, AxisSettings, Beta, CanvasDrawSettings, Capability, Chart, ChartSettings, CohenKappa, CompileGRRObject, Correlation, Count, CreateBarChart, CreateBoxPlot, CreateBoxandWhiskerGraph, CreateCapabilityHistogram, CreateCapabilityPlot, CreateCategoricalBarGraph, CreateContinuousBarGraph, CreateControlChart, CreateLayeredChart, CreateMatrixChart, CreateScatterPlot, CreateSplitGraph, CreateStackedChart, CreateSummaryChart, CustomDraw, CustomGridline, CustomGridlines, DataCollection, DataSet, DescriptiveStatistics, Distributions, ERF, EngFont, EquivalenceTesting, EvaluateMeasurementProcess, FixedPointIteration, FleissKappa, G1, G1Graphs, GRR, GRRByPartChart, GRRData, GRRGraphs, GRROperator, GRRPartxOperatorChart, GRRReplication, GageEvaluation, GageLinearityAndBias, Gamma, GeneralizedLinearModel, GetFunctionValues, GetHistogramDataset, GoodnessOfFit, HypothesisTesting, IndividualDistributionIdentification, IndividualDistributionPlots, IndividualValuePlot, LastObservationsChart, Line, MANOVA, Margin, MaximumLikelihoodParameters, Mean, Median, MedianMovingRange, MovingRange, MovingRangeChart, MultipleLinearRegression, NestedGRR, NestedGRRByPartChart, NewtonRaphson, NonlinearRootFinding, NonparametricTesting, ObjectToArray, OutlierTesting, ParameterizedCDF, PlochhammerSymbol, Point, Probability, Product, QQPlot, QQPlotChart, Quartile, Rang, SerializeData, SimpleLinearRegression, Specification, StDev, StandardError, Sum, SumSq, SumSquaredDeviation, TrimmedMean, UnbiasingConstant, Unique, Variance, VarianceComponents };

package/dist/index.js

@@ -119,6 +119,7 @@ __export(QESuite_exports, {
   SumSquaredDeviation: () => SumSquaredDeviation,
   TrimmedMean: () => TrimmedMean,
   UnbiasingConstant: () => UnbiasingConstant,
+  Unique: () => Unique,
   Variance: () => Variance,
   VarianceComponents: () => VarianceComponents
 });
@@ -5586,7 +5587,7 @@ var ANOVA = {
           let rp = ANOVA.PostHoc.Duncan.Rp(treatmentNumber, DFE);
           let Rp = rp * S_x;
           returnArray.push({
-            treatment: infoArray[i].treatment,
+            Treatment: infoArray[i].treatment,
             Mean: infoArray[i].Mean,
             Rp,
             D: infoArray[i].Mean - Rp,
@@ -5711,8 +5712,8 @@ var ANOVA = {
             let diff = infoArray[i].Mean - infoArray[rev].Mean;
             let t = diff / Math.sqrt(MSE * (1 / infoArray[i].N + 1 / infoArray[rev].N));
             returnArray.push({
-              name: `${infoArray[i].treatment} - ${infoArray[rev].treatment}`,
-              meanDifference: diff,
+              Name: `${infoArray[i].treatment} - ${infoArray[rev].treatment}`,
+              MeanDifference: diff,
               T: t,
               LSD: tCrit * Math.sqrt(MSE * (1 / infoArray[i].N + 1 / infoArray[rev].N)),
               p: 2 * Distributions.T.RightTail(t, DFE)
@@ -5738,7 +5739,7 @@ var ANOVA = {
           let mean = Mean(data[key]);
           let tCrit = q_prime * Math.sqrt(MSE / n);
           returnArray.push({
-            treatment: key,
+            Treatment: key,
             N: n,
             Mean: mean,
             T: tCrit,
@@ -6628,6 +6629,73 @@ var ANOVA = {
   }
 };
 var MANOVA = {
+  WilksLambda(responses, factors, interactions = []) {
+    let rKeys = Object.keys(responses);
+    let fKeys = Object.keys(factors);
+    let coefficientCount = Sum(fKeys.map((k) => {
+      return Unique(factors[k]).length;
+    })) - 1 + interactions.length;
+    let dfError = responses[rKeys[0]].length - coefficientCount - 1;
+    let varTerms = [...fKeys];
+    let p = rKeys.length;
+    let v = dfError;
+    let returner = [];
+    fKeys.map((fk) => {
+      let factor = factors[fk];
+      let q = Unique(factor).length - 1;
+      let E = MANOVA.ErrorMatrix(responses, factor);
+      let H = MANOVA.HypothesisMatrix(responses, factor);
+      let lambda = Matrix.Determinant(E) / Matrix.Determinant(Matrix.Sum(H, E));
+      let m = 0.5 * (Math.abs(p - 1) - 1);
+      let s = Math.min(p, 1);
+      let u = 0.25 * (p * q - 2);
+      let n = 0.5 * (v - p - 1);
+      let r = v - 0.5 * (p - q + 1);
+      let t = Math.pow(p, 2) + Math.pow(q, 2) - 5 > 0 ? Math.sqrt((Math.pow(p, 2) * Math.pow(q, 2) - 4) / (Math.pow(p, 2) + Math.pow(q, 2) - 5)) : 1;
+      let f = (1 - Math.pow(lambda, 1 / t)) / Math.pow(lambda, 1 / t) * (r * t - 2 * u) / (p * q);
+      let denDF = r * t - 2 * u;
+      let numDF = p * q;
+      returner.push({
+        DenominatorDF: denDF,
+        F: f,
+        Factor: fk,
+        Lambda: lambda,
+        NumeratorDF: numDF,
+        p: Distributions.F.RightTail(f, denDF, numDF)
+      });
+    });
+    interactions.forEach((int) => {
+      let factor = [];
+      for (let i = 0; i < responses[rKeys[0]].length; i++) {
+        factor.push(int.map((inter) => {
+          return factors[inter][i];
+        }).join("*"));
+      }
+      console.log(factor);
+      let q = Unique(factor).length - 1;
+      let E = MANOVA.ErrorMatrix(responses, factor);
+      let H = MANOVA.HypothesisMatrix(responses, factor);
+      let lambda = Matrix.Determinant(E) / Matrix.Determinant(Matrix.Sum(H, E));
+      let m = 0.5 * (Math.abs(p - 1) - 1);
+      let s = Math.min(p, 1);
+      let u = 0.25 * (p * q - 2);
+      let n = 0.5 * (v - p - 1);
+      let r = v - 0.5 * (p - q + 1);
+      let t = Math.pow(p, 2) + Math.pow(q, 2) - 5 > 0 ? Math.sqrt((Math.pow(p, 2) * Math.pow(q, 2) - 4) / (Math.pow(p, 2) + Math.pow(q, 2) - 5)) : 1;
+      let f = (1 - Math.pow(lambda, 1 / t)) / Math.pow(lambda, 1 / t) * (r * t - 2 * u) / (p * q);
+      let denDF = r * t - 2 * u;
+      let numDF = p * q;
+      returner.push({
+        DenominatorDF: denDF,
+        F: f,
+        Factor: int.join("*"),
+        Lambda: lambda,
+        NumeratorDF: numDF,
+        p: Distributions.F.RightTail(f, denDF, numDF)
+      });
+    });
+    return returner;
+  },
   Report(responses, factors, interactions = []) {
     let returner = [];
     let rKeys = Object.keys(responses);
@@ -6647,7 +6715,6 @@ var MANOVA = {
       let u = 0.25 * (p * q - 2);
       let m = 0.5 * (Math.abs(p - q) - 1);
       let n = 0.5 * (v - p - 1);
-      console.log({ m, n, p, q, r, s, t, u, v });
       let TotalMatrix = Matrix.New(p, p);
       let ErrorMatrix = Matrix.New(p, p);
       Matrix.Fill(ErrorMatrix, 0);
@@ -6674,6 +6741,7 @@ var MANOVA = {
         let factorSSCP = Matrix.Multiply(Matrix.Transpose(rDevMatrix), rDevMatrix);
         ErrorMatrix = Matrix.Sum(ErrorMatrix, factorSSCP);
       });
+      console.log({ DeviationMatrix, nMatrix });
       let devXn = Matrix.Scale(DeviationMatrix, nMatrix);
       let HypothesisMatrix = Matrix.Multiply(Matrix.Transpose(DeviationMatrix), devXn);
       console.log({ HypothesisMatrix, ErrorMatrix });
@@ -6736,7 +6804,49 @@ var MANOVA = {
     });
     return returner;
   },
+  ErrorMatrix(responses, factor) {
+    let factorLevels = Unique(factor);
+    let errorMatrix = Matrix.New(Object.keys(responses).length, Object.keys(responses).length);
+    Matrix.Fill(errorMatrix, 0);
+    factorLevels.forEach((lv, lvi) => {
+      let nResponses = factor.filter((f) => f == lv).length;
+      let filteredResponses = {};
+      let responseMeans = {};
+      Object.keys(responses).forEach((k) => {
+        let subgroup = responses[k].filter((r, ri) => factor[ri] == lv);
+        filteredResponses[k] = subgroup;
+        responseMeans[k] = Mean(subgroup);
+      });
+      let treatmentMatrix = Matrix.New(nResponses, Object.keys(responses).length);
+      let deviationMatrix = Matrix.New(nResponses, Object.keys(responses).length);
+      for (let i = 0; i < nResponses; i++) {
+        Object.keys(responses).forEach((k, ki) => {
+          treatmentMatrix[i][ki] = filteredResponses[k][i];
+          deviationMatrix[i][ki] = filteredResponses[k][i] - responseMeans[k];
+        });
+      }
+      let SSCP = Matrix.Multiply(Matrix.Transpose(deviationMatrix), treatmentMatrix);
+      errorMatrix = Matrix.Sum(errorMatrix, SSCP);
+    });
+    return errorMatrix;
+  },
   HypothesisMatrix(responses, factor) {
+    let factorGroups = Unique(factor);
+    let responseMeans = {};
+    let deviationMatrix = Matrix.New(factorGroups.length, Object.keys(responses).length);
+    let nMatrix = Matrix.New(factorGroups.length, Object.keys(responses).length);
+    Object.keys(responses).forEach((k, ki) => {
+      let responseMean = Mean(responses[k]);
+      responseMeans[k] = responseMean;
+      factorGroups.forEach((g, gi) => {
+        let subgroup = responses[k].filter((d, di) => factor[di] == g);
+        deviationMatrix[gi][ki] = Mean(subgroup) - responseMean;
+        nMatrix[gi][ki] = subgroup.length;
+      });
+    });
+    let devXn = Matrix.Scale(deviationMatrix, nMatrix);
+    let hypothesisMatrix = Matrix.Multiply(Matrix.Transpose(deviationMatrix), devXn);
+    return hypothesisMatrix;
   }
 };
 function AverageMovingRange(data, w = 2) {
@@ -8211,6 +8321,54 @@ var OutlierTesting = {
       rejectNull: p <= alpha,
       chart: outlierPlot
     };
+  },
+  /**
+   * Performs Dixon's Test to test for and identify (if one exists) an outlier. H_0: data has no outliers | H_A: there is one outlier
+   * @param data The array of numbers to evaluate
+   * @param tails OPTIONAL: Whether to perform single or two-tailed Grubbs' Test. DEFAULT = 0: two-tailed | -1: lower tailed (checks if smallest value is an outlier) | 1: upper tailed (checks if largest value is an outlier)
+   * @param alpha OPTIONAL: The significance value to evaluate outliers. DEFAULT = 0.05
+  */
+  Dixon(data, i = 1, j = 0, alpha = 0.05, name) {
+    let y = [...data].sort((a, b) => a - b);
+    let n = y.length;
+    let r_ij = Math.max((y[i] - y[0]) / (y[n - 1 - j] - y[0]), (y[n - 1] - y[n - 1 - i]) / (y[n - 1] - y[j]));
+    let C = Probability.Factorialize(n) * Math.pow(2 * Math.PI, -3 / 2) / (Probability.Factorialize(j) * Probability.Factorialize(n - i - j - 2) * Probability.Factorialize(i - 1));
+    const j_xvr = (x, v, r) => {
+    };
+    const t = (r, v) => {
+      return (1 + Math.pow(r, 2)) * Math.pow(v, 2) / 2;
+    };
+    const u2 = (x) => {
+      return 3 * Math.pow(x, 2) / 2;
+    };
+    let p = NaN;
+    let outliers = [];
+    if (p <= alpha) {
+      if ((y[i] - y[0]) / (y[n - 1 - j] - y[0]) > (y[n - 1] - y[n - 1 - i]) / (y[n - 1] - y[j])) {
+        outliers = [{ i: data.findIndex((d) => d == y[0]), value: y[0] }];
+      } else {
+        outliers = [{ i: data.findIndex((d) => d == y[n - 1]), value: y[n - 1] }];
+      }
+    }
+    let plotData = data.map((d, i2) => {
+      if (outliers.findIndex((o) => o.i == i2) >= 0) {
+        return new Point(d, 1, new CanvasDrawSettings("red", "red", 1, 1));
+      } else {
+        return new Point(d, 1, new CanvasDrawSettings("#2183b2", "#2183b2"));
+      }
+    });
+    let plotSettings = new ChartSettings();
+    plotSettings.axis.y = new AxisSettings(0, 2);
+    plotSettings.axis.y.showAutoGridlines = false;
+    plotSettings.axis.y.customGridlines = [new CustomGridline(1)];
+    let title = name ? `Outlier Plot of ${name}` : `Outlier Plot`;
+    let outlierPlot = CreateScatterPlot([new DataSet(plotData)], title, plotSettings);
+    return {
+      test: i == 1 && j == 0 ? "Dixon's Q" : `Dixon's r${i}${j}`,
+      r: r_ij,
+      p,
+      chart: outlierPlot
+    };
   }
 };
 var Probability = {
@@ -8540,7 +8698,7 @@ var UnbiasingConstant = {
     return 3.4873 + 0.0250141 * N - 9823e-8 * Math.pow(N, 2);
   },
   D2prime(g, m) {
-    let d2PrimeTable = [];
+    let d2PrimeTable = [[{ "df": 1, "d2": 1.41421 }, { "df": 2, "d2": 1.91155 }, { "df": 2.9, "d2": 2.23887 }, { "df": 3.8, "d2": 2.48124 }, { "df": 4.7, "d2": 2.67253 }, { "df": 5.5, "d2": 2.82981 }, { "df": 6.3, "d2": 2.96288 }, { "df": 7, "d2": 3.07794 }, { "df": 7.7, "d2": 3.17905 }, { "df": 8.3, "d2": 3.26909 }, { "df": 9, "d2": 3.35016 }, { "df": 9.6, "d2": 3.42378 }, { "df": 10.2, "d2": 3.49116 }, { "df": 10.8, "d2": 3.55333 }, { "df": 11.3, "d2": 3.61071 }, { "df": 11.9, "d2": 3.66422 }, { "df": 12.4, "d2": 3.71424 }, { "df": 12.9, "d2": 3.76118 }, { "df": 13.4, "d2": 3.80537 }], [{ "df": 1.9, "d2": 1.27931 }, { "df": 3.8, "d2": 1.80538 }, { "df": 5.7, "d2": 2.15069 }, { "df": 7.5, "d2": 2.40484 }, { "df": 9.2, "d2": 2.60438 }, { "df": 10.8, "d2": 2.76779 }, { "df": 12.3, "d2": 2.90562 }, { "df": 13.8, "d2": 3.02446 }, { "df": 15.1, "d2": 3.12869 }, { "df": 16.5, "d2": 3.22134 }, { "df": 17.8, "d2": 3.30463 }, { "df": 19, "d2": 3.38017 }, { "df": 20.2, "d2": 3.44922 }, { "df": 21.3, "d2": 3.51287 }, { "df": 22.4, "d2": 3.57156 }, { "df": 23.5, "d2": 3.62625 }, { "df": 24.5, "d2": 3.67734 }, { "df": 25.5, "d2": 3.72524 }, { "df": 26.5, "d2": 3.77032 }], [{ "df": 2.8, "d2": 1.23105 }, { "df": 5.7, "d2": 1.76858 }, { "df": 8.4, "d2": 2.12049 }, { "df": 11.1, "d2": 2.37883 }, { "df": 13.6, "d2": 2.58127 }, { "df": 16, "d2": 2.74681 }, { "df": 18.3, "d2": 2.88628 }, { "df": 20.5, "d2": 3.00643 }, { "df": 22.6, "d2": 3.11173 }, { "df": 24.6, "d2": 3.20526 }, { "df": 26.5, "d2": 3.28931 }, { "df": 28.4, "d2": 3.3655 }, { "df": 30.1, "d2": 3.43512 }, { "df": 31.9, "d2": 3.49927 }, { "df": 33.5, "d2": 3.55842 }, { "df": 35.1, "d2": 3.61351 }, { "df": 36.7, "d2": 3.66495 }, { "df": 38.2, "d2": 3.71319 }, { "df": 39.7, "d2": 3.75857 }], [{ "df": 3.7, "d2": 1.20621 }, { "df": 7.5, "d2": 1.74989 }, { "df": 11.2, "d2": 2.10522 }, { "df": 14.7, "d2": 2.36571 }, { "df": 18.1, "d2": 2.56964 }, { "df": 21.3, "d2": 2.73626 }, { "df": 24.4, "d2": 2.87656 }, { "df": 27.3, "d2": 2.99737 }, { "df": 30.1, "d2": 3.10321 }, { "df": 32.7, "d2": 3.1972 }, { "df": 35.3, "d2": 3.28163 }, { "df": 37.7, "d2": 3.35815 }, { "df": 40.1, "d2": 3.42805 }, { "df": 42.4, "d2": 3.49246 }, { "df": 44.6, "d2": 3.55183 }, { "df": 46.7, "d2": 3.60712 }, { "df": 48.8, "d2": 3.65875 }, { "df": 50.8, "d2": 3.70715 }, { "df": 52.8, "d2": 3.75268 }], [{ "df": 4.6, "d2": 1.19105 }, { "df": 9.3, "d2": 1.73857 }, { "df": 13.9, "d2": 2.09601 }, { "df": 18.4, "d2": 2.35781 }, { "df": 22.6, "d2": 2.56263 }, { "df": 26.6, "d2": 2.72991 }, { "df": 30.4, "d2": 2.87071 }, { "df": 34, "d2": 2.99192 }, { "df": 37.5, "d2": 3.09808 }, { "df": 40.8, "d2": 3.19235 }, { "df": 44, "d2": 3.27701 }, { "df": 47.1, "d2": 3.35372 }, { "df": 50.1, "d2": 3.42381 }, { "df": 52.9, "d2": 3.48836 }, { "df": 55.7, "d2": 3.54787 }, { "df": 58.4, "d2": 3.60328 }, { "df": 61, "d2": 3.65502 }, { "df": 63.5, "d2": 3.70352 }, { "df": 65.9, "d2": 3.74914 }], [{ "df": 5.5, "d2": 1.18083 }, { "df": 11.1, "d2": 1.73099 }, { "df": 16.7, "d2": 2.08985 }, { "df": 22, "d2": 2.35253 }, { "df": 27, "d2": 2.55795 }, { "df": 31.8, "d2": 2.72567 }, { "df": 36.4, "d2": 2.8668 }, { "df": 40.8, "d2": 2.98829 }, { "df": 45, "d2": 3.09467 }, { "df": 49, "d2": 3.18911 }, { "df": 52.8, "d2": 3.27392 }, { "df": 56.5, "d2": 3.35077 }, { "df": 60.1, "d2": 3.42097 }, { "df": 63.5, "d2": 3.48563 }, { "df": 66.8, "d2": 3.54522 }, { "df": 70, "d2": 3.60072 }, { "df": 73.1, "d2": 3.65253 }, { "df": 76.1, "d2": 3.70109 }, { "df": 79.1, "d2": 3.74678 }], [{ "df": 6.4, "d2": 1.17348 }, { "df": 12.9, "d2": 1.72555 }, { "df": 19.4, "d2": 2.08543 }, { "df": 25.6, "d2": 2.34875 }, { "df": 31.5, "d2": 2.5546 }, { "df": 37.1, "d2": 2.72263 }, { "df": 42.5, "d2": 2.86401 }, { "df": 47.6, "d2": 2.98568 }, { "df": 52.4, "d2": 3.09222 }, { "df": 57.1, "d2": 3.18679 }, { "df": 61.6, "d2": 3.27172 }, { "df": 65.9, "d2": 3.34866 }, { "df": 70, "d2": 3.41894 }, { "df": 74, "d2": 3.48368 }, { "df": 77.9, "d2": 3.54333 }, { "df": 81.6, "d2": 3.59888 }, { "df": 85.3, "d2": 3.65075 }, { "df": 88.8, "d2": 3.69936 }, { "df": 92.2, "d2": 3.74509 }], [{ "df": 7.2, "d2": 1.16794 }, { "df": 14.8, "d2": 1.72147 }, { "df": 22.1, "d2": 2.08212 }, { "df": 29.2, "d2": 2.34591 }, { "df": 36, "d2": 2.55208 }, { "df": 42.4, "d2": 2.72036 }, { "df": 48.5, "d2": 2.86192 }, { "df": 54.3, "d2": 2.98373 }, { "df": 59.9, "d2": 3.09039 }, { "df": 65.2, "d2": 3.18506 }, { "df": 70.3, "d2": 3.27006 }, { "df": 75.2, "d2": 3.34708 }, { "df": 80, "d2": 3.41742 }, { "df": 84.6, "d2": 3.48221 }, { "df": 89, "d2": 3.54192 }, { "df": 93.3, "d2": 3.59751 }, { "df": 97.4, "d2": 3.64941 }, { "df": 101.4, "d2": 3.69806 }, { "df": 105.3, "d2": 3.74382 }], [{ "df": 8.1, "d2": 1.16361 }, { "df": 16.6, "d2": 1.71828 }, { "df": 24.9, "d2": 2.07953 }, { "df": 32.9, "d2": 2.3437 }, { "df": 40.4, "d2": 2.55013 }, { "df": 47.7, "d2": 2.71858 }, { "df": 54.5, "d2": 2.86028 }, { "df": 61.1, "d2": 2.98221 }, { "df": 67.3, "d2": 3.08896 }, { "df": 73.3, "d2": 3.1837 }, { "df": 79.1, "d2": 3.26878 }, { "df": 84.6, "d2": 3.34585 }, { "df": 90, "d2": 3.41624 }, { "df": 95.1, "d2": 3.48107 }, { "df": 100.1, "d2": 3.54081 }, { "df": 104.9, "d2": 3.59644 }, { "df": 109.5, "d2": 3.64838 }, { "df": 114.1, "d2": 3.69705 }, { "df": 118.5, "d2": 3.74284 }], [{ "df": 9, "d2": 1.16014 }, { "df": 18.4, "d2": 1.71573 }, { "df": 27.6, "d2": 2.07746 }, { "df": 36.5, "d2": 2.34192 }, { "df": 44.9, "d2": 2.54856 }, { "df": 52.9, "d2": 2.71717 }, { "df": 60.6, "d2": 2.85898 }, { "df": 67.8, "d2": 2.981 }, { "df": 74.8, "d2": 3.08781 }, { "df": 81.5, "d2": 3.18262 }, { "df": 87.9, "d2": 3.26775 }, { "df": 94, "d2": 3.34486 }, { "df": 99.9, "d2": 3.41529 }, { "df": 105.6, "d2": 3.48016 }, { "df": 111.2, "d2": 3.53993 }, { "df": 116.5, "d2": 3.59559 }, { "df": 121.7, "d2": 3.64755 }, { "df": 126.7, "d2": 3.69625 }, { "df": 131.6, "d2": 3.74205 }], [{ "df": 9.9, "d2": 1.15729 }, { "df": 20.2, "d2": 1.71363 }, { "df": 30.4, "d2": 2.07577 }, { "df": 40.1, "d2": 2.34048 }, { "df": 49.4, "d2": 2.54728 }, { "df": 58.2, "d2": 2.716 }, { "df": 66.6, "d2": 2.85791 }, { "df": 74.6, "d2": 2.98 }, { "df": 82.2, "d2": 3.08688 }, { "df": 89.6, "d2": 3.18174 }, { "df": 96.6, "d2": 3.2669 }, { "df": 103.4, "d2": 3.34406 }, { "df": 109.9, "d2": 3.41452 }, { "df": 116.2, "d2": 3.47941 }, { "df": 122.3, "d2": 3.53921 }, { "df": 128.1, "d2": 3.59489 }, { "df": 133.8, "d2": 3.64687 }, { "df": 139.4, "d2": 3.69558 }, { "df": 144.7, "d2": 3.74141 }], [{ "df": 10.7, "d2": 1.1549 }, { "df": 22, "d2": 1.71189 }, { "df": 33.1, "d2": 2.07436 }, { "df": 43.7, "d2": 2.33927 }, { "df": 53.8, "d2": 2.54621 }, { "df": 63.5, "d2": 2.71504 }, { "df": 72.6, "d2": 2.85702 }, { "df": 81.3, "d2": 2.97917 }, { "df": 89.7, "d2": 3.0861 }, { "df": 97.7, "d2": 3.181 }, { "df": 105.4, "d2": 3.2662 }, { "df": 112.7, "d2": 3.34339 }, { "df": 119.9, "d2": 3.41387 }, { "df": 126.7, "d2": 3.47879 }, { "df": 133.3, "d2": 3.53861 }, { "df": 139.8, "d2": 3.5943 }, { "df": 146, "d2": 3.6463 }, { "df": 152, "d2": 3.69503 }, { "df": 157.9, "d2": 3.74087 }], [{ "df": 11.6, "d2": 1.15289 }, { "df": 23.8, "d2": 1.71041 }, { "df": 35.8, "d2": 2.07316 }, { "df": 47.3, "d2": 2.33824 }, { "df": 58.3, "d2": 2.5453 }, { "df": 68.7, "d2": 2.71422 }, { "df": 78.6, "d2": 2.85627 }, { "df": 88.1, "d2": 2.97847 }, { "df": 97.1, "d2": 3.08544 }, { "df": 105.8, "d2": 3.18037 }, { "df": 114.1, "d2": 3.26561 }, { "df": 122.1, "d2": 3.34282 }, { "df": 129.8, "d2": 3.41333 }, { "df": 137.3, "d2": 3.47826 }, { "df": 144.4, "d2": 3.5381 }, { "df": 151.4, "d2": 3.59381 }, { "df": 158.1, "d2": 3.64582 }, { "df": 164.7, "d2": 3.69457 }, { "df": 171, "d2": 3.74041 }], [{ "df": 12.5, "d2": 1.15115 }, { "df": 25.7, "d2": 1.70914 }, { "df": 38.6, "d2": 2.07213 }, { "df": 51, "d2": 2.33737 }, { "df": 62.8, "d2": 2.54452 }, { "df": 74, "d2": 2.71351 }, { "df": 84.7, "d2": 2.85562 }, { "df": 94.9, "d2": 2.97787 }, { "df": 104.6, "d2": 3.08487 }, { "df": 113.9, "d2": 3.17984 }, { "df": 122.9, "d2": 3.2651 }, { "df": 131.5, "d2": 3.34233 }, { "df": 139.8, "d2": 3.41286 }, { "df": 147.8, "d2": 3.47781 }, { "df": 155.5, "d2": 3.53766 }, { "df": 163, "d2": 3.59339 }, { "df": 170.3, "d2": 3.64541 }, { "df": 177.3, "d2": 3.69417 }, { "df": 184.2, "d2": 3.74002 }], [{ "df": 13.4, "d2": 1.14965 }, { "df": 27.5, "d2": 1.70804 }, { "df": 41.3, "d2": 2.07125 }, { "df": 54.6, "d2": 2.33661 }, { "df": 67.2, "d2": 2.54385 }, { "df": 79.3, "d2": 2.7129 }, { "df": 90.7, "d2": 2.85506 }, { "df": 101.6, "d2": 2.97735 }, { "df": 112.1, "d2": 3.08438 }, { "df": 122.1, "d2": 3.17938 }, { "df": 131.7, "d2": 3.26465 }, { "df": 140.9, "d2": 3.34191 }, { "df": 149.8, "d2": 3.41245 }, { "df": 158.3, "d2": 3.47742 }, { "df": 166.6, "d2": 3.53728 }, { "df": 174.6, "d2": 3.59302 }, { "df": 182.4, "d2": 3.64505 }, { "df": 190, "d2": 3.69382 }, { "df": 197.3, "d2": 3.73969 }], [{ "df": 14.3, "d2": 1.14833 }, { "df": 29.3, "d2": 1.70708 }, { "df": 44.1, "d2": 2.07047 }, { "df": 58.2, "d2": 2.33594 }, { "df": 71.7, "d2": 2.54326 }, { "df": 84.5, "d2": 2.71237 }, { "df": 96.7, "d2": 2.85457 }, { "df": 108.4, "d2": 2.97689 }, { "df": 119.5, "d2": 3.08395 }, { "df": 130.2, "d2": 3.17897 }, { "df": 140.4, "d2": 3.26427 }, { "df": 150.2, "d2": 3.34154 }, { "df": 159.7, "d2": 3.4121 }, { "df": 168.9, "d2": 3.47707 }, { "df": 177.7, "d2": 3.53695 }, { "df": 186.3, "d2": 3.5927 }, { "df": 194.6, "d2": 3.64474 }, { "df": 202.6, "d2": 3.69351 }, { "df": 210.4, "d2": 3.73939 }], [{ "df": 15.1, "d2": 1.14717 }, { "df": 31.1, "d2": 1.70623 }, { "df": 46.8, "d2": 2.06978 }, { "df": 61.8, "d2": 2.33535 }, { "df": 76.2, "d2": 2.54274 }, { "df": 89.8, "d2": 2.7119 }, { "df": 102.8, "d2": 2.85413 }, { "df": 115.1, "d2": 2.97649 }, { "df": 127, "d2": 3.08358 }, { "df": 138.3, "d2": 3.17861 }, { "df": 149.2, "d2": 3.26393 }, { "df": 159.6, "d2": 3.34121 }, { "df": 169.7, "d2": 3.41178 }, { "df": 179.4, "d2": 3.47677 }, { "df": 188.8, "d2": 3.53666 }, { "df": 197.9, "d2": 3.59242 }, { "df": 206.7, "d2": 3.64447 }, { "df": 215.2, "d2": 3.69325 }, { "df": 223.6, "d2": 3.73913 }], [{ "df": 16, "d2": 1.14613 }, { "df": 32.9, "d2": 1.70547 }, { "df": 49.5, "d2": 2.06917 }, { "df": 65.5, "d2": 2.33483 }, { "df": 80.6, "d2": 2.54228 }, { "df": 95.1, "d2": 2.71148 }, { "df": 108.8, "d2": 2.85375 }, { "df": 121.9, "d2": 2.97613 }, { "df": 134.4, "d2": 3.08324 }, { "df": 146.4, "d2": 3.17829 }, { "df": 157.9, "d2": 3.26362 }, { "df": 169, "d2": 3.34092 }, { "df": 179.7, "d2": 3.4115 }, { "df": 190, "d2": 3.4765 }, { "df": 199.9, "d2": 3.5364 }, { "df": 209.5, "d2": 3.59216 }, { "df": 218.8, "d2": 3.64422 }, { "df": 227.9, "d2": 3.69301 }, { "df": 236.7, "d2": 3.7389 }], [{ "df": 16.9, "d2": 1.1452 }, { "df": 34.7, "d2": 1.7048 }, { "df": 52.3, "d2": 2.06862 }, { "df": 69.1, "d2": 2.33436 }, { "df": 85.1, "d2": 2.54187 }, { "df": 100.3, "d2": 2.71111 }, { "df": 114.8, "d2": 2.85341 }, { "df": 128.7, "d2": 2.97581 }, { "df": 141.9, "d2": 3.08294 }, { "df": 154.5, "d2": 3.17801 }, { "df": 166.7, "d2": 3.26335 }, { "df": 178.4, "d2": 3.34066 }, { "df": 189.6, "d2": 3.41125 }, { "df": 200.5, "d2": 3.47626 }, { "df": 211, "d2": 3.53617 }, { "df": 221.1, "d2": 3.59194 }, { "df": 231, "d2": 3.644 }, { "df": 240.5, "d2": 3.6928 }, { "df": 249.8, "d2": 3.73869 }], [{ "df": 17.8, "d2": 1.14437 }, { "df": 36.5, "d2": 1.70419 }, { "df": 55, "d2": 2.06813 }, { "df": 72.7, "d2": 2.33394 }, { "df": 89.6, "d2": 2.54149 }, { "df": 105.6, "d2": 2.71077 }, { "df": 120.9, "d2": 2.8531 }, { "df": 135.4, "d2": 2.97552 }, { "df": 149.3, "d2": 3.08267 }, { "df": 162.7, "d2": 3.17775 }, { "df": 175.5, "d2": 3.26311 }, { "df": 187.8, "d2": 3.34042 }, { "df": 199.6, "d2": 3.41103 }, { "df": 211, "d2": 3.47605 }, { "df": 222.1, "d2": 3.53596 }, { "df": 232.8, "d2": 3.59174 }, { "df": 243.1, "d2": 3.6438 }, { "df": 253.2, "d2": 3.6926 }, { "df": 263, "d2": 3.7385 }]];
     if (g > 20) {
       let d2s = [1.12838, 0.69257, 2.05875, 2.32593, 2.53441, 2.70436, 2.8472, 2.97003, 3.07751, 3.17287, 3.25845, 3.33598, 3.40676, 3.47193, 3.53198, 3.58788, 3.64006, 3.68896, 3.735];
       let cds = [0.876, 1.815, 2.7378, 3.623, 4.4658, 5.2673, 6.0305, 6.7582, 7.4539, 8.1207, 8.7602, 9.3751, 9.9679, 10.5396, 11.0913, 11.6259, 12.144, 12.6468, 13.1362];
@@ -10170,6 +10328,8 @@ function QQPlotChart(data, distribution = "Normal") {
   ];
   chartSettings.axis.y.customGridlines = customGridlines;
   chartSettings.axis.y.showAutoGridlines = false;
+  chartSettings.width += 100;
+  chartSettings.margins.right += 100;
   if (distribution.toLowerCase().includes("exponential")) {
     chartSettings.axis.x.max = 100;
     chartSettings.axis.x.min = 0.1;
@@ -10255,7 +10415,7 @@ function QQPlotChart(data, distribution = "Normal") {
         { key: "SW(R):", value: SWR.W },
         { key: "p-value:", value: SWR.p },
         { key: "", value: "" }
-      ], 1075, 50, 175, 100, "", 30);
+      ], chartSettings.width - chartSettings.margins.right + 50, 50, chartSettings.margins.right - 50, 100, "", 30, new EngFont("black", "Calibri", "28px"));
     }
   } else {
     let dist = distribution.replace("-", "").replace("arameter", "").replace(/\s/g, "");
@@ -10265,7 +10425,7 @@ function QQPlotChart(data, distribution = "Normal") {
     let AD = GoodnessOfFit.AndersonDarling.Test(data, dist);
     let ctx = chart.getContext("2d");
     if (ctx && (AD.AD || AD.AD === Infinity || Number.isNaN(AD.AD)) && (AD.p || AD.p === Infinity || Number.isNaN(AD.p))) {
-      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }], 1075, 150, 175, 100, "", 30);
+      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }], chartSettings.width - chartSettings.margins.right + 50, 175, chartSettings.margins.right - 50, 100, "", 30, new EngFont("black", "Calibri", "28px"));
     }
   }
   return chart;
@@ -10338,7 +10498,7 @@ function FactorCode(level, factorLevels) {
   let centered = zeroBased - factorLevels.length / 2 + 0.5;
   return centered;
 }
-function GeneralizedLinearModel(dependent, factors, covariates = {}, interactions = []) {
+function GeneralizedLinearModel(dependent, factors, covariates = {}, interactions = [], simple) {
   let refactored = { ...covariates };
   let unfactored = { ...covariates, ...factors };
   let refactoredWLevels = { ...covariates };
@@ -11403,8 +11563,8 @@ function BiasVsReferencePlot(biasAndLinearity, parts, referenceValues, responses
   let biasInfo = [];
   if (biasAndLinearity.GageLinearity.Linearity) {
     modelInfo.push(["Linearity", RoundDigits(biasAndLinearity.GageLinearity.Linearity, 6), "%Linearity", RoundDigits(biasAndLinearity.GageLinearity.PercentLinearity * 100, 2) + "%"]);
-    biasInfo = [["Reference Value", "Bias", "%Bias"], ...biasAndLinearity.GageBias.map((d) => {
-      return [RoundDigits(d.ReferenceValue, 6, true), RoundDigits(d.Bias, 6), RoundDigits(d.PercentBias * 100, 2)];
+    biasInfo = [["Reference Value", "Bias", "%Bias", "P"], ...biasAndLinearity.GageBias.map((d) => {
+      return [RoundDigits(d.ReferenceValue, 6, true), RoundDigits(d.Bias, 6), RoundDigits(d.PercentBias * 100, 2), RoundDigits(d.p, 3)];
     })];
     let miniX = chartSettings.width - chartSettings.margins.right + 75;
     let miniW = chartSettings.margins.right - 125;
@@ -11639,24 +11799,23 @@ function GageLinearityAndBias(parts, referenceValues, responses, processVariatio
     let bias = Sum(biases) / n2;
     let avgMeasure = Mean(measurements);
     let t = NaN;
+    let { d2, df } = UnbiasingConstant.D2prime(n2, Object.keys(byParts).length);
     if (useStDevMethod) {
       let \u03C3 = Math.sqrt(SumSq(measurements.map((d) => {
         return d - avgMeasure;
       })) / (n2 - 1));
       t = bias / (\u03C3 / Math.sqrt(n2));
     } else {
-      let \u03C3 = (Math.max(...biases) - Math.min(...biases)) / UnbiasingConstant.D2(n2);
-      console.log({ \u03C3, d2: UnbiasingConstant.D2(n2) });
+      let \u03C3 = (Math.max(...biases) - Math.min(...biases)) / d2;
       t = bias / (\u03C3 / Math.sqrt(n2));
     }
-    let df = n2 - 1;
-    let pValue = 2 * (1 - Distributions.T.cdf(Math.abs(t), df));
-    console.log({ t, pValue });
+    console.log(t, df);
+    let pValue = Distributions.T.cdf(-Math.abs(t), df) + Distributions.T.RightTail(Math.abs(t), df);
     partBiases.push({
       ReferenceValue: ref,
       Bias: bias,
-      PercentBias: processVariation ? Math.abs(bias / processVariation) : NaN
-      // p: pValue
+      PercentBias: processVariation ? Math.abs(bias / processVariation) : NaN,
+      p: pValue
     });
     biasAsY.push(bias);
     refAsX.push(ref);
@@ -11664,12 +11823,13 @@ function GageLinearityAndBias(parts, referenceValues, responses, processVariatio
   let avgBias = Mean(partBiases.map((b) => {
     return b.Bias;
   }));
-  let avgT = avgBias / (Math.max(...biasAsY) - Math.min(...biasAsY)) / UnbiasingConstant.D2(responses.length) / Math.sqrt(responses.length);
+  let { d2: avgD2, df: avgDF } = UnbiasingConstant.D2prime(4, 5);
+  let avgT = avgBias / (Math.max(...biasAsY) - Math.min(...biasAsY)) / avgD2 / Math.sqrt(responses.length);
   let avgGroup = {
     ReferenceValue: "Average",
     Bias: avgBias,
     PercentBias: processVariation ? Math.abs(avgBias / processVariation) : NaN,
-    p: Distributions.T.cdf(Math.abs(avgT), responses.length - 1)
+    p: Distributions.T.cdf(Math.abs(avgT), avgDF)
   };
   let gageBiases = [avgGroup, ...partBiases];
   let regression = SimpleLinearRegression(biasAsY, refAsX);
@@ -12689,6 +12849,7 @@ var Matrix = {
   SumSquaredDeviation,
   TrimmedMean,
   UnbiasingConstant,
+  Unique,
   Variance,
   VarianceComponents
 });

package/dist/index.mjs

@@ -5462,7 +5462,7 @@ var ANOVA = {
           let rp = ANOVA.PostHoc.Duncan.Rp(treatmentNumber, DFE);
           let Rp = rp * S_x;
           returnArray.push({
-            treatment: infoArray[i].treatment,
+            Treatment: infoArray[i].treatment,
             Mean: infoArray[i].Mean,
             Rp,
             D: infoArray[i].Mean - Rp,
@@ -5587,8 +5587,8 @@ var ANOVA = {
             let diff = infoArray[i].Mean - infoArray[rev].Mean;
             let t = diff / Math.sqrt(MSE * (1 / infoArray[i].N + 1 / infoArray[rev].N));
             returnArray.push({
-              name: `${infoArray[i].treatment} - ${infoArray[rev].treatment}`,
-              meanDifference: diff,
+              Name: `${infoArray[i].treatment} - ${infoArray[rev].treatment}`,
+              MeanDifference: diff,
               T: t,
               LSD: tCrit * Math.sqrt(MSE * (1 / infoArray[i].N + 1 / infoArray[rev].N)),
               p: 2 * Distributions.T.RightTail(t, DFE)
@@ -5614,7 +5614,7 @@ var ANOVA = {
           let mean = Mean(data[key]);
           let tCrit = q_prime * Math.sqrt(MSE / n);
           returnArray.push({
-            treatment: key,
+            Treatment: key,
             N: n,
             Mean: mean,
             T: tCrit,
@@ -6504,6 +6504,73 @@ var ANOVA = {
   }
 };
 var MANOVA = {
+  WilksLambda(responses, factors, interactions = []) {
+    let rKeys = Object.keys(responses);
+    let fKeys = Object.keys(factors);
+    let coefficientCount = Sum(fKeys.map((k) => {
+      return Unique(factors[k]).length;
+    })) - 1 + interactions.length;
+    let dfError = responses[rKeys[0]].length - coefficientCount - 1;
+    let varTerms = [...fKeys];
+    let p = rKeys.length;
+    let v = dfError;
+    let returner = [];
+    fKeys.map((fk) => {
+      let factor = factors[fk];
+      let q = Unique(factor).length - 1;
+      let E = MANOVA.ErrorMatrix(responses, factor);
+      let H = MANOVA.HypothesisMatrix(responses, factor);
+      let lambda = Matrix.Determinant(E) / Matrix.Determinant(Matrix.Sum(H, E));
+      let m = 0.5 * (Math.abs(p - 1) - 1);
+      let s = Math.min(p, 1);
+      let u = 0.25 * (p * q - 2);
+      let n = 0.5 * (v - p - 1);
+      let r = v - 0.5 * (p - q + 1);
+      let t = Math.pow(p, 2) + Math.pow(q, 2) - 5 > 0 ? Math.sqrt((Math.pow(p, 2) * Math.pow(q, 2) - 4) / (Math.pow(p, 2) + Math.pow(q, 2) - 5)) : 1;
+      let f = (1 - Math.pow(lambda, 1 / t)) / Math.pow(lambda, 1 / t) * (r * t - 2 * u) / (p * q);
+      let denDF = r * t - 2 * u;
+      let numDF = p * q;
+      returner.push({
+        DenominatorDF: denDF,
+        F: f,
+        Factor: fk,
+        Lambda: lambda,
+        NumeratorDF: numDF,
+        p: Distributions.F.RightTail(f, denDF, numDF)
+      });
+    });
+    interactions.forEach((int) => {
+      let factor = [];
+      for (let i = 0; i < responses[rKeys[0]].length; i++) {
+        factor.push(int.map((inter) => {
+          return factors[inter][i];
+        }).join("*"));
+      }
+      console.log(factor);
+      let q = Unique(factor).length - 1;
+      let E = MANOVA.ErrorMatrix(responses, factor);
+      let H = MANOVA.HypothesisMatrix(responses, factor);
+      let lambda = Matrix.Determinant(E) / Matrix.Determinant(Matrix.Sum(H, E));
+      let m = 0.5 * (Math.abs(p - 1) - 1);
+      let s = Math.min(p, 1);
+      let u = 0.25 * (p * q - 2);
+      let n = 0.5 * (v - p - 1);
+      let r = v - 0.5 * (p - q + 1);
+      let t = Math.pow(p, 2) + Math.pow(q, 2) - 5 > 0 ? Math.sqrt((Math.pow(p, 2) * Math.pow(q, 2) - 4) / (Math.pow(p, 2) + Math.pow(q, 2) - 5)) : 1;
+      let f = (1 - Math.pow(lambda, 1 / t)) / Math.pow(lambda, 1 / t) * (r * t - 2 * u) / (p * q);
+      let denDF = r * t - 2 * u;
+      let numDF = p * q;
+      returner.push({
+        DenominatorDF: denDF,
+        F: f,
+        Factor: int.join("*"),
+        Lambda: lambda,
+        NumeratorDF: numDF,
+        p: Distributions.F.RightTail(f, denDF, numDF)
+      });
+    });
+    return returner;
+  },
   Report(responses, factors, interactions = []) {
     let returner = [];
     let rKeys = Object.keys(responses);
@@ -6523,7 +6590,6 @@ var MANOVA = {
       let u = 0.25 * (p * q - 2);
       let m = 0.5 * (Math.abs(p - q) - 1);
       let n = 0.5 * (v - p - 1);
-      console.log({ m, n, p, q, r, s, t, u, v });
       let TotalMatrix = Matrix.New(p, p);
       let ErrorMatrix = Matrix.New(p, p);
       Matrix.Fill(ErrorMatrix, 0);
@@ -6550,6 +6616,7 @@ var MANOVA = {
         let factorSSCP = Matrix.Multiply(Matrix.Transpose(rDevMatrix), rDevMatrix);
         ErrorMatrix = Matrix.Sum(ErrorMatrix, factorSSCP);
       });
+      console.log({ DeviationMatrix, nMatrix });
       let devXn = Matrix.Scale(DeviationMatrix, nMatrix);
       let HypothesisMatrix = Matrix.Multiply(Matrix.Transpose(DeviationMatrix), devXn);
       console.log({ HypothesisMatrix, ErrorMatrix });
@@ -6612,7 +6679,49 @@ var MANOVA = {
     });
     return returner;
   },
+  ErrorMatrix(responses, factor) {
+    let factorLevels = Unique(factor);
+    let errorMatrix = Matrix.New(Object.keys(responses).length, Object.keys(responses).length);
+    Matrix.Fill(errorMatrix, 0);
+    factorLevels.forEach((lv, lvi) => {
+      let nResponses = factor.filter((f) => f == lv).length;
+      let filteredResponses = {};
+      let responseMeans = {};
+      Object.keys(responses).forEach((k) => {
+        let subgroup = responses[k].filter((r, ri) => factor[ri] == lv);
+        filteredResponses[k] = subgroup;
+        responseMeans[k] = Mean(subgroup);
+      });
+      let treatmentMatrix = Matrix.New(nResponses, Object.keys(responses).length);
+      let deviationMatrix = Matrix.New(nResponses, Object.keys(responses).length);
+      for (let i = 0; i < nResponses; i++) {
+        Object.keys(responses).forEach((k, ki) => {
+          treatmentMatrix[i][ki] = filteredResponses[k][i];
+          deviationMatrix[i][ki] = filteredResponses[k][i] - responseMeans[k];
+        });
+      }
+      let SSCP = Matrix.Multiply(Matrix.Transpose(deviationMatrix), treatmentMatrix);
+      errorMatrix = Matrix.Sum(errorMatrix, SSCP);
+    });
+    return errorMatrix;
+  },
   HypothesisMatrix(responses, factor) {
+    let factorGroups = Unique(factor);
+    let responseMeans = {};
+    let deviationMatrix = Matrix.New(factorGroups.length, Object.keys(responses).length);
+    let nMatrix = Matrix.New(factorGroups.length, Object.keys(responses).length);
+    Object.keys(responses).forEach((k, ki) => {
+      let responseMean = Mean(responses[k]);
+      responseMeans[k] = responseMean;
+      factorGroups.forEach((g, gi) => {
+        let subgroup = responses[k].filter((d, di) => factor[di] == g);
+        deviationMatrix[gi][ki] = Mean(subgroup) - responseMean;
+        nMatrix[gi][ki] = subgroup.length;
+      });
+    });
+    let devXn = Matrix.Scale(deviationMatrix, nMatrix);
+    let hypothesisMatrix = Matrix.Multiply(Matrix.Transpose(deviationMatrix), devXn);
+    return hypothesisMatrix;
   }
 };
 function AverageMovingRange(data, w = 2) {
@@ -8087,6 +8196,54 @@ var OutlierTesting = {
       rejectNull: p <= alpha,
       chart: outlierPlot
     };
+  },
+  /**
+   * Performs Dixon's Test to test for and identify (if one exists) an outlier. H_0: data has no outliers | H_A: there is one outlier
+   * @param data The array of numbers to evaluate
+   * @param tails OPTIONAL: Whether to perform single or two-tailed Grubbs' Test. DEFAULT = 0: two-tailed | -1: lower tailed (checks if smallest value is an outlier) | 1: upper tailed (checks if largest value is an outlier)
+   * @param alpha OPTIONAL: The significance value to evaluate outliers. DEFAULT = 0.05
+  */
+  Dixon(data, i = 1, j = 0, alpha = 0.05, name) {
+    let y = [...data].sort((a, b) => a - b);
+    let n = y.length;
+    let r_ij = Math.max((y[i] - y[0]) / (y[n - 1 - j] - y[0]), (y[n - 1] - y[n - 1 - i]) / (y[n - 1] - y[j]));
+    let C = Probability.Factorialize(n) * Math.pow(2 * Math.PI, -3 / 2) / (Probability.Factorialize(j) * Probability.Factorialize(n - i - j - 2) * Probability.Factorialize(i - 1));
+    const j_xvr = (x, v, r) => {
+    };
+    const t = (r, v) => {
+      return (1 + Math.pow(r, 2)) * Math.pow(v, 2) / 2;
+    };
+    const u2 = (x) => {
+      return 3 * Math.pow(x, 2) / 2;
+    };
+    let p = NaN;
+    let outliers = [];
+    if (p <= alpha) {
+      if ((y[i] - y[0]) / (y[n - 1 - j] - y[0]) > (y[n - 1] - y[n - 1 - i]) / (y[n - 1] - y[j])) {
+        outliers = [{ i: data.findIndex((d) => d == y[0]), value: y[0] }];
+      } else {
+        outliers = [{ i: data.findIndex((d) => d == y[n - 1]), value: y[n - 1] }];
+      }
+    }
+    let plotData = data.map((d, i2) => {
+      if (outliers.findIndex((o) => o.i == i2) >= 0) {
+        return new Point(d, 1, new CanvasDrawSettings("red", "red", 1, 1));
+      } else {
+        return new Point(d, 1, new CanvasDrawSettings("#2183b2", "#2183b2"));
+      }
+    });
+    let plotSettings = new ChartSettings();
+    plotSettings.axis.y = new AxisSettings(0, 2);
+    plotSettings.axis.y.showAutoGridlines = false;
+    plotSettings.axis.y.customGridlines = [new CustomGridline(1)];
+    let title = name ? `Outlier Plot of ${name}` : `Outlier Plot`;
+    let outlierPlot = CreateScatterPlot([new DataSet(plotData)], title, plotSettings);
+    return {
+      test: i == 1 && j == 0 ? "Dixon's Q" : `Dixon's r${i}${j}`,
+      r: r_ij,
+      p,
+      chart: outlierPlot
+    };
   }
 };
 var Probability = {
@@ -8416,7 +8573,7 @@ var UnbiasingConstant = {
     return 3.4873 + 0.0250141 * N - 9823e-8 * Math.pow(N, 2);
   },
   D2prime(g, m) {
-    let d2PrimeTable = [];
+    let d2PrimeTable = [[{ "df": 1, "d2": 1.41421 }, { "df": 2, "d2": 1.91155 }, { "df": 2.9, "d2": 2.23887 }, { "df": 3.8, "d2": 2.48124 }, { "df": 4.7, "d2": 2.67253 }, { "df": 5.5, "d2": 2.82981 }, { "df": 6.3, "d2": 2.96288 }, { "df": 7, "d2": 3.07794 }, { "df": 7.7, "d2": 3.17905 }, { "df": 8.3, "d2": 3.26909 }, { "df": 9, "d2": 3.35016 }, { "df": 9.6, "d2": 3.42378 }, { "df": 10.2, "d2": 3.49116 }, { "df": 10.8, "d2": 3.55333 }, { "df": 11.3, "d2": 3.61071 }, { "df": 11.9, "d2": 3.66422 }, { "df": 12.4, "d2": 3.71424 }, { "df": 12.9, "d2": 3.76118 }, { "df": 13.4, "d2": 3.80537 }], [{ "df": 1.9, "d2": 1.27931 }, { "df": 3.8, "d2": 1.80538 }, { "df": 5.7, "d2": 2.15069 }, { "df": 7.5, "d2": 2.40484 }, { "df": 9.2, "d2": 2.60438 }, { "df": 10.8, "d2": 2.76779 }, { "df": 12.3, "d2": 2.90562 }, { "df": 13.8, "d2": 3.02446 }, { "df": 15.1, "d2": 3.12869 }, { "df": 16.5, "d2": 3.22134 }, { "df": 17.8, "d2": 3.30463 }, { "df": 19, "d2": 3.38017 }, { "df": 20.2, "d2": 3.44922 }, { "df": 21.3, "d2": 3.51287 }, { "df": 22.4, "d2": 3.57156 }, { "df": 23.5, "d2": 3.62625 }, { "df": 24.5, "d2": 3.67734 }, { "df": 25.5, "d2": 3.72524 }, { "df": 26.5, "d2": 3.77032 }], [{ "df": 2.8, "d2": 1.23105 }, { "df": 5.7, "d2": 1.76858 }, { "df": 8.4, "d2": 2.12049 }, { "df": 11.1, "d2": 2.37883 }, { "df": 13.6, "d2": 2.58127 }, { "df": 16, "d2": 2.74681 }, { "df": 18.3, "d2": 2.88628 }, { "df": 20.5, "d2": 3.00643 }, { "df": 22.6, "d2": 3.11173 }, { "df": 24.6, "d2": 3.20526 }, { "df": 26.5, "d2": 3.28931 }, { "df": 28.4, "d2": 3.3655 }, { "df": 30.1, "d2": 3.43512 }, { "df": 31.9, "d2": 3.49927 }, { "df": 33.5, "d2": 3.55842 }, { "df": 35.1, "d2": 3.61351 }, { "df": 36.7, "d2": 3.66495 }, { "df": 38.2, "d2": 3.71319 }, { "df": 39.7, "d2": 3.75857 }], [{ "df": 3.7, "d2": 1.20621 }, { "df": 7.5, "d2": 1.74989 }, { "df": 11.2, "d2": 2.10522 }, { "df": 14.7, "d2": 2.36571 }, { "df": 18.1, "d2": 2.56964 }, { "df": 21.3, "d2": 2.73626 }, { "df": 24.4, "d2": 2.87656 }, { "df": 27.3, "d2": 2.99737 }, { "df": 30.1, "d2": 3.10321 }, { "df": 32.7, "d2": 3.1972 }, { "df": 35.3, "d2": 3.28163 }, { "df": 37.7, "d2": 3.35815 }, { "df": 40.1, "d2": 3.42805 }, { "df": 42.4, "d2": 3.49246 }, { "df": 44.6, "d2": 3.55183 }, { "df": 46.7, "d2": 3.60712 }, { "df": 48.8, "d2": 3.65875 }, { "df": 50.8, "d2": 3.70715 }, { "df": 52.8, "d2": 3.75268 }], [{ "df": 4.6, "d2": 1.19105 }, { "df": 9.3, "d2": 1.73857 }, { "df": 13.9, "d2": 2.09601 }, { "df": 18.4, "d2": 2.35781 }, { "df": 22.6, "d2": 2.56263 }, { "df": 26.6, "d2": 2.72991 }, { "df": 30.4, "d2": 2.87071 }, { "df": 34, "d2": 2.99192 }, { "df": 37.5, "d2": 3.09808 }, { "df": 40.8, "d2": 3.19235 }, { "df": 44, "d2": 3.27701 }, { "df": 47.1, "d2": 3.35372 }, { "df": 50.1, "d2": 3.42381 }, { "df": 52.9, "d2": 3.48836 }, { "df": 55.7, "d2": 3.54787 }, { "df": 58.4, "d2": 3.60328 }, { "df": 61, "d2": 3.65502 }, { "df": 63.5, "d2": 3.70352 }, { "df": 65.9, "d2": 3.74914 }], [{ "df": 5.5, "d2": 1.18083 }, { "df": 11.1, "d2": 1.73099 }, { "df": 16.7, "d2": 2.08985 }, { "df": 22, "d2": 2.35253 }, { "df": 27, "d2": 2.55795 }, { "df": 31.8, "d2": 2.72567 }, { "df": 36.4, "d2": 2.8668 }, { "df": 40.8, "d2": 2.98829 }, { "df": 45, "d2": 3.09467 }, { "df": 49, "d2": 3.18911 }, { "df": 52.8, "d2": 3.27392 }, { "df": 56.5, "d2": 3.35077 }, { "df": 60.1, "d2": 3.42097 }, { "df": 63.5, "d2": 3.48563 }, { "df": 66.8, "d2": 3.54522 }, { "df": 70, "d2": 3.60072 }, { "df": 73.1, "d2": 3.65253 }, { "df": 76.1, "d2": 3.70109 }, { "df": 79.1, "d2": 3.74678 }], [{ "df": 6.4, "d2": 1.17348 }, { "df": 12.9, "d2": 1.72555 }, { "df": 19.4, "d2": 2.08543 }, { "df": 25.6, "d2": 2.34875 }, { "df": 31.5, "d2": 2.5546 }, { "df": 37.1, "d2": 2.72263 }, { "df": 42.5, "d2": 2.86401 }, { "df": 47.6, "d2": 2.98568 }, { "df": 52.4, "d2": 3.09222 }, { "df": 57.1, "d2": 3.18679 }, { "df": 61.6, "d2": 3.27172 }, { "df": 65.9, "d2": 3.34866 }, { "df": 70, "d2": 3.41894 }, { "df": 74, "d2": 3.48368 }, { "df": 77.9, "d2": 3.54333 }, { "df": 81.6, "d2": 3.59888 }, { "df": 85.3, "d2": 3.65075 }, { "df": 88.8, "d2": 3.69936 }, { "df": 92.2, "d2": 3.74509 }], [{ "df": 7.2, "d2": 1.16794 }, { "df": 14.8, "d2": 1.72147 }, { "df": 22.1, "d2": 2.08212 }, { "df": 29.2, "d2": 2.34591 }, { "df": 36, "d2": 2.55208 }, { "df": 42.4, "d2": 2.72036 }, { "df": 48.5, "d2": 2.86192 }, { "df": 54.3, "d2": 2.98373 }, { "df": 59.9, "d2": 3.09039 }, { "df": 65.2, "d2": 3.18506 }, { "df": 70.3, "d2": 3.27006 }, { "df": 75.2, "d2": 3.34708 }, { "df": 80, "d2": 3.41742 }, { "df": 84.6, "d2": 3.48221 }, { "df": 89, "d2": 3.54192 }, { "df": 93.3, "d2": 3.59751 }, { "df": 97.4, "d2": 3.64941 }, { "df": 101.4, "d2": 3.69806 }, { "df": 105.3, "d2": 3.74382 }], [{ "df": 8.1, "d2": 1.16361 }, { "df": 16.6, "d2": 1.71828 }, { "df": 24.9, "d2": 2.07953 }, { "df": 32.9, "d2": 2.3437 }, { "df": 40.4, "d2": 2.55013 }, { "df": 47.7, "d2": 2.71858 }, { "df": 54.5, "d2": 2.86028 }, { "df": 61.1, "d2": 2.98221 }, { "df": 67.3, "d2": 3.08896 }, { "df": 73.3, "d2": 3.1837 }, { "df": 79.1, "d2": 3.26878 }, { "df": 84.6, "d2": 3.34585 }, { "df": 90, "d2": 3.41624 }, { "df": 95.1, "d2": 3.48107 }, { "df": 100.1, "d2": 3.54081 }, { "df": 104.9, "d2": 3.59644 }, { "df": 109.5, "d2": 3.64838 }, { "df": 114.1, "d2": 3.69705 }, { "df": 118.5, "d2": 3.74284 }], [{ "df": 9, "d2": 1.16014 }, { "df": 18.4, "d2": 1.71573 }, { "df": 27.6, "d2": 2.07746 }, { "df": 36.5, "d2": 2.34192 }, { "df": 44.9, "d2": 2.54856 }, { "df": 52.9, "d2": 2.71717 }, { "df": 60.6, "d2": 2.85898 }, { "df": 67.8, "d2": 2.981 }, { "df": 74.8, "d2": 3.08781 }, { "df": 81.5, "d2": 3.18262 }, { "df": 87.9, "d2": 3.26775 }, { "df": 94, "d2": 3.34486 }, { "df": 99.9, "d2": 3.41529 }, { "df": 105.6, "d2": 3.48016 }, { "df": 111.2, "d2": 3.53993 }, { "df": 116.5, "d2": 3.59559 }, { "df": 121.7, "d2": 3.64755 }, { "df": 126.7, "d2": 3.69625 }, { "df": 131.6, "d2": 3.74205 }], [{ "df": 9.9, "d2": 1.15729 }, { "df": 20.2, "d2": 1.71363 }, { "df": 30.4, "d2": 2.07577 }, { "df": 40.1, "d2": 2.34048 }, { "df": 49.4, "d2": 2.54728 }, { "df": 58.2, "d2": 2.716 }, { "df": 66.6, "d2": 2.85791 }, { "df": 74.6, "d2": 2.98 }, { "df": 82.2, "d2": 3.08688 }, { "df": 89.6, "d2": 3.18174 }, { "df": 96.6, "d2": 3.2669 }, { "df": 103.4, "d2": 3.34406 }, { "df": 109.9, "d2": 3.41452 }, { "df": 116.2, "d2": 3.47941 }, { "df": 122.3, "d2": 3.53921 }, { "df": 128.1, "d2": 3.59489 }, { "df": 133.8, "d2": 3.64687 }, { "df": 139.4, "d2": 3.69558 }, { "df": 144.7, "d2": 3.74141 }], [{ "df": 10.7, "d2": 1.1549 }, { "df": 22, "d2": 1.71189 }, { "df": 33.1, "d2": 2.07436 }, { "df": 43.7, "d2": 2.33927 }, { "df": 53.8, "d2": 2.54621 }, { "df": 63.5, "d2": 2.71504 }, { "df": 72.6, "d2": 2.85702 }, { "df": 81.3, "d2": 2.97917 }, { "df": 89.7, "d2": 3.0861 }, { "df": 97.7, "d2": 3.181 }, { "df": 105.4, "d2": 3.2662 }, { "df": 112.7, "d2": 3.34339 }, { "df": 119.9, "d2": 3.41387 }, { "df": 126.7, "d2": 3.47879 }, { "df": 133.3, "d2": 3.53861 }, { "df": 139.8, "d2": 3.5943 }, { "df": 146, "d2": 3.6463 }, { "df": 152, "d2": 3.69503 }, { "df": 157.9, "d2": 3.74087 }], [{ "df": 11.6, "d2": 1.15289 }, { "df": 23.8, "d2": 1.71041 }, { "df": 35.8, "d2": 2.07316 }, { "df": 47.3, "d2": 2.33824 }, { "df": 58.3, "d2": 2.5453 }, { "df": 68.7, "d2": 2.71422 }, { "df": 78.6, "d2": 2.85627 }, { "df": 88.1, "d2": 2.97847 }, { "df": 97.1, "d2": 3.08544 }, { "df": 105.8, "d2": 3.18037 }, { "df": 114.1, "d2": 3.26561 }, { "df": 122.1, "d2": 3.34282 }, { "df": 129.8, "d2": 3.41333 }, { "df": 137.3, "d2": 3.47826 }, { "df": 144.4, "d2": 3.5381 }, { "df": 151.4, "d2": 3.59381 }, { "df": 158.1, "d2": 3.64582 }, { "df": 164.7, "d2": 3.69457 }, { "df": 171, "d2": 3.74041 }], [{ "df": 12.5, "d2": 1.15115 }, { "df": 25.7, "d2": 1.70914 }, { "df": 38.6, "d2": 2.07213 }, { "df": 51, "d2": 2.33737 }, { "df": 62.8, "d2": 2.54452 }, { "df": 74, "d2": 2.71351 }, { "df": 84.7, "d2": 2.85562 }, { "df": 94.9, "d2": 2.97787 }, { "df": 104.6, "d2": 3.08487 }, { "df": 113.9, "d2": 3.17984 }, { "df": 122.9, "d2": 3.2651 }, { "df": 131.5, "d2": 3.34233 }, { "df": 139.8, "d2": 3.41286 }, { "df": 147.8, "d2": 3.47781 }, { "df": 155.5, "d2": 3.53766 }, { "df": 163, "d2": 3.59339 }, { "df": 170.3, "d2": 3.64541 }, { "df": 177.3, "d2": 3.69417 }, { "df": 184.2, "d2": 3.74002 }], [{ "df": 13.4, "d2": 1.14965 }, { "df": 27.5, "d2": 1.70804 }, { "df": 41.3, "d2": 2.07125 }, { "df": 54.6, "d2": 2.33661 }, { "df": 67.2, "d2": 2.54385 }, { "df": 79.3, "d2": 2.7129 }, { "df": 90.7, "d2": 2.85506 }, { "df": 101.6, "d2": 2.97735 }, { "df": 112.1, "d2": 3.08438 }, { "df": 122.1, "d2": 3.17938 }, { "df": 131.7, "d2": 3.26465 }, { "df": 140.9, "d2": 3.34191 }, { "df": 149.8, "d2": 3.41245 }, { "df": 158.3, "d2": 3.47742 }, { "df": 166.6, "d2": 3.53728 }, { "df": 174.6, "d2": 3.59302 }, { "df": 182.4, "d2": 3.64505 }, { "df": 190, "d2": 3.69382 }, { "df": 197.3, "d2": 3.73969 }], [{ "df": 14.3, "d2": 1.14833 }, { "df": 29.3, "d2": 1.70708 }, { "df": 44.1, "d2": 2.07047 }, { "df": 58.2, "d2": 2.33594 }, { "df": 71.7, "d2": 2.54326 }, { "df": 84.5, "d2": 2.71237 }, { "df": 96.7, "d2": 2.85457 }, { "df": 108.4, "d2": 2.97689 }, { "df": 119.5, "d2": 3.08395 }, { "df": 130.2, "d2": 3.17897 }, { "df": 140.4, "d2": 3.26427 }, { "df": 150.2, "d2": 3.34154 }, { "df": 159.7, "d2": 3.4121 }, { "df": 168.9, "d2": 3.47707 }, { "df": 177.7, "d2": 3.53695 }, { "df": 186.3, "d2": 3.5927 }, { "df": 194.6, "d2": 3.64474 }, { "df": 202.6, "d2": 3.69351 }, { "df": 210.4, "d2": 3.73939 }], [{ "df": 15.1, "d2": 1.14717 }, { "df": 31.1, "d2": 1.70623 }, { "df": 46.8, "d2": 2.06978 }, { "df": 61.8, "d2": 2.33535 }, { "df": 76.2, "d2": 2.54274 }, { "df": 89.8, "d2": 2.7119 }, { "df": 102.8, "d2": 2.85413 }, { "df": 115.1, "d2": 2.97649 }, { "df": 127, "d2": 3.08358 }, { "df": 138.3, "d2": 3.17861 }, { "df": 149.2, "d2": 3.26393 }, { "df": 159.6, "d2": 3.34121 }, { "df": 169.7, "d2": 3.41178 }, { "df": 179.4, "d2": 3.47677 }, { "df": 188.8, "d2": 3.53666 }, { "df": 197.9, "d2": 3.59242 }, { "df": 206.7, "d2": 3.64447 }, { "df": 215.2, "d2": 3.69325 }, { "df": 223.6, "d2": 3.73913 }], [{ "df": 16, "d2": 1.14613 }, { "df": 32.9, "d2": 1.70547 }, { "df": 49.5, "d2": 2.06917 }, { "df": 65.5, "d2": 2.33483 }, { "df": 80.6, "d2": 2.54228 }, { "df": 95.1, "d2": 2.71148 }, { "df": 108.8, "d2": 2.85375 }, { "df": 121.9, "d2": 2.97613 }, { "df": 134.4, "d2": 3.08324 }, { "df": 146.4, "d2": 3.17829 }, { "df": 157.9, "d2": 3.26362 }, { "df": 169, "d2": 3.34092 }, { "df": 179.7, "d2": 3.4115 }, { "df": 190, "d2": 3.4765 }, { "df": 199.9, "d2": 3.5364 }, { "df": 209.5, "d2": 3.59216 }, { "df": 218.8, "d2": 3.64422 }, { "df": 227.9, "d2": 3.69301 }, { "df": 236.7, "d2": 3.7389 }], [{ "df": 16.9, "d2": 1.1452 }, { "df": 34.7, "d2": 1.7048 }, { "df": 52.3, "d2": 2.06862 }, { "df": 69.1, "d2": 2.33436 }, { "df": 85.1, "d2": 2.54187 }, { "df": 100.3, "d2": 2.71111 }, { "df": 114.8, "d2": 2.85341 }, { "df": 128.7, "d2": 2.97581 }, { "df": 141.9, "d2": 3.08294 }, { "df": 154.5, "d2": 3.17801 }, { "df": 166.7, "d2": 3.26335 }, { "df": 178.4, "d2": 3.34066 }, { "df": 189.6, "d2": 3.41125 }, { "df": 200.5, "d2": 3.47626 }, { "df": 211, "d2": 3.53617 }, { "df": 221.1, "d2": 3.59194 }, { "df": 231, "d2": 3.644 }, { "df": 240.5, "d2": 3.6928 }, { "df": 249.8, "d2": 3.73869 }], [{ "df": 17.8, "d2": 1.14437 }, { "df": 36.5, "d2": 1.70419 }, { "df": 55, "d2": 2.06813 }, { "df": 72.7, "d2": 2.33394 }, { "df": 89.6, "d2": 2.54149 }, { "df": 105.6, "d2": 2.71077 }, { "df": 120.9, "d2": 2.8531 }, { "df": 135.4, "d2": 2.97552 }, { "df": 149.3, "d2": 3.08267 }, { "df": 162.7, "d2": 3.17775 }, { "df": 175.5, "d2": 3.26311 }, { "df": 187.8, "d2": 3.34042 }, { "df": 199.6, "d2": 3.41103 }, { "df": 211, "d2": 3.47605 }, { "df": 222.1, "d2": 3.53596 }, { "df": 232.8, "d2": 3.59174 }, { "df": 243.1, "d2": 3.6438 }, { "df": 253.2, "d2": 3.6926 }, { "df": 263, "d2": 3.7385 }]];
     if (g > 20) {
       let d2s = [1.12838, 0.69257, 2.05875, 2.32593, 2.53441, 2.70436, 2.8472, 2.97003, 3.07751, 3.17287, 3.25845, 3.33598, 3.40676, 3.47193, 3.53198, 3.58788, 3.64006, 3.68896, 3.735];
       let cds = [0.876, 1.815, 2.7378, 3.623, 4.4658, 5.2673, 6.0305, 6.7582, 7.4539, 8.1207, 8.7602, 9.3751, 9.9679, 10.5396, 11.0913, 11.6259, 12.144, 12.6468, 13.1362];
@@ -10046,6 +10203,8 @@ function QQPlotChart(data, distribution = "Normal") {
   ];
   chartSettings.axis.y.customGridlines = customGridlines;
   chartSettings.axis.y.showAutoGridlines = false;
+  chartSettings.width += 100;
+  chartSettings.margins.right += 100;
   if (distribution.toLowerCase().includes("exponential")) {
     chartSettings.axis.x.max = 100;
     chartSettings.axis.x.min = 0.1;
@@ -10131,7 +10290,7 @@ function QQPlotChart(data, distribution = "Normal") {
         { key: "SW(R):", value: SWR.W },
         { key: "p-value:", value: SWR.p },
         { key: "", value: "" }
-      ], 1075, 50, 175, 100, "", 30);
+      ], chartSettings.width - chartSettings.margins.right + 50, 50, chartSettings.margins.right - 50, 100, "", 30, new EngFont("black", "Calibri", "28px"));
     }
   } else {
     let dist = distribution.replace("-", "").replace("arameter", "").replace(/\s/g, "");
@@ -10141,7 +10300,7 @@ function QQPlotChart(data, distribution = "Normal") {
     let AD = GoodnessOfFit.AndersonDarling.Test(data, dist);
     let ctx = chart.getContext("2d");
     if (ctx && (AD.AD || AD.AD === Infinity || Number.isNaN(AD.AD)) && (AD.p || AD.p === Infinity || Number.isNaN(AD.p))) {
-      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }], 1075, 150, 175, 100, "", 30);
+      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }], chartSettings.width - chartSettings.margins.right + 50, 175, chartSettings.margins.right - 50, 100, "", 30, new EngFont("black", "Calibri", "28px"));
     }
   }
   return chart;
@@ -10214,7 +10373,7 @@ function FactorCode(level, factorLevels) {
   let centered = zeroBased - factorLevels.length / 2 + 0.5;
   return centered;
 }
-function GeneralizedLinearModel(dependent, factors, covariates = {}, interactions = []) {
+function GeneralizedLinearModel(dependent, factors, covariates = {}, interactions = [], simple) {
   let refactored = { ...covariates };
   let unfactored = { ...covariates, ...factors };
   let refactoredWLevels = { ...covariates };
@@ -11279,8 +11438,8 @@ function BiasVsReferencePlot(biasAndLinearity, parts, referenceValues, responses
   let biasInfo = [];
   if (biasAndLinearity.GageLinearity.Linearity) {
     modelInfo.push(["Linearity", RoundDigits(biasAndLinearity.GageLinearity.Linearity, 6), "%Linearity", RoundDigits(biasAndLinearity.GageLinearity.PercentLinearity * 100, 2) + "%"]);
-    biasInfo = [["Reference Value", "Bias", "%Bias"], ...biasAndLinearity.GageBias.map((d) => {
-      return [RoundDigits(d.ReferenceValue, 6, true), RoundDigits(d.Bias, 6), RoundDigits(d.PercentBias * 100, 2)];
+    biasInfo = [["Reference Value", "Bias", "%Bias", "P"], ...biasAndLinearity.GageBias.map((d) => {
+      return [RoundDigits(d.ReferenceValue, 6, true), RoundDigits(d.Bias, 6), RoundDigits(d.PercentBias * 100, 2), RoundDigits(d.p, 3)];
     })];
     let miniX = chartSettings.width - chartSettings.margins.right + 75;
     let miniW = chartSettings.margins.right - 125;
@@ -11515,24 +11674,23 @@ function GageLinearityAndBias(parts, referenceValues, responses, processVariatio
     let bias = Sum(biases) / n2;
     let avgMeasure = Mean(measurements);
     let t = NaN;
+    let { d2, df } = UnbiasingConstant.D2prime(n2, Object.keys(byParts).length);
     if (useStDevMethod) {
       let \u03C3 = Math.sqrt(SumSq(measurements.map((d) => {
         return d - avgMeasure;
       })) / (n2 - 1));
       t = bias / (\u03C3 / Math.sqrt(n2));
     } else {
-      let \u03C3 = (Math.max(...biases) - Math.min(...biases)) / UnbiasingConstant.D2(n2);
-      console.log({ \u03C3, d2: UnbiasingConstant.D2(n2) });
+      let \u03C3 = (Math.max(...biases) - Math.min(...biases)) / d2;
       t = bias / (\u03C3 / Math.sqrt(n2));
     }
-    let df = n2 - 1;
-    let pValue = 2 * (1 - Distributions.T.cdf(Math.abs(t), df));
-    console.log({ t, pValue });
+    console.log(t, df);
+    let pValue = Distributions.T.cdf(-Math.abs(t), df) + Distributions.T.RightTail(Math.abs(t), df);
     partBiases.push({
       ReferenceValue: ref,
       Bias: bias,
-      PercentBias: processVariation ? Math.abs(bias / processVariation) : NaN
-      // p: pValue
+      PercentBias: processVariation ? Math.abs(bias / processVariation) : NaN,
+      p: pValue
     });
     biasAsY.push(bias);
     refAsX.push(ref);
@@ -11540,12 +11698,13 @@ function GageLinearityAndBias(parts, referenceValues, responses, processVariatio
   let avgBias = Mean(partBiases.map((b) => {
     return b.Bias;
   }));
-  let avgT = avgBias / (Math.max(...biasAsY) - Math.min(...biasAsY)) / UnbiasingConstant.D2(responses.length) / Math.sqrt(responses.length);
+  let { d2: avgD2, df: avgDF } = UnbiasingConstant.D2prime(4, 5);
+  let avgT = avgBias / (Math.max(...biasAsY) - Math.min(...biasAsY)) / avgD2 / Math.sqrt(responses.length);
   let avgGroup = {
     ReferenceValue: "Average",
     Bias: avgBias,
     PercentBias: processVariation ? Math.abs(avgBias / processVariation) : NaN,
-    p: Distributions.T.cdf(Math.abs(avgT), responses.length - 1)
+    p: Distributions.T.cdf(Math.abs(avgT), avgDF)
   };
   let gageBiases = [avgGroup, ...partBiases];
   let regression = SimpleLinearRegression(biasAsY, refAsX);
@@ -12564,6 +12723,7 @@ export {
   SumSquaredDeviation,
   TrimmedMean,
   UnbiasingConstant,
+  Unique,
   Variance,
   VarianceComponents
 };

package/index.ts

@@ -1,4 +1,4 @@
-/** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ /** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  */
 // Classes
 
@@ -5858,7 +5858,7 @@ export const ANOVA = {
           let Rp = rp*S_x;
 
           returnArray.push({
-            treatment: infoArray[i].treatment,
+            Treatment: infoArray[i].treatment,
             Mean: infoArray[i].Mean,
             Rp: Rp,
             D: infoArray[i].Mean - Rp,
@@ -5983,8 +5983,8 @@ export const ANOVA = {
             let diff = infoArray[i].Mean - infoArray[rev].Mean;
             let t = diff/Math.sqrt(MSE*(1/infoArray[i].N + 1/infoArray[rev].N))
             returnArray.push({
-              name: `${infoArray[i].treatment} - ${infoArray[rev].treatment}`,
-              meanDifference: diff,
+              Name: `${infoArray[i].treatment} - ${infoArray[rev].treatment}`,
+              MeanDifference: diff,
               T: t,
               LSD: tCrit*Math.sqrt(MSE*(1/infoArray[i].N + 1/infoArray[rev].N)),
               p: 2*Distributions.T.RightTail(t, DFE)
@@ -6008,7 +6008,7 @@ export const ANOVA = {
           let mean = Mean(data[key]);
           let tCrit = q_prime*Math.sqrt(MSE/n)
           returnArray.push({
-            treatment: key,
+            Treatment: key,
             N: n,
             Mean: mean,
             T: tCrit,
@@ -6968,6 +6968,83 @@ export const ANOVA = {
 }
 
 export const MANOVA = {
+  WilksLambda(responses: any, factors: any, interactions: any[] = []){
+    let rKeys = Object.keys(responses);
+    let fKeys = Object.keys(factors);
+    let coefficientCount = Sum(fKeys.map(k => {return Unique(factors[k]).length})) - 1 + interactions.length;
+    let dfError = responses[rKeys[0]].length - coefficientCount - 1;
+
+    let varTerms = [...fKeys];
+
+    let p = rKeys.length;
+    let v = dfError;
+
+    
+    let returner: any = [];
+    fKeys.map((fk) => {
+      let factor = factors[fk]      
+      let q = Unique(factor).length - 1;
+      let E = MANOVA.ErrorMatrix(responses, factor);
+      let H = MANOVA.HypothesisMatrix(responses, factor);
+      let lambda = Matrix.Determinant(E) / Matrix.Determinant(Matrix.Sum(H, E))
+
+      let m  = 0.5 * (Math.abs(p - 1) - 1)
+      let s = Math.min(p, 1);
+      let u = 0.25*(p*q - 2)
+
+      let n = 0.5 * (v - p - 1)
+      let r = v - 0.5*(p - q + 1)
+      let t = (Math.pow(p, 2) + Math.pow(q, 2) - 5 > 0) ? Math.sqrt((Math.pow(p,2) * Math.pow(q,2) - 4)/(Math.pow(p,2) + Math.pow(q,2) - 5)) : 1
+
+      let f = (1 - Math.pow(lambda, 1/t)) / Math.pow(lambda, 1/t) * (r*t - 2*u) / (p * q)
+      let denDF = (r*t - 2*u)
+      let numDF = p*q
+
+      returner.push({
+        DenominatorDF: denDF,
+        F: f,
+        Factor: fk,
+        Lambda: lambda,
+        NumeratorDF: numDF,
+        p: Distributions.F.RightTail(f, denDF, numDF),
+      })
+    })
+
+    interactions.forEach(int => {
+      let factor: any = [];
+      for(let i = 0; i < responses[rKeys[0]].length; i++){
+        factor.push(int.map(inter => {return factors[inter][i]}).join("*"))
+      }
+      console.log(factor);
+  
+      let q = Unique(factor).length - 1;
+      let E = MANOVA.ErrorMatrix(responses, factor);
+      let H = MANOVA.HypothesisMatrix(responses, factor);
+      let lambda = Matrix.Determinant(E) / Matrix.Determinant(Matrix.Sum(H, E))
+      
+      let m  = 0.5 * (Math.abs(p - 1) - 1)
+      let s = Math.min(p, 1);
+      let u = 0.25*(p*q - 2)
+
+      let n = 0.5 * (v - p - 1)
+      let r = v - 0.5*(p - q + 1)
+      let t = (Math.pow(p, 2) + Math.pow(q, 2) - 5 > 0) ? Math.sqrt((Math.pow(p,2) * Math.pow(q,2) - 4)/(Math.pow(p,2) + Math.pow(q,2) - 5)) : 1
+
+      let f = (1 - Math.pow(lambda, 1/t)) / Math.pow(lambda, 1/t) * (r*t - 2*u) / (p * q)
+      let denDF = (r*t - 2*u)
+      let numDF = p*q
+
+      returner.push({
+        DenominatorDF: denDF,
+        F: f,
+        Factor: int.join("*"),
+        Lambda: lambda,
+        NumeratorDF: numDF,
+        p: Distributions.F.RightTail(f, denDF, numDF),
+      })
+    })
+    return returner
+  },
   Report(responses: any, factors: any, interactions: any[] = []){
     let returner: any = []
     let rKeys = Object.keys(responses);
@@ -6986,7 +7063,6 @@ export const MANOVA = {
       let u = 0.25 * (p * q - 2);
       let m = 0.5 * (Math.abs(p - q) - 1)
       let n = 0.5 * (v - p - 1)
-      console.log({m, n, p, q, r, s, t, u, v})
 
       let TotalMatrix = Matrix.New(p, p); //Total SSCP
       let ErrorMatrix: any = Matrix.New(p, p); //Within
@@ -7017,6 +7093,7 @@ export const MANOVA = {
         ErrorMatrix = Matrix.Sum(ErrorMatrix, factorSSCP); 
       })
 
+      console.log({DeviationMatrix, nMatrix})
       let devXn = Matrix.Scale(DeviationMatrix, nMatrix);
       let HypothesisMatrix = Matrix.Multiply(Matrix.Transpose(DeviationMatrix), devXn)//, responses[rKeys[0]].length/factorGroups.length)
       console.log({HypothesisMatrix, ErrorMatrix})
@@ -7098,8 +7175,50 @@ export const MANOVA = {
     // })
     return returner
   },
+  ErrorMatrix(responses: any, factor: any){
+    let factorLevels = Unique(factor);
+    let errorMatrix: any = Matrix.New(Object.keys(responses).length, Object.keys(responses).length);
+    Matrix.Fill(errorMatrix, 0);
+    factorLevels.forEach((lv, lvi) => {
+      let nResponses = factor.filter(f => f == lv).length
+      let filteredResponses = {}
+      let responseMeans = {}
+      Object.keys(responses).forEach(k => {
+        let subgroup = responses[k].filter((r, ri) => factor[ri] == lv);
+        filteredResponses[k] = subgroup;
+        responseMeans[k] = Mean(subgroup);
+      })
+      let treatmentMatrix = Matrix.New(nResponses, Object.keys(responses).length);
+      let deviationMatrix = Matrix.New(nResponses, Object.keys(responses).length);
+      for(let i = 0; i < nResponses; i++){
+        Object.keys(responses).forEach((k, ki) => {
+          treatmentMatrix[i][ki] = filteredResponses[k][i]
+          deviationMatrix[i][ki] = filteredResponses[k][i] - responseMeans[k]
+        })
+      }
+      let SSCP = Matrix.Multiply(Matrix.Transpose(deviationMatrix), treatmentMatrix);
+      errorMatrix = Matrix.Sum(errorMatrix, SSCP)
+    })
+    return errorMatrix
+  },
   HypothesisMatrix(responses: any, factor: any){
-    
+    let factorGroups = Unique(factor)
+    let responseMeans: any =  {}
+    let deviationMatrix: any = Matrix.New(factorGroups.length, Object.keys(responses).length)
+    let nMatrix: any  = Matrix.New(factorGroups.length, Object.keys(responses).length)
+    Object.keys(responses).forEach((k, ki) => {
+      let responseMean =  Mean(responses[k]);
+      responseMeans[k] = responseMean
+      factorGroups.forEach((g, gi) => {
+        let subgroup = responses[k].filter((d, di) => factor[di] == g);
+        deviationMatrix[gi][ki] = (Mean(subgroup) - responseMean)
+        nMatrix[gi][ki] = subgroup.length
+      })
+    })
+    let devXn = Matrix.Scale(deviationMatrix, nMatrix);
+    let hypothesisMatrix = Matrix.Multiply(Matrix.Transpose(deviationMatrix), devXn);
+
+    return hypothesisMatrix
   }
 }
 
@@ -8626,6 +8745,49 @@ export const OutlierTesting = {
       rejectNull: p <= alpha,
       chart: outlierPlot
     }
+  },
+   /**
+   * Performs Dixon's Test to test for and identify (if one exists) an outlier. H_0: data has no outliers | H_A: there is one outlier
+   * @param data The array of numbers to evaluate
+   * @param tails OPTIONAL: Whether to perform single or two-tailed Grubbs' Test. DEFAULT = 0: two-tailed | -1: lower tailed (checks if smallest value is an outlier) | 1: upper tailed (checks if largest value is an outlier)
+   * @param alpha OPTIONAL: The significance value to evaluate outliers. DEFAULT = 0.05
+  */
+  Dixon(data: any[], i: number = 1, j: number = 0, alpha: number = 0.05, name?: string){
+    let y = [...data].sort((a, b) => a-b);
+    let n = y.length;
+
+    let r_ij = Math.max((y[i] - y[0])/(y[n - 1 - j] - y[0]), (y[n-1] - y[n - 1 - i])/(y[n-1] - y[j]))
+
+    let C = (Probability.Factorialize(n)*Math.pow(2*Math.PI, -3/2))/(Probability.Factorialize(j) * Probability.Factorialize(n - i - j - 2) * Probability.Factorialize(i - 1))
+    const j_xvr = (x, v, r) => {}
+    const t = (r, v) => {return (1 + Math.pow(r, 2))*Math.pow(v, 2)/2} 
+    const u2 = (x) => {return 3*Math.pow(x, 2)/2}
+
+    let p = NaN;
+    let outliers:any = [];
+    if(p <= alpha){
+      // Add outliers
+      if((y[i] - y[0])/(y[n - 1 - j] - y[0]) > (y[n-1] - y[n - 1 - i])/(y[n-1] - y[j])){
+        outliers = [{i: data.findIndex(d => d == y[0]), value: y[0]}]
+      }else{
+        outliers = [{i: data.findIndex(d => d == y[n-1]), value: y[n-1]}]
+      }
+    }
+    
+    let plotData = data.map((d, i) => {if(outliers.findIndex(o => o.i == i) >= 0){return new Point(d, 1, new CanvasDrawSettings("red","red",1,1))}else{return new Point(d, 1, new CanvasDrawSettings('#2183b2', '#2183b2'))}})
+    let plotSettings = new ChartSettings();
+    plotSettings.axis.y = new AxisSettings(0,2);
+    plotSettings.axis.y.showAutoGridlines = false;
+    plotSettings.axis.y.customGridlines = [new CustomGridline(1)];
+    let title = name ? `Outlier Plot of ${name}` : `Outlier Plot`
+    let outlierPlot = CreateScatterPlot([new DataSet(plotData)], title, plotSettings)
+
+    return {
+      test: (i == 1 && j == 0) ? "Dixon's Q" : `Dixon's r${i}${j}`,
+      r: r_ij,
+      p: p,
+      chart: outlierPlot
+    }
   }
 }
 
@@ -9000,7 +9162,7 @@ export const UnbiasingConstant = {
     return 3.4873 + 0.0250141 * N - 0.00009823 * Math.pow(N,2);
   },
   D2prime(g: number, m: number){
-    let d2PrimeTable:any = [];
+    let d2PrimeTable:any = [[{"df":1,"d2":1.41421},{"df":2,"d2":1.91155},{"df":2.9,"d2":2.23887},{"df":3.8,"d2":2.48124},{"df":4.7,"d2":2.67253},{"df":5.5,"d2":2.82981},{"df":6.3,"d2":2.96288},{"df":7,"d2":3.07794},{"df":7.7,"d2":3.17905},{"df":8.3,"d2":3.26909},{"df":9,"d2":3.35016},{"df":9.6,"d2":3.42378},{"df":10.2,"d2":3.49116},{"df":10.8,"d2":3.55333},{"df":11.3,"d2":3.61071},{"df":11.9,"d2":3.66422},{"df":12.4,"d2":3.71424},{"df":12.9,"d2":3.76118},{"df":13.4,"d2":3.80537}],[{"df":1.9,"d2":1.27931},{"df":3.8,"d2":1.80538},{"df":5.7,"d2":2.15069},{"df":7.5,"d2":2.40484},{"df":9.2,"d2":2.60438},{"df":10.8,"d2":2.76779},{"df":12.3,"d2":2.90562},{"df":13.8,"d2":3.02446},{"df":15.1,"d2":3.12869},{"df":16.5,"d2":3.22134},{"df":17.8,"d2":3.30463},{"df":19,"d2":3.38017},{"df":20.2,"d2":3.44922},{"df":21.3,"d2":3.51287},{"df":22.4,"d2":3.57156},{"df":23.5,"d2":3.62625},{"df":24.5,"d2":3.67734},{"df":25.5,"d2":3.72524},{"df":26.5,"d2":3.77032}],[{"df":2.8,"d2":1.23105},{"df":5.7,"d2":1.76858},{"df":8.4,"d2":2.12049},{"df":11.1,"d2":2.37883},{"df":13.6,"d2":2.58127},{"df":16,"d2":2.74681},{"df":18.3,"d2":2.88628},{"df":20.5,"d2":3.00643},{"df":22.6,"d2":3.11173},{"df":24.6,"d2":3.20526},{"df":26.5,"d2":3.28931},{"df":28.4,"d2":3.3655},{"df":30.1,"d2":3.43512},{"df":31.9,"d2":3.49927},{"df":33.5,"d2":3.55842},{"df":35.1,"d2":3.61351},{"df":36.7,"d2":3.66495},{"df":38.2,"d2":3.71319},{"df":39.7,"d2":3.75857}],[{"df":3.7,"d2":1.20621},{"df":7.5,"d2":1.74989},{"df":11.2,"d2":2.10522},{"df":14.7,"d2":2.36571},{"df":18.1,"d2":2.56964},{"df":21.3,"d2":2.73626},{"df":24.4,"d2":2.87656},{"df":27.3,"d2":2.99737},{"df":30.1,"d2":3.10321},{"df":32.7,"d2":3.1972},{"df":35.3,"d2":3.28163},{"df":37.7,"d2":3.35815},{"df":40.1,"d2":3.42805},{"df":42.4,"d2":3.49246},{"df":44.6,"d2":3.55183},{"df":46.7,"d2":3.60712},{"df":48.8,"d2":3.65875},{"df":50.8,"d2":3.70715},{"df":52.8,"d2":3.75268}],[{"df":4.6,"d2":1.19105},{"df":9.3,"d2":1.73857},{"df":13.9,"d2":2.09601},{"df":18.4,"d2":2.35781},{"df":22.6,"d2":2.56263},{"df":26.6,"d2":2.72991},{"df":30.4,"d2":2.87071},{"df":34,"d2":2.99192},{"df":37.5,"d2":3.09808},{"df":40.8,"d2":3.19235},{"df":44,"d2":3.27701},{"df":47.1,"d2":3.35372},{"df":50.1,"d2":3.42381},{"df":52.9,"d2":3.48836},{"df":55.7,"d2":3.54787},{"df":58.4,"d2":3.60328},{"df":61,"d2":3.65502},{"df":63.5,"d2":3.70352},{"df":65.9,"d2":3.74914}],[{"df":5.5,"d2":1.18083},{"df":11.1,"d2":1.73099},{"df":16.7,"d2":2.08985},{"df":22,"d2":2.35253},{"df":27,"d2":2.55795},{"df":31.8,"d2":2.72567},{"df":36.4,"d2":2.8668},{"df":40.8,"d2":2.98829},{"df":45,"d2":3.09467},{"df":49,"d2":3.18911},{"df":52.8,"d2":3.27392},{"df":56.5,"d2":3.35077},{"df":60.1,"d2":3.42097},{"df":63.5,"d2":3.48563},{"df":66.8,"d2":3.54522},{"df":70,"d2":3.60072},{"df":73.1,"d2":3.65253},{"df":76.1,"d2":3.70109},{"df":79.1,"d2":3.74678}],[{"df":6.4,"d2":1.17348},{"df":12.9,"d2":1.72555},{"df":19.4,"d2":2.08543},{"df":25.6,"d2":2.34875},{"df":31.5,"d2":2.5546},{"df":37.1,"d2":2.72263},{"df":42.5,"d2":2.86401},{"df":47.6,"d2":2.98568},{"df":52.4,"d2":3.09222},{"df":57.1,"d2":3.18679},{"df":61.6,"d2":3.27172},{"df":65.9,"d2":3.34866},{"df":70,"d2":3.41894},{"df":74,"d2":3.48368},{"df":77.9,"d2":3.54333},{"df":81.6,"d2":3.59888},{"df":85.3,"d2":3.65075},{"df":88.8,"d2":3.69936},{"df":92.2,"d2":3.74509}],[{"df":7.2,"d2":1.16794},{"df":14.8,"d2":1.72147},{"df":22.1,"d2":2.08212},{"df":29.2,"d2":2.34591},{"df":36,"d2":2.55208},{"df":42.4,"d2":2.72036},{"df":48.5,"d2":2.86192},{"df":54.3,"d2":2.98373},{"df":59.9,"d2":3.09039},{"df":65.2,"d2":3.18506},{"df":70.3,"d2":3.27006},{"df":75.2,"d2":3.34708},{"df":80,"d2":3.41742},{"df":84.6,"d2":3.48221},{"df":89,"d2":3.54192},{"df":93.3,"d2":3.59751},{"df":97.4,"d2":3.64941},{"df":101.4,"d2":3.69806},{"df":105.3,"d2":3.74382}],[{"df":8.1,"d2":1.16361},{"df":16.6,"d2":1.71828},{"df":24.9,"d2":2.07953},{"df":32.9,"d2":2.3437},{"df":40.4,"d2":2.55013},{"df":47.7,"d2":2.71858},{"df":54.5,"d2":2.86028},{"df":61.1,"d2":2.98221},{"df":67.3,"d2":3.08896},{"df":73.3,"d2":3.1837},{"df":79.1,"d2":3.26878},{"df":84.6,"d2":3.34585},{"df":90,"d2":3.41624},{"df":95.1,"d2":3.48107},{"df":100.1,"d2":3.54081},{"df":104.9,"d2":3.59644},{"df":109.5,"d2":3.64838},{"df":114.1,"d2":3.69705},{"df":118.5,"d2":3.74284}],[{"df":9,"d2":1.16014},{"df":18.4,"d2":1.71573},{"df":27.6,"d2":2.07746},{"df":36.5,"d2":2.34192},{"df":44.9,"d2":2.54856},{"df":52.9,"d2":2.71717},{"df":60.6,"d2":2.85898},{"df":67.8,"d2":2.981},{"df":74.8,"d2":3.08781},{"df":81.5,"d2":3.18262},{"df":87.9,"d2":3.26775},{"df":94,"d2":3.34486},{"df":99.9,"d2":3.41529},{"df":105.6,"d2":3.48016},{"df":111.2,"d2":3.53993},{"df":116.5,"d2":3.59559},{"df":121.7,"d2":3.64755},{"df":126.7,"d2":3.69625},{"df":131.6,"d2":3.74205}],[{"df":9.9,"d2":1.15729},{"df":20.2,"d2":1.71363},{"df":30.4,"d2":2.07577},{"df":40.1,"d2":2.34048},{"df":49.4,"d2":2.54728},{"df":58.2,"d2":2.716},{"df":66.6,"d2":2.85791},{"df":74.6,"d2":2.98},{"df":82.2,"d2":3.08688},{"df":89.6,"d2":3.18174},{"df":96.6,"d2":3.2669},{"df":103.4,"d2":3.34406},{"df":109.9,"d2":3.41452},{"df":116.2,"d2":3.47941},{"df":122.3,"d2":3.53921},{"df":128.1,"d2":3.59489},{"df":133.8,"d2":3.64687},{"df":139.4,"d2":3.69558},{"df":144.7,"d2":3.74141}],[{"df":10.7,"d2":1.1549},{"df":22,"d2":1.71189},{"df":33.1,"d2":2.07436},{"df":43.7,"d2":2.33927},{"df":53.8,"d2":2.54621},{"df":63.5,"d2":2.71504},{"df":72.6,"d2":2.85702},{"df":81.3,"d2":2.97917},{"df":89.7,"d2":3.0861},{"df":97.7,"d2":3.181},{"df":105.4,"d2":3.2662},{"df":112.7,"d2":3.34339},{"df":119.9,"d2":3.41387},{"df":126.7,"d2":3.47879},{"df":133.3,"d2":3.53861},{"df":139.8,"d2":3.5943},{"df":146,"d2":3.6463},{"df":152,"d2":3.69503},{"df":157.9,"d2":3.74087}],[{"df":11.6,"d2":1.15289},{"df":23.8,"d2":1.71041},{"df":35.8,"d2":2.07316},{"df":47.3,"d2":2.33824},{"df":58.3,"d2":2.5453},{"df":68.7,"d2":2.71422},{"df":78.6,"d2":2.85627},{"df":88.1,"d2":2.97847},{"df":97.1,"d2":3.08544},{"df":105.8,"d2":3.18037},{"df":114.1,"d2":3.26561},{"df":122.1,"d2":3.34282},{"df":129.8,"d2":3.41333},{"df":137.3,"d2":3.47826},{"df":144.4,"d2":3.5381},{"df":151.4,"d2":3.59381},{"df":158.1,"d2":3.64582},{"df":164.7,"d2":3.69457},{"df":171,"d2":3.74041}],[{"df":12.5,"d2":1.15115},{"df":25.7,"d2":1.70914},{"df":38.6,"d2":2.07213},{"df":51,"d2":2.33737},{"df":62.8,"d2":2.54452},{"df":74,"d2":2.71351},{"df":84.7,"d2":2.85562},{"df":94.9,"d2":2.97787},{"df":104.6,"d2":3.08487},{"df":113.9,"d2":3.17984},{"df":122.9,"d2":3.2651},{"df":131.5,"d2":3.34233},{"df":139.8,"d2":3.41286},{"df":147.8,"d2":3.47781},{"df":155.5,"d2":3.53766},{"df":163,"d2":3.59339},{"df":170.3,"d2":3.64541},{"df":177.3,"d2":3.69417},{"df":184.2,"d2":3.74002}],[{"df":13.4,"d2":1.14965},{"df":27.5,"d2":1.70804},{"df":41.3,"d2":2.07125},{"df":54.6,"d2":2.33661},{"df":67.2,"d2":2.54385},{"df":79.3,"d2":2.7129},{"df":90.7,"d2":2.85506},{"df":101.6,"d2":2.97735},{"df":112.1,"d2":3.08438},{"df":122.1,"d2":3.17938},{"df":131.7,"d2":3.26465},{"df":140.9,"d2":3.34191},{"df":149.8,"d2":3.41245},{"df":158.3,"d2":3.47742},{"df":166.6,"d2":3.53728},{"df":174.6,"d2":3.59302},{"df":182.4,"d2":3.64505},{"df":190,"d2":3.69382},{"df":197.3,"d2":3.73969}],[{"df":14.3,"d2":1.14833},{"df":29.3,"d2":1.70708},{"df":44.1,"d2":2.07047},{"df":58.2,"d2":2.33594},{"df":71.7,"d2":2.54326},{"df":84.5,"d2":2.71237},{"df":96.7,"d2":2.85457},{"df":108.4,"d2":2.97689},{"df":119.5,"d2":3.08395},{"df":130.2,"d2":3.17897},{"df":140.4,"d2":3.26427},{"df":150.2,"d2":3.34154},{"df":159.7,"d2":3.4121},{"df":168.9,"d2":3.47707},{"df":177.7,"d2":3.53695},{"df":186.3,"d2":3.5927},{"df":194.6,"d2":3.64474},{"df":202.6,"d2":3.69351},{"df":210.4,"d2":3.73939}],[{"df":15.1,"d2":1.14717},{"df":31.1,"d2":1.70623},{"df":46.8,"d2":2.06978},{"df":61.8,"d2":2.33535},{"df":76.2,"d2":2.54274},{"df":89.8,"d2":2.7119},{"df":102.8,"d2":2.85413},{"df":115.1,"d2":2.97649},{"df":127,"d2":3.08358},{"df":138.3,"d2":3.17861},{"df":149.2,"d2":3.26393},{"df":159.6,"d2":3.34121},{"df":169.7,"d2":3.41178},{"df":179.4,"d2":3.47677},{"df":188.8,"d2":3.53666},{"df":197.9,"d2":3.59242},{"df":206.7,"d2":3.64447},{"df":215.2,"d2":3.69325},{"df":223.6,"d2":3.73913}],[{"df":16,"d2":1.14613},{"df":32.9,"d2":1.70547},{"df":49.5,"d2":2.06917},{"df":65.5,"d2":2.33483},{"df":80.6,"d2":2.54228},{"df":95.1,"d2":2.71148},{"df":108.8,"d2":2.85375},{"df":121.9,"d2":2.97613},{"df":134.4,"d2":3.08324},{"df":146.4,"d2":3.17829},{"df":157.9,"d2":3.26362},{"df":169,"d2":3.34092},{"df":179.7,"d2":3.4115},{"df":190,"d2":3.4765},{"df":199.9,"d2":3.5364},{"df":209.5,"d2":3.59216},{"df":218.8,"d2":3.64422},{"df":227.9,"d2":3.69301},{"df":236.7,"d2":3.7389}],[{"df":16.9,"d2":1.1452},{"df":34.7,"d2":1.7048},{"df":52.3,"d2":2.06862},{"df":69.1,"d2":2.33436},{"df":85.1,"d2":2.54187},{"df":100.3,"d2":2.71111},{"df":114.8,"d2":2.85341},{"df":128.7,"d2":2.97581},{"df":141.9,"d2":3.08294},{"df":154.5,"d2":3.17801},{"df":166.7,"d2":3.26335},{"df":178.4,"d2":3.34066},{"df":189.6,"d2":3.41125},{"df":200.5,"d2":3.47626},{"df":211,"d2":3.53617},{"df":221.1,"d2":3.59194},{"df":231,"d2":3.644},{"df":240.5,"d2":3.6928},{"df":249.8,"d2":3.73869}],[{"df":17.8,"d2":1.14437},{"df":36.5,"d2":1.70419},{"df":55,"d2":2.06813},{"df":72.7,"d2":2.33394},{"df":89.6,"d2":2.54149},{"df":105.6,"d2":2.71077},{"df":120.9,"d2":2.8531},{"df":135.4,"d2":2.97552},{"df":149.3,"d2":3.08267},{"df":162.7,"d2":3.17775},{"df":175.5,"d2":3.26311},{"df":187.8,"d2":3.34042},{"df":199.6,"d2":3.41103},{"df":211,"d2":3.47605},{"df":222.1,"d2":3.53596},{"df":232.8,"d2":3.59174},{"df":243.1,"d2":3.6438},{"df":253.2,"d2":3.6926},{"df":263,"d2":3.7385}]]
     if(g > 20){
       // https://dastmardi.ir/wp-content/uploads/2017/07/MSA_Reference_Manual_4th_Edition.pdf#page=212&zoom=100,0,97
       let d2s = [1.12838, .69257, 2.05875, 2.32593, 2.53441, 2.70436, 2.8472, 2.97003, 3.07751, 3.17287, 3.25845, 3.33598, 3.40676, 3.47193, 3.53198, 3.58788, 3.64006, 3.68896, 3.735]
@@ -10934,6 +11096,8 @@ export function QQPlotChart(data: number[], distribution: string = 'Normal'){
 
   chartSettings.axis.y.customGridlines = customGridlines;
   chartSettings.axis.y.showAutoGridlines = false;
+  chartSettings.width += 100;
+  chartSettings.margins.right += 100;
 
   if(distribution.toLowerCase().includes('exponential')){
       chartSettings.axis.x.max = 100;
@@ -11018,7 +11182,7 @@ export function QQPlotChart(data: number[], distribution: string = 'Normal'){
         { 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, 50, 175, 100, "", 30);
+        { key: "SW(R):", value: SWR.W }, { key: "p-value:", value: SWR.p }, { key: "", value: "" },], chartSettings.width - chartSettings.margins.right + 50, 50, chartSettings.margins.right - 50, 100, "", 30, new EngFont("black","Calibri","28px"));
     }
   }else{
     let dist = distribution.replace("-", "").replace("arameter", "").replace(/\s/g, "");
@@ -11031,7 +11195,7 @@ export function QQPlotChart(data: number[], distribution: string = 'Normal'){
     let ctx: any = chart.getContext("2d");
     
     if(ctx && (AD.AD || AD.AD === Infinity || Number.isNaN(AD.AD)) && (AD.p || AD.p === Infinity || Number.isNaN(AD.p))){
-      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }], 1075, 150, 175, 100, "", 30);
+      AddInfoTable(ctx, [{ key: "AD:", value: AD.AD }, { key: "p-value:", value: AD.p }],  chartSettings.width - chartSettings.margins.right + 50, 175, chartSettings.margins.right - 50, 100, "", 30, new EngFont("black","Calibri","28px"));
     }
 
   }
@@ -11129,8 +11293,7 @@ function FactorCode(level: any, factorLevels: string[]){
   return centered
 }
 
-
-export function GeneralizedLinearModel(dependent: number[], factors: any, covariates: any = {}, interactions: any[][] = []){
+export function GeneralizedLinearModel(dependent: number[], factors: any, covariates: any = {}, interactions: any[][] = [], simple?: boolean){
   let refactored = {...covariates};
   let unfactored = {...covariates, ...factors}
   let refactoredWLevels = {...covariates}
@@ -11258,7 +11421,6 @@ export function MultipleLinearRegression(dependent: number[], factors: any, inte
       row[fi] = fObj[f][ri]
     })
   })
-  
   let inv_xTx = Matrix.Inverse(Matrix.Multiply(Matrix.Transpose(x), x))
   let xTy = Matrix.Multiply(Matrix.Transpose(x), y);
   let betaEstimators = Matrix.Multiply(inv_xTx, xTy);
@@ -12225,7 +12387,7 @@ function BiasVsReferencePlot(biasAndLinearity: any, parts: any[], referenceValue
   let biasInfo: any = []
   if(biasAndLinearity.GageLinearity.Linearity){
     modelInfo.push(["Linearity", RoundDigits(biasAndLinearity.GageLinearity.Linearity, 6),"%Linearity",RoundDigits(biasAndLinearity.GageLinearity.PercentLinearity*100, 2)+"%"])
-    biasInfo = [["Reference Value", "Bias", "%Bias"], ...biasAndLinearity.GageBias.map(d => {return [RoundDigits(d.ReferenceValue, 6, true), RoundDigits(d.Bias, 6), RoundDigits(d.PercentBias*100, 2)]})]
+    biasInfo = [["Reference Value", "Bias", "%Bias", "P"], ...biasAndLinearity.GageBias.map(d => {return [RoundDigits(d.ReferenceValue, 6, true), RoundDigits(d.Bias, 6), RoundDigits(d.PercentBias*100, 2), RoundDigits(d.p, 3)]})]
   
     // Draw %Process Variation Chart
     let miniX = chartSettings.width - chartSettings.margins.right + 75
@@ -12492,24 +12654,22 @@ export function GageLinearityAndBias(parts: any[], referenceValues: number[], re
     let avgMeasure = Mean(measurements);
 
     let t = NaN;
+    let {d2: d2, df: df} = UnbiasingConstant.D2prime(n, Object.keys(byParts).length)
     if(useStDevMethod){
       let σ = Math.sqrt(SumSq(measurements.map(d => {return d - avgMeasure}))/(n - 1))
       t = bias / (σ/Math.sqrt(n))
     }else{
-      let σ = (Math.max(...biases) - Math.min(...biases)) / UnbiasingConstant.D2(n);
-      console.log({σ, d2: UnbiasingConstant.D2(n)})
+      let σ = (Math.max(...biases) - Math.min(...biases)) / d2;
       t = bias / (σ/Math.sqrt(n))
     }
-
-    let df = n - 1;
-    let pValue = 2 * (1 - Distributions.T.cdf(Math.abs(t), df))
-    console.log({t, pValue})
+    console.log(t, df)
+    let pValue = Distributions.T.cdf(-Math.abs(t), df) + Distributions.T.RightTail(Math.abs(t), df)
 
     partBiases.push({
       ReferenceValue: ref,
       Bias: bias,
       PercentBias: processVariation ? Math.abs(bias/processVariation) : NaN,
-      // p: pValue
+      p: pValue
     })
 
     biasAsY.push(bias)
@@ -12517,12 +12677,13 @@ export function GageLinearityAndBias(parts: any[], referenceValues: number[], re
   })
 
   let avgBias = Mean(partBiases.map(b => {return b.Bias}))
-  let avgT = avgBias/(Math.max(...biasAsY) - Math.min(...biasAsY))/(UnbiasingConstant.D2(responses.length))/Math.sqrt(responses.length)
+  let {d2: avgD2, df: avgDF} = UnbiasingConstant.D2prime(4, 5)
+  let avgT = avgBias/(Math.max(...biasAsY) - Math.min(...biasAsY))/(avgD2)/Math.sqrt(responses.length)
   let avgGroup = {
     ReferenceValue: "Average",
     Bias: avgBias,
     PercentBias: processVariation ? Math.abs(avgBias/processVariation) : NaN,
-    p: Distributions.T.cdf(Math.abs(avgT), responses.length - 1)
+    p: Distributions.T.cdf(Math.abs(avgT), avgDF)
   }
   let gageBiases = [avgGroup, ...partBiases]
 
@@ -13331,7 +13492,7 @@ function AllIdentical(values: any[]){
   }
   return identical
 }
-function Unique(values: any[]){
+export function Unique(values: any[]){
   let unique: any[] = [];
   values.forEach(v => {
     if(unique.indexOf(v) < 0){
@@ -13545,4 +13706,5 @@ const Matrix = {
   },
 }
 // End Matrix Functions
-/** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ */
+/** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ *///  Validation
+ 

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "qesuite",
-  "version": "1.0.58",
+  "version": "1.0.59",
   "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_59.md

@@ -0,0 +1,3 @@
+- Changed Font Size for normality test results in QQPlotChart()
+- Adjusted Output Property names of Post Hoc tests to align with Capitalization Format
+- Added an Export for the Unique Function