@danielsimonjr/mathts-functions 0.68.3 → 0.68.4

package/dist/algebra/symbolicEqual.d.ts

@@ -9,7 +9,7 @@ interface SimplifyOptions {
 }
 interface SymbolicEqualDependencies {
     parse: (expr: string) => MathNode;
-    simplify: (node: MathNode, rules: unknown[], options: SimplifyOptions) => MathNode;
+    simplify: (node: MathNode, scope: Map<string, unknown>, options: SimplifyOptions) => MathNode;
     typed: TypedFunction;
     OperatorNode: OperatorNodeConstructor;
 }

package/dist/algebra/symbolicEqual.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"symbolicEqual.d.ts","sourceRoot":"","sources":["../../src/algebra/symbolicEqual.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,kBAAkB,CAAC;AACjD,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAO/D,UAAU,uBAAuB;IAC/B,KAAK,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,GAAG,QAAQ,CAAC;CAC1D;AAED,UAAU,eAAe;IACvB,OAAO,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;IAClC,CAAC,GAAG,EAAE,MAAM,GAAG,OAAO,CAAC;CACxB;AAED,UAAU,yBAAyB;IACjC,KAAK,EAAE,CAAC,IAAI,EAAE,MAAM,KAAK,QAAQ,CAAC;IAClC,QAAQ,EAAE,CAAC,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,OAAO,EAAE,EAAE,OAAO,EAAE,eAAe,KAAK,QAAQ,CAAC;IACnF,KAAK,EAAE,aAAa,CAAC;IACrB,YAAY,EAAE,uBAAuB,CAAC;CACvC;AAKD,eAAO,MAAM,mBAAmB,yFAwD/B,CAAC"}
+{"version":3,"file":"symbolicEqual.d.ts","sourceRoot":"","sources":["../../src/algebra/symbolicEqual.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,kBAAkB,CAAC;AACjD,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAO/D,UAAU,uBAAuB;IAC/B,KAAK,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,GAAG,QAAQ,CAAC;CAC1D;AAED,UAAU,eAAe;IACvB,OAAO,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;IAClC,CAAC,GAAG,EAAE,MAAM,GAAG,OAAO,CAAC;CACxB;AAED,UAAU,yBAAyB;IACjC,KAAK,EAAE,CAAC,IAAI,EAAE,MAAM,KAAK,QAAQ,CAAC;IAClC,QAAQ,EAAE,CAAC,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,GAAG,CAAC,MAAM,EAAE,OAAO,CAAC,EAAE,OAAO,EAAE,eAAe,KAAK,QAAQ,CAAC;IAC9F,KAAK,EAAE,aAAa,CAAC;IACrB,YAAY,EAAE,uBAAuB,CAAC;CACvC;AAKD,eAAO,MAAM,mBAAmB,yFA0D/B,CAAC"}

package/dist/distribution-functions.d.ts

@@ -27,4 +27,11 @@ export declare const laplaceQuantile: (p: number, mu?: number, b?: number) => nu
 export declare const logisticPDF: (x: number, mu?: number, s?: number) => number;
 export declare const logisticCDF: (x: number, mu?: number, s?: number) => number;
 export declare const logisticQuantile: (p: number, mu?: number, s?: number) => number;
+/**
+ * Two-sided p-value of a Student-t statistic `t` with `df` degrees of freedom,
+ * computed DIRECTLY in the tail as the regularized incomplete beta
+ * `I_{df/(df+t²)}(df/2, 1/2)`. `2·(1 − studentTCDF(|t|, df))` rounds a far
+ * tail to 0, and a steep regression's slope test IS that tail.
+ */
+export declare function studentTTwoSidedP(t: number, df: number): number;
 //# sourceMappingURL=distribution-functions.d.ts.map

package/dist/distribution-functions.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"distribution-functions.d.ts","sourceRoot":"","sources":["../src/distribution-functions.ts"],"names":[],"mappings":"AAmBA,2DAA2D;AAC3D,eAAO,MAAM,cAAc,MAAO,MAAM,kCAAsB,MAC/B,CAAC;AAEhC,+DAA+D;AAC/D,eAAO,MAAM,WAAW,MAAO,MAAM,MAAM,MAAM,KAAG,MAA6B,CAAC;AAClF,eAAO,MAAM,gBAAgB,MAAO,MAAM,MAAM,MAAM,KAAG,MAAkC,CAAC;AAE5F,qEAAqE;AACrE,eAAO,MAAM,aAAa,MAAO,MAAM,MAAM,MAAM,KAAG,MAA+B,CAAC;AACtF,eAAO,MAAM,kBAAkB,MAAO,MAAM,MAAM,MAAM,KAAG,MAAoC,CAAC;AAEhG,wEAAwE;AACxE,eAAO,MAAM,IAAI,MAAO,MAAM,MAAM,MAAM,MAAM,MAAM,KAAG,MAAgC,CAAC;AAC1F,eAAO,MAAM,SAAS,MAAO,MAAM,MAAM,MAAM,MAAM,MAAM,KAAG,MAAqC,CAAC;AAEpG,oDAAoD;AACpD,eAAO,MAAM,QAAQ,MAAO,MAAM,SAAS,MAAM,oBAAa,MAAoC,CAAC;AACnG,eAAO,MAAM,aAAa,MAAO,MAAM,SAAS,MAAM,oBAAa,MAClC,CAAC;AAElC,0DAA0D;AAC1D,eAAO,MAAM,OAAO,MAAO,MAAM,KAAK,MAAM,KAAK,MAAM,KAAG,MAA4B,CAAC;AACvF,eAAO,MAAM,YAAY,MAAO,MAAM,KAAK,MAAM,KAAK,MAAM,KAAG,MAAiC,CAAC;AAIjG,4EAA4E;AAC5E,eAAO,MAAM,SAAS,MAAO,MAAM,kCAAsB,MACF,CAAC;AACxD,eAAO,MAAM,SAAS,MAAO,MAAM,kCAAsB,MACZ,CAAC;AAC9C,eAAO,MAAM,cAAc,MAAO,MAAM,kCAAsB,MAClB,CAAC;AAE7C,oFAAoF;AACpF,eAAO,MAAM,UAAU,MAAO,MAAM,8BAAkB,MACX,CAAC;AAC5C,eAAO,MAAM,UAAU,MAAO,MAAM,8BAAkB,MACqB,CAAC;AAC5E,eAAO,MAAM,eAAe,MAAO,MAAM,8BAAkB,MACQ,CAAC;AAEpE,gEAAgE;AAChE,eAAO,MAAM,WAAW,MAAO,MAAM,8BAAkB,MAGtD,CAAC;AACF,eAAO,MAAM,WAAW,MAAO,MAAM,8BAAkB,MAA2C,CAAC;AACnG,eAAO,MAAM,gBAAgB,MAAO,MAAM,8BAAkB,MAC5B,CAAC"}
+{"version":3,"file":"distribution-functions.d.ts","sourceRoot":"","sources":["../src/distribution-functions.ts"],"names":[],"mappings":"AAoBA,2DAA2D;AAC3D,eAAO,MAAM,cAAc,MAAO,MAAM,kCAAsB,MAC/B,CAAC;AAEhC,+DAA+D;AAC/D,eAAO,MAAM,WAAW,MAAO,MAAM,MAAM,MAAM,KAAG,MAA6B,CAAC;AAClF,eAAO,MAAM,gBAAgB,MAAO,MAAM,MAAM,MAAM,KAAG,MAAkC,CAAC;AAE5F,qEAAqE;AACrE,eAAO,MAAM,aAAa,MAAO,MAAM,MAAM,MAAM,KAAG,MAA+B,CAAC;AACtF,eAAO,MAAM,kBAAkB,MAAO,MAAM,MAAM,MAAM,KAAG,MAAoC,CAAC;AAEhG,wEAAwE;AACxE,eAAO,MAAM,IAAI,MAAO,MAAM,MAAM,MAAM,MAAM,MAAM,KAAG,MAAgC,CAAC;AAC1F,eAAO,MAAM,SAAS,MAAO,MAAM,MAAM,MAAM,MAAM,MAAM,KAAG,MAAqC,CAAC;AAEpG,oDAAoD;AACpD,eAAO,MAAM,QAAQ,MAAO,MAAM,SAAS,MAAM,oBAAa,MAAoC,CAAC;AACnG,eAAO,MAAM,aAAa,MAAO,MAAM,SAAS,MAAM,oBAAa,MAClC,CAAC;AAElC,0DAA0D;AAC1D,eAAO,MAAM,OAAO,MAAO,MAAM,KAAK,MAAM,KAAK,MAAM,KAAG,MAA4B,CAAC;AACvF,eAAO,MAAM,YAAY,MAAO,MAAM,KAAK,MAAM,KAAK,MAAM,KAAG,MAAiC,CAAC;AAIjG,4EAA4E;AAC5E,eAAO,MAAM,SAAS,MAAO,MAAM,kCAAsB,MACF,CAAC;AACxD,eAAO,MAAM,SAAS,MAAO,MAAM,kCAAsB,MACZ,CAAC;AAC9C,eAAO,MAAM,cAAc,MAAO,MAAM,kCAAsB,MAClB,CAAC;AAE7C,oFAAoF;AACpF,eAAO,MAAM,UAAU,MAAO,MAAM,8BAAkB,MACX,CAAC;AAC5C,eAAO,MAAM,UAAU,MAAO,MAAM,8BAAkB,MACqB,CAAC;AAC5E,eAAO,MAAM,eAAe,MAAO,MAAM,8BAAkB,MACQ,CAAC;AAEpE,gEAAgE;AAChE,eAAO,MAAM,WAAW,MAAO,MAAM,8BAAkB,MAGtD,CAAC;AACF,eAAO,MAAM,WAAW,MAAO,MAAM,8BAAkB,MAA2C,CAAC;AACnG,eAAO,MAAM,gBAAgB,MAAO,MAAM,8BAAkB,MAC5B,CAAC;AAEjC;;;;;GAKG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM,CAM/D"}

package/dist/index.js

@@ -18401,6 +18401,13 @@ var logisticPDF = (x, mu = 0, s = 1) => {
 };
 var logisticCDF = (x, mu = 0, s = 1) => 1 / (1 + Math.exp(-(x - mu) / s));
 var logisticQuantile = (p, mu = 0, s = 1) => mu + s * Math.log(p / (1 - p));
+function studentTTwoSidedP(t, df) {
+  if (!(df > 0)) return NaN;
+  if (Number.isNaN(t)) return NaN;
+  if (!Number.isFinite(t)) return 0;
+  if (t === 0) return 1;
+  return betaincScalar(df / 2, 0.5, df / (df + t * t));
+}
 
 // src/stats/inference-extra.ts
 function chi2Contingency(table, opts = {}) {
@@ -42290,7 +42297,7 @@ var createSymbolicEqual = /* @__PURE__ */ factory(
   ({ parse: _parse2, simplify: simplify2, typed: typed3, OperatorNode }) => {
     function _symbolicEqual(e1, e2, options = {}) {
       const diff2 = new OperatorNode("-", "subtract", [e1, e2]);
-      const simplified = simplify2(diff2, [], options);
+      const simplified = simplify2(diff2, /* @__PURE__ */ new Map(), options);
       return isConstantNode(simplified) && !simplified.value;
     }
     return typed3(name248, {
@@ -48753,6 +48760,15 @@ function matVec(a, v) {
   return a.map((row2) => row2.reduce((sum3, val, j) => sum3 + val * v[j], 0));
 }
 function invert(a) {
+  const scale4 = a.map((row2, i) => Math.sqrt(Math.abs(row2[i])));
+  if (scale4.some((s) => !(s > 0) || !Number.isFinite(s))) {
+    throw new Error("ols: X^T X is singular (predictors are collinear)");
+  }
+  const equilibrated = a.map((row2, i) => row2.map((v, j) => v / (scale4[i] * scale4[j])));
+  const inv3 = invertUnscaled(equilibrated);
+  return inv3.map((row2, i) => row2.map((v, j) => v / (scale4[i] * scale4[j])));
+}
+function invertUnscaled(a) {
   const n = a.length;
   const M = a.map((row2, i) => [...row2, ...Array.from({ length: n }, (_, j) => i === j ? 1 : 0)]);
   for (let col = 0; col < n; col++) {
@@ -48764,7 +48780,7 @@ function invert(a) {
         pivotRow = row2;
       }
     }
-    if (maxVal < 1e-15) throw new Error("ols: X^T X is singular (predictors are collinear)");
+    if (maxVal < 1e-12) throw new Error("ols: X^T X is singular (predictors are collinear)");
     if (pivotRow !== col) [M[col], M[pivotRow]] = [M[pivotRow], M[col]];
     const pivot = M[col][col];
     for (let j = 0; j < 2 * n; j++) M[col][j] /= pivot;
@@ -48801,7 +48817,7 @@ function ols(X, y, opts) {
   const sigma2 = rss / dfResid;
   const stderr = XtXinv.map((row2, j) => Math.sqrt(sigma2 * row2[j]));
   const tValues = coefficients.map((b, j) => b / stderr[j]);
-  const pValues = tValues.map((t) => 2 * (1 - studentTCDF(Math.abs(t), dfResid)));
+  const pValues = tValues.map((t) => studentTTwoSidedP(t, dfResid));
   const r2 = tss === 0 ? 1 : 1 - rss / tss;
   const adjR2 = 1 - (1 - r2) * (n - 1) / dfResid;
   const fStat = p > 1 ? (tss - rss) / (p - 1) / (rss / dfResid) : NaN;
@@ -48953,7 +48969,7 @@ function linearRegression(x, y) {
   const df = n - 2;
   const stdErr = Math.sqrt((syy - slope * sxy) / df / sxx);
   const t = slope / stdErr;
-  const pValue = 2 * (1 - studentTCDF(Math.abs(t), df));
+  const pValue = studentTTwoSidedP(t, df);
   return { slope, intercept, rValue, rSquared: rValue * rValue, pValue, stdErr };
 }
 

package/dist/ml/ols.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"ols.d.ts","sourceRoot":"","sources":["../../src/ml/ols.ts"],"names":[],"mappings":"AAUA,yBAAyB;AACzB,MAAM,WAAW,UAAU;IACzB,+CAA+C;IAC/C,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB;AAED,wEAAwE;AACxE,MAAM,WAAW,SAAS;IACxB,0DAA0D;IAC1D,YAAY,EAAE,MAAM,EAAE,CAAC;IACvB,0CAA0C;IAC1C,MAAM,EAAE,MAAM,EAAE,CAAC;IACjB,8DAA8D;IAC9D,OAAO,EAAE,MAAM,EAAE,CAAC;IAClB,4DAA4D;IAC5D,OAAO,EAAE,MAAM,EAAE,CAAC;IAClB,oCAAoC;IACpC,EAAE,EAAE,MAAM,CAAC;IACX,gDAAgD;IAChD,KAAK,EAAE,MAAM,CAAC;IACd,kEAAkE;IAClE,KAAK,EAAE,MAAM,CAAC;IACd,wBAAwB;IACxB,SAAS,EAAE,MAAM,EAAE,CAAC;CACrB;AAqED;;;;;;;;;;;;GAYG;AACH,wBAAgB,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,UAAU,GAAG,SAAS,CAmC5E"}
+{"version":3,"file":"ols.d.ts","sourceRoot":"","sources":["../../src/ml/ols.ts"],"names":[],"mappings":"AAUA,yBAAyB;AACzB,MAAM,WAAW,UAAU;IACzB,+CAA+C;IAC/C,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB;AAED,wEAAwE;AACxE,MAAM,WAAW,SAAS;IACxB,0DAA0D;IAC1D,YAAY,EAAE,MAAM,EAAE,CAAC;IACvB,0CAA0C;IAC1C,MAAM,EAAE,MAAM,EAAE,CAAC;IACjB,8DAA8D;IAC9D,OAAO,EAAE,MAAM,EAAE,CAAC;IAClB,4DAA4D;IAC5D,OAAO,EAAE,MAAM,EAAE,CAAC;IAClB,oCAAoC;IACpC,EAAE,EAAE,MAAM,CAAC;IACX,gDAAgD;IAChD,KAAK,EAAE,MAAM,CAAC;IACd,kEAAkE;IAClE,KAAK,EAAE,MAAM,CAAC;IACd,wBAAwB;IACxB,SAAS,EAAE,MAAM,EAAE,CAAC;CACrB;AAsFD;;;;;;;;;;;;GAYG;AACH,wBAAgB,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,UAAU,GAAG,SAAS,CAmC5E"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-functions",
-  "version": "0.68.3",
+  "version": "0.68.4",
   "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
   "author": "Daniel Simon Jr.",
   "license": "MIT",
@@ -34,7 +34,7 @@
     "build:prod": "tsup src/index.ts --format esm --clean --minify --treeshake && tsc -p tsconfig.dts.json && node scripts/copy-wasm.mjs"
   },
   "dependencies": {
-    "@danielsimonjr/mathts-core": "^0.18.0",
+    "@danielsimonjr/mathts-core": "^0.18.1",
     "@danielsimonjr/mathts-expression": "^0.10.2",
     "@danielsimonjr/mathts-gpu": "^0.2.1",
     "@danielsimonjr/mathts-matrix": "^0.7.8",