qesuite 1.0.41 → 1.0.43

package/dist/index.d.mts

@@ -606,9 +606,11 @@ declare const Distributions: {
         Lilliefors: {
             Test(data: number[]): {
                 KS: number;
-                p: number;
+                p: void;
             };
-            p(KS: number, N: number): number;
+            p(KS: number, N: number): void;
+            pApprox(KS: number, N: number): number;
+            pChart(KS: number, N: number): number | "<0.01";
         };
         RyanJoiner: {
             Test(data: number[]): {
@@ -791,8 +793,6 @@ declare function IndividualDistributionPlots(data: number[]): {
     boxCox: HTMLCanvasElement;
     exponential: HTMLCanvasElement;
     exponential2p: HTMLCanvasElement;
-    gamma: HTMLCanvasElement;
-    gamma3p: HTMLCanvasElement;
     largestExtremeValue: HTMLCanvasElement;
     logistic: HTMLCanvasElement;
     loglogistic: HTMLCanvasElement;

package/dist/index.d.ts

@@ -606,9 +606,11 @@ declare const Distributions: {
         Lilliefors: {
             Test(data: number[]): {
                 KS: number;
-                p: number;
+                p: void;
             };
-            p(KS: number, N: number): number;
+            p(KS: number, N: number): void;
+            pApprox(KS: number, N: number): number;
+            pChart(KS: number, N: number): number | "<0.01";
         };
         RyanJoiner: {
             Test(data: number[]): {
@@ -791,8 +793,6 @@ declare function IndividualDistributionPlots(data: number[]): {
     boxCox: HTMLCanvasElement;
     exponential: HTMLCanvasElement;
     exponential2p: HTMLCanvasElement;
-    gamma: HTMLCanvasElement;
-    gamma3p: HTMLCanvasElement;
     largestExtremeValue: HTMLCanvasElement;
     logistic: HTMLCanvasElement;
     loglogistic: HTMLCanvasElement;

package/dist/index.js

@@ -1566,6 +1566,19 @@ var Distributions = {
         };
       },
       p(KS, N) {
+        let p = Distributions.Normal.Lilliefors.pChart(KS, N);
+        if (p === 1) {
+          p = Distributions.Normal.Lilliefors.pApprox(KS, N);
+        }
+      },
+      pApprox(KS, N) {
+        let b0 = 0.37872256037043;
+        let b1 = 1.3074818507879;
+        let b2 = 0.08861783849346;
+        let A = -(b1 + N) + Math.sqrt(Math.pow(b1 + N, 2) - 4 * b2 * (b0 - Math.pow(KS, -2)));
+        return -0.37782822932809 + 1.67819837908004 * A - 3.02959249450445 * Math.pow(A, 2) + 2.80015798142101 * Math.pow(A, 3) - 1.39874347510845 * Math.pow(A, 4) + 0.40466213484419 * Math.pow(A, 5) - 0.06353440854207 * Math.pow(A, 6) + 0.00287462087623 * Math.pow(A, 7) + 6965001311e-13 * Math.pow(A, 8) - 11872227037e-14 * Math.pow(A, 9) + 575586834e-14 * Math.pow(A, 10);
+      },
+      pChart(KS, N) {
         let pSelections = [];
         let pChart = [
           [0.3027, 0.3216, 0.3456, 0.3754, 0.4129],
@@ -1631,9 +1644,9 @@ var Distributions = {
         if (selectedValue === 0) {
           return 1;
         }
-        ;
         if (selectedValue === 5) {
-          return 0;
+          console.log("KS too LOW");
+          return "<0.01";
         }
         ;
         let p1 = pValues[selectedValue - 1];
@@ -2279,8 +2292,8 @@ function IndividualDistributionPlots(data) {
     boxCox: QQPlotChart(boxCoxData, "Normal BoxCox"),
     exponential: QQPlotChart(data, "Exponential"),
     exponential2p: QQPlotChart(data, "2-Parameter Exponential"),
-    gamma: QQPlotChart(data, "Gamma"),
-    gamma3p: QQPlotChart(data, "3-Parameter Gamma"),
+    // gamma: QQPlotChart(data, 'Gamma'),
+    // gamma3p: QQPlotChart(data, '3-Parameter Gamma'),
     largestExtremeValue: QQPlotChart(data, "Largest Extreme Value"),
     logistic: QQPlotChart(data, "Logistic"),
     loglogistic: QQPlotChart(data, "2-Parameter LogLogistic"),
@@ -2295,7 +2308,9 @@ function IndividualDistributionPlots(data) {
 }
 var QQPlot = {
   p(data, method = "Hazen") {
-    let orderedData = [...data].sort();
+    let orderedData = [...data].sort(function(a, b) {
+      return a - b;
+    });
     let PValues = [];
     let N = orderedData.length;
     for (let i = 1; i <= N; i++) {
@@ -3582,8 +3597,6 @@ function QQPlotChart(data, distribution = "Normal") {
   let dataset = new DataSet(qqData[0].map((d, i) => {
     return new Point(QQDistributionModificationX(d, distribution, data), QQDistributionModification(qqData[1][i], distribution));
   }), "Data", new CanvasDrawSettings(GetThemeColor(0), GetThemeColor(0), 1, 0, 10));
-  let mean = Mean(data);
-  let std = StDev.S(data);
   let chartSettings = new ChartSettings();
   let customGridlines = [
     new CustomGridline(QQDistributionModification(0.999, distribution), "99.9", void 0, new EngFont("black", void 0, "22px")),

package/dist/index.mjs

@@ -1479,6 +1479,19 @@ var Distributions = {
         };
       },
       p(KS, N) {
+        let p = Distributions.Normal.Lilliefors.pChart(KS, N);
+        if (p === 1) {
+          p = Distributions.Normal.Lilliefors.pApprox(KS, N);
+        }
+      },
+      pApprox(KS, N) {
+        let b0 = 0.37872256037043;
+        let b1 = 1.3074818507879;
+        let b2 = 0.08861783849346;
+        let A = -(b1 + N) + Math.sqrt(Math.pow(b1 + N, 2) - 4 * b2 * (b0 - Math.pow(KS, -2)));
+        return -0.37782822932809 + 1.67819837908004 * A - 3.02959249450445 * Math.pow(A, 2) + 2.80015798142101 * Math.pow(A, 3) - 1.39874347510845 * Math.pow(A, 4) + 0.40466213484419 * Math.pow(A, 5) - 0.06353440854207 * Math.pow(A, 6) + 0.00287462087623 * Math.pow(A, 7) + 6965001311e-13 * Math.pow(A, 8) - 11872227037e-14 * Math.pow(A, 9) + 575586834e-14 * Math.pow(A, 10);
+      },
+      pChart(KS, N) {
         let pSelections = [];
         let pChart = [
           [0.3027, 0.3216, 0.3456, 0.3754, 0.4129],
@@ -1544,9 +1557,9 @@ var Distributions = {
         if (selectedValue === 0) {
           return 1;
         }
-        ;
         if (selectedValue === 5) {
-          return 0;
+          console.log("KS too LOW");
+          return "<0.01";
         }
         ;
         let p1 = pValues[selectedValue - 1];
@@ -2192,8 +2205,8 @@ function IndividualDistributionPlots(data) {
     boxCox: QQPlotChart(boxCoxData, "Normal BoxCox"),
     exponential: QQPlotChart(data, "Exponential"),
     exponential2p: QQPlotChart(data, "2-Parameter Exponential"),
-    gamma: QQPlotChart(data, "Gamma"),
-    gamma3p: QQPlotChart(data, "3-Parameter Gamma"),
+    // gamma: QQPlotChart(data, 'Gamma'),
+    // gamma3p: QQPlotChart(data, '3-Parameter Gamma'),
     largestExtremeValue: QQPlotChart(data, "Largest Extreme Value"),
     logistic: QQPlotChart(data, "Logistic"),
     loglogistic: QQPlotChart(data, "2-Parameter LogLogistic"),
@@ -2208,7 +2221,9 @@ function IndividualDistributionPlots(data) {
 }
 var QQPlot = {
   p(data, method = "Hazen") {
-    let orderedData = [...data].sort();
+    let orderedData = [...data].sort(function(a, b) {
+      return a - b;
+    });
     let PValues = [];
     let N = orderedData.length;
     for (let i = 1; i <= N; i++) {
@@ -3495,8 +3510,6 @@ function QQPlotChart(data, distribution = "Normal") {
   let dataset = new DataSet(qqData[0].map((d, i) => {
     return new Point(QQDistributionModificationX(d, distribution, data), QQDistributionModification(qqData[1][i], distribution));
   }), "Data", new CanvasDrawSettings(GetThemeColor(0), GetThemeColor(0), 1, 0, 10));
-  let mean = Mean(data);
-  let std = StDev.S(data);
   let chartSettings = new ChartSettings();
   let customGridlines = [
     new CustomGridline(QQDistributionModification(0.999, distribution), "99.9", void 0, new EngFont("black", void 0, "22px")),

package/index.ts

@@ -1648,6 +1648,19 @@ export const Distributions = {
         }
       },
       p(KS: number, N: number){
+        let p = Distributions.Normal.Lilliefors.pChart(KS, N);
+        if(p === 1){
+          p = Distributions.Normal.Lilliefors.pApprox(KS, N)
+        }
+      },
+      pApprox(KS: number, N: number){
+        let b0 = 0.37872256037043;
+        let b1 = 1.30748185078790;
+        let b2 = 0.08861783849346;
+        let A = (-(b1+N)+Math.sqrt(Math.pow(b1+N, 2)-4*b2*(b0-Math.pow(KS, -2))))
+        return -.37782822932809+1.67819837908004*A-3.02959249450445*Math.pow(A,2)+2.80015798142101*Math.pow(A,3)-1.39874347510845*Math.pow(A,4)+0.40466213484419*Math.pow(A,5)-0.06353440854207*Math.pow(A,6) +0.00287462087623*Math.pow(A,7)+0.00069650013110*Math.pow(A,8)-0.00011872227037*Math.pow(A,9)+0.00000575586834*Math.pow(A,10)
+      },
+      pChart(KS: number, N: number){
         let pSelections: number[] = [];
         let pChart = [
           [0.3027,0.3216,0.3456,0.3754,0.4129],
@@ -1710,8 +1723,13 @@ export const Distributions = {
             selectedValue += 1
           }
         })
-        if(selectedValue === 0){return 1};
-        if(selectedValue === 5){return 0};
+        if(selectedValue === 0){
+          return 1
+        }
+        if(selectedValue === 5){
+          console.log("KS too LOW")
+          return "<0.01"
+        };
         let p1 = pValues[selectedValue-1];
         let p2 = pValues[selectedValue];
         let KS1 = pSelections[selectedValue-1];
@@ -2466,8 +2484,8 @@ export function IndividualDistributionPlots(data: number[]){
     boxCox: QQPlotChart(boxCoxData, "Normal BoxCox"),
     exponential: QQPlotChart(data, 'Exponential'),
     exponential2p: QQPlotChart(data, '2-Parameter Exponential'),
-    gamma: QQPlotChart(data, 'Gamma'),
-    gamma3p: QQPlotChart(data, '3-Parameter Gamma'),
+    // gamma: QQPlotChart(data, 'Gamma'),
+    // gamma3p: QQPlotChart(data, '3-Parameter Gamma'),
     largestExtremeValue: QQPlotChart(data, 'Largest Extreme Value'),
     logistic: QQPlotChart(data, 'Logistic'),
     loglogistic: QQPlotChart(data, '2-Parameter LogLogistic'),
@@ -2484,7 +2502,7 @@ export function IndividualDistributionPlots(data: number[]){
 
 export const QQPlot = {
   p(data: number[], method: string = 'Hazen'){
-    let orderedData = [...data].sort();
+    let orderedData = [...data].sort(function(a,b){return a-b});
     let PValues: number[] = [];
     let N = orderedData.length
     for(let i = 1; i <= N; i++){
@@ -3373,6 +3391,7 @@ function PlotPoint(ctx: any, value: Point, dataSetDisplay: CanvasDrawSettings, c
       
       let x = ConvertToCanvasPt(value.x, chartSettings.axis.x.min, chartSettings.axis.x.max, chartSettings.width - chartSettings.margins.left - chartSettings.margins.right - padding.left - padding.right, chartSettings.axis.x);
       let y = ConvertToCanvasPt(value.y, chartSettings.axis.y.max, chartSettings.axis.y.min, chartSettings.height - chartSettings.margins.top - chartSettings.margins.bottom - padding.top - padding.bottom, chartSettings.axis.y);       
+      
       if(y >= 0 && y <= chartSettings.height - chartSettings.margins.bottom - chartSettings.margins.top && x >= 0 && x <= chartSettings.width - chartSettings.margins.right - chartSettings.margins.left){     
           DrawPolygon(ctx, shape, x + chartSettings.margins.left + padding.left, y + chartSettings.margins.top + padding.top, shapeRadius);
       }
@@ -4075,8 +4094,6 @@ export function QQPlotChart(data: number[], distribution: string = 'Normal'){
   let qqData = QQPlot.p([...data]);
 
   let dataset = new DataSet(qqData[0].map((d, i) => {return new Point(QQDistributionModificationX(d, distribution, data), QQDistributionModification(qqData[1][i], distribution))}), 'Data', new CanvasDrawSettings(GetThemeColor(0), GetThemeColor(0), 1, 0, 10));
-  let mean = Mean(data);
-  let std = StDev.S(data);
   
   // Create Chart Settings
   let chartSettings = new ChartSettings();

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "qesuite",
-  "version": "1.0.41",
+  "version": "1.0.43",
   "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_42.md

@@ -0,0 +1,3 @@
+- Added approximation method for evaluating p-value in Lilliefors test
+
+- take a look at distribution identification make sure only graphs that are supported with AD and MLE are generated.

package/versions/1_0_43.md

@@ -0,0 +1,2 @@
+- separated both methods of p-value calculation in Lilliefors test.
+- fixed bug in QQPlot that sorted negative values descending