ppef 1.0.0 → 1.1.0

package/dist/robustness/perturbations.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"perturbations.d.ts","sourceRoot":"","sources":["../../src/robustness/perturbations.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAKH,OAAO,KAAK,EAAE,YAAY,EAAE,kBAAkB,EAAE,MAAM,0BAA0B,CAAC;AAEjF;;;GAGG;AACH,eAAO,MAAM,uBAAuB,EAAE,YAgDrC,CAAC;AAEF;;;GAGG;AACH,eAAO,MAAM,qBAAqB,EAAE,YA4BnC,CAAC;AAEF;;;GAGG;AACH,eAAO,MAAM,uBAAuB,EAAE,YA+BrC,CAAC;AAEF;;;GAGG;AACH,eAAO,MAAM,uBAAuB,EAAE,YA+BrC,CAAC;AAEF;;GAEG;AACH,eAAO,MAAM,aAAa,EAAE,YAAY,EAKvC,CAAC;AAEF;;;;;GAKG;AACH,eAAO,MAAM,eAAe,GAAI,IAAI,MAAM,KAAG,YAAY,GAAG,SACrB,CAAC;AAExC;;;;;GAKG;AACH,eAAO,MAAM,kBAAkB,GAAI,QAAQ,kBAAkB,KAAG,YAU/D,CAAC"}

package/dist/robustness/perturbations.js

@@ -0,0 +1,184 @@
+/**
+ * Perturbation Implementations
+ *
+ * Concrete perturbations for robustness analysis.
+ * Each perturbation modifies evaluation cases in controlled ways.
+ */
+import { createHash } from "node:crypto";
+/**
+ * Edge removal perturbation.
+ * Randomly removes a fraction of edges from the graph.
+ */
+export const edgeRemovalPerturbation = {
+    id: "edge-removal",
+    name: "Edge Removal",
+    description: "Randomly remove a fraction of edges",
+    type: "structural",
+    intensity: 0.1, // Default: remove 10% of edges
+    apply(evaluationCase, seed = 42) {
+        const intensity = this.intensity ?? 0.1;
+        // Note: rng would be used when actually perturbing the graph at load time
+        // const rng = new SeededRandom(seed);
+        // Create perturbed case ID
+        const perturbedId = createHash("sha256")
+            .update(`${evaluationCase.caseId}-edge-removal-${intensity}-${seed}`)
+            .digest("hex")
+            .slice(0, 16);
+        // Mark artefacts as perturbed (actual perturbation happens at load time)
+        const perturbedArtefacts = (evaluationCase.inputs.artefacts ?? []).map((a) => ({
+            ...a,
+            metadata: {
+                ...a.metadata,
+                perturbation: "edge-removal",
+                perturbationIntensity: intensity,
+                perturbationSeed: seed,
+            },
+        }));
+        return {
+            ...evaluationCase,
+            caseId: perturbedId,
+            name: `${evaluationCase.name ?? "Case"} (${Math.round(intensity * 100)}% edges removed)`,
+            inputs: {
+                ...evaluationCase.inputs,
+                summary: {
+                    ...evaluationCase.inputs.summary,
+                    perturbation: "edge-removal",
+                    perturbationIntensity: intensity,
+                    perturbationSeed: seed,
+                },
+                artefacts: perturbedArtefacts,
+            },
+            tags: [...(evaluationCase.tags ?? []), "perturbed", "edge-removal"],
+        };
+    },
+};
+/**
+ * Seed shift perturbation.
+ * Shifts seed nodes to their neighbors.
+ */
+export const seedShiftPerturbation = {
+    id: "seed-shift",
+    name: "Seed Shift",
+    description: "Shift seed nodes to random neighbors",
+    type: "seed",
+    intensity: 1, // Default: shift all seeds
+    apply: (evaluationCase, seed = 42) => {
+        const perturbedId = createHash("sha256")
+            .update(`${evaluationCase.caseId}-seed-shift-${seed}`)
+            .digest("hex")
+            .slice(0, 16);
+        return {
+            ...evaluationCase,
+            caseId: perturbedId,
+            name: `${evaluationCase.name ?? "Case"} (seeds shifted)`,
+            inputs: {
+                ...evaluationCase.inputs,
+                summary: {
+                    ...evaluationCase.inputs.summary,
+                    perturbation: "seed-shift",
+                    perturbationSeed: seed,
+                },
+            },
+            tags: [...(evaluationCase.tags ?? []), "perturbed", "seed-shift"],
+        };
+    },
+};
+/**
+ * Node removal perturbation.
+ * Randomly removes non-seed nodes from the graph.
+ */
+export const nodeRemovalPerturbation = {
+    id: "node-removal",
+    name: "Node Removal",
+    description: "Randomly remove non-seed nodes",
+    type: "structural",
+    intensity: 0.05, // Default: remove 5% of nodes
+    apply(evaluationCase, seed = 42) {
+        const intensity = this.intensity ?? 0.05;
+        const perturbedId = createHash("sha256")
+            .update(`${evaluationCase.caseId}-node-removal-${intensity}-${seed}`)
+            .digest("hex")
+            .slice(0, 16);
+        return {
+            ...evaluationCase,
+            caseId: perturbedId,
+            name: `${evaluationCase.name ?? "Case"} (${Math.round(intensity * 100)}% nodes removed)`,
+            inputs: {
+                ...evaluationCase.inputs,
+                summary: {
+                    ...evaluationCase.inputs.summary,
+                    perturbation: "node-removal",
+                    perturbationIntensity: intensity,
+                    perturbationSeed: seed,
+                },
+            },
+            tags: [...(evaluationCase.tags ?? []), "perturbed", "node-removal"],
+        };
+    },
+};
+/**
+ * Weight noise perturbation.
+ * Adds Gaussian noise to edge weights.
+ */
+export const weightNoisePerturbation = {
+    id: "weight-noise",
+    name: "Weight Noise",
+    description: "Add Gaussian noise to edge weights",
+    type: "noise",
+    intensity: 0.1, // Default: 10% noise (std dev as fraction of weight)
+    apply(evaluationCase, seed = 42) {
+        const intensity = this.intensity ?? 0.1;
+        const perturbedId = createHash("sha256")
+            .update(`${evaluationCase.caseId}-weight-noise-${intensity}-${seed}`)
+            .digest("hex")
+            .slice(0, 16);
+        return {
+            ...evaluationCase,
+            caseId: perturbedId,
+            name: `${evaluationCase.name ?? "Case"} (${Math.round(intensity * 100)}% weight noise)`,
+            inputs: {
+                ...evaluationCase.inputs,
+                summary: {
+                    ...evaluationCase.inputs.summary,
+                    perturbation: "weight-noise",
+                    perturbationIntensity: intensity,
+                    perturbationSeed: seed,
+                },
+            },
+            tags: [...(evaluationCase.tags ?? []), "perturbed", "weight-noise"],
+        };
+    },
+};
+/**
+ * All built-in perturbations.
+ */
+export const PERTURBATIONS = [
+    edgeRemovalPerturbation,
+    seedShiftPerturbation,
+    nodeRemovalPerturbation,
+    weightNoisePerturbation,
+];
+/**
+ * Get a perturbation by ID.
+ *
+ * @param id - Perturbation identifier
+ * @returns Perturbation or undefined
+ */
+export const getPerturbation = (id) => PERTURBATIONS.find((p) => p.id === id);
+/**
+ * Create a perturbation with custom intensity.
+ *
+ * @param config - Perturbation configuration
+ * @returns Perturbation with configured intensity
+ */
+export const createPerturbation = (config) => {
+    const basePerturbation = getPerturbation(config.type);
+    if (!basePerturbation) {
+        throw new Error(`Unknown perturbation type: ${config.type}`);
+    }
+    return {
+        ...basePerturbation,
+        intensity: config.intensity,
+    };
+};
+//# sourceMappingURL=perturbations.js.map

package/dist/robustness/perturbations.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"perturbations.js","sourceRoot":"","sources":["../../src/robustness/perturbations.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,OAAO,EAAE,UAAU,EAAE,MAAM,aAAa,CAAC;AAKzC;;;GAGG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAiB;IACpD,EAAE,EAAE,cAAc;IAClB,IAAI,EAAE,cAAc;IACpB,WAAW,EAAE,qCAAqC;IAClD,IAAI,EAAE,YAAY;IAClB,SAAS,EAAE,GAAG,EAAE,+BAA+B;IAE/C,KAAK,CAAC,cAA8B,EAAE,IAAI,GAAG,EAAE;QAC9C,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,IAAI,GAAG,CAAC;QACxC,0EAA0E;QAC1E,sCAAsC;QAEtC,2BAA2B;QAC3B,MAAM,WAAW,GAAG,UAAU,CAAC,QAAQ,CAAC;aACtC,MAAM,CAAC,GAAG,cAAc,CAAC,MAAM,iBAAiB,SAAS,IAAI,IAAI,EAAE,CAAC;aACpE,MAAM,CAAC,KAAK,CAAC;aACb,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAEf,yEAAyE;QACzE,MAAM,kBAAkB,GAAwB,CAAC,cAAc,CAAC,MAAM,CAAC,SAAS,IAAI,EAAE,CAAC,CAAC,GAAG,CAC1F,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;YACP,GAAG,CAAC;YACJ,QAAQ,EAAE;gBACT,GAAG,CAAC,CAAC,QAAQ;gBACb,YAAY,EAAE,cAAc;gBAC5B,qBAAqB,EAAE,SAAS;gBAChC,gBAAgB,EAAE,IAAI;aACtB;SACD,CAAC,CACF,CAAC;QAEF,OAAO;YACN,GAAG,cAAc;YACjB,MAAM,EAAE,WAAW;YACnB,IAAI,EAAE,GAAG,cAAc,CAAC,IAAI,IAAI,MAAM,KAAK,IAAI,CAAC,KAAK,CAAC,SAAS,GAAG,GAAG,CAAC,kBAAkB;YACxF,MAAM,EAAE;gBACP,GAAG,cAAc,CAAC,MAAM;gBACxB,OAAO,EAAE;oBACR,GAAG,cAAc,CAAC,MAAM,CAAC,OAAO;oBAChC,YAAY,EAAE,cAAc;oBAC5B,qBAAqB,EAAE,SAAS;oBAChC,gBAAgB,EAAE,IAAI;iBACtB;gBACD,SAAS,EAAE,kBAAkB;aAC7B;YACD,IAAI,EAAE,CAAC,GAAG,CAAC,cAAc,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE,WAAW,EAAE,cAAc,CAAC;SACnE,CAAC;IACH,CAAC;CACD,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,qBAAqB,GAAiB;IAClD,EAAE,EAAE,YAAY;IAChB,IAAI,EAAE,YAAY;IAClB,WAAW,EAAE,sCAAsC;IACnD,IAAI,EAAE,MAAM;IACZ,SAAS,EAAE,CAAC,EAAE,2BAA2B;IAEzC,KAAK,EAAE,CAAC,cAA8B,EAAE,IAAI,GAAG,EAAE,EAAkB,EAAE;QACpE,MAAM,WAAW,GAAG,UAAU,CAAC,QAAQ,CAAC;aACtC,MAAM,CAAC,GAAG,cAAc,CAAC,MAAM,eAAe,IAAI,EAAE,CAAC;aACrD,MAAM,CAAC,KAAK,CAAC;aACb,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAEf,OAAO;YACN,GAAG,cAAc;YACjB,MAAM,EAAE,WAAW;YACnB,IAAI,EAAE,GAAG,cAAc,CAAC,IAAI,IAAI,MAAM,kBAAkB;YACxD,MAAM,EAAE;gBACP,GAAG,cAAc,CAAC,MAAM;gBACxB,OAAO,EAAE;oBACR,GAAG,cAAc,CAAC,MAAM,CAAC,OAAO;oBAChC,YAAY,EAAE,YAAY;oBAC1B,gBAAgB,EAAE,IAAI;iBACtB;aACD;YACD,IAAI,EAAE,CAAC,GAAG,CAAC,cAAc,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE,WAAW,EAAE,YAAY,CAAC;SACjE,CAAC;IACH,CAAC;CACD,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAiB;IACpD,EAAE,EAAE,cAAc;IAClB,IAAI,EAAE,cAAc;IACpB,WAAW,EAAE,gCAAgC;IAC7C,IAAI,EAAE,YAAY;IAClB,SAAS,EAAE,IAAI,EAAE,8BAA8B;IAE/C,KAAK,CAAC,cAA8B,EAAE,IAAI,GAAG,EAAE;QAC9C,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,IAAI,IAAI,CAAC;QAEzC,MAAM,WAAW,GAAG,UAAU,CAAC,QAAQ,CAAC;aACtC,MAAM,CAAC,GAAG,cAAc,CAAC,MAAM,iBAAiB,SAAS,IAAI,IAAI,EAAE,CAAC;aACpE,MAAM,CAAC,KAAK,CAAC;aACb,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAEf,OAAO;YACN,GAAG,cAAc;YACjB,MAAM,EAAE,WAAW;YACnB,IAAI,EAAE,GAAG,cAAc,CAAC,IAAI,IAAI,MAAM,KAAK,IAAI,CAAC,KAAK,CAAC,SAAS,GAAG,GAAG,CAAC,kBAAkB;YACxF,MAAM,EAAE;gBACP,GAAG,cAAc,CAAC,MAAM;gBACxB,OAAO,EAAE;oBACR,GAAG,cAAc,CAAC,MAAM,CAAC,OAAO;oBAChC,YAAY,EAAE,cAAc;oBAC5B,qBAAqB,EAAE,SAAS;oBAChC,gBAAgB,EAAE,IAAI;iBACtB;aACD;YACD,IAAI,EAAE,CAAC,GAAG,CAAC,cAAc,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE,WAAW,EAAE,cAAc,CAAC;SACnE,CAAC;IACH,CAAC;CACD,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAiB;IACpD,EAAE,EAAE,cAAc;IAClB,IAAI,EAAE,cAAc;IACpB,WAAW,EAAE,oCAAoC;IACjD,IAAI,EAAE,OAAO;IACb,SAAS,EAAE,GAAG,EAAE,qDAAqD;IAErE,KAAK,CAAC,cAA8B,EAAE,IAAI,GAAG,EAAE;QAC9C,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,IAAI,GAAG,CAAC;QAExC,MAAM,WAAW,GAAG,UAAU,CAAC,QAAQ,CAAC;aACtC,MAAM,CAAC,GAAG,cAAc,CAAC,MAAM,iBAAiB,SAAS,IAAI,IAAI,EAAE,CAAC;aACpE,MAAM,CAAC,KAAK,CAAC;aACb,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAEf,OAAO;YACN,GAAG,cAAc;YACjB,MAAM,EAAE,WAAW;YACnB,IAAI,EAAE,GAAG,cAAc,CAAC,IAAI,IAAI,MAAM,KAAK,IAAI,CAAC,KAAK,CAAC,SAAS,GAAG,GAAG,CAAC,iBAAiB;YACvF,MAAM,EAAE;gBACP,GAAG,cAAc,CAAC,MAAM;gBACxB,OAAO,EAAE;oBACR,GAAG,cAAc,CAAC,MAAM,CAAC,OAAO;oBAChC,YAAY,EAAE,cAAc;oBAC5B,qBAAqB,EAAE,SAAS;oBAChC,gBAAgB,EAAE,IAAI;iBACtB;aACD;YACD,IAAI,EAAE,CAAC,GAAG,CAAC,cAAc,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE,WAAW,EAAE,cAAc,CAAC;SACnE,CAAC;IACH,CAAC;CACD,CAAC;AAEF;;GAEG;AACH,MAAM,CAAC,MAAM,aAAa,GAAmB;IAC5C,uBAAuB;IACvB,qBAAqB;IACrB,uBAAuB;IACvB,uBAAuB;CACvB,CAAC;AAEF;;;;;GAKG;AACH,MAAM,CAAC,MAAM,eAAe,GAAG,CAAC,EAAU,EAA4B,EAAE,CACvE,aAAa,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC;AAExC;;;;;GAKG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAAC,MAA0B,EAAgB,EAAE;IAC9E,MAAM,gBAAgB,GAAG,eAAe,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IACtD,IAAI,CAAC,gBAAgB,EAAE,CAAC;QACvB,MAAM,IAAI,KAAK,CAAC,8BAA8B,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC;IAC9D,CAAC;IAED,OAAO;QACN,GAAG,gBAAgB;QACnB,SAAS,EAAE,MAAM,CAAC,SAAS;KAC3B,CAAC;AACH,CAAC,CAAC"}

package/dist/statistical/__tests__/mann-whitney-u.unit.test.d.ts

@@ -0,0 +1,7 @@
+/**
+ * Unit tests for Mann-Whitney U test and related statistical functions.
+ *
+ * Tests the non-parametric statistical test for comparing two independent samples.
+ */
+export {};
+//# sourceMappingURL=mann-whitney-u.unit.test.d.ts.map

package/dist/statistical/__tests__/mann-whitney-u.unit.test.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"mann-whitney-u.unit.test.d.ts","sourceRoot":"","sources":["../../../src/statistical/__tests__/mann-whitney-u.unit.test.ts"],"names":[],"mappings":"AAAA;;;;GAIG"}

package/dist/statistical/__tests__/mann-whitney-u.unit.test.js

@@ -0,0 +1,185 @@
+/**
+ * Unit tests for Mann-Whitney U test and related statistical functions.
+ *
+ * Tests the non-parametric statistical test for comparing two independent samples.
+ */
+import { describe, it } from "node:test";
+import { strict as assert } from "node:assert";
+import { normalCDF, mannWhitneyUTest, cohensD, confidenceInterval } from "../mann-whitney-u.js";
+describe("normalCDF", () => {
+    it("should return 0.5 for z = 0", () => {
+        const result = normalCDF(0);
+        assert.ok(Math.abs(result - 0.5) < 0.001);
+    });
+    it("should approach 1 for large positive z", () => {
+        const result = normalCDF(5);
+        assert.ok(result > 0.99);
+    });
+    it("should approach 0 for large negative z", () => {
+        const result = normalCDF(-5);
+        assert.ok(result < 0.01);
+    });
+    it("should be symmetric around 0", () => {
+        const z = 1.5;
+        const positive = normalCDF(z);
+        const negative = normalCDF(-z);
+        assert.ok(Math.abs(positive + negative - 1) < 0.001);
+    });
+    it("should handle negative z values", () => {
+        const result = normalCDF(-1);
+        assert.ok(result < 0.5);
+        assert.ok(result > 0.1);
+    });
+});
+describe("mannWhitneyUTest", () => {
+    it("should calculate U statistic for distinct samples", () => {
+        const sampleA = [1, 2, 3];
+        const sampleB = [4, 5, 6];
+        const result = mannWhitneyUTest(sampleA, sampleB);
+        assert.ok(typeof result.u === "number");
+        assert.ok(result.u >= 0);
+        assert.ok(typeof result.pValue === "number");
+        assert.ok(typeof result.significant === "boolean");
+    });
+    it("should detect significant difference between distinct samples", () => {
+        const sampleA = [1, 1, 1, 1, 1];
+        const sampleB = [100, 100, 100, 100, 100];
+        const result = mannWhitneyUTest(sampleA, sampleB);
+        assert.ok(result.pValue < 0.05);
+        assert.strictEqual(result.significant, true);
+    });
+    it("should not detect difference in similar samples", () => {
+        const sampleA = [5, 5, 5, 5, 5];
+        const sampleB = [5, 5, 5, 5, 5];
+        const result = mannWhitneyUTest(sampleA, sampleB);
+        // When samples are identical, p-value should be 1
+        assert.ok(result.pValue >= 0.95);
+        assert.strictEqual(result.significant, false);
+    });
+    it("should handle overlapping samples", () => {
+        const sampleA = [1, 2, 3, 4, 5];
+        const sampleB = [3, 4, 5, 6, 7];
+        const result = mannWhitneyUTest(sampleA, sampleB);
+        assert.ok(result.pValue >= 0 && result.pValue <= 1);
+        assert.ok(typeof result.significant === "boolean");
+    });
+    it("should handle samples with ties", () => {
+        const sampleA = [1, 2, 2, 3];
+        const sampleB = [2, 3, 4, 5];
+        const result = mannWhitneyUTest(sampleA, sampleB);
+        assert.ok(typeof result.u === "number");
+        assert.ok(result.pValue >= 0 && result.pValue <= 1);
+    });
+    it("should handle single-element samples", () => {
+        const sampleA = [5];
+        const sampleB = [10];
+        const result = mannWhitneyUTest(sampleA, sampleB);
+        assert.ok(typeof result.u === "number");
+        assert.ok(typeof result.pValue === "number");
+    });
+    it("should handle empty sample A", () => {
+        const sampleA = [];
+        const sampleB = [1, 2, 3];
+        const result = mannWhitneyUTest(sampleA, sampleB);
+        assert.strictEqual(result.u, 0);
+    });
+    it("should handle empty sample B", () => {
+        const sampleA = [1, 2, 3];
+        const sampleB = [];
+        const result = mannWhitneyUTest(sampleA, sampleB);
+        assert.strictEqual(result.u, 0);
+    });
+});
+describe("cohensD", () => {
+    it("should calculate effect size for distinct samples", () => {
+        const sampleA = [1, 2, 3, 4, 5];
+        const sampleB = [6, 7, 8, 9, 10];
+        const result = cohensD(sampleA, sampleB);
+        assert.ok(result > 0);
+        assert.ok(typeof result === "number");
+    });
+    it("should be zero for identical samples", () => {
+        const sampleA = [5, 5, 5, 5, 5];
+        const sampleB = [5, 5, 5, 5, 5];
+        const result = cohensD(sampleA, sampleB);
+        assert.strictEqual(result, 0);
+    });
+    it("should be symmetric", () => {
+        const sampleA = [1, 2, 3, 4, 5];
+        const sampleB = [10, 20, 30, 40, 50];
+        const d1 = cohensD(sampleA, sampleB);
+        const d2 = cohensD(sampleB, sampleA);
+        assert.ok(Math.abs(d1 - d2) < 0.001);
+    });
+    it("should handle large effect sizes", () => {
+        const sampleA = [1, 2, 1, 2, 1]; // Mean ~1.4, small variance
+        const sampleB = [100, 101, 100, 101, 100]; // Mean ~100.4, small variance
+        const result = cohensD(sampleA, sampleB);
+        assert.ok(result > 0.8); // Large effect due to separation of means
+    });
+    it("should handle small effect sizes", () => {
+        const sampleA = [10, 11, 12, 13, 14];
+        const sampleB = [12, 13, 14, 15, 16];
+        const result = cohensD(sampleA, sampleB);
+        assert.ok(result >= 0);
+        assert.ok(result < 2); // Should be reasonable
+    });
+    it("should handle samples with variance", () => {
+        const sampleA = [1, 2, 3, 4, 5];
+        const sampleB = [3, 4, 5, 6, 7];
+        const result = cohensD(sampleA, sampleB);
+        assert.ok(result >= 0);
+    });
+    it("should handle samples with same mean but different variance", () => {
+        const sampleA = [0, 10, 0, 10, 0, 10];
+        const sampleB = [4, 6, 4, 6, 4, 6];
+        // Both have mean 5
+        const meanA = sampleA.reduce((a, b) => a + b, 0) / sampleA.length;
+        const meanB = sampleB.reduce((a, b) => a + b, 0) / sampleB.length;
+        const result = cohensD(sampleA, sampleB);
+        assert.strictEqual(meanA, 5);
+        assert.strictEqual(meanB, 5);
+        assert.strictEqual(result, 0);
+    });
+});
+describe("confidenceInterval", () => {
+    it("should calculate interval for sample values", () => {
+        const values = [1, 2, 3, 4, 5];
+        const result = confidenceInterval(values);
+        assert.ok(typeof result.lower === "number");
+        assert.ok(typeof result.upper === "number");
+        assert.ok(result.lower < result.upper);
+    });
+    it("should be centered around mean", () => {
+        const values = [10, 20, 30, 40, 50];
+        const mean = values.reduce((a, b) => a + b, 0) / values.length;
+        const result = confidenceInterval(values);
+        const midpoint = (result.lower + result.upper) / 2;
+        assert.ok(Math.abs(midpoint - mean) < 0.1);
+    });
+    it("should handle identical values", () => {
+        const values = [5, 5, 5, 5, 5];
+        const result = confidenceInterval(values);
+        // With no variance, lower and upper should be same
+        assert.ok(Math.abs(result.lower - 5) < 0.01);
+        assert.ok(Math.abs(result.upper - 5) < 0.01);
+    });
+    it("should handle two values", () => {
+        const values = [0, 10];
+        const result = confidenceInterval(values);
+        assert.ok(result.lower < result.upper);
+        assert.ok(result.lower < 5);
+        assert.ok(result.upper > 5);
+    });
+    it("should have reasonable interval width for low variance", () => {
+        const values = [50, 51, 50, 51, 50, 51];
+        const mean = values.reduce((a, b) => a + b, 0) / values.length;
+        const result = confidenceInterval(values);
+        // Low variance should give narrow interval
+        const width = result.upper - result.lower;
+        assert.ok(width < 2); // Should be narrow
+        assert.ok(result.lower <= mean);
+        assert.ok(result.upper >= mean);
+    });
+});
+//# sourceMappingURL=mann-whitney-u.unit.test.js.map

package/dist/statistical/__tests__/mann-whitney-u.unit.test.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mann-whitney-u.unit.test.js","sourceRoot":"","sources":["../../../src/statistical/__tests__/mann-whitney-u.unit.test.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,OAAO,EAAE,QAAQ,EAAE,EAAE,EAAE,MAAM,WAAW,CAAC;AACzC,OAAO,EAAE,MAAM,IAAI,MAAM,EAAE,MAAM,aAAa,CAAC;AAE/C,OAAO,EAAE,SAAS,EAAE,gBAAgB,EAAE,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAEhG,QAAQ,CAAC,WAAW,EAAE,GAAG,EAAE;IAC1B,EAAE,CAAC,6BAA6B,EAAE,GAAG,EAAE;QACtC,MAAM,MAAM,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC5B,MAAM,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;IAC3C,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,wCAAwC,EAAE,GAAG,EAAE;QACjD,MAAM,MAAM,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC5B,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,IAAI,CAAC,CAAC;IAC1B,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,wCAAwC,EAAE,GAAG,EAAE;QACjD,MAAM,MAAM,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,IAAI,CAAC,CAAC;IAC1B,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACvC,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QAC9B,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;IACtD,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,iCAAiC,EAAE,GAAG,EAAE;QAC1C,MAAM,MAAM,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;QACxB,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC;AACJ,CAAC,CAAC,CAAC;AAEH,QAAQ,CAAC,kBAAkB,EAAE,GAAG,EAAE;IACjC,EAAE,CAAC,mDAAmD,EAAE,GAAG,EAAE;QAC5D,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1B,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE1B,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAElD,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,CAAC,KAAK,QAAQ,CAAC,CAAC;QACxC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;QACzB,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,MAAM,KAAK,QAAQ,CAAC,CAAC;QAC7C,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,WAAW,KAAK,SAAS,CAAC,CAAC;IACpD,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,+DAA+D,EAAE,GAAG,EAAE;QACxE,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAE1C,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAElD,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,MAAM,GAAG,IAAI,CAAC,CAAC;QAChC,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,WAAW,EAAE,IAAI,CAAC,CAAC;IAC9C,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,iDAAiD,EAAE,GAAG,EAAE;QAC1D,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEhC,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAElD,kDAAkD;QAClD,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,MAAM,IAAI,IAAI,CAAC,CAAC;QACjC,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC;IAC/C,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,mCAAmC,EAAE,GAAG,EAAE;QAC5C,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEhC,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAElD,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,MAAM,IAAI,CAAC,IAAI,MAAM,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC;QACpD,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,WAAW,KAAK,SAAS,CAAC,CAAC;IACpD,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,iCAAiC,EAAE,GAAG,EAAE;QAC1C,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7B,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE7B,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAElD,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,CAAC,KAAK,QAAQ,CAAC,CAAC;QACxC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,MAAM,IAAI,CAAC,IAAI,MAAM,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC;IACrD,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,sCAAsC,EAAE,GAAG,EAAE;QAC/C,MAAM,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,OAAO,GAAG,CAAC,EAAE,CAAC,CAAC;QAErB,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAElD,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,CAAC,KAAK,QAAQ,CAAC,CAAC;QACxC,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,MAAM,KAAK,QAAQ,CAAC,CAAC;IAC9C,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACvC,MAAM,OAAO,GAAa,EAAE,CAAC;QAC7B,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE1B,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAElD,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACvC,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1B,MAAM,OAAO,GAAa,EAAE,CAAC;QAE7B,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAElD,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC;AACJ,CAAC,CAAC,CAAC;AAEH,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE;IACxB,EAAE,CAAC,mDAAmD,EAAE,GAAG,EAAE;QAC5D,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAEjC,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAEzC,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACtB,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,KAAK,QAAQ,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,sCAAsC,EAAE,GAAG,EAAE;QAC/C,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEhC,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAEzC,MAAM,CAAC,WAAW,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,qBAAqB,EAAE,GAAG,EAAE;QAC9B,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAErC,MAAM,EAAE,GAAG,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QACrC,MAAM,EAAE,GAAG,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAErC,MAAM,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC;IACtC,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,kCAAkC,EAAE,GAAG,EAAE;QAC3C,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,4BAA4B;QAC7D,MAAM,OAAO,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,8BAA8B;QAEzE,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAEzC,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,0CAA0C;IACpE,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,kCAAkC,EAAE,GAAG,EAAE;QAC3C,MAAM,OAAO,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACrC,MAAM,OAAO,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAErC,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAEzC,MAAM,CAAC,EAAE,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC;QACvB,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,uBAAuB;IAC/C,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,qCAAqC,EAAE,GAAG,EAAE;QAC9C,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEhC,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAEzC,MAAM,CAAC,EAAE,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC;IACxB,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,6DAA6D,EAAE,GAAG,EAAE;QACtE,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACtC,MAAM,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEnC,mBAAmB;QACnB,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC;QAClE,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC;QAElE,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAEzC,MAAM,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;QAC7B,MAAM,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;QAC7B,MAAM,CAAC,WAAW,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;AACJ,CAAC,CAAC,CAAC;AAEH,QAAQ,CAAC,oBAAoB,EAAE,GAAG,EAAE;IACnC,EAAE,CAAC,6CAA6C,EAAE,GAAG,EAAE;QACtD,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE/B,MAAM,MAAM,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;QAE1C,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,KAAK,KAAK,QAAQ,CAAC,CAAC;QAC5C,MAAM,CAAC,EAAE,CAAC,OAAO,MAAM,CAAC,KAAK,KAAK,QAAQ,CAAC,CAAC;QAC5C,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;IACxC,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,gCAAgC,EAAE,GAAG,EAAE;QACzC,MAAM,MAAM,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACpC,MAAM,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QAE/D,MAAM,MAAM,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;QAE1C,MAAM,QAAQ,GAAG,CAAC,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;QACnD,MAAM,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC;IAC5C,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,gCAAgC,EAAE,GAAG,EAAE;QACzC,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE/B,MAAM,MAAM,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;QAE1C,mDAAmD;QACnD,MAAM,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;QAC7C,MAAM,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IAC9C,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,0BAA0B,EAAE,GAAG,EAAE;QACnC,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAEvB,MAAM,MAAM,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;QAE1C,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;QACvC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC5B,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IAC7B,CAAC,CAAC,CAAC;IAEH,EAAE,CAAC,wDAAwD,EAAE,GAAG,EAAE;QACjE,MAAM,MAAM,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QAE/D,MAAM,MAAM,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;QAE1C,2CAA2C;QAC3C,MAAM,KAAK,GAAG,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC;QAC1C,MAAM,CAAC,EAAE,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,mBAAmB;QACzC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,KAAK,IAAI,IAAI,CAAC,CAAC;QAChC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,KAAK,IAAI,IAAI,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC;AACJ,CAAC,CAAC,CAAC"}

package/dist/statistical/index.d.ts

@@ -0,0 +1,8 @@
+/**
+ * Statistical Test Utilities
+ *
+ * Collection of statistical functions for hypothesis testing,
+ * effect size calculation, and confidence intervals.
+ */
+export { mannWhitneyUTest, normalCDF, cohensD, confidenceInterval } from "./mann-whitney-u.js";
+//# sourceMappingURL=index.d.ts.map

package/dist/statistical/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/statistical/index.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,OAAO,EAAE,gBAAgB,EAAE,SAAS,EAAE,OAAO,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC"}

package/dist/statistical/index.js

@@ -0,0 +1,8 @@
+/**
+ * Statistical Test Utilities
+ *
+ * Collection of statistical functions for hypothesis testing,
+ * effect size calculation, and confidence intervals.
+ */
+export { mannWhitneyUTest, normalCDF, cohensD, confidenceInterval } from "./mann-whitney-u.js";
+//# sourceMappingURL=index.js.map

package/dist/statistical/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/statistical/index.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,OAAO,EAAE,gBAAgB,EAAE,SAAS,EAAE,OAAO,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC"}

package/dist/statistical/mann-whitney-u.d.ts

@@ -0,0 +1,62 @@
+/**
+ * Mann-Whitney U Test
+ *
+ * Statistical test for comparing two independent samples.
+ * Tests whether two populations have the same distribution.
+ *
+ * H0: Both populations have the same distribution
+ * H1: Populations have different distributions
+ *
+ * Returns p-value (smaller = more significant difference)
+ */
+/**
+ * Standard normal cumulative distribution function.
+ * Uses the Abramowitz and Stegun approximation.
+ *
+ * @param z - Z-score
+ * @returns Cumulative probability from -infinity to z
+ */
+export declare const normalCDF: (z: number) => number;
+/**
+ * Mann-Whitney U test for comparing two independent samples.
+ *
+ * Non-parametric test that does not assume normal distribution.
+ * Tests whether two populations have the same distribution.
+ *
+ * @param sampleA - First sample array
+ * @param sampleB - Second sample array
+ * @returns Object containing U statistic, p-value, and significance flag
+ */
+export declare const mannWhitneyUTest: (sampleA: number[], sampleB: number[]) => {
+    u: number;
+    pValue: number;
+    significant: boolean;
+};
+/**
+ * Calculate Cohen's d effect size.
+ *
+ * Measures the standardized difference between two means.
+ *
+ * Interpretation:
+ * - 0.2: Small effect
+ * - 0.5: Medium effect
+ * - 0.8: Large effect
+ *
+ * @param sampleA - First sample array
+ * @param sampleB - Second sample array
+ * @returns Cohen's d effect size
+ */
+export declare const cohensD: (sampleA: number[], sampleB: number[]) => number;
+/**
+ * Calculate confidence interval for a mean.
+ *
+ * Uses t-distribution approximation (1.96 for 95% CI with large samples).
+ *
+ * @param values - Sample values
+ * @returns Object with lower and upper bounds
+ */
+export declare const confidenceInterval: (values: number[]) => {
+    lower: number;
+    upper: number;
+};
+//# sourceMappingURL=mann-whitney-u.d.ts.map

package/dist/statistical/mann-whitney-u.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"mann-whitney-u.d.ts","sourceRoot":"","sources":["../../src/statistical/mann-whitney-u.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,GAAI,GAAG,MAAM,KAAG,MAarC,CAAC;AAEF;;;;;;;;;GASG;AACH,eAAO,MAAM,gBAAgB,GAC5B,SAAS,MAAM,EAAE,EACjB,SAAS,MAAM,EAAE,KACf;IACF,CAAC,EAAE,MAAM,CAAC;IACV,MAAM,EAAE,MAAM,CAAC;IACf,WAAW,EAAE,OAAO,CAAC;CAgDrB,CAAC;AAEF;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,OAAO,GAAI,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,KAAG,MAa9D,CAAC;AAEF;;;;;;;GAOG;AACH,eAAO,MAAM,kBAAkB,GAAI,QAAQ,MAAM,EAAE,KAAG;IAAE,KAAK,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAYnF,CAAC"}

package/dist/statistical/mann-whitney-u.js

@@ -0,0 +1,127 @@
+/**
+ * Mann-Whitney U Test
+ *
+ * Statistical test for comparing two independent samples.
+ * Tests whether two populations have the same distribution.
+ *
+ * H0: Both populations have the same distribution
+ * H1: Populations have different distributions
+ *
+ * Returns p-value (smaller = more significant difference)
+ */
+/**
+ * Standard normal cumulative distribution function.
+ * Uses the Abramowitz and Stegun approximation.
+ *
+ * @param z - Z-score
+ * @returns Cumulative probability from -infinity to z
+ */
+export const normalCDF = (z) => {
+    const sign = z < 0 ? -1 : 1;
+    z = Math.abs(z) / Math.sqrt(2);
+    const a1 = 0.254_829_592;
+    const a2 = -0.284_496_736;
+    const a3 = 1.421_413_741;
+    const a4 = -1.453_152_027;
+    const a5 = 1.061_405_429;
+    const p = 0.327_591_1;
+    const t = 1 / (1 + p * z);
+    const y = 1 - ((((a5 * t + a4) * t + a3) * t + a2) * t + a1) * t * Math.exp(-z * z);
+    return 0.5 * (1 + sign * y);
+};
+/**
+ * Mann-Whitney U test for comparing two independent samples.
+ *
+ * Non-parametric test that does not assume normal distribution.
+ * Tests whether two populations have the same distribution.
+ *
+ * @param sampleA - First sample array
+ * @param sampleB - Second sample array
+ * @returns Object containing U statistic, p-value, and significance flag
+ */
+export const mannWhitneyUTest = (sampleA, sampleB) => {
+    // Rank all values combined
+    const combined = [...sampleA, ...sampleB];
+    const sorted = [...combined].sort((a, b) => a - b);
+    // Assign ranks (handle ties)
+    const ranks = new Map();
+    for (const [index, value] of sorted.entries()) {
+        if (!ranks.has(value)) {
+            ranks.set(value, []);
+        }
+        const positions = ranks.get(value);
+        if (positions) {
+            positions.push(index + 1);
+        }
+    }
+    // Average rank for tied values
+    const avgRanks = new Map();
+    for (const [value, positions] of ranks) {
+        avgRanks.set(value, positions.reduce((a, b) => a + b, 0) / positions.length);
+    }
+    // Sum ranks for each sample
+    const rankSumA = sampleA.reduce((sum, value) => sum + (avgRanks.get(value) ?? 0), 0);
+    const rankSumB = sampleB.reduce((sum, value) => sum + (avgRanks.get(value) ?? 0), 0);
+    // Calculate U statistics
+    const n1 = sampleA.length;
+    const n2 = sampleB.length;
+    const u1 = rankSumA - (n1 * (n1 + 1)) / 2;
+    const u2 = rankSumB - (n2 * (n2 + 1)) / 2;
+    const u = Math.min(u1, u2);
+    // Calculate z-score for large samples
+    const meanU = (n1 * n2) / 2;
+    const stdU = Math.sqrt((n1 * n2 * (n1 + n2 + 1)) / 12);
+    const z = stdU > 0 ? (u - meanU) / stdU : 0;
+    // Two-tailed p-value from z-score (approximation)
+    const pValue = 2 * (1 - normalCDF(Math.abs(z)));
+    return {
+        u,
+        pValue,
+        significant: pValue < 0.05, // 95% confidence level
+    };
+};
+/**
+ * Calculate Cohen's d effect size.
+ *
+ * Measures the standardized difference between two means.
+ *
+ * Interpretation:
+ * - 0.2: Small effect
+ * - 0.5: Medium effect
+ * - 0.8: Large effect
+ *
+ * @param sampleA - First sample array
+ * @param sampleB - Second sample array
+ * @returns Cohen's d effect size
+ */
+export const cohensD = (sampleA, sampleB) => {
+    const n1 = sampleA.length;
+    const n2 = sampleB.length;
+    const mean1 = sampleA.reduce((a, b) => a + b, 0) / n1;
+    const mean2 = sampleB.reduce((a, b) => a + b, 0) / n2;
+    const variable1 = sampleA.reduce((sum, value) => sum + (value - mean1) ** 2, 0) / (n1 - 1);
+    const variable2 = sampleB.reduce((sum, value) => sum + (value - mean2) ** 2, 0) / (n2 - 1);
+    const pooledStd = Math.sqrt(((n1 - 1) * variable1 + (n2 - 1) * variable2) / (n1 + n2 - 2));
+    return pooledStd > 0 ? Math.abs(mean1 - mean2) / pooledStd : 0;
+};
+/**
+ * Calculate confidence interval for a mean.
+ *
+ * Uses t-distribution approximation (1.96 for 95% CI with large samples).
+ *
+ * @param values - Sample values
+ * @returns Object with lower and upper bounds
+ */
+export const confidenceInterval = (values) => {
+    const n = values.length;
+    const mean = values.reduce((a, b) => a + b, 0) / n;
+    const std = Math.sqrt(values.reduce((sum, value) => sum + (value - mean) ** 2, 0) / (n - 1));
+    const se = std / Math.sqrt(n);
+    const t = 1.96; // Approximation for large samples (95% CI)
+    const margin = t * se;
+    return {
+        lower: mean - margin,
+        upper: mean + margin,
+    };
+};
+//# sourceMappingURL=mann-whitney-u.js.map

package/dist/statistical/mann-whitney-u.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mann-whitney-u.js","sourceRoot":"","sources":["../../src/statistical/mann-whitney-u.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,CAAS,EAAU,EAAE;IAC9C,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5B,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC/B,MAAM,EAAE,GAAG,aAAa,CAAC;IACzB,MAAM,EAAE,GAAG,CAAC,aAAa,CAAC;IAC1B,MAAM,EAAE,GAAG,aAAa,CAAC;IACzB,MAAM,EAAE,GAAG,CAAC,aAAa,CAAC;IAC1B,MAAM,EAAE,GAAG,aAAa,CAAC;IACzB,MAAM,CAAC,GAAG,WAAW,CAAC;IAEtB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACpF,OAAO,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;AAC7B,CAAC,CAAC;AAEF;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAC/B,OAAiB,EACjB,OAAiB,EAKhB,EAAE;IACH,2BAA2B;IAC3B,MAAM,QAAQ,GAAG,CAAC,GAAG,OAAO,EAAE,GAAG,OAAO,CAAC,CAAC;IAC1C,MAAM,MAAM,GAAG,CAAC,GAAG,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAEnD,6BAA6B;IAC7B,MAAM,KAAK,GAAG,IAAI,GAAG,EAAoB,CAAC;IAC1C,KAAK,MAAM,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE,EAAE,CAAC;QAC/C,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC;YACvB,KAAK,CAAC,GAAG,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACtB,CAAC;QACD,MAAM,SAAS,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACnC,IAAI,SAAS,EAAE,CAAC;YACf,SAAS,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC3B,CAAC;IACF,CAAC;IAED,+BAA+B;IAC/B,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAkB,CAAC;IAC3C,KAAK,MAAM,CAAC,KAAK,EAAE,SAAS,CAAC,IAAI,KAAK,EAAE,CAAC;QACxC,QAAQ,CAAC,GAAG,CAAC,KAAK,EAAE,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;IAC9E,CAAC;IAED,4BAA4B;IAC5B,MAAM,QAAQ,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrF,MAAM,QAAQ,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAErF,yBAAyB;IACzB,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,CAAC;IAC1B,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,CAAC;IAC1B,MAAM,EAAE,GAAG,QAAQ,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC1C,MAAM,EAAE,GAAG,QAAQ,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC1C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAE3B,sCAAsC;IACtC,MAAM,KAAK,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;IAC5B,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACvD,MAAM,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAE5C,kDAAkD;IAClD,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAEhD,OAAO;QACN,CAAC;QACD,MAAM;QACN,WAAW,EAAE,MAAM,GAAG,IAAI,EAAE,uBAAuB;KACnD,CAAC;AACH,CAAC,CAAC;AAEF;;;;;;;;;;;;;GAaG;AACH,MAAM,CAAC,MAAM,OAAO,GAAG,CAAC,OAAiB,EAAE,OAAiB,EAAU,EAAE;IACvE,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,CAAC;IAC1B,MAAM,EAAE,GAAG,OAAO,CAAC,MAAM,CAAC;IAE1B,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;IACtD,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;IAEtD,MAAM,SAAS,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IAC3F,MAAM,SAAS,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IAE3F,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,SAAS,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IAE3F,OAAO,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;AAChE,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAAC,MAAgB,EAAoC,EAAE;IACxF,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,MAAM,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;IACnD,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC7F,MAAM,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,2CAA2C;IAE3D,MAAM,MAAM,GAAG,CAAC,GAAG,EAAE,CAAC;IACtB,OAAO;QACN,KAAK,EAAE,IAAI,GAAG,MAAM;QACpB,KAAK,EAAE,IAAI,GAAG,MAAM;KACpB,CAAC;AACH,CAAC,CAAC"}