pragmastat 5.0.0 → 5.2.0

package/README.md

@@ -3,10 +3,10 @@
 Install from npm:
 
 ```bash
-npm i pragmastat@5.0.0
+npm i pragmastat@5.2.0
 ```
 
-Source code: https://github.com/AndreyAkinshin/pragmastat/tree/v5.0.0/ts
+Source code: https://github.com/AndreyAkinshin/pragmastat/tree/v5.2.0/ts
 
 Pragmastat on npm: https://www.npmjs.com/package/pragmastat
 

package/dist/assumptions.d.ts

@@ -0,0 +1,33 @@
+/**
+ * Assumption validation framework for Pragmastat estimators.
+ *
+ * Assumption IDs (canonical priority order):
+ *   1. validity - non-empty input with finite defined real values
+ *   2. positivity - values must be strictly positive
+ *   3. sparity - sample must be non tie-dominant (Spread > 0)
+ *
+ * When multiple assumptions are violated, the violation with highest priority
+ * is reported. For two-sample functions, subject X is checked before Y.
+ */
+export declare const AssumptionId: {
+    readonly VALIDITY: "validity";
+    readonly POSITIVITY: "positivity";
+    readonly SPARITY: "sparity";
+};
+export type AssumptionId = (typeof AssumptionId)[keyof typeof AssumptionId];
+export type Subject = 'x' | 'y';
+export interface Violation {
+    id: AssumptionId;
+    subject: Subject;
+}
+export declare class AssumptionError extends Error {
+    readonly violation: Violation;
+    constructor(violation: Violation);
+    static validity(_functionName: string, subject: Subject): AssumptionError;
+    static positivity(_functionName: string, subject: Subject): AssumptionError;
+    static sparity(_functionName: string, subject: Subject): AssumptionError;
+}
+export declare function checkValidity(values: number[], subject: Subject, functionName: string): void;
+export declare function checkPositivity(values: number[], subject: Subject, functionName: string): void;
+export declare function checkSparity(values: number[], subject: Subject, functionName: string): void;
+//# sourceMappingURL=assumptions.d.ts.map

package/dist/assumptions.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"assumptions.d.ts","sourceRoot":"","sources":["../src/assumptions.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;GAUG;AAEH,eAAO,MAAM,YAAY;;;;CAIf,CAAC;AAEX,MAAM,MAAM,YAAY,GAAG,CAAC,OAAO,YAAY,CAAC,CAAC,MAAM,OAAO,YAAY,CAAC,CAAC;AAE5E,MAAM,MAAM,OAAO,GAAG,GAAG,GAAG,GAAG,CAAC;AAEhC,MAAM,WAAW,SAAS;IACxB,EAAE,EAAE,YAAY,CAAC;IACjB,OAAO,EAAE,OAAO,CAAC;CAClB;AAED,qBAAa,eAAgB,SAAQ,KAAK;IACxC,QAAQ,CAAC,SAAS,EAAE,SAAS,CAAC;gBAElB,SAAS,EAAE,SAAS;IAMhC,MAAM,CAAC,QAAQ,CAAC,aAAa,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,GAAG,eAAe;IAIzE,MAAM,CAAC,UAAU,CAAC,aAAa,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,GAAG,eAAe;IAI3E,MAAM,CAAC,OAAO,CAAC,aAAa,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,GAAG,eAAe;CAGzE;AAED,wBAAgB,aAAa,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,YAAY,EAAE,MAAM,GAAG,IAAI,CAO5F;AAED,wBAAgB,eAAe,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,YAAY,EAAE,MAAM,GAAG,IAAI,CAI9F;AAED,wBAAgB,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,YAAY,EAAE,MAAM,GAAG,IAAI,CAQ3F"}

package/dist/assumptions.js

@@ -0,0 +1,63 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.AssumptionError = exports.AssumptionId = void 0;
+exports.checkValidity = checkValidity;
+exports.checkPositivity = checkPositivity;
+exports.checkSparity = checkSparity;
+const fastSpread_1 = require("./fastSpread");
+/**
+ * Assumption validation framework for Pragmastat estimators.
+ *
+ * Assumption IDs (canonical priority order):
+ *   1. validity - non-empty input with finite defined real values
+ *   2. positivity - values must be strictly positive
+ *   3. sparity - sample must be non tie-dominant (Spread > 0)
+ *
+ * When multiple assumptions are violated, the violation with highest priority
+ * is reported. For two-sample functions, subject X is checked before Y.
+ */
+exports.AssumptionId = {
+    VALIDITY: 'validity',
+    POSITIVITY: 'positivity',
+    SPARITY: 'sparity',
+};
+class AssumptionError extends Error {
+    constructor(violation) {
+        super(`${violation.id}(${violation.subject})`);
+        this.name = 'AssumptionError';
+        this.violation = violation;
+    }
+    static validity(_functionName, subject) {
+        return new AssumptionError({ id: exports.AssumptionId.VALIDITY, subject });
+    }
+    static positivity(_functionName, subject) {
+        return new AssumptionError({ id: exports.AssumptionId.POSITIVITY, subject });
+    }
+    static sparity(_functionName, subject) {
+        return new AssumptionError({ id: exports.AssumptionId.SPARITY, subject });
+    }
+}
+exports.AssumptionError = AssumptionError;
+function checkValidity(values, subject, functionName) {
+    if (values.length === 0) {
+        throw AssumptionError.validity(functionName, subject);
+    }
+    if (!values.every((v) => Number.isFinite(v))) {
+        throw AssumptionError.validity(functionName, subject);
+    }
+}
+function checkPositivity(values, subject, functionName) {
+    if (values.some((v) => v <= 0)) {
+        throw AssumptionError.positivity(functionName, subject);
+    }
+}
+function checkSparity(values, subject, functionName) {
+    if (values.length < 2) {
+        throw AssumptionError.sparity(functionName, subject);
+    }
+    const spread = (0, fastSpread_1.fastSpread)(values);
+    if (spread <= 0) {
+        throw AssumptionError.sparity(functionName, subject);
+    }
+}
+//# sourceMappingURL=assumptions.js.map

package/dist/assumptions.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"assumptions.js","sourceRoot":"","sources":["../src/assumptions.ts"],"names":[],"mappings":";;;AAmDA,sCAOC;AAED,0CAIC;AAED,oCAQC;AA1ED,6CAA0C;AAE1C;;;;;;;;;;GAUG;AAEU,QAAA,YAAY,GAAG;IAC1B,QAAQ,EAAE,UAAU;IACpB,UAAU,EAAE,YAAY;IACxB,OAAO,EAAE,SAAS;CACV,CAAC;AAWX,MAAa,eAAgB,SAAQ,KAAK;IAGxC,YAAY,SAAoB;QAC9B,KAAK,CAAC,GAAG,SAAS,CAAC,EAAE,IAAI,SAAS,CAAC,OAAO,GAAG,CAAC,CAAC;QAC/C,IAAI,CAAC,IAAI,GAAG,iBAAiB,CAAC;QAC9B,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;IAC7B,CAAC;IAED,MAAM,CAAC,QAAQ,CAAC,aAAqB,EAAE,OAAgB;QACrD,OAAO,IAAI,eAAe,CAAC,EAAE,EAAE,EAAE,oBAAY,CAAC,QAAQ,EAAE,OAAO,EAAE,CAAC,CAAC;IACrE,CAAC;IAED,MAAM,CAAC,UAAU,CAAC,aAAqB,EAAE,OAAgB;QACvD,OAAO,IAAI,eAAe,CAAC,EAAE,EAAE,EAAE,oBAAY,CAAC,UAAU,EAAE,OAAO,EAAE,CAAC,CAAC;IACvE,CAAC;IAED,MAAM,CAAC,OAAO,CAAC,aAAqB,EAAE,OAAgB;QACpD,OAAO,IAAI,eAAe,CAAC,EAAE,EAAE,EAAE,oBAAY,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC,CAAC;IACpE,CAAC;CACF;AApBD,0CAoBC;AAED,SAAgB,aAAa,CAAC,MAAgB,EAAE,OAAgB,EAAE,YAAoB;IACpF,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACxB,MAAM,eAAe,CAAC,QAAQ,CAAC,YAAY,EAAE,OAAO,CAAC,CAAC;IACxD,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QAC7C,MAAM,eAAe,CAAC,QAAQ,CAAC,YAAY,EAAE,OAAO,CAAC,CAAC;IACxD,CAAC;AACH,CAAC;AAED,SAAgB,eAAe,CAAC,MAAgB,EAAE,OAAgB,EAAE,YAAoB;IACtF,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;QAC/B,MAAM,eAAe,CAAC,UAAU,CAAC,YAAY,EAAE,OAAO,CAAC,CAAC;IAC1D,CAAC;AACH,CAAC;AAED,SAAgB,YAAY,CAAC,MAAgB,EAAE,OAAgB,EAAE,YAAoB;IACnF,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACtB,MAAM,eAAe,CAAC,OAAO,CAAC,YAAY,EAAE,OAAO,CAAC,CAAC;IACvD,CAAC;IACD,MAAM,MAAM,GAAG,IAAA,uBAAU,EAAC,MAAM,CAAC,CAAC;IAClC,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;QAChB,MAAM,eAAe,CAAC,OAAO,CAAC,YAAY,EAAE,OAAO,CAAC,CAAC;IACvD,CAAC;AACH,CAAC"}

package/dist/estimators.d.ts

@@ -11,14 +11,24 @@ export declare function center(x: number[]): number;
 /**
  * Calculate the Spread - median of all pairwise absolute differences |x[i] - x[j]|
  * Uses fast O(n log n) algorithm.
+ *
+ * Assumptions:
+ *   sparity(x) - sample must be non tie-dominant (Spread > 0)
+ *
  * @param x Array of sample values
  * @returns The spread estimate
+ * @throws AssumptionError if sample is empty, contains NaN/Inf, or is tie-dominant
  */
 export declare function spread(x: number[]): number;
 /**
  * Calculate the RelSpread - ratio of Spread to absolute Center
+ *
+ * Assumptions:
+ *   positivity(x) - all values must be strictly positive (ensures Center > 0)
+ *
  * @param x Array of sample values
  * @returns The relative spread estimate
+ * @throws AssumptionError if sample is empty, contains NaN/Inf, or contains non-positive values
  */
 export declare function relSpread(x: number[]): number;
 /**
@@ -31,23 +41,41 @@ export declare function relSpread(x: number[]): number;
 export declare function shift(x: number[], y: number[]): number;
 /**
  * Calculate the Ratio - median of all pairwise ratios (x[i] / y[j])
+ *
+ * Assumptions:
+ *   positivity(x) - all values in x must be strictly positive
+ *   positivity(y) - all values in y must be strictly positive
+ *
  * @param x First sample
  * @param y Second sample
  * @returns The ratio estimate
+ * @throws AssumptionError if either sample is empty, contains NaN/Inf, or contains non-positive values
  */
 export declare function ratio(x: number[], y: number[]): number;
 /**
  * Calculate the AvgSpread - weighted average of spreads: (n*Spread(x) + m*Spread(y))/(n+m)
+ *
+ * Assumptions:
+ *   sparity(x) - first sample must be non tie-dominant (Spread > 0)
+ *   sparity(y) - second sample must be non tie-dominant (Spread > 0)
+ *
  * @param x First sample
  * @param y Second sample
  * @returns The combined spread estimate
+ * @throws AssumptionError if either sample is empty, contains NaN/Inf, or is tie-dominant
  */
 export declare function avgSpread(x: number[], y: number[]): number;
 /**
  * Calculate the Disparity - Shift / AvgSpread
+ *
+ * Assumptions:
+ *   sparity(x) - first sample must be non tie-dominant (Spread > 0)
+ *   sparity(y) - second sample must be non tie-dominant (Spread > 0)
+ *
  * @param x First sample
  * @param y Second sample
  * @returns The disparity estimate
+ * @throws AssumptionError if either sample is empty, contains NaN/Inf, or is tie-dominant
  */
 export declare function disparity(x: number[], y: number[]): number;
 /**
@@ -67,6 +95,7 @@ export interface Bounds {
  * @param y Second sample
  * @param misrate Misclassification rate (probability that true shift falls outside bounds)
  * @returns An object containing the lower and upper bounds
+ * @throws AssumptionError if either sample is empty or contains NaN/Inf
  */
 export declare function shiftBounds(x: number[], y: number[], misrate: number): Bounds;
 //# sourceMappingURL=estimators.d.ts.map

package/dist/estimators.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"estimators.d.ts","sourceRoot":"","sources":["../src/estimators.ts"],"names":[],"mappings":"AAAA;;GAEG;AAQH;;;;;GAKG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAK1C;AAED;;;;;GAKG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAK1C;AAED;;;;GAIG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAa7C;AAED;;;;;;GAMG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAStD;AAED;;;;;GAKG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAqBtD;AAED;;;;;GAKG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAa1D;AAED;;;;;GAKG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAa1D;AAED;;GAEG;AACH,MAAM,WAAW,MAAM;IACrB,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;CACf;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,CAoC7E"}
+{"version":3,"file":"estimators.d.ts","sourceRoot":"","sources":["../src/estimators.ts"],"names":[],"mappings":"AAAA;;GAEG;AASH;;;;;GAKG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAI1C;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAM1C;AAED;;;;;;;;;GASG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAW7C;AAED;;;;;;GAMG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAMtD;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAoBtD;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAkB1D;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,CAqB1D;AAED;;GAEG;AACH,MAAM,WAAW,MAAM;IACrB,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;CACf;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM,CAqC7E"}

package/dist/estimators.js

@@ -16,6 +16,7 @@ const fastCenter_1 = require("./fastCenter");
 const fastSpread_1 = require("./fastSpread");
 const fastShift_1 = require("./fastShift");
 const pairwiseMargin_1 = require("./pairwiseMargin");
+const assumptions_1 = require("./assumptions");
 /**
  * Calculate the Center - median of all pairwise averages (x[i] + x[j])/2
  * Uses fast O(n log n) algorithm.
@@ -23,37 +24,48 @@ const pairwiseMargin_1 = require("./pairwiseMargin");
  * @returns The center estimate
  */
 function center(x) {
-    if (x.length === 0) {
-        throw new Error('Input array cannot be empty');
-    }
+    // Check validity (priority 0)
+    (0, assumptions_1.checkValidity)(x, 'x', 'Center');
     return (0, fastCenter_1.fastCenter)(x);
 }
 /**
  * Calculate the Spread - median of all pairwise absolute differences |x[i] - x[j]|
  * Uses fast O(n log n) algorithm.
+ *
+ * Assumptions:
+ *   sparity(x) - sample must be non tie-dominant (Spread > 0)
+ *
  * @param x Array of sample values
  * @returns The spread estimate
+ * @throws AssumptionError if sample is empty, contains NaN/Inf, or is tie-dominant
  */
 function spread(x) {
-    if (x.length === 0) {
-        throw new Error('Input array cannot be empty');
-    }
+    // Check validity (priority 0)
+    (0, assumptions_1.checkValidity)(x, 'x', 'Spread');
+    // Check sparity (priority 2)
+    (0, assumptions_1.checkSparity)(x, 'x', 'Spread');
     return (0, fastSpread_1.fastSpread)(x);
 }
 /**
  * Calculate the RelSpread - ratio of Spread to absolute Center
+ *
+ * Assumptions:
+ *   positivity(x) - all values must be strictly positive (ensures Center > 0)
+ *
  * @param x Array of sample values
  * @returns The relative spread estimate
+ * @throws AssumptionError if sample is empty, contains NaN/Inf, or contains non-positive values
  */
 function relSpread(x) {
-    if (x.length === 0) {
-        throw new Error('Input array cannot be empty');
-    }
-    const s = spread(x);
-    const c = center(x);
-    if (c === 0) {
-        throw new Error('RelSpread is undefined when Center equals zero');
-    }
+    // Check validity (priority 0)
+    (0, assumptions_1.checkValidity)(x, 'x', 'RelSpread');
+    // Check positivity (priority 1)
+    (0, assumptions_1.checkPositivity)(x, 'x', 'RelSpread');
+    // Calculate center (we know x is valid, center should succeed)
+    const c = (0, fastCenter_1.fastCenter)(x);
+    // Calculate spread (using internal implementation since we already validated)
+    const s = (0, fastSpread_1.fastSpread)(x);
+    // center is guaranteed positive because all values are positive
     return s / Math.abs(c);
 }
 /**
@@ -64,29 +76,34 @@ function relSpread(x) {
  * @returns The shift estimate
  */
 function shift(x, y) {
-    const nx = x.length;
-    const ny = y.length;
-    if (nx === 0 || ny === 0) {
-        throw new Error('Input arrays cannot be empty');
-    }
+    // Check validity (priority 0)
+    (0, assumptions_1.checkValidity)(x, 'x', 'Shift');
+    (0, assumptions_1.checkValidity)(y, 'y', 'Shift');
     return (0, fastShift_1.fastShift)(x, y, [0.5], false)[0];
 }
 /**
  * Calculate the Ratio - median of all pairwise ratios (x[i] / y[j])
+ *
+ * Assumptions:
+ *   positivity(x) - all values in x must be strictly positive
+ *   positivity(y) - all values in y must be strictly positive
+ *
  * @param x First sample
  * @param y Second sample
  * @returns The ratio estimate
+ * @throws AssumptionError if either sample is empty, contains NaN/Inf, or contains non-positive values
  */
 function ratio(x, y) {
+    // Check validity for x (priority 0, subject x)
+    (0, assumptions_1.checkValidity)(x, 'x', 'Ratio');
+    // Check validity for y (priority 0, subject y)
+    (0, assumptions_1.checkValidity)(y, 'y', 'Ratio');
+    // Check positivity for x (priority 1, subject x)
+    (0, assumptions_1.checkPositivity)(x, 'x', 'Ratio');
+    // Check positivity for y (priority 1, subject y)
+    (0, assumptions_1.checkPositivity)(y, 'y', 'Ratio');
     const nx = x.length;
     const ny = y.length;
-    if (nx === 0 || ny === 0) {
-        throw new Error('Input arrays cannot be empty');
-    }
-    // Check that all y values are strictly positive
-    if (y.some((val) => val <= 0)) {
-        throw new Error('All values in y must be strictly positive');
-    }
     const pairwiseRatios = [];
     for (let i = 0; i < nx; i++) {
         for (let j = 0; j < ny; j++) {
@@ -97,37 +114,62 @@ function ratio(x, y) {
 }
 /**
  * Calculate the AvgSpread - weighted average of spreads: (n*Spread(x) + m*Spread(y))/(n+m)
+ *
+ * Assumptions:
+ *   sparity(x) - first sample must be non tie-dominant (Spread > 0)
+ *   sparity(y) - second sample must be non tie-dominant (Spread > 0)
+ *
  * @param x First sample
  * @param y Second sample
  * @returns The combined spread estimate
+ * @throws AssumptionError if either sample is empty, contains NaN/Inf, or is tie-dominant
  */
 function avgSpread(x, y) {
+    // Check validity for x (priority 0, subject x)
+    (0, assumptions_1.checkValidity)(x, 'x', 'AvgSpread');
+    // Check validity for y (priority 0, subject y)
+    (0, assumptions_1.checkValidity)(y, 'y', 'AvgSpread');
+    // Check sparity for x (priority 2, subject x)
+    (0, assumptions_1.checkSparity)(x, 'x', 'AvgSpread');
+    // Check sparity for y (priority 2, subject y)
+    (0, assumptions_1.checkSparity)(y, 'y', 'AvgSpread');
     const nx = x.length;
     const ny = y.length;
-    if (nx === 0 || ny === 0) {
-        throw new Error('Input arrays cannot be empty');
-    }
-    // Calculate weighted average of individual spreads
-    const spreadX = spread(x);
-    const spreadY = spread(y);
+    // Calculate spreads (using internal implementation since we already validated)
+    const spreadX = (0, fastSpread_1.fastSpread)(x);
+    const spreadY = (0, fastSpread_1.fastSpread)(y);
     return (nx * spreadX + ny * spreadY) / (nx + ny);
 }
 /**
  * Calculate the Disparity - Shift / AvgSpread
+ *
+ * Assumptions:
+ *   sparity(x) - first sample must be non tie-dominant (Spread > 0)
+ *   sparity(y) - second sample must be non tie-dominant (Spread > 0)
+ *
  * @param x First sample
  * @param y Second sample
  * @returns The disparity estimate
+ * @throws AssumptionError if either sample is empty, contains NaN/Inf, or is tie-dominant
  */
 function disparity(x, y) {
-    if (x.length === 0 || y.length === 0) {
-        throw new Error('Input arrays cannot be empty');
-    }
-    const shiftVal = shift(x, y);
-    const combinedSpread = avgSpread(x, y);
-    if (combinedSpread === 0) {
-        return Infinity;
-    }
-    return shiftVal / combinedSpread;
+    // Check validity for x (priority 0, subject x)
+    (0, assumptions_1.checkValidity)(x, 'x', 'Disparity');
+    // Check validity for y (priority 0, subject y)
+    (0, assumptions_1.checkValidity)(y, 'y', 'Disparity');
+    // Check sparity for x (priority 2, subject x)
+    (0, assumptions_1.checkSparity)(x, 'x', 'Disparity');
+    // Check sparity for y (priority 2, subject y)
+    (0, assumptions_1.checkSparity)(y, 'y', 'Disparity');
+    const nx = x.length;
+    const ny = y.length;
+    // Calculate shift (we know inputs are valid)
+    const shiftVal = (0, fastShift_1.fastShift)(x, y, [0.5], false)[0];
+    // Calculate avg_spread (using internal implementation since we already validated)
+    const spreadX = (0, fastSpread_1.fastSpread)(x);
+    const spreadY = (0, fastSpread_1.fastSpread)(y);
+    const avgSpreadVal = (nx * spreadX + ny * spreadY) / (nx + ny);
+    return shiftVal / avgSpreadVal;
 }
 /**
  * Provides bounds on the Shift estimator with specified misclassification rate (ShiftBounds)
@@ -139,11 +181,13 @@ function disparity(x, y) {
  * @param y Second sample
  * @param misrate Misclassification rate (probability that true shift falls outside bounds)
  * @returns An object containing the lower and upper bounds
+ * @throws AssumptionError if either sample is empty or contains NaN/Inf
  */
 function shiftBounds(x, y, misrate) {
-    if (x.length === 0 || y.length === 0) {
-        throw new Error('Input arrays cannot be empty');
-    }
+    // Check validity for x
+    (0, assumptions_1.checkValidity)(x, 'x', 'ShiftBounds');
+    // Check validity for y
+    (0, assumptions_1.checkValidity)(y, 'y', 'ShiftBounds');
     const n = x.length;
     const m = y.length;
     // Sort both arrays

package/dist/estimators.js.map

@@ -1 +1 @@
-{"version":3,"file":"estimators.js","sourceRoot":"","sources":["../src/estimators.ts"],"names":[],"mappings":";AAAA;;GAEG;;AAcH,wBAKC;AAQD,wBAKC;AAOD,8BAaC;AASD,sBASC;AAQD,sBAqBC;AAQD,8BAaC;AAQD,8BAaC;AAqBD,kCAoCC;AApMD,mCAAiC;AACjC,6CAA0C;AAC1C,6CAA0C;AAC1C,2CAAwC;AACxC,qDAAkD;AAElD;;;;;GAKG;AACH,SAAgB,MAAM,CAAC,CAAW;IAChC,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,CAAC;IACjD,CAAC;IACD,OAAO,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;AACvB,CAAC;AAED;;;;;GAKG;AACH,SAAgB,MAAM,CAAC,CAAW;IAChC,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,CAAC;IACjD,CAAC;IACD,OAAO,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;AACvB,CAAC;AAED;;;;GAIG;AACH,SAAgB,SAAS,CAAC,CAAW;IACnC,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,CAAC;IACjD,CAAC;IAED,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IACpB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAEpB,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,MAAM,IAAI,KAAK,CAAC,gDAAgD,CAAC,CAAC;IACpE,CAAC;IAED,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACzB,CAAC;AAED;;;;;;GAMG;AACH,SAAgB,KAAK,CAAC,CAAW,EAAE,CAAW;IAC5C,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IAEpB,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QACzB,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IAClD,CAAC;IAED,OAAO,IAAA,qBAAS,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1C,CAAC;AAED;;;;;GAKG;AACH,SAAgB,KAAK,CAAC,CAAW,EAAE,CAAW;IAC5C,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IAEpB,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QACzB,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IAClD,CAAC;IAED,gDAAgD;IAChD,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC;QAC9B,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;IAC/D,CAAC;IAED,MAAM,cAAc,GAAa,EAAE,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,cAAc,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACnC,CAAC;IACH,CAAC;IAED,OAAO,IAAA,cAAM,EAAC,cAAc,CAAC,CAAC;AAChC,CAAC;AAED;;;;;GAKG;AACH,SAAgB,SAAS,CAAC,CAAW,EAAE,CAAW;IAChD,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IAEpB,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QACzB,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IAClD,CAAC;IAED,mDAAmD;IACnD,MAAM,OAAO,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAC1B,MAAM,OAAO,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAE1B,OAAO,CAAC,EAAE,GAAG,OAAO,GAAG,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;AACnD,CAAC;AAED;;;;;GAKG;AACH,SAAgB,SAAS,CAAC,CAAW,EAAE,CAAW;IAChD,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACrC,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IAClD,CAAC;IAED,MAAM,QAAQ,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7B,MAAM,cAAc,GAAG,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAEvC,IAAI,cAAc,KAAK,CAAC,EAAE,CAAC;QACzB,OAAO,QAAQ,CAAC;IAClB,CAAC;IAED,OAAO,QAAQ,GAAG,cAAc,CAAC;AACnC,CAAC;AAUD;;;;;;;;;;GAUG;AACH,SAAgB,WAAW,CAAC,CAAW,EAAE,CAAW,EAAE,OAAe;IACnE,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACrC,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IAClD,CAAC;IAED,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IAEnB,mBAAmB;IACnB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACxC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAExC,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;IAEpB,6FAA6F;IAC7F,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAC5B,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IACxC,CAAC;IAED,MAAM,MAAM,GAAG,IAAA,+BAAc,EAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;IAC7C,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACjF,MAAM,KAAK,GAAG,UAAU,CAAC;IACzB,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,GAAG,UAAU,CAAC;IAEtC,6BAA6B;IAC7B,MAAM,WAAW,GAAG,KAAK,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,KAAK,GAAG,WAAW,CAAC;IAClC,MAAM,MAAM,GAAG,MAAM,GAAG,WAAW,CAAC;IAEpC,6DAA6D;IAC7D,MAAM,CAAC,IAAI,EAAE,KAAK,CAAC,GAAG,IAAA,qBAAS,EAAC,EAAE,EAAE,EAAE,EAAE,CAAC,KAAK,EAAE,MAAM,CAAC,EAAE,IAAI,CAAC,CAAC;IAC/D,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IAEpC,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;AAC1B,CAAC"}
+{"version":3,"file":"estimators.js","sourceRoot":"","sources":["../src/estimators.ts"],"names":[],"mappings":";AAAA;;GAEG;;AAeH,wBAIC;AAaD,wBAMC;AAYD,8BAWC;AASD,sBAMC;AAcD,sBAoBC;AAcD,8BAkBC;AAcD,8BAqBC;AAsBD,kCAqCC;AA1OD,mCAAiC;AACjC,6CAA0C;AAC1C,6CAA0C;AAC1C,2CAAwC;AACxC,qDAAkD;AAClD,+CAA6E;AAE7E;;;;;GAKG;AACH,SAAgB,MAAM,CAAC,CAAW;IAChC,8BAA8B;IAC9B,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IAChC,OAAO,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;AACvB,CAAC;AAED;;;;;;;;;;GAUG;AACH,SAAgB,MAAM,CAAC,CAAW;IAChC,8BAA8B;IAC9B,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IAChC,6BAA6B;IAC7B,IAAA,0BAAY,EAAC,CAAC,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IAC/B,OAAO,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;AACvB,CAAC;AAED;;;;;;;;;GASG;AACH,SAAgB,SAAS,CAAC,CAAW;IACnC,8BAA8B;IAC9B,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IACnC,gCAAgC;IAChC,IAAA,6BAAe,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IACrC,+DAA+D;IAC/D,MAAM,CAAC,GAAG,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;IACxB,8EAA8E;IAC9E,MAAM,CAAC,GAAG,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;IACxB,gEAAgE;IAChE,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACzB,CAAC;AAED;;;;;;GAMG;AACH,SAAgB,KAAK,CAAC,CAAW,EAAE,CAAW;IAC5C,8BAA8B;IAC9B,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;IAC/B,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;IAE/B,OAAO,IAAA,qBAAS,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1C,CAAC;AAED;;;;;;;;;;;GAWG;AACH,SAAgB,KAAK,CAAC,CAAW,EAAE,CAAW;IAC5C,+CAA+C;IAC/C,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;IAC/B,+CAA+C;IAC/C,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;IAC/B,iDAAiD;IACjD,IAAA,6BAAe,EAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;IACjC,iDAAiD;IACjD,IAAA,6BAAe,EAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;IAEjC,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IACpB,MAAM,cAAc,GAAa,EAAE,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,cAAc,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACnC,CAAC;IACH,CAAC;IAED,OAAO,IAAA,cAAM,EAAC,cAAc,CAAC,CAAC;AAChC,CAAC;AAED;;;;;;;;;;;GAWG;AACH,SAAgB,SAAS,CAAC,CAAW,EAAE,CAAW;IAChD,+CAA+C;IAC/C,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IACnC,+CAA+C;IAC/C,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IACnC,8CAA8C;IAC9C,IAAA,0BAAY,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IAClC,8CAA8C;IAC9C,IAAA,0BAAY,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IAElC,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IAEpB,+EAA+E;IAC/E,MAAM,OAAO,GAAG,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,OAAO,GAAG,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;IAE9B,OAAO,CAAC,EAAE,GAAG,OAAO,GAAG,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;AACnD,CAAC;AAED;;;;;;;;;;;GAWG;AACH,SAAgB,SAAS,CAAC,CAAW,EAAE,CAAW;IAChD,+CAA+C;IAC/C,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IACnC,+CAA+C;IAC/C,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IACnC,8CAA8C;IAC9C,IAAA,0BAAY,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IAClC,8CAA8C;IAC9C,IAAA,0BAAY,EAAC,CAAC,EAAE,GAAG,EAAE,WAAW,CAAC,CAAC;IAElC,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IAEpB,6CAA6C;IAC7C,MAAM,QAAQ,GAAG,IAAA,qBAAS,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAClD,kFAAkF;IAClF,MAAM,OAAO,GAAG,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,OAAO,GAAG,IAAA,uBAAU,EAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,YAAY,GAAG,CAAC,EAAE,GAAG,OAAO,GAAG,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IAE/D,OAAO,QAAQ,GAAG,YAAY,CAAC;AACjC,CAAC;AAUD;;;;;;;;;;;GAWG;AACH,SAAgB,WAAW,CAAC,CAAW,EAAE,CAAW,EAAE,OAAe;IACnE,uBAAuB;IACvB,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,aAAa,CAAC,CAAC;IACrC,uBAAuB;IACvB,IAAA,2BAAa,EAAC,CAAC,EAAE,GAAG,EAAE,aAAa,CAAC,CAAC;IAErC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IAEnB,mBAAmB;IACnB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACxC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAExC,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;IAEpB,6FAA6F;IAC7F,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAC5B,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IACxC,CAAC;IAED,MAAM,MAAM,GAAG,IAAA,+BAAc,EAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;IAC7C,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACjF,MAAM,KAAK,GAAG,UAAU,CAAC;IACzB,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,GAAG,UAAU,CAAC;IAEtC,6BAA6B;IAC7B,MAAM,WAAW,GAAG,KAAK,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,KAAK,GAAG,WAAW,CAAC;IAClC,MAAM,MAAM,GAAG,MAAM,GAAG,WAAW,CAAC;IAEpC,6DAA6D;IAC7D,MAAM,CAAC,IAAI,EAAE,KAAK,CAAC,GAAG,IAAA,qBAAS,EAAC,EAAE,EAAE,EAAE,EAAE,CAAC,KAAK,EAAE,MAAM,CAAC,EAAE,IAAI,CAAC,CAAC;IAC/D,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IAEpC,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC;AAC1B,CAAC"}

package/dist/index.d.ts

@@ -3,6 +3,7 @@
  *
  * A collection of robust statistical estimators for real-world data analysis.
  */
+export { AssumptionId, AssumptionError, type Subject, type Violation } from './assumptions';
 export { center, spread, relSpread, shift, ratio, avgSpread, disparity, shiftBounds, type Bounds, } from './estimators';
 export { median } from './utils';
 export { pairwiseMargin } from './pairwiseMargin';

package/dist/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,OAAO,EACL,MAAM,EACN,MAAM,EACN,SAAS,EACT,KAAK,EACL,KAAK,EACL,SAAS,EACT,SAAS,EACT,WAAW,EACX,KAAK,MAAM,GACZ,MAAM,cAAc,CAAC;AAEtB,OAAO,EAAE,MAAM,EAAE,MAAM,SAAS,CAAC;AAEjC,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAElD,OAAO,EAAE,GAAG,EAAE,MAAM,OAAO,CAAC;AAE5B,OAAO,EAAE,KAAK,YAAY,EAAE,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,GAAG,EAAE,KAAK,EAAE,MAAM,uBAAuB,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAGH,OAAO,EAAE,YAAY,EAAE,eAAe,EAAE,KAAK,OAAO,EAAE,KAAK,SAAS,EAAE,MAAM,eAAe,CAAC;AAG5F,OAAO,EACL,MAAM,EACN,MAAM,EACN,SAAS,EACT,KAAK,EACL,KAAK,EACL,SAAS,EACT,SAAS,EACT,WAAW,EACX,KAAK,MAAM,GACZ,MAAM,cAAc,CAAC;AAEtB,OAAO,EAAE,MAAM,EAAE,MAAM,SAAS,CAAC;AAEjC,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAElD,OAAO,EAAE,GAAG,EAAE,MAAM,OAAO,CAAC;AAE5B,OAAO,EAAE,KAAK,YAAY,EAAE,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,GAAG,EAAE,KAAK,EAAE,MAAM,uBAAuB,CAAC"}

package/dist/index.js

@@ -5,7 +5,12 @@
  * A collection of robust statistical estimators for real-world data analysis.
  */
 Object.defineProperty(exports, "__esModule", { value: true });
-exports.Power = exports.Exp = exports.Multiplic = exports.Additive = exports.Uniform = exports.Rng = exports.pairwiseMargin = exports.median = exports.shiftBounds = exports.disparity = exports.avgSpread = exports.ratio = exports.shift = exports.relSpread = exports.spread = exports.center = void 0;
+exports.Power = exports.Exp = exports.Multiplic = exports.Additive = exports.Uniform = exports.Rng = exports.pairwiseMargin = exports.median = exports.shiftBounds = exports.disparity = exports.avgSpread = exports.ratio = exports.shift = exports.relSpread = exports.spread = exports.center = exports.AssumptionError = exports.AssumptionId = void 0;
+// Assumptions
+var assumptions_1 = require("./assumptions");
+Object.defineProperty(exports, "AssumptionId", { enumerable: true, get: function () { return assumptions_1.AssumptionId; } });
+Object.defineProperty(exports, "AssumptionError", { enumerable: true, get: function () { return assumptions_1.AssumptionError; } });
+// Estimators
 var estimators_1 = require("./estimators");
 Object.defineProperty(exports, "center", { enumerable: true, get: function () { return estimators_1.center; } });
 Object.defineProperty(exports, "spread", { enumerable: true, get: function () { return estimators_1.spread; } });

package/dist/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":";AAAA;;;;GAIG;;;AAEH,2CAUsB;AATpB,oGAAA,MAAM,OAAA;AACN,oGAAA,MAAM,OAAA;AACN,uGAAA,SAAS,OAAA;AACT,mGAAA,KAAK,OAAA;AACL,mGAAA,KAAK,OAAA;AACL,uGAAA,SAAS,OAAA;AACT,uGAAA,SAAS,OAAA;AACT,yGAAA,WAAW,OAAA;AAIb,iCAAiC;AAAxB,+FAAA,MAAM,OAAA;AAEf,mDAAkD;AAAzC,gHAAA,cAAc,OAAA;AAEvB,6BAA4B;AAAnB,0FAAA,GAAG,OAAA;AAEZ,+CAAoG;AAAxE,gGAAA,OAAO,OAAA;AAAE,iGAAA,QAAQ,OAAA;AAAE,kGAAA,SAAS,OAAA;AAAE,4FAAA,GAAG,OAAA;AAAE,8FAAA,KAAK,OAAA"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":";AAAA;;;;GAIG;;;AAEH,cAAc;AACd,6CAA4F;AAAnF,2GAAA,YAAY,OAAA;AAAE,8GAAA,eAAe,OAAA;AAEtC,aAAa;AACb,2CAUsB;AATpB,oGAAA,MAAM,OAAA;AACN,oGAAA,MAAM,OAAA;AACN,uGAAA,SAAS,OAAA;AACT,mGAAA,KAAK,OAAA;AACL,mGAAA,KAAK,OAAA;AACL,uGAAA,SAAS,OAAA;AACT,uGAAA,SAAS,OAAA;AACT,yGAAA,WAAW,OAAA;AAIb,iCAAiC;AAAxB,+FAAA,MAAM,OAAA;AAEf,mDAAkD;AAAzC,gHAAA,cAAc,OAAA;AAEvB,6BAA4B;AAAnB,0FAAA,GAAG,OAAA;AAEZ,+CAAoG;AAAxE,gGAAA,OAAO,OAAA;AAAE,iGAAA,QAAQ,OAAA;AAAE,kGAAA,SAAS,OAAA;AAAE,4FAAA,GAAG,OAAA;AAAE,8FAAA,KAAK,OAAA"}

package/dist/rng.d.ts

@@ -47,6 +47,14 @@ export declare class Rng {
      * @returns Random value in [0, 1).
      */
     uniform(): number;
+    /**
+     * Generate a uniform random float in [min, max).
+     *
+     * @param min - Minimum value (inclusive).
+     * @param max - Maximum value (exclusive).
+     * @returns Random value in [min, max). Returns min if min >= max.
+     */
+    uniformRange(min: number, max: number): number;
     /**
      * Generate a uniform random integer in [min, max).
      *
@@ -59,6 +67,24 @@ export declare class Rng {
      * @returns Random integer in [min, max). Returns min if min >= max.
      */
     uniformInt(min: number, max: number): number;
+    /**
+     * Generate a uniform random BigInt in [min, max).
+     *
+     * Uses modulo reduction which introduces slight bias for ranges that don't
+     * evenly divide 2^64. This bias is negligible for statistical simulations
+     * but not suitable for cryptographic applications.
+     *
+     * @param min - Minimum value (inclusive).
+     * @param max - Maximum value (exclusive).
+     * @returns Random BigInt in [min, max). Returns min if min >= max.
+     */
+    uniformBigInt(min: bigint, max: bigint): bigint;
+    /**
+     * Generate a uniform random boolean with P(true) = 0.5.
+     *
+     * @returns Random boolean value.
+     */
+    uniformBool(): boolean;
     /**
      * Return a shuffled copy of the input array.
      *

package/dist/rng.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"rng.d.ts","sourceRoot":"","sources":["../src/rng.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAIH;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,qBAAa,GAAG;IACd,OAAO,CAAC,KAAK,CAAqB;IAElC;;;;;;OAMG;gBACS,IAAI,CAAC,EAAE,MAAM,GAAG,MAAM;IAelC;;;;;;OAMG;IACH,OAAO,IAAI,MAAM;IAIjB;;;;;;;;;;OAUG;IACH,UAAU,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;IAK5C;;;;;;;;OAQG;IACH,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,EAAE;IAavB;;;;;;;;;;OAUG;IACH,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE;CAuBlC"}
+{"version":3,"file":"rng.d.ts","sourceRoot":"","sources":["../src/rng.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAIH;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,qBAAa,GAAG;IACd,OAAO,CAAC,KAAK,CAAqB;IAElC;;;;;;OAMG;gBACS,IAAI,CAAC,EAAE,MAAM,GAAG,MAAM;IAmBlC;;;;;;OAMG;IACH,OAAO,IAAI,MAAM;IAIjB;;;;;;OAMG;IACH,YAAY,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;IAQ9C;;;;;;;;;;OAUG;IACH,UAAU,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;IAK5C;;;;;;;;;;OAUG;IACH,aAAa,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;IAQ/C;;;;OAIG;IACH,WAAW,IAAI,OAAO;IAItB;;;;;;;;OAQG;IACH,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,EAAE;IAavB;;;;;;;;;;OAUG;IACH,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE;CAuBlC"}

package/dist/rng.js

@@ -55,6 +55,9 @@ class Rng {
         }
         this.inner = new xoshiro256_1.Xoshiro256PlusPlus(seedBigInt);
     }
+    // ========================================================================
+    // Floating Point Methods
+    // ========================================================================
     /**
      * Generate a uniform random float in [0, 1).
      *
@@ -65,6 +68,19 @@ class Rng {
     uniform() {
         return this.inner.uniform();
     }
+    /**
+     * Generate a uniform random float in [min, max).
+     *
+     * @param min - Minimum value (inclusive).
+     * @param max - Maximum value (exclusive).
+     * @returns Random value in [min, max). Returns min if min >= max.
+     */
+    uniformRange(min, max) {
+        return this.inner.uniformRange(min, max);
+    }
+    // ========================================================================
+    // Integer Methods
+    // ========================================================================
     /**
      * Generate a uniform random integer in [min, max).
      *
@@ -80,6 +96,31 @@ class Rng {
         const result = this.inner.uniformInt(BigInt(min), BigInt(max));
         return Number(result);
     }
+    /**
+     * Generate a uniform random BigInt in [min, max).
+     *
+     * Uses modulo reduction which introduces slight bias for ranges that don't
+     * evenly divide 2^64. This bias is negligible for statistical simulations
+     * but not suitable for cryptographic applications.
+     *
+     * @param min - Minimum value (inclusive).
+     * @param max - Maximum value (exclusive).
+     * @returns Random BigInt in [min, max). Returns min if min >= max.
+     */
+    uniformBigInt(min, max) {
+        return this.inner.uniformInt(min, max);
+    }
+    // ========================================================================
+    // Boolean Methods
+    // ========================================================================
+    /**
+     * Generate a uniform random boolean with P(true) = 0.5.
+     *
+     * @returns Random boolean value.
+     */
+    uniformBool() {
+        return this.inner.uniformBool();
+    }
     /**
      * Return a shuffled copy of the input array.
      *

package/dist/rng.js.map

@@ -1 +1 @@
-{"version":3,"file":"rng.js","sourceRoot":"","sources":["../src/rng.ts"],"names":[],"mappings":";AAAA;;;;;GAKG;;;AAEH,6CAA6D;AAE7D;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,MAAa,GAAG;IAGd;;;;;;OAMG;IACH,YAAY,IAAsB;QAChC,IAAI,UAAkB,CAAC;QAEvB,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,UAAU,GAAG,MAAM,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC;QAClC,CAAC;aAAM,IAAI,OAAO,IAAI,KAAK,QAAQ,EAAE,CAAC;YACpC,UAAU,GAAG,IAAA,sBAAS,EAAC,IAAI,CAAC,CAAC;QAC/B,CAAC;aAAM,CAAC;YACN,2EAA2E;YAC3E,UAAU,GAAG,MAAM,CAAC,MAAM,CAAC,EAAE,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC;QAC/C,CAAC;QAED,IAAI,CAAC,KAAK,GAAG,IAAI,+BAAkB,CAAC,UAAU,CAAC,CAAC;IAClD,CAAC;IAED;;;;;;OAMG;IACH,OAAO;QACL,OAAO,IAAI,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC;IAC9B,CAAC;IAED;;;;;;;;;;OAUG;IACH,UAAU,CAAC,GAAW,EAAE,GAAW;QACjC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/D,OAAO,MAAM,CAAC,MAAM,CAAC,CAAC;IACxB,CAAC;IAED;;;;;;;;OAQG;IACH,OAAO,CAAI,CAAM;QACf,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QAExB,mCAAmC;QACnC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACpC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,CAAC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;;;;;;;;;OAUG;IACH,MAAM,CAAI,CAAM,EAAE,CAAS;QACzB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;QACpD,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;YACX,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC;QAChB,CAAC;QAED,MAAM,MAAM,GAAQ,EAAE,CAAC;QACvB,IAAI,SAAS,GAAG,CAAC,CAAC;QAElB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,MAAM,SAAS,GAAG,CAAC,GAAG,CAAC,CAAC;YACxB,4DAA4D;YAC5D,IAAI,IAAI,CAAC,OAAO,EAAE,GAAG,SAAS,GAAG,SAAS,EAAE,CAAC;gBAC3C,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClB,SAAS,EAAE,CAAC;YACd,CAAC;QACH,CAAC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AA5GD,kBA4GC"}
+{"version":3,"file":"rng.js","sourceRoot":"","sources":["../src/rng.ts"],"names":[],"mappings":";AAAA;;;;;GAKG;;;AAEH,6CAA6D;AAE7D;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,MAAa,GAAG;IAGd;;;;;;OAMG;IACH,YAAY,IAAsB;QAChC,IAAI,UAAkB,CAAC;QAEvB,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,UAAU,GAAG,MAAM,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC;QAClC,CAAC;aAAM,IAAI,OAAO,IAAI,KAAK,QAAQ,EAAE,CAAC;YACpC,UAAU,GAAG,IAAA,sBAAS,EAAC,IAAI,CAAC,CAAC;QAC/B,CAAC;aAAM,CAAC;YACN,2EAA2E;YAC3E,UAAU,GAAG,MAAM,CAAC,MAAM,CAAC,EAAE,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC;QAC/C,CAAC;QAED,IAAI,CAAC,KAAK,GAAG,IAAI,+BAAkB,CAAC,UAAU,CAAC,CAAC;IAClD,CAAC;IAED,2EAA2E;IAC3E,yBAAyB;IACzB,2EAA2E;IAE3E;;;;;;OAMG;IACH,OAAO;QACL,OAAO,IAAI,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC;IAC9B,CAAC;IAED;;;;;;OAMG;IACH,YAAY,CAAC,GAAW,EAAE,GAAW;QACnC,OAAO,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAC3C,CAAC;IAED,2EAA2E;IAC3E,kBAAkB;IAClB,2EAA2E;IAE3E;;;;;;;;;;OAUG;IACH,UAAU,CAAC,GAAW,EAAE,GAAW;QACjC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/D,OAAO,MAAM,CAAC,MAAM,CAAC,CAAC;IACxB,CAAC;IAED;;;;;;;;;;OAUG;IACH,aAAa,CAAC,GAAW,EAAE,GAAW;QACpC,OAAO,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACzC,CAAC;IAED,2EAA2E;IAC3E,kBAAkB;IAClB,2EAA2E;IAE3E;;;;OAIG;IACH,WAAW;QACT,OAAO,IAAI,CAAC,KAAK,CAAC,WAAW,EAAE,CAAC;IAClC,CAAC;IAED;;;;;;;;OAQG;IACH,OAAO,CAAI,CAAM;QACf,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QAExB,mCAAmC;QACnC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACpC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,CAAC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;;;;;;;;;OAUG;IACH,MAAM,CAAI,CAAM,EAAE,CAAS;QACzB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;QACpD,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;YACX,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC;QAChB,CAAC;QAED,MAAM,MAAM,GAAQ,EAAE,CAAC;QACvB,IAAI,SAAS,GAAG,CAAC,CAAC;QAElB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,MAAM,SAAS,GAAG,CAAC,GAAG,CAAC,CAAC;YACxB,4DAA4D;YAC5D,IAAI,IAAI,CAAC,OAAO,EAAE,GAAG,SAAS,GAAG,SAAS,EAAE,CAAC;gBAC3C,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClB,SAAS,EAAE,CAAC;YACd,CAAC;QACH,CAAC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AA3JD,kBA2JC"}

package/dist/utils.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"utils.d.ts","sourceRoot":"","sources":["../src/utils.ts"],"names":[],"mappings":"AAAA;;GAEG;AAEH;;;;GAIG;AACH,wBAAgB,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAa/C"}
+{"version":3,"file":"utils.d.ts","sourceRoot":"","sources":["../src/utils.ts"],"names":[],"mappings":"AAAA;;GAEG;AAIH;;;;GAIG;AACH,wBAAgB,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAY/C"}

package/dist/utils.js

@@ -4,15 +4,15 @@
  */
 Object.defineProperty(exports, "__esModule", { value: true });
 exports.median = median;
+const assumptions_1 = require("./assumptions");
 /**
  * Calculate the median of an array of numbers
  * @param values Array of numbers
  * @returns The median value
  */
 function median(values) {
-    if (values.length === 0) {
-        throw new Error('Input array cannot be empty');
-    }
+    // Check validity (priority 0)
+    (0, assumptions_1.checkValidity)(values, 'x', 'Median');
     const sorted = [...values].sort((a, b) => a - b);
     const mid = Math.floor(sorted.length / 2);
     if (sorted.length % 2 === 0) {

package/dist/utils.js.map

@@ -1 +1 @@
-{"version":3,"file":"utils.js","sourceRoot":"","sources":["../src/utils.ts"],"names":[],"mappings":";AAAA;;GAEG;;AAOH,wBAaC;AAlBD;;;;GAIG;AACH,SAAgB,MAAM,CAAC,MAAgB;IACrC,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,CAAC;IACjD,CAAC;IAED,MAAM,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACjD,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAE1C,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC;QAC5B,OAAO,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IAC7C,CAAC;SAAM,CAAC;QACN,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC;IACrB,CAAC;AACH,CAAC"}
+{"version":3,"file":"utils.js","sourceRoot":"","sources":["../src/utils.ts"],"names":[],"mappings":";AAAA;;GAEG;;AASH,wBAYC;AAnBD,+CAA8C;AAE9C;;;;GAIG;AACH,SAAgB,MAAM,CAAC,MAAgB;IACrC,8BAA8B;IAC9B,IAAA,2BAAa,EAAC,MAAM,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC;IAErC,MAAM,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACjD,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAE1C,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC;QAC5B,OAAO,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IAC7C,CAAC;SAAM,CAAC;QACN,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC;IACrB,CAAC;AACH,CAAC"}

package/dist/xoshiro256.d.ts

@@ -16,11 +16,13 @@ export declare class Xoshiro256PlusPlus {
     constructor(seed: bigint);
     nextU64(): bigint;
     uniform(): number;
+    uniformRange(min: number, max: number): number;
     /**
      * Generate a uniform integer in [min, max).
      * @throws RangeError if max - min exceeds i64 range.
      */
     uniformInt(min: bigint, max: bigint): bigint;
+    uniformBool(): boolean;
 }
 /**
  * Compute FNV-1a 64-bit hash of a string

package/dist/xoshiro256.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"xoshiro256.d.ts","sourceRoot":"","sources":["../src/xoshiro256.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AA2BH;;;;;GAKG;AACH,qBAAa,kBAAkB;IAC7B,OAAO,CAAC,KAAK,CAAmC;gBAEpC,IAAI,EAAE,MAAM;IAKxB,OAAO,IAAI,MAAM;IAiBjB,OAAO,IAAI,MAAM;IAMjB;;;OAGG;IACH,UAAU,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;CAW7C;AAMD;;GAEG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAW3C"}
+{"version":3,"file":"xoshiro256.d.ts","sourceRoot":"","sources":["../src/xoshiro256.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AA2BH;;;;;GAKG;AACH,qBAAa,kBAAkB;IAC7B,OAAO,CAAC,KAAK,CAAmC;gBAEpC,IAAI,EAAE,MAAM;IAKxB,OAAO,IAAI,MAAM;IAqBjB,OAAO,IAAI,MAAM;IAMjB,YAAY,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;IAS9C;;;OAGG;IACH,UAAU,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;IAgB5C,WAAW,IAAI,OAAO;CAGvB;AAMD;;GAEG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAW3C"}

package/dist/xoshiro256.js

@@ -51,11 +51,22 @@ class Xoshiro256PlusPlus {
         s[3] = rotl(s[3], 45n);
         return result;
     }
+    // ========================================================================
+    // Floating Point Methods
+    // ========================================================================
     uniform() {
         // Use upper 53 bits for maximum precision in float64
         const u64 = this.nextU64();
         return Number(u64 >> 11n) * (1.0 / Number(1n << 53n));
     }
+    uniformRange(min, max) {
+        if (min >= max)
+            return min;
+        return min + (max - min) * this.uniform();
+    }
+    // ========================================================================
+    // Integer Methods
+    // ========================================================================
     /**
      * Generate a uniform integer in [min, max).
      * @throws RangeError if max - min exceeds i64 range.
@@ -71,6 +82,12 @@ class Xoshiro256PlusPlus {
         }
         return min + (this.nextU64() % range);
     }
+    // ========================================================================
+    // Boolean Methods
+    // ========================================================================
+    uniformBool() {
+        return this.uniform() < 0.5;
+    }
 }
 exports.Xoshiro256PlusPlus = Xoshiro256PlusPlus;
 // FNV-1a hash constants

package/dist/xoshiro256.js.map

@@ -1 +1 @@
-{"version":3,"file":"xoshiro256.js","sourceRoot":"","sources":["../src/xoshiro256.ts"],"names":[],"mappings":";AAAA;;;;;;GAMG;;;AAwFH,8BAWC;AAjGD,MAAM,OAAO,GAAG,mBAAmB,CAAC;AAEpC,SAAS,IAAI,CAAC,CAAS,EAAE,CAAS;IAChC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;AACjD,CAAC;AAED;;GAEG;AACH,MAAM,UAAU;IAGd,YAAY,IAAY;QACtB,IAAI,CAAC,KAAK,GAAG,IAAI,GAAG,OAAO,CAAC;IAC9B,CAAC;IAED,IAAI;QACF,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,mBAAmB,CAAC,GAAG,OAAO,CAAC;QAC1D,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACnB,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,mBAAmB,CAAC,GAAG,OAAO,CAAC;QACvD,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,mBAAmB,CAAC,GAAG,OAAO,CAAC;QACvD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC;IACpC,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAa,kBAAkB;IAG7B,YAAY,IAAY;QACtB,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC;QAChC,IAAI,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC;IAC5D,CAAC;IAED,OAAO;QACL,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACrB,MAAM,MAAM,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;QAErE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,OAAO,CAAC;QAElC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAEb,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACV,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAEvB,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,OAAO;QACL,qDAAqD;QACrD,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,EAAE,CAAC;QAC3B,OAAO,MAAM,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,EAAE,IAAI,GAAG,CAAC,CAAC,CAAC;IACxD,CAAC;IAED;;;OAGG;IACH,UAAU,CAAC,GAAW,EAAE,GAAW;QACjC,IAAI,GAAG,IAAI,GAAG,EAAE,CAAC;YACf,OAAO,GAAG,CAAC;QACb,CAAC;QACD,MAAM,KAAK,GAAG,GAAG,GAAG,GAAG,CAAC;QACxB,8DAA8D;QAC9D,IAAI,KAAK,GAAG,mBAAmB,EAAE,CAAC;YAChC,MAAM,IAAI,UAAU,CAAC,qDAAqD,CAAC,CAAC;QAC9E,CAAC;QACD,OAAO,GAAG,GAAG,CAAC,IAAI,CAAC,OAAO,EAAE,GAAG,KAAK,CAAC,CAAC;IACxC,CAAC;CACF;AA9CD,gDA8CC;AAED,wBAAwB;AACxB,MAAM,gBAAgB,GAAG,mBAAmB,CAAC;AAC7C,MAAM,SAAS,GAAG,mBAAmB,CAAC;AAEtC;;GAEG;AACH,SAAgB,SAAS,CAAC,CAAS;IACjC,IAAI,IAAI,GAAG,gBAAgB,CAAC;IAC5B,MAAM,OAAO,GAAG,IAAI,WAAW,EAAE,CAAC;IAClC,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAEhC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC;QACrB,IAAI,GAAG,CAAC,IAAI,GAAG,SAAS,CAAC,GAAG,OAAO,CAAC;IACtC,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC"}
+{"version":3,"file":"xoshiro256.js","sourceRoot":"","sources":["../src/xoshiro256.ts"],"names":[],"mappings":";AAAA;;;;;;GAMG;;;AA6GH,8BAWC;AAtHD,MAAM,OAAO,GAAG,mBAAmB,CAAC;AAEpC,SAAS,IAAI,CAAC,CAAS,EAAE,CAAS;IAChC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;AACjD,CAAC;AAED;;GAEG;AACH,MAAM,UAAU;IAGd,YAAY,IAAY;QACtB,IAAI,CAAC,KAAK,GAAG,IAAI,GAAG,OAAO,CAAC;IAC9B,CAAC;IAED,IAAI;QACF,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,mBAAmB,CAAC,GAAG,OAAO,CAAC;QAC1D,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACnB,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,mBAAmB,CAAC,GAAG,OAAO,CAAC;QACvD,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,mBAAmB,CAAC,GAAG,OAAO,CAAC;QACvD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC;IACpC,CAAC;CACF;AAED;;;;;GAKG;AACH,MAAa,kBAAkB;IAG7B,YAAY,IAAY;QACtB,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC;QAChC,IAAI,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC;IAC5D,CAAC;IAED,OAAO;QACL,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACrB,MAAM,MAAM,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;QAErE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,OAAO,CAAC;QAElC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAEb,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACV,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAEvB,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,2EAA2E;IAC3E,yBAAyB;IACzB,2EAA2E;IAE3E,OAAO;QACL,qDAAqD;QACrD,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,EAAE,CAAC;QAC3B,OAAO,MAAM,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,EAAE,IAAI,GAAG,CAAC,CAAC,CAAC;IACxD,CAAC;IAED,YAAY,CAAC,GAAW,EAAE,GAAW;QACnC,IAAI,GAAG,IAAI,GAAG;YAAE,OAAO,GAAG,CAAC;QAC3B,OAAO,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,EAAE,CAAC;IAC5C,CAAC;IAED,2EAA2E;IAC3E,kBAAkB;IAClB,2EAA2E;IAE3E;;;OAGG;IACH,UAAU,CAAC,GAAW,EAAE,GAAW;QACjC,IAAI,GAAG,IAAI,GAAG,EAAE,CAAC;YACf,OAAO,GAAG,CAAC;QACb,CAAC;QACD,MAAM,KAAK,GAAG,GAAG,GAAG,GAAG,CAAC;QACxB,8DAA8D;QAC9D,IAAI,KAAK,GAAG,mBAAmB,EAAE,CAAC;YAChC,MAAM,IAAI,UAAU,CAAC,qDAAqD,CAAC,CAAC;QAC9E,CAAC;QACD,OAAO,GAAG,GAAG,CAAC,IAAI,CAAC,OAAO,EAAE,GAAG,KAAK,CAAC,CAAC;IACxC,CAAC;IAED,2EAA2E;IAC3E,kBAAkB;IAClB,2EAA2E;IAE3E,WAAW;QACT,OAAO,IAAI,CAAC,OAAO,EAAE,GAAG,GAAG,CAAC;IAC9B,CAAC;CACF;AAnED,gDAmEC;AAED,wBAAwB;AACxB,MAAM,gBAAgB,GAAG,mBAAmB,CAAC;AAC7C,MAAM,SAAS,GAAG,mBAAmB,CAAC;AAEtC;;GAEG;AACH,SAAgB,SAAS,CAAC,CAAS;IACjC,IAAI,IAAI,GAAG,gBAAgB,CAAC;IAC5B,MAAM,OAAO,GAAG,IAAI,WAAW,EAAE,CAAC;IAClC,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAEhC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC;QACrB,IAAI,GAAG,CAAC,IAAI,GAAG,SAAS,CAAC,GAAG,OAAO,CAAC;IACtC,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC"}

package/package.json

@@ -1,20 +1,9 @@
 {
   "name": "pragmastat",
-  "version": "5.0.0",
+  "version": "5.2.0",
   "description": "Pragmastat: Pragmatic Statistical Toolkit",
   "main": "dist/index.js",
   "types": "dist/index.d.ts",
-  "scripts": {
-    "build": "tsc",
-    "test": "jest",
-    "test:watch": "jest --watch",
-    "test:coverage": "jest --coverage",
-    "lint": "eslint src tests --ext .ts",
-    "format": "prettier --write \"src/**/*.ts\" \"tests/**/*.ts\"",
-    "format:check": "prettier --check \"src/**/*.ts\" \"tests/**/*.ts\"",
-    "clean": "rm -rf dist coverage",
-    "prepublishOnly": "npm run clean && npm run build && npm run test"
-  },
   "keywords": [
     "statistics"
   ],
@@ -41,5 +30,15 @@
     "url": "https://github.com/AndreyAkinshin/pragmastat.git"
   },
   "homepage": "https://pragmastat.dev",
-  "doi": "10.5281/zenodo.17236778"
+  "doi": "10.5281/zenodo.17236778",
+  "scripts": {
+    "build": "tsc",
+    "test": "jest",
+    "test:watch": "jest --watch",
+    "test:coverage": "jest --coverage",
+    "lint": "eslint src tests --ext .ts",
+    "format": "prettier --write \"src/**/*.ts\" \"tests/**/*.ts\"",
+    "format:check": "prettier --check \"src/**/*.ts\" \"tests/**/*.ts\"",
+    "clean": "rm -rf dist coverage"
+  }
 }