qesuite 1.0.49 → 1.0.51

package/dist/index.d.mts

@@ -1019,6 +1019,9 @@ declare const ANOVA: {
             p: number;
         };
     };
+    Multifactor: {
+        MainEffectsPlot(factors: any, response: any): HTMLCanvasElement;
+    };
 };
 declare function SumSquaredDeviation(data: number[], referenceValue?: any): number;
 declare const Probability: {

package/dist/index.d.ts

@@ -1019,6 +1019,9 @@ declare const ANOVA: {
             p: number;
         };
     };
+    Multifactor: {
+        MainEffectsPlot(factors: any, response: any): HTMLCanvasElement;
+    };
 };
 declare function SumSquaredDeviation(data: number[], referenceValue?: any): number;
 declare const Probability: {

package/dist/index.js

@@ -2766,15 +2766,17 @@ var ANOVA = {
       let scaleDF = Math.sqrt((rCount - 1) / (rCount * tCount));
       let UDL = grandMean + pooledSTD * hAlpha * scaleDF;
       let LDL = grandMean - pooledSTD * hAlpha * scaleDF;
+      let chartSettings = new ChartSettings();
       returnArray.forEach((t, x) => {
         let canvasDraw = t.Mean < UDL && t.Mean > LDL ? new CanvasDrawSettings("#2183b2", "#2183b2", void 0, 0) : new CanvasDrawSettings("#F00", "#F00", void 0, 1);
         let newDataSet = new DataSet([
           new Point(x + 1, t.Mean, canvasDraw),
           new Point(x + 1, grandMean, new CanvasDrawSettings("#000", "#000", 0, 0, 0))
         ], t.treatment, void 0, new Line("straight", new CanvasDrawSettings("#000", "#000", 1)));
+        chartSettings.axis.x.customGridlines.push(new CustomGridline(x + 1, t.treatment, void 0, void 0, void 0, void 0, true));
         dataSets.push(newDataSet);
       });
-      let chartSettings = new ChartSettings();
+      chartSettings.axis.x.showAutoGridlines = false;
       chartSettings.axis.y.customGridlines = [new CustomGridline(grandMean, `x\u0305\u0305=${Math.round(grandMean * 1e3) / 1e3}`, new CanvasDrawSettings(GetThemeColor(2), GetThemeColor(2)), void 0, true), new CustomGridline(LDL, `LDL=${Math.round(LDL * 1e3) / 1e3}`, new CanvasDrawSettings("#F00", "#F00"), void 0, true, void 0, true), new CustomGridline(UDL, `UDL=${Math.round(UDL * 1e3) / 1e3}`, new CanvasDrawSettings("#F00", "#F00"), void 0, true, void 0, true)];
       let chart = CreateScatterPlot(dataSets, name ? `Analysis of Means for ${name}` : "Analysis of Means", chartSettings);
       return {
@@ -2852,8 +2854,8 @@ var ANOVA = {
           mean,
           stdev: StDev.S(data[key]),
           SE: standardError,
-          CI_Low: mean - tInv * standardError,
-          CI_High: mean + tInv * standardError
+          CI_Low: mean - Math.abs(tInv * standardError),
+          CI_High: mean + Math.abs(tInv * standardError)
         };
         summaryTable.ConfidenceIntervals.push(CI_table);
       }
@@ -3676,6 +3678,44 @@ var ANOVA = {
         p: Distributions.F.TwoTail(F, factor1.length - 1, factor2.length - 1)
       };
     }
+  },
+  Multifactor: {
+    MainEffectsPlot(factors, response) {
+      let responseName = Object.keys(response)[0];
+      let factorKeys = Object.keys(factors);
+      let individualCharts = [];
+      factorKeys.forEach((k) => {
+        let arr = factors[k];
+        let uniqueLevels = arr.filter((value, index, array) => array.indexOf(value) === index);
+        let dataSet = new DataSet([]);
+        dataSet.line = new Line("straight");
+        let chartSettings = new ChartSettings();
+        chartSettings.axis.x.showAutoGridlines = false;
+        chartSettings.axis.x.min = 0.5;
+        chartSettings.axis.x.max = uniqueLevels.length + 0.5;
+        chartSettings.axis.x.autoSize = false;
+        uniqueLevels.forEach((level, levelI) => {
+          let mean = 0;
+          let N = 0;
+          arr.forEach((v, i) => {
+            if (v === level) {
+              mean += response[responseName][i];
+              N += 1;
+            }
+          });
+          mean = mean / N;
+          dataSet.values.push(new Point(levelI + 1, mean, new CanvasDrawSettings("#2183b2", "#000", void 0, void 0, 5)));
+          chartSettings.axis.x.customGridlines.push(new CustomGridline(levelI + 1, level, void 0, void 0, void 0, true, false));
+        });
+        chartSettings.width += 1e3;
+        individualCharts.push(new Chart([dataSet], "scatter", k, chartSettings));
+      });
+      let chartSet = new ChartSettings();
+      chartSet.width += 1e3;
+      chartSet.axis.y.tickOnly = true;
+      let graph = CreateSplitGraph(individualCharts, `Main Effects Plot of ${Object.keys(response)[0]}`, chartSet);
+      return graph;
+    }
   }
 };
 function SumSquaredDeviation(data, referenceValue) {
@@ -3851,7 +3891,7 @@ function AddAxisGridlines(ctx, chartSettings, axis) {
       }
       let x = ConvertToCanvasPt(unadjustedValue, chartSettings.axis.x.min, chartSettings.axis.x.max, chartSettings.width - chartSettings.margins.left - chartSettings.margins.right - (((_v = chartSettings.axis.x.padding) == null ? void 0 : _v.left) ? (_w = chartSettings.axis.x.padding) == null ? void 0 : _w.left : 0) - (((_x = chartSettings.axis.x.padding) == null ? void 0 : _x.right) ? (_y = chartSettings.axis.x.padding) == null ? void 0 : _y.right : 0), chartSettings.axis.x) + chartSettings.margins.left + (((_z = chartSettings.axis.x.padding) == null ? void 0 : _z.left) ? (_A = chartSettings.axis.x.padding) == null ? void 0 : _A.left : 0);
       ctx.beginPath();
-      ctx.moveTo(x, chartSettings.axis.x.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yT);
+      ctx.moveTo(x, chartSettings.axis.x.tickOnly ? chartSettings.height - chartSettings.margins.bottom + 10 : yT);
       ctx.lineTo(x, yB);
       if (((_B = chartSettings.axis.x.labels) == null ? void 0 : _B.length) === void 0) {
         chartSettings.axis.x.labels = [];
@@ -3873,9 +3913,9 @@ function AddAxisGridlines(ctx, chartSettings, axis) {
     if (gridline.dashed) {
       ctx.setLineDash([20, 15]);
     }
-    ctx.moveTo(x, gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yT);
+    ctx.moveTo(x, gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom + 10 : yT);
     ctx.lineTo(x, gridline.opposite && !gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yB);
-    ctx.fillText(gridline.label || RoundTO(gridline.value, 3), x, gridline.opposite ? chartSettings.margins.top - 10 : chartSettings.height - chartSettings.margins.bottom + 40);
+    ctx.fillText(gridline.label || RoundTO(gridline.value, 3), x, gridline.opposite ? chartSettings.margins.top - 10 : chartSettings.height - chartSettings.margins.bottom + 32);
     ctx.stroke();
     ctx.setLineDash([]);
     ctx.closePath();
@@ -4603,8 +4643,6 @@ function CreateSplitGraph(charts, title, chartSettings = new ChartSettings()) {
       ctx.fillText(chart.title, chart.chartSettings.margins.left, chart.chartSettings.margins.top - 10);
       ctx == null ? void 0 : ctx.closePath();
     }
-    chart.chartSettings.axis.x.min = NaN;
-    chart.chartSettings.axis.x.max = NaN;
     AutoAxis(chart.data, chart.chartSettings);
     AddAxisGridlines(ctx, chart.chartSettings, "x");
     PlotData(ctx, chart.data, chart.chartSettings);
@@ -5347,8 +5385,6 @@ function GRRRChart(data) {
       gl.display.strokeColor = "rgb(230,230,230)";
       setSettings.axis.x.customGridlines.push(gl);
     }
-    setSettings.axis.x.min = 1;
-    setSettings.axis.x.min = Npart;
     setSettings.axis.x.showAutoGridlines = false;
     newDataSet.line.type = "straight";
     newDataSet.line.display.strokeWidth = 2;
@@ -5432,8 +5468,6 @@ function GRRXbarChart(data) {
       gl.display.strokeColor = "rgb(230,230,230)";
       setSettings.axis.x.customGridlines.push(gl);
     }
-    setSettings.axis.x.min = 1;
-    setSettings.axis.x.min = Npart;
     setSettings.axis.x.showAutoGridlines = false;
     newDataSet.line.type = "straight";
     newDataSet.line.display.strokeWidth = 2;

package/dist/index.mjs

@@ -2674,15 +2674,17 @@ var ANOVA = {
       let scaleDF = Math.sqrt((rCount - 1) / (rCount * tCount));
       let UDL = grandMean + pooledSTD * hAlpha * scaleDF;
       let LDL = grandMean - pooledSTD * hAlpha * scaleDF;
+      let chartSettings = new ChartSettings();
       returnArray.forEach((t, x) => {
         let canvasDraw = t.Mean < UDL && t.Mean > LDL ? new CanvasDrawSettings("#2183b2", "#2183b2", void 0, 0) : new CanvasDrawSettings("#F00", "#F00", void 0, 1);
         let newDataSet = new DataSet([
           new Point(x + 1, t.Mean, canvasDraw),
           new Point(x + 1, grandMean, new CanvasDrawSettings("#000", "#000", 0, 0, 0))
         ], t.treatment, void 0, new Line("straight", new CanvasDrawSettings("#000", "#000", 1)));
+        chartSettings.axis.x.customGridlines.push(new CustomGridline(x + 1, t.treatment, void 0, void 0, void 0, void 0, true));
         dataSets.push(newDataSet);
       });
-      let chartSettings = new ChartSettings();
+      chartSettings.axis.x.showAutoGridlines = false;
       chartSettings.axis.y.customGridlines = [new CustomGridline(grandMean, `x\u0305\u0305=${Math.round(grandMean * 1e3) / 1e3}`, new CanvasDrawSettings(GetThemeColor(2), GetThemeColor(2)), void 0, true), new CustomGridline(LDL, `LDL=${Math.round(LDL * 1e3) / 1e3}`, new CanvasDrawSettings("#F00", "#F00"), void 0, true, void 0, true), new CustomGridline(UDL, `UDL=${Math.round(UDL * 1e3) / 1e3}`, new CanvasDrawSettings("#F00", "#F00"), void 0, true, void 0, true)];
       let chart = CreateScatterPlot(dataSets, name ? `Analysis of Means for ${name}` : "Analysis of Means", chartSettings);
       return {
@@ -2760,8 +2762,8 @@ var ANOVA = {
           mean,
           stdev: StDev.S(data[key]),
           SE: standardError,
-          CI_Low: mean - tInv * standardError,
-          CI_High: mean + tInv * standardError
+          CI_Low: mean - Math.abs(tInv * standardError),
+          CI_High: mean + Math.abs(tInv * standardError)
         };
         summaryTable.ConfidenceIntervals.push(CI_table);
       }
@@ -3584,6 +3586,44 @@ var ANOVA = {
         p: Distributions.F.TwoTail(F, factor1.length - 1, factor2.length - 1)
       };
     }
+  },
+  Multifactor: {
+    MainEffectsPlot(factors, response) {
+      let responseName = Object.keys(response)[0];
+      let factorKeys = Object.keys(factors);
+      let individualCharts = [];
+      factorKeys.forEach((k) => {
+        let arr = factors[k];
+        let uniqueLevels = arr.filter((value, index, array) => array.indexOf(value) === index);
+        let dataSet = new DataSet([]);
+        dataSet.line = new Line("straight");
+        let chartSettings = new ChartSettings();
+        chartSettings.axis.x.showAutoGridlines = false;
+        chartSettings.axis.x.min = 0.5;
+        chartSettings.axis.x.max = uniqueLevels.length + 0.5;
+        chartSettings.axis.x.autoSize = false;
+        uniqueLevels.forEach((level, levelI) => {
+          let mean = 0;
+          let N = 0;
+          arr.forEach((v, i) => {
+            if (v === level) {
+              mean += response[responseName][i];
+              N += 1;
+            }
+          });
+          mean = mean / N;
+          dataSet.values.push(new Point(levelI + 1, mean, new CanvasDrawSettings("#2183b2", "#000", void 0, void 0, 5)));
+          chartSettings.axis.x.customGridlines.push(new CustomGridline(levelI + 1, level, void 0, void 0, void 0, true, false));
+        });
+        chartSettings.width += 1e3;
+        individualCharts.push(new Chart([dataSet], "scatter", k, chartSettings));
+      });
+      let chartSet = new ChartSettings();
+      chartSet.width += 1e3;
+      chartSet.axis.y.tickOnly = true;
+      let graph = CreateSplitGraph(individualCharts, `Main Effects Plot of ${Object.keys(response)[0]}`, chartSet);
+      return graph;
+    }
   }
 };
 function SumSquaredDeviation(data, referenceValue) {
@@ -3759,7 +3799,7 @@ function AddAxisGridlines(ctx, chartSettings, axis) {
       }
       let x = ConvertToCanvasPt(unadjustedValue, chartSettings.axis.x.min, chartSettings.axis.x.max, chartSettings.width - chartSettings.margins.left - chartSettings.margins.right - (((_v = chartSettings.axis.x.padding) == null ? void 0 : _v.left) ? (_w = chartSettings.axis.x.padding) == null ? void 0 : _w.left : 0) - (((_x = chartSettings.axis.x.padding) == null ? void 0 : _x.right) ? (_y = chartSettings.axis.x.padding) == null ? void 0 : _y.right : 0), chartSettings.axis.x) + chartSettings.margins.left + (((_z = chartSettings.axis.x.padding) == null ? void 0 : _z.left) ? (_A = chartSettings.axis.x.padding) == null ? void 0 : _A.left : 0);
       ctx.beginPath();
-      ctx.moveTo(x, chartSettings.axis.x.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yT);
+      ctx.moveTo(x, chartSettings.axis.x.tickOnly ? chartSettings.height - chartSettings.margins.bottom + 10 : yT);
       ctx.lineTo(x, yB);
       if (((_B = chartSettings.axis.x.labels) == null ? void 0 : _B.length) === void 0) {
         chartSettings.axis.x.labels = [];
@@ -3781,9 +3821,9 @@ function AddAxisGridlines(ctx, chartSettings, axis) {
     if (gridline.dashed) {
       ctx.setLineDash([20, 15]);
     }
-    ctx.moveTo(x, gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yT);
+    ctx.moveTo(x, gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom + 10 : yT);
     ctx.lineTo(x, gridline.opposite && !gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yB);
-    ctx.fillText(gridline.label || RoundTO(gridline.value, 3), x, gridline.opposite ? chartSettings.margins.top - 10 : chartSettings.height - chartSettings.margins.bottom + 40);
+    ctx.fillText(gridline.label || RoundTO(gridline.value, 3), x, gridline.opposite ? chartSettings.margins.top - 10 : chartSettings.height - chartSettings.margins.bottom + 32);
     ctx.stroke();
     ctx.setLineDash([]);
     ctx.closePath();
@@ -4511,8 +4551,6 @@ function CreateSplitGraph(charts, title, chartSettings = new ChartSettings()) {
       ctx.fillText(chart.title, chart.chartSettings.margins.left, chart.chartSettings.margins.top - 10);
       ctx == null ? void 0 : ctx.closePath();
     }
-    chart.chartSettings.axis.x.min = NaN;
-    chart.chartSettings.axis.x.max = NaN;
     AutoAxis(chart.data, chart.chartSettings);
     AddAxisGridlines(ctx, chart.chartSettings, "x");
     PlotData(ctx, chart.data, chart.chartSettings);
@@ -5255,8 +5293,6 @@ function GRRRChart(data) {
       gl.display.strokeColor = "rgb(230,230,230)";
       setSettings.axis.x.customGridlines.push(gl);
     }
-    setSettings.axis.x.min = 1;
-    setSettings.axis.x.min = Npart;
     setSettings.axis.x.showAutoGridlines = false;
     newDataSet.line.type = "straight";
     newDataSet.line.display.strokeWidth = 2;
@@ -5340,8 +5376,6 @@ function GRRXbarChart(data) {
       gl.display.strokeColor = "rgb(230,230,230)";
       setSettings.axis.x.customGridlines.push(gl);
     }
-    setSettings.axis.x.min = 1;
-    setSettings.axis.x.min = Npart;
     setSettings.axis.x.showAutoGridlines = false;
     newDataSet.line.type = "straight";
     newDataSet.line.display.strokeWidth = 2;

package/index.ts

@@ -977,6 +977,28 @@ export const Gamma = {
   **/
   Trigamma(x: number){
     return 1 / x + 1 / (2 * (x ^ 2)) + 1 / (6 * (x ^ 3)) - 1 / (30 * (x ^ 5)) + 1 / (42 * (x ^ 7)) - 1 / (30 * (x ^ 9)) + 5 / (66 * (x ^ 11)) - 691 / (2730 * (x ^ 13)) + 15 / (6 * (x ^ 15))
+  },
+  /**
+  * Evaluates the lower incomplete gamma function for the specified value.
+  * @customfunction
+  * @param n 
+  * @param x
+  **/
+  LowerIncomplete(n: number, x: number){
+    return Gamma.fn(n) - Gamma.UpperIncomplete(n, x);
+  },
+  /**
+  * Evaluates the upper incomplete gamma function for the specified value.
+  * @customfunction
+  * @param n 
+  * @param x
+  **/
+  UpperIncomplete(n: number, x: number){
+    let expSum = 0;
+    for(let k = 0; k < n; k++){
+      expSum += Math.pow(x, k) / Probability.Factorialize(k);
+    }
+    return Probability.Factorialize(n-1) * Math.exp(-x) * expSum
   }
 }
 
@@ -1134,11 +1156,44 @@ export const Distributions = {
     }  
   },
   ChiSq: {
+    pdf(x: number, df: number){
+      let num = Math.pow(x, df/2-1) * Math.exp(-x/2)
+      let den = Math.pow(2, df/2) * Gamma.fn(df/2);
+      return num/den
+    },
     cdf(x: number, df: number){
       if(df === 2){
         return 1 - Math.exp(-x/2);
       }
-      return NaN
+      if(df%2){
+
+        let test = 0;
+        let stepCount = 10000000;
+        let stepSize = x/stepCount;
+  
+        for(let i = 0; i < stepCount + 1; i++){
+          test += stepSize * Distributions.ChiSq.pdf(i*stepSize, df);
+        }
+        return test
+        
+        //this approximation is best used when x < df
+        if(x < df){
+          let y = Math.sqrt(df/2)*Math.log(df/x);
+          const aFunc = (c1, c2, c3, c4) => {
+            return c1 + c2 * Math.pow(2/df, .5) + c3 * (2/df) + c4 * Math.pow(2/df,1/5)
+          }
+          let A = aFunc(.5, .1323, .0036, -.0038)
+          let a1 = aFunc(.1968, -.0452, -.0128, -.0168);
+          let a2 = aFunc(0.1152,-0.099,0.0539,-0.0168);
+          let a3 = aFunc(0.0004,0.0442,-0.0866,0.0398);
+          let a4 = aFunc(0.0195,-0.0629,0.0708,0.0269);
+          
+          return A/Math.pow(1 + a1*y + a2*Math.pow(y,2) + a3*Math.pow(y,3) + a4*Math.pow(y,4),4)
+        }
+      }
+      let num = Gamma.LowerIncomplete(df/2, x/2)
+      let den = Gamma.fn(df/2)
+      return num/den;
     },
     RightTail(x: number, df: number){
       return 1 - Distributions.ChiSq.cdf(x, df);
@@ -2839,6 +2894,59 @@ export function Sum(data: any[]){
   return sum;
 }
 
+function SumSquareByGroups(groups: any, response: any) {
+  let groupKeys = Object.keys(groups);
+  let groupSets: any[] = [];
+  let treatmentCount = 0;
+  groupKeys.forEach((group) => {
+    let set = new Set(groups[group]);
+    groupSets.push(set)
+    treatmentCount += set.size;
+  });
+
+  let placeFactors: any[] = [];
+  groupSets.forEach((s, sI) => {
+    if(sI === groupSets.length - 1){
+      placeFactors.push(1);
+    }else{
+      placeFactors.push(groupSets[sI+1].size)
+    }
+  })
+
+  let returnObject = {};
+  let groupNObject = {};
+  response[Object.keys(response)[0]].forEach((value, index) => {
+    let indexer: any[] = [];
+    groupKeys.forEach((k, kI) => {
+      let currentKey = groups[k][index]
+      groupSets[kI].forEach((s,sI) => {
+        if(s === currentKey){
+          indexer.push(sI)
+        }
+      })
+    })
+    returnObject[JSON.stringify(indexer)] = 0;
+    groupNObject[JSON.stringify(indexer)] = 0;
+  })
+  response[Object.keys(response)[0]].forEach((value, index) => {
+    let indexer: any = [];
+    groupKeys.forEach((k, kI) => {
+      let currentKey = groups[k][index]
+      groupSets[kI].forEach((s,sI) => {
+        if(s === currentKey){
+          indexer.push(sI)
+        }
+      })
+    })
+    returnObject[JSON.stringify(indexer)] += value;
+    groupNObject[JSON.stringify(indexer)] += 1;
+  })
+  let sumsq = 0;
+  Object.keys(returnObject).forEach(k => {
+    sumsq += Math.pow(returnObject[k], 2)/groupNObject[k]
+  })
+  return sumsq;
+}
 /**
 * Calculates the mean of the array supplied.
 * @customfunction
@@ -3018,17 +3126,18 @@ export const ANOVA = {
       let UDL = grandMean + pooledSTD*hAlpha*scaleDF
       let LDL = grandMean - pooledSTD*hAlpha*scaleDF
 
+      let chartSettings = new ChartSettings();
       returnArray.forEach((t,x) => {
         let canvasDraw: CanvasDrawSettings = (t.Mean < UDL && t.Mean > LDL) ? new CanvasDrawSettings('#2183b2', '#2183b2', undefined, 0) : new CanvasDrawSettings('#F00', '#F00',undefined,1)
         let newDataSet = new DataSet([
           new Point(x+1, t.Mean, canvasDraw),
           new Point(x+1, grandMean, new CanvasDrawSettings('#000', '#000', 0, 0, 0))
         ],t.treatment,undefined, new Line("straight",new CanvasDrawSettings('#000','#000',1)));
-
+        chartSettings.axis.x.customGridlines.push(new CustomGridline(x+1, t.treatment, undefined, undefined, undefined, undefined, true))
         dataSets.push(newDataSet);
       })
 
-      let chartSettings = new ChartSettings();
+      chartSettings.axis.x.showAutoGridlines = false;
       chartSettings.axis.y.customGridlines = [new CustomGridline(grandMean, `x̅̅=${Math.round(grandMean*1000)/1000}`, new CanvasDrawSettings(GetThemeColor(2), GetThemeColor(2)),undefined, true),new CustomGridline(LDL, `LDL=${Math.round(LDL*1000)/1000}`, new CanvasDrawSettings('#F00', '#F00'),undefined,true,undefined,true), new CustomGridline(UDL, `UDL=${Math.round(UDL*1000)/1000}`, new CanvasDrawSettings('#F00', '#F00'),undefined,true,undefined,true)]
       let chart = CreateScatterPlot(dataSets, name ? `Analysis of Means for ${name}` : "Analysis of Means", chartSettings);
 
@@ -3115,8 +3224,8 @@ export const ANOVA = {
           mean: mean,
           stdev: StDev.S(data[key]),
           SE: standardError,
-          CI_Low: mean - tInv*standardError,
-          CI_High: mean + tInv*standardError
+          CI_Low: mean - Math.abs(tInv*standardError),
+          CI_High: mean + Math.abs(tInv*standardError)
         }
         summaryTable.ConfidenceIntervals.push(CI_table);
       }
@@ -3977,6 +4086,184 @@ export const ANOVA = {
         p: Distributions.F.TwoTail(F, factor1.length - 1, factor2.length - 1)
       }
     }
+  },
+  Multifactor: {
+    MainEffectsPlot(factors: any, response: any){
+      let responseName = Object.keys(response)[0]
+      let factorKeys = Object.keys(factors)
+      let individualCharts: Chart[] = [];
+      factorKeys.forEach(k => {
+        let arr = factors[k];
+        let uniqueLevels = arr.filter((value, index, array) => array.indexOf(value) === index);
+        
+        let dataSet = new DataSet([]);
+        dataSet.line = new Line('straight')
+
+        let chartSettings = new ChartSettings();
+        chartSettings.axis.x.showAutoGridlines = false
+        chartSettings.axis.x.min = .5;
+        chartSettings.axis.x.max = uniqueLevels.length + .5
+        chartSettings.axis.x.autoSize = false
+
+        uniqueLevels.forEach((level, levelI) => {
+          let mean = 0;
+          let N = 0;
+          arr.forEach((v,i) => {
+            if(v === level){
+              mean += response[responseName][i]
+              N += 1
+            }
+          })
+          mean = mean/N;
+          dataSet.values.push(new Point(levelI+1, mean, new CanvasDrawSettings('#2183b2','#000',undefined,undefined,5)));
+          chartSettings.axis.x.customGridlines.push(new CustomGridline(levelI + 1, level,undefined, undefined,undefined, true, false))
+        })
+        chartSettings.width += 1000
+
+        individualCharts.push(new Chart([dataSet], "scatter", k, chartSettings))
+      })
+      let chartSet = new ChartSettings()
+      chartSet.width += 1000
+      chartSet.axis.y.tickOnly = true
+      let graph = CreateSplitGraph(individualCharts, `Main Effects Plot of ${Object.keys(response)[0]}`, chartSet)
+      console.log(graph)
+      return graph
+    },
+    Table(factors: any, response: any){
+      let anovaTable = {};
+      let factorKeys = Object.keys(factors);
+      let responseKey = Object.keys(response)[0];
+      let Ntotal = response[responseKey].length;
+      let DFE = 1;
+      let correctionFactor = Math.pow(Sum(response[responseKey]),2)/Ntotal;
+      // Individual Variance
+      factorKeys.forEach(f => {
+        let fSet = new Set(factors[f])
+        let Nf = fSet.size;
+        let SS_f = 0;
+
+        fSet.forEach(level => {
+          let levelTotal = 0;
+          let levelN = 0;
+          factors[f].forEach((v,j) => {
+            if(v === level){
+              levelTotal += response[responseKey][j];
+              levelN += 1;
+            }
+          })
+          SS_f += Math.pow(levelTotal, 2) / levelN
+        })
+        SS_f = SS_f - correctionFactor;
+
+        anovaTable[f] = {
+          DF: Nf-1,
+          SS: SS_f,
+          MS: SS_f/(Nf-1),
+          F: undefined,
+          p: undefined
+        }
+        DFE = DFE*Nf;
+      })
+      DFE = Ntotal - DFE;
+
+      // Interaction Effects
+      for(let i = 0; i < factorKeys.length - 1; i++){
+        let f_i_Set = new Set(factors[factorKeys[i]])
+        for(let j = i+1; j < factorKeys.length; j++){
+          let f_j_Set = new Set(factors[factorKeys[j]])          
+          let SS_subtotal = 0;
+          f_i_Set.forEach(fi => {
+            f_j_Set.forEach(fj => {   
+              let groupSum = 0;
+              let groupN = 0;
+              response[responseKey].forEach((val, index) => {
+                if(factors[factorKeys[i]][index] === fi && factors[factorKeys[j]][index] === fj){
+                  groupSum += val
+                  groupN += 1;
+                }
+              })
+              SS_subtotal += Math.pow(groupSum, 2)/groupN;
+            })
+          })
+          SS_subtotal = SS_subtotal-correctionFactor;
+
+          let iRow = anovaTable[factorKeys[i]]
+          let jRow = anovaTable[factorKeys[j]]
+          let fi_fj_DF = iRow.DF*jRow.DF
+          let fi_fj_SS = SS_subtotal - iRow.SS - jRow.SS;
+          anovaTable[`${factorKeys[i]}*${factorKeys[j]}`] = {DF: fi_fj_DF, SS: fi_fj_SS, MS: fi_fj_SS/fi_fj_DF, F: undefined, p: undefined}
+        }
+      }
+      
+      // Total Interaction Effects
+      let interactionSS_subtotal = SumSquareByGroups(factors, response) - correctionFactor;
+      let fullInteractionDF = 1;
+      let SSE = 0;
+      
+      factorKeys.forEach(source => {
+        fullInteractionDF = fullInteractionDF * anovaTable[source].DF;
+      })
+      Object.keys(anovaTable).forEach((source) => {
+        interactionSS_subtotal = interactionSS_subtotal - anovaTable[source].SS;
+        SSE += anovaTable[source].SS;
+      });
+      let fullInteractionName = factorKeys.join("*");
+      anovaTable[fullInteractionName] = {DF: fullInteractionDF, SS: interactionSS_subtotal, MS: interactionSS_subtotal/fullInteractionDF, F: undefined, p: undefined}
+      
+      SSE += interactionSS_subtotal;
+      let SStotal = SumSquaredDeviation(response[responseKey], 0) - correctionFactor;
+      SSE = SStotal - SSE;
+      let MSE = SSE/DFE;
+
+      Object.keys(anovaTable).forEach(source => {
+        let sourceF = anovaTable[source].MS/MSE;
+        anovaTable[source].F = sourceF;
+        anovaTable[source].p = Distributions.F.RightTail(sourceF, anovaTable[source].DF, DFE);
+      })
+
+      anovaTable["Error"] = {DF: DFE, SS: SSE, MS: MSE}
+      anovaTable["Total"] = {DF: Ntotal - 1, SS: SStotal}
+
+      return anovaTable;
+    },
+    GLM(factors: any, response: any){
+      let factorKeys = Object.keys(factors);
+      let anovaTable = {}
+      let responseKey = Object.keys(response)[0];
+      let Ntotal = response[responseKey].length;
+      let DFE = 1;
+      let responseData = response[responseKey];
+
+      factorKeys.forEach(k => {
+        let SS = 0;
+        let levels = new Set(factors[k]);
+        let groupedData = GroupData(factors[k], responseData);
+        
+        let descriptiveStats = {};
+        let mean_0;
+        Object.keys(groupedData).forEach((k, k_i) => {
+          let mean = Mean(groupedData[k]);
+          descriptiveStats[k].mean = mean;
+          descriptiveStats[k].std = StDev.S(groupedData[k]);
+          descriptiveStats[k].N = groupedData[k].length;
+          if(k_i === 0){
+            mean_0 = mean;
+            groupedData[k].forEach(d => {
+              SS += mean;
+            })
+          }else{
+            groupedData[k].forEach(d => {
+              SS += mean - Math.pow(mean - d, 2)
+            })
+          }
+        })
+        console.log(k, SS);
+      })
+
+
+
+      return anovaTable;
+    }
   }
 }
 
@@ -4003,6 +4290,22 @@ export const Probability = {
     }
   }
 }
+
+function GroupData(groupingArray, responseArray){
+  if(groupingArray.length != responseArray.length){return {forEachGroup(){}}}
+  let returnObj = {}
+
+  let levels = new Set(groupingArray);
+
+  levels.forEach((l: any) => {
+    returnObj[l] = [];
+  })
+
+  for(let i = 0; i < groupingArray.length; i++){
+    returnObj[groupingArray[i]].push(responseArray[i])
+  }
+  return returnObj;
+}
 // End General Statistics
 // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 /** ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -4209,7 +4512,7 @@ function AddAxisGridlines(ctx:  any, chartSettings: ChartSettings, axis: string)
               let x = ConvertToCanvasPt(unadjustedValue, chartSettings.axis.x.min, chartSettings.axis.x.max, chartSettings.width - chartSettings.margins.left - chartSettings.margins.right - (chartSettings.axis.x.padding?.left ? chartSettings.axis.x.padding?.left : 0) - (chartSettings.axis.x.padding?.right ? chartSettings.axis.x.padding?.right : 0), chartSettings.axis.x) + chartSettings.margins.left + (chartSettings.axis.x.padding?.left ? chartSettings.axis.x.padding?.left : 0);
               
               ctx.beginPath();
-                  ctx.moveTo(x, chartSettings.axis.x.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yT);
+                  ctx.moveTo(x, chartSettings.axis.x.tickOnly ? chartSettings.height - chartSettings.margins.bottom + 10 : yT);
                   ctx.lineTo(x, yB);
                   if(chartSettings.axis.x.labels?.length === undefined){
                       chartSettings.axis.x.labels = []
@@ -4233,10 +4536,10 @@ function AddAxisGridlines(ctx:  any, chartSettings: ChartSettings, axis: string)
           
           ctx.beginPath();
               if(gridline.dashed){ctx.setLineDash([20, 15])}
-              ctx.moveTo(x, gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yT);
+              ctx.moveTo(x, gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom + 10 : yT);
               ctx.lineTo(x, gridline.opposite && !gridline.tickOnly ? chartSettings.height - chartSettings.margins.bottom : yB);
               
-              ctx.fillText(gridline.label || RoundTO(gridline.value, 3), x, gridline.opposite ? chartSettings.margins.top - 10 : chartSettings.height - chartSettings.margins.bottom + 40);
+              ctx.fillText(gridline.label || RoundTO(gridline.value, 3), x, gridline.opposite ? chartSettings.margins.top - 10 : chartSettings.height - chartSettings.margins.bottom + 32);
               ctx.stroke();
               ctx.setLineDash([])
           ctx.closePath();
@@ -5110,8 +5413,8 @@ export function CreateSplitGraph(charts: Chart[], title: string, chartSettings:
       }
 
       // Automate X-axis
-      chart.chartSettings.axis.x.min = NaN;
-      chart.chartSettings.axis.x.max = NaN;
+      // chart.chartSettings.axis.x.min = NaN;
+      // chart.chartSettings.axis.x.max = NaN;
       AutoAxis(chart.data, chart.chartSettings);
       AddAxisGridlines(ctx, chart.chartSettings, 'x');
 
@@ -6029,8 +6332,6 @@ function GRRRChart(data: GRRData){
           gl.display.strokeColor = 'rgb(230,230,230)'
           setSettings.axis.x.customGridlines.push(gl)
       }
-      setSettings.axis.x.min = 1
-      setSettings.axis.x.min = Npart;
       setSettings.axis.x.showAutoGridlines = false;
       newDataSet.line.type = 'straight';
       newDataSet.line.display.strokeWidth = 2;
@@ -6123,8 +6424,6 @@ function GRRXbarChart(data: GRRData){
           gl.display.strokeColor = 'rgb(230,230,230)'
           setSettings.axis.x.customGridlines.push(gl)
       }
-      setSettings.axis.x.min = 1
-      setSettings.axis.x.min = Npart;
       setSettings.axis.x.showAutoGridlines = false;
       newDataSet.line.type = 'straight';
       newDataSet.line.display.strokeWidth = 2;

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "qesuite",
-  "version": "1.0.49",
+  "version": "1.0.51",
   "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_50.md

@@ -0,0 +1,14 @@
+Added Features:
+    - ANOVA
+        - Multifactor
+            - MainEffectsPlot
+
+Modified Features
+    - Split graph
+        - removed lines that set x.min and x.max to NaN to fix custom gridline bug
+    - GRRRChart & GGRXBarChart
+        - removed lines of code that set x.min and x.max to fix NaN custom gridlining bug (see previous note^)
+    - AnalysisOfMeans.OneWay() Chart
+        - changed x-axis to match categorical treatment groupings
+    - ANOVA.OneWay.Table
+        - fixed CI_Low/CI_High flipped bug in Confidence Intervals by using absolute value of Tcrit

package/versions/1_0_51.md

@@ -0,0 +1 @@
+- Added approximations for ChiSq Distribution CDF