@datagrok-libraries/statistics 0.1.3 → 0.1.6
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/.eslintrc.json +28 -28
- package/README.md +2 -2
- package/package.json +36 -32
- package/src/correlation-coefficient.d.ts +15 -0
- package/src/correlation-coefficient.d.ts.map +1 -0
- package/src/correlation-coefficient.js +100 -0
- package/src/multiple-tests.d.ts +34 -34
- package/src/multiple-tests.d.ts.map +1 -1
- package/src/multiple-tests.js +202 -206
- package/src/tests.d.ts +11 -10
- package/src/tests.d.ts.map +1 -1
- package/src/tests.js +77 -77
- package/tsconfig.json +71 -71
- package/src/multiple-tests.ts +0 -234
- package/src/tests.ts +0 -101
package/.eslintrc.json
CHANGED
|
@@ -1,29 +1,29 @@
|
|
|
1
|
-
{
|
|
2
|
-
"env": {
|
|
3
|
-
"browser": true,
|
|
4
|
-
"es2021": true
|
|
5
|
-
},
|
|
6
|
-
"extends": [
|
|
7
|
-
"google"
|
|
8
|
-
],
|
|
9
|
-
"parserOptions": {
|
|
10
|
-
"ecmaVersion": 12,
|
|
11
|
-
"sourceType": "module"
|
|
12
|
-
},
|
|
13
|
-
"rules": {
|
|
14
|
-
"indent": [
|
|
15
|
-
"error",
|
|
16
|
-
2
|
|
17
|
-
],
|
|
18
|
-
"max-len": [
|
|
19
|
-
"error",
|
|
20
|
-
120
|
|
21
|
-
],
|
|
22
|
-
"require-jsdoc": "off",
|
|
23
|
-
"spaced-comment": "off"
|
|
24
|
-
},
|
|
25
|
-
"parser": "@typescript-eslint/parser",
|
|
26
|
-
"plugins": [
|
|
27
|
-
"@typescript-eslint"
|
|
28
|
-
]
|
|
1
|
+
{
|
|
2
|
+
"env": {
|
|
3
|
+
"browser": true,
|
|
4
|
+
"es2021": true
|
|
5
|
+
},
|
|
6
|
+
"extends": [
|
|
7
|
+
"google"
|
|
8
|
+
],
|
|
9
|
+
"parserOptions": {
|
|
10
|
+
"ecmaVersion": 12,
|
|
11
|
+
"sourceType": "module"
|
|
12
|
+
},
|
|
13
|
+
"rules": {
|
|
14
|
+
"indent": [
|
|
15
|
+
"error",
|
|
16
|
+
2
|
|
17
|
+
],
|
|
18
|
+
"max-len": [
|
|
19
|
+
"error",
|
|
20
|
+
120
|
|
21
|
+
],
|
|
22
|
+
"require-jsdoc": "off",
|
|
23
|
+
"spaced-comment": "off"
|
|
24
|
+
},
|
|
25
|
+
"parser": "@typescript-eslint/parser",
|
|
26
|
+
"plugins": [
|
|
27
|
+
"@typescript-eslint"
|
|
28
|
+
]
|
|
29
29
|
}
|
package/README.md
CHANGED
|
@@ -1,3 +1,3 @@
|
|
|
1
|
-
# statistics
|
|
2
|
-
|
|
1
|
+
# statistics
|
|
2
|
+
|
|
3
3
|
statistics is a [package](https://datagrok.ai/help/develop/develop#packages) for the [Datagrok](https://datagrok.ai) platform.
|
package/package.json
CHANGED
|
@@ -1,32 +1,36 @@
|
|
|
1
|
-
{
|
|
2
|
-
|
|
3
|
-
|
|
4
|
-
|
|
5
|
-
|
|
6
|
-
|
|
7
|
-
|
|
8
|
-
|
|
9
|
-
|
|
10
|
-
|
|
11
|
-
|
|
12
|
-
|
|
13
|
-
|
|
14
|
-
|
|
15
|
-
|
|
16
|
-
|
|
17
|
-
|
|
18
|
-
|
|
19
|
-
|
|
20
|
-
|
|
21
|
-
|
|
22
|
-
|
|
23
|
-
|
|
24
|
-
|
|
25
|
-
|
|
26
|
-
|
|
27
|
-
|
|
28
|
-
|
|
29
|
-
|
|
30
|
-
|
|
31
|
-
|
|
32
|
-
|
|
1
|
+
{
|
|
2
|
+
"name": "@datagrok-libraries/statistics",
|
|
3
|
+
"publishConfig": {
|
|
4
|
+
"access": "public"
|
|
5
|
+
},
|
|
6
|
+
"friendlyName": "statistics",
|
|
7
|
+
"version": "0.1.6",
|
|
8
|
+
"description": "",
|
|
9
|
+
"dependencies": {
|
|
10
|
+
"datagrok-api": ">0.95.4",
|
|
11
|
+
"@datagrok-libraries/utils": ">=0.0.22",
|
|
12
|
+
"cash-dom": "latest",
|
|
13
|
+
"dayjs": "latest",
|
|
14
|
+
"jstat": "latest"
|
|
15
|
+
},
|
|
16
|
+
"devDependencies": {
|
|
17
|
+
"ts-loader": "^9.2.6",
|
|
18
|
+
"typescript": "^4.2.3",
|
|
19
|
+
"eslint": "latest",
|
|
20
|
+
"eslint-config-google": "latest",
|
|
21
|
+
"@typescript-eslint/eslint-plugin": "latest",
|
|
22
|
+
"@typescript-eslint/parser": "latest"
|
|
23
|
+
},
|
|
24
|
+
"scripts": {
|
|
25
|
+
"link": "npm link",
|
|
26
|
+
"link-api": "npm link datagrok-api",
|
|
27
|
+
"link-utils": "npm link @datagrok-libraries/utils",
|
|
28
|
+
"link-all": "npm link datagrok-api @datagrok-libraries/utils",
|
|
29
|
+
"debug-statistics": "grok publish --rebuild",
|
|
30
|
+
"release-statistics": "grok publish --rebuild --release",
|
|
31
|
+
"build-statistics": "tsc",
|
|
32
|
+
"build": "tsc",
|
|
33
|
+
"lint": "eslint src --ext .ts",
|
|
34
|
+
"lint-fix": "eslint src --ext .ts --fix"
|
|
35
|
+
}
|
|
36
|
+
}
|
|
@@ -0,0 +1,15 @@
|
|
|
1
|
+
import { Options } from '../../utils/src/type-declarations';
|
|
2
|
+
/**
|
|
3
|
+
* Calculates Kendall's tau, a correlation measure for ordinal data, and an associated p-value.
|
|
4
|
+
* Returns: Kendall's tau, two-tailed p-value.
|
|
5
|
+
* Derived from older SciPy: http://web.mit.edu/6.863/spring2011/packages/scipy_src/scipy/stats/stats.py
|
|
6
|
+
*
|
|
7
|
+
* @export
|
|
8
|
+
* @param {number[]} x The first array.
|
|
9
|
+
* @param {number[]} y The second array.
|
|
10
|
+
* @return {Options} The result.
|
|
11
|
+
* @link https://github.com/pdfernhout/narrafirma/blob/
|
|
12
|
+
* c9c122d577a4b8868ce603badabbb7d10f45740c/webapp/source/statistics/kendallsTau.ts#L8
|
|
13
|
+
*/
|
|
14
|
+
export declare function kendallsTau(x: number[], y: number[]): Options;
|
|
15
|
+
//# sourceMappingURL=correlation-coefficient.d.ts.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"correlation-coefficient.d.ts","sourceRoot":"","sources":["correlation-coefficient.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,OAAO,EAAC,MAAM,mCAAmC,CAAC;AAwD1D;;;;;;;;;;;GAWG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,OAAO,CAiC7D"}
|
|
@@ -0,0 +1,100 @@
|
|
|
1
|
+
//import {Options} from '@datagrok-libraries/utils/src/type-declarations';
|
|
2
|
+
/**
|
|
3
|
+
* Returns the error function erf(x).
|
|
4
|
+
*
|
|
5
|
+
* @param {number} x An argument.
|
|
6
|
+
* @return {number} The result.
|
|
7
|
+
* @link https://github.com/jstat/jstat/blob/65ce096a99f753d6a22482e5e74accbfc1c33767/dist/jstat.js#L1562
|
|
8
|
+
*/
|
|
9
|
+
function erf(x) {
|
|
10
|
+
const cof = [-1.3026537197817094, 6.4196979235649026e-1, 1.9476473204185836e-2,
|
|
11
|
+
-9.561514786808631e-3, -9.46595344482036e-4, 3.66839497852761e-4,
|
|
12
|
+
4.2523324806907e-5, -2.0278578112534e-5, -1.624290004647e-6,
|
|
13
|
+
1.303655835580e-6, 1.5626441722e-8, -8.5238095915e-8,
|
|
14
|
+
6.529054439e-9, 5.059343495e-9, -9.91364156e-10,
|
|
15
|
+
-2.27365122e-10, 9.6467911e-11, 2.394038e-12,
|
|
16
|
+
-6.886027e-12, 8.94487e-13, 3.13092e-13,
|
|
17
|
+
-1.12708e-13, 3.81e-16, 7.106e-15,
|
|
18
|
+
-1.523e-15, -9.4e-17, 1.21e-16,
|
|
19
|
+
-2.8e-17];
|
|
20
|
+
let j = cof.length - 1;
|
|
21
|
+
let isneg = false;
|
|
22
|
+
let d = 0;
|
|
23
|
+
let dd = 0;
|
|
24
|
+
let t = 0;
|
|
25
|
+
let ty = 0;
|
|
26
|
+
let tmp = 0;
|
|
27
|
+
let res = 0;
|
|
28
|
+
if (x < 0) {
|
|
29
|
+
x = -x;
|
|
30
|
+
isneg = true;
|
|
31
|
+
}
|
|
32
|
+
t = 2 / (2 + x);
|
|
33
|
+
ty = 4 * t - 2;
|
|
34
|
+
for (; j > 0; j--) {
|
|
35
|
+
tmp = d;
|
|
36
|
+
d = ty * d - dd + cof[j];
|
|
37
|
+
dd = tmp;
|
|
38
|
+
}
|
|
39
|
+
res = t * Math.exp(-x * x + 0.5 * (cof[0] + ty * d) - dd);
|
|
40
|
+
return isneg ? res - 1 : 1 - res;
|
|
41
|
+
}
|
|
42
|
+
/**
|
|
43
|
+
* Returns the complmentary error function erfc(x)
|
|
44
|
+
*
|
|
45
|
+
* @param {number} x An argument.
|
|
46
|
+
* @return {number} The result.
|
|
47
|
+
* @link https://github.com/jstat/jstat/blob/65ce096a99f753d6a22482e5e74accbfc1c33767/dist/jstat.js#L1599
|
|
48
|
+
*/
|
|
49
|
+
function erfc(x) {
|
|
50
|
+
return 1 - erf(x);
|
|
51
|
+
}
|
|
52
|
+
/**
|
|
53
|
+
* Calculates Kendall's tau, a correlation measure for ordinal data, and an associated p-value.
|
|
54
|
+
* Returns: Kendall's tau, two-tailed p-value.
|
|
55
|
+
* Derived from older SciPy: http://web.mit.edu/6.863/spring2011/packages/scipy_src/scipy/stats/stats.py
|
|
56
|
+
*
|
|
57
|
+
* @export
|
|
58
|
+
* @param {number[]} x The first array.
|
|
59
|
+
* @param {number[]} y The second array.
|
|
60
|
+
* @return {Options} The result.
|
|
61
|
+
* @link https://github.com/pdfernhout/narrafirma/blob/
|
|
62
|
+
* c9c122d577a4b8868ce603badabbb7d10f45740c/webapp/source/statistics/kendallsTau.ts#L8
|
|
63
|
+
*/
|
|
64
|
+
export function kendallsTau(x, y) {
|
|
65
|
+
let n1 = 0;
|
|
66
|
+
let n2 = 0;
|
|
67
|
+
let iss = 0;
|
|
68
|
+
for (let j = 0; j < x.length - 1; j++) {
|
|
69
|
+
for (let k = j + 1; k < y.length; k++) {
|
|
70
|
+
const a1 = x[j] - x[k];
|
|
71
|
+
const a2 = y[j] - y[k];
|
|
72
|
+
const aa = a1 * a2;
|
|
73
|
+
if (aa) {
|
|
74
|
+
// neither array has a tie
|
|
75
|
+
n1 = n1 + 1;
|
|
76
|
+
n2 = n2 + 1;
|
|
77
|
+
if (aa > 0) {
|
|
78
|
+
iss = iss + 1;
|
|
79
|
+
}
|
|
80
|
+
else {
|
|
81
|
+
iss = iss - 1;
|
|
82
|
+
}
|
|
83
|
+
}
|
|
84
|
+
else {
|
|
85
|
+
if (a1) {
|
|
86
|
+
n1 = n1 + 1;
|
|
87
|
+
}
|
|
88
|
+
if (a2) {
|
|
89
|
+
n2 = n2 + 1;
|
|
90
|
+
}
|
|
91
|
+
}
|
|
92
|
+
}
|
|
93
|
+
}
|
|
94
|
+
const tau = iss / Math.sqrt(n1 * n2);
|
|
95
|
+
const svar = (4.0 * x.length + 10.0) / (9.0 * x.length * (x.length - 1));
|
|
96
|
+
const z = tau / Math.sqrt(svar);
|
|
97
|
+
const prob = erfc(Math.abs(z) / 1.4142136);
|
|
98
|
+
return { test: tau, z: z, prob: prob };
|
|
99
|
+
}
|
|
100
|
+
//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoiY29ycmVsYXRpb24tY29lZmZpY2llbnQuanMiLCJzb3VyY2VSb290IjoiIiwic291cmNlcyI6WyJjb3JyZWxhdGlvbi1jb2VmZmljaWVudC50cyJdLCJuYW1lcyI6W10sIm1hcHBpbmdzIjoiQUFDQSwwRUFBMEU7QUFFMUU7Ozs7OztHQU1HO0FBQ0gsU0FBUyxHQUFHLENBQUMsQ0FBUztJQUNwQixNQUFNLEdBQUcsR0FBRyxDQUFDLENBQUMsa0JBQWtCLEVBQUUscUJBQXFCLEVBQUUscUJBQXFCO1FBQzVFLENBQUMsb0JBQW9CLEVBQUUsQ0FBQyxtQkFBbUIsRUFBRSxtQkFBbUI7UUFDaEUsa0JBQWtCLEVBQUUsQ0FBQyxrQkFBa0IsRUFBRSxDQUFDLGlCQUFpQjtRQUMzRCxpQkFBaUIsRUFBRSxlQUFlLEVBQUUsQ0FBQyxlQUFlO1FBQ3BELGNBQWMsRUFBRSxjQUFjLEVBQUUsQ0FBQyxjQUFjO1FBQy9DLENBQUMsY0FBYyxFQUFFLGFBQWEsRUFBRSxZQUFZO1FBQzVDLENBQUMsWUFBWSxFQUFFLFdBQVcsRUFBRSxXQUFXO1FBQ3ZDLENBQUMsV0FBVyxFQUFFLFFBQVEsRUFBRSxTQUFTO1FBQ2pDLENBQUMsU0FBUyxFQUFFLENBQUMsT0FBTyxFQUFFLFFBQVE7UUFDOUIsQ0FBQyxPQUFPLENBQUMsQ0FBQztJQUNaLElBQUksQ0FBQyxHQUFHLEdBQUcsQ0FBQyxNQUFNLEdBQUcsQ0FBQyxDQUFDO0lBQ3ZCLElBQUksS0FBSyxHQUFHLEtBQUssQ0FBQztJQUNsQixJQUFJLENBQUMsR0FBRyxDQUFDLENBQUM7SUFDVixJQUFJLEVBQUUsR0FBRyxDQUFDLENBQUM7SUFDWCxJQUFJLENBQUMsR0FBVyxDQUFDLENBQUM7SUFBQyxJQUFJLEVBQUUsR0FBVyxDQUFDLENBQUM7SUFBQyxJQUFJLEdBQUcsR0FBVyxDQUFDLENBQUM7SUFBQyxJQUFJLEdBQUcsR0FBVyxDQUFDLENBQUM7SUFFaEYsSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFO1FBQ1QsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1FBQ1AsS0FBSyxHQUFHLElBQUksQ0FBQztLQUNkO0lBRUQsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQztJQUNoQixFQUFFLEdBQUcsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUM7SUFFZixPQUFPLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxFQUFFLEVBQUU7UUFDakIsR0FBRyxHQUFHLENBQUMsQ0FBQztRQUNSLENBQUMsR0FBRyxFQUFFLEdBQUcsQ0FBQyxHQUFHLEVBQUUsR0FBRyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDekIsRUFBRSxHQUFHLEdBQUcsQ0FBQztLQUNWO0lBRUQsR0FBRyxHQUFHLENBQUMsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsR0FBRyxHQUFHLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsRUFBRSxHQUFHLENBQUMsQ0FBQyxHQUFHLEVBQUUsQ0FBQyxDQUFDO0lBQzFELE9BQU8sS0FBSyxDQUFDLENBQUMsQ0FBQyxHQUFHLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsR0FBRyxDQUFDO0FBQ25DLENBQUM7QUFFRDs7Ozs7O0dBTUc7QUFDSCxTQUFTLElBQUksQ0FBQyxDQUFTO0lBQ3JCLE9BQU8sQ0FBQyxHQUFHLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztBQUNwQixDQUFDO0FBRUQ7Ozs7Ozs7Ozs7O0dBV0c7QUFDSCxNQUFNLFVBQVUsV0FBVyxDQUFDLENBQVcsRUFBRSxDQUFXO0lBQ2xELElBQUksRUFBRSxHQUFHLENBQUMsQ0FBQztJQUNYLElBQUksRUFBRSxHQUFHLENBQUMsQ0FBQztJQUNYLElBQUksR0FBRyxHQUFHLENBQUMsQ0FBQztJQUNaLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxDQUFDLENBQUMsTUFBTSxHQUFHLENBQUMsRUFBRSxDQUFDLEVBQUUsRUFBRTtRQUNyQyxLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxHQUFHLENBQUMsQ0FBQyxNQUFNLEVBQUUsQ0FBQyxFQUFFLEVBQUU7WUFDckMsTUFBTSxFQUFFLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztZQUN2QixNQUFNLEVBQUUsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO1lBQ3ZCLE1BQU0sRUFBRSxHQUFHLEVBQUUsR0FBRyxFQUFFLENBQUM7WUFDbkIsSUFBSSxFQUFFLEVBQUU7Z0JBQ04sMEJBQTBCO2dCQUMxQixFQUFFLEdBQUcsRUFBRSxHQUFHLENBQUMsQ0FBQztnQkFDWixFQUFFLEdBQUcsRUFBRSxHQUFHLENBQUMsQ0FBQztnQkFDWixJQUFJLEVBQUUsR0FBRyxDQUFDLEVBQUU7b0JBQ1YsR0FBRyxHQUFHLEdBQUcsR0FBRyxDQUFDLENBQUM7aUJBQ2Y7cUJBQU07b0JBQ0wsR0FBRyxHQUFHLEdBQUcsR0FBRyxDQUFDLENBQUM7aUJBQ2Y7YUFDRjtpQkFBTTtnQkFDTCxJQUFJLEVBQUUsRUFBRTtvQkFDTixFQUFFLEdBQUcsRUFBRSxHQUFHLENBQUMsQ0FBQztpQkFDYjtnQkFDRCxJQUFJLEVBQUUsRUFBRTtvQkFDTixFQUFFLEdBQUcsRUFBRSxHQUFHLENBQUMsQ0FBQztpQkFDYjthQUNGO1NBQ0Y7S0FDRjtJQUNELE1BQU0sR0FBRyxHQUFHLEdBQUcsR0FBRyxJQUFJLENBQUMsSUFBSSxDQUFDLEVBQUUsR0FBRyxFQUFFLENBQUMsQ0FBQztJQUNyQyxNQUFNLElBQUksR0FBRyxDQUFDLEdBQUcsR0FBRyxDQUFDLENBQUMsTUFBTSxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsR0FBRyxHQUFHLENBQUMsQ0FBQyxNQUFNLEdBQUcsQ0FBQyxDQUFDLENBQUMsTUFBTSxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUM7SUFDekUsTUFBTSxDQUFDLEdBQUcsR0FBRyxHQUFHLElBQUksQ0FBQyxJQUFJLENBQUMsSUFBSSxDQUFDLENBQUM7SUFDaEMsTUFBTSxJQUFJLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsU0FBUyxDQUFDLENBQUM7SUFDM0MsT0FBTyxFQUFDLElBQUksRUFBRSxHQUFHLEVBQUUsQ0FBQyxFQUFFLENBQUMsRUFBRSxJQUFJLEVBQUUsSUFBSSxFQUFDLENBQUM7QUFDdkMsQ0FBQyIsInNvdXJjZXNDb250ZW50IjpbImltcG9ydCB7T3B0aW9uc30gZnJvbSAnLi4vLi4vdXRpbHMvc3JjL3R5cGUtZGVjbGFyYXRpb25zJztcbi8vaW1wb3J0IHtPcHRpb25zfSBmcm9tICdAZGF0YWdyb2stbGlicmFyaWVzL3V0aWxzL3NyYy90eXBlLWRlY2xhcmF0aW9ucyc7XG5cbi8qKlxuICogUmV0dXJucyB0aGUgZXJyb3IgZnVuY3Rpb24gZXJmKHgpLlxuICpcbiAqIEBwYXJhbSB7bnVtYmVyfSB4IEFuIGFyZ3VtZW50LlxuICogQHJldHVybiB7bnVtYmVyfSBUaGUgcmVzdWx0LlxuICogQGxpbmsgaHR0cHM6Ly9naXRodWIuY29tL2pzdGF0L2pzdGF0L2Jsb2IvNjVjZTA5NmE5OWY3NTNkNmEyMjQ4MmU1ZTc0YWNjYmZjMWMzMzc2Ny9kaXN0L2pzdGF0LmpzI0wxNTYyXG4gKi9cbmZ1bmN0aW9uIGVyZih4OiBudW1iZXIpOiBudW1iZXIge1xuICBjb25zdCBjb2YgPSBbLTEuMzAyNjUzNzE5NzgxNzA5NCwgNi40MTk2OTc5MjM1NjQ5MDI2ZS0xLCAxLjk0NzY0NzMyMDQxODU4MzZlLTIsXG4gICAgLTkuNTYxNTE0Nzg2ODA4NjMxZS0zLCAtOS40NjU5NTM0NDQ4MjAzNmUtNCwgMy42NjgzOTQ5Nzg1Mjc2MWUtNCxcbiAgICA0LjI1MjMzMjQ4MDY5MDdlLTUsIC0yLjAyNzg1NzgxMTI1MzRlLTUsIC0xLjYyNDI5MDAwNDY0N2UtNixcbiAgICAxLjMwMzY1NTgzNTU4MGUtNiwgMS41NjI2NDQxNzIyZS04LCAtOC41MjM4MDk1OTE1ZS04LFxuICAgIDYuNTI5MDU0NDM5ZS05LCA1LjA1OTM0MzQ5NWUtOSwgLTkuOTEzNjQxNTZlLTEwLFxuICAgIC0yLjI3MzY1MTIyZS0xMCwgOS42NDY3OTExZS0xMSwgMi4zOTQwMzhlLTEyLFxuICAgIC02Ljg4NjAyN2UtMTIsIDguOTQ0ODdlLTEzLCAzLjEzMDkyZS0xMyxcbiAgICAtMS4xMjcwOGUtMTMsIDMuODFlLTE2LCA3LjEwNmUtMTUsXG4gICAgLTEuNTIzZS0xNSwgLTkuNGUtMTcsIDEuMjFlLTE2LFxuICAgIC0yLjhlLTE3XTtcbiAgbGV0IGogPSBjb2YubGVuZ3RoIC0gMTtcbiAgbGV0IGlzbmVnID0gZmFsc2U7XG4gIGxldCBkID0gMDtcbiAgbGV0IGRkID0gMDtcbiAgbGV0IHQ6IG51bWJlciA9IDA7IGxldCB0eTogbnVtYmVyID0gMDsgbGV0IHRtcDogbnVtYmVyID0gMDsgbGV0IHJlczogbnVtYmVyID0gMDtcblxuICBpZiAoeCA8IDApIHtcbiAgICB4ID0gLXg7XG4gICAgaXNuZWcgPSB0cnVlO1xuICB9XG5cbiAgdCA9IDIgLyAoMiArIHgpO1xuICB0eSA9IDQgKiB0IC0gMjtcblxuICBmb3IgKDsgaiA+IDA7IGotLSkge1xuICAgIHRtcCA9IGQ7XG4gICAgZCA9IHR5ICogZCAtIGRkICsgY29mW2pdO1xuICAgIGRkID0gdG1wO1xuICB9XG5cbiAgcmVzID0gdCAqIE1hdGguZXhwKC14ICogeCArIDAuNSAqIChjb2ZbMF0gKyB0eSAqIGQpIC0gZGQpO1xuICByZXR1cm4gaXNuZWcgPyByZXMgLSAxIDogMSAtIHJlcztcbn1cblxuLyoqXG4gKiBSZXR1cm5zIHRoZSBjb21wbG1lbnRhcnkgZXJyb3IgZnVuY3Rpb24gZXJmYyh4KVxuICpcbiAqIEBwYXJhbSB7bnVtYmVyfSB4IEFuIGFyZ3VtZW50LlxuICogQHJldHVybiB7bnVtYmVyfSBUaGUgcmVzdWx0LlxuICogQGxpbmsgaHR0cHM6Ly9naXRodWIuY29tL2pzdGF0L2pzdGF0L2Jsb2IvNjVjZTA5NmE5OWY3NTNkNmEyMjQ4MmU1ZTc0YWNjYmZjMWMzMzc2Ny9kaXN0L2pzdGF0LmpzI0wxNTk5XG4gKi9cbmZ1bmN0aW9uIGVyZmMoeDogbnVtYmVyKTogbnVtYmVyIHtcbiAgcmV0dXJuIDEgLSBlcmYoeCk7XG59XG5cbi8qKlxuICogQ2FsY3VsYXRlcyBLZW5kYWxsJ3MgdGF1LCBhIGNvcnJlbGF0aW9uIG1lYXN1cmUgZm9yIG9yZGluYWwgZGF0YSwgYW5kIGFuIGFzc29jaWF0ZWQgcC12YWx1ZS5cbiAqIFJldHVybnM6IEtlbmRhbGwncyB0YXUsIHR3by10YWlsZWQgcC12YWx1ZS5cbiAqIERlcml2ZWQgZnJvbSBvbGRlciBTY2lQeTogaHR0cDovL3dlYi5taXQuZWR1LzYuODYzL3NwcmluZzIwMTEvcGFja2FnZXMvc2NpcHlfc3JjL3NjaXB5L3N0YXRzL3N0YXRzLnB5XG4gKlxuICogQGV4cG9ydFxuICogQHBhcmFtIHtudW1iZXJbXX0geCBUaGUgZmlyc3QgYXJyYXkuXG4gKiBAcGFyYW0ge251bWJlcltdfSB5IFRoZSBzZWNvbmQgYXJyYXkuXG4gKiBAcmV0dXJuIHtPcHRpb25zfSBUaGUgcmVzdWx0LlxuICogQGxpbmsgaHR0cHM6Ly9naXRodWIuY29tL3BkZmVybmhvdXQvbmFycmFmaXJtYS9ibG9iL1xuICogYzljMTIyZDU3N2E0Yjg4NjhjZTYwM2JhZGFiYmI3ZDEwZjQ1NzQwYy93ZWJhcHAvc291cmNlL3N0YXRpc3RpY3Mva2VuZGFsbHNUYXUudHMjTDhcbiAqL1xuZXhwb3J0IGZ1bmN0aW9uIGtlbmRhbGxzVGF1KHg6IG51bWJlcltdLCB5OiBudW1iZXJbXSk6IE9wdGlvbnMge1xuICBsZXQgbjEgPSAwO1xuICBsZXQgbjIgPSAwO1xuICBsZXQgaXNzID0gMDtcbiAgZm9yIChsZXQgaiA9IDA7IGogPCB4Lmxlbmd0aCAtIDE7IGorKykge1xuICAgIGZvciAobGV0IGsgPSBqICsgMTsgayA8IHkubGVuZ3RoOyBrKyspIHtcbiAgICAgIGNvbnN0IGExID0geFtqXSAtIHhba107XG4gICAgICBjb25zdCBhMiA9IHlbal0gLSB5W2tdO1xuICAgICAgY29uc3QgYWEgPSBhMSAqIGEyO1xuICAgICAgaWYgKGFhKSB7XG4gICAgICAgIC8vIG5laXRoZXIgYXJyYXkgaGFzIGEgdGllXG4gICAgICAgIG4xID0gbjEgKyAxO1xuICAgICAgICBuMiA9IG4yICsgMTtcbiAgICAgICAgaWYgKGFhID4gMCkge1xuICAgICAgICAgIGlzcyA9IGlzcyArIDE7XG4gICAgICAgIH0gZWxzZSB7XG4gICAgICAgICAgaXNzID0gaXNzIC0gMTtcbiAgICAgICAgfVxuICAgICAgfSBlbHNlIHtcbiAgICAgICAgaWYgKGExKSB7XG4gICAgICAgICAgbjEgPSBuMSArIDE7XG4gICAgICAgIH1cbiAgICAgICAgaWYgKGEyKSB7XG4gICAgICAgICAgbjIgPSBuMiArIDE7XG4gICAgICAgIH1cbiAgICAgIH1cbiAgICB9XG4gIH1cbiAgY29uc3QgdGF1ID0gaXNzIC8gTWF0aC5zcXJ0KG4xICogbjIpO1xuICBjb25zdCBzdmFyID0gKDQuMCAqIHgubGVuZ3RoICsgMTAuMCkgLyAoOS4wICogeC5sZW5ndGggKiAoeC5sZW5ndGggLSAxKSk7XG4gIGNvbnN0IHogPSB0YXUgLyBNYXRoLnNxcnQoc3Zhcik7XG4gIGNvbnN0IHByb2IgPSBlcmZjKE1hdGguYWJzKHopIC8gMS40MTQyMTM2KTtcbiAgcmV0dXJuIHt0ZXN0OiB0YXUsIHo6IHosIHByb2I6IHByb2J9O1xufVxuIl19
|
package/src/multiple-tests.d.ts
CHANGED
|
@@ -1,35 +1,35 @@
|
|
|
1
|
-
/**
|
|
2
|
-
* pvalue correction for false discovery rate
|
|
3
|
-
*
|
|
4
|
-
* @export
|
|
5
|
-
* @param {Float32Array} pvals Set of p-values of the individual tests.
|
|
6
|
-
* @param {number} [alpha=0.05] Family-wise error rate. Defaults to 0.05.
|
|
7
|
-
* @param {string} [method='n'] {'i', 'indep', 'p', 'poscorr', 'n', 'negcorr'}, optional
|
|
8
|
-
* Which method to use for FDR correction.
|
|
9
|
-
* ``{'i', 'indep', 'p', 'poscorr'}`` all refer to ``fdr_bh``
|
|
10
|
-
* (Benjamini/Hochberg for independent or positively
|
|
11
|
-
* correlated tests). ``{'n', 'negcorr'}`` both refer to ``fdr_by``
|
|
12
|
-
* (Benjamini/Yekutieli for general or negatively correlated tests).
|
|
13
|
-
* Defaults to 'n'.
|
|
14
|
-
* @param {boolean} [isSorted=false] If False (default), the p_values will be sorted, but the corrected
|
|
15
|
-
* pvalues are in the original order. If True, then it assumed that the
|
|
16
|
-
* pvalues are already sorted in ascending order.
|
|
17
|
-
* @return {[Array<boolean>, Float32Array]} rejected : ndarray, bool
|
|
18
|
-
* True if a hypothesis is rejected, False if not
|
|
19
|
-
* pvalue-corrected : ndarray
|
|
20
|
-
* pvalues adjusted for multiple hypothesis testing to limit FDR
|
|
21
|
-
* @see
|
|
22
|
-
* If there is prior information on the fraction of true hypothesis, then alpha
|
|
23
|
-
* should be set to ``alpha * m/m_0`` where m is the number of tests,
|
|
24
|
-
* given by the p-values, and m_0 is an estimate of the true hypothesis.
|
|
25
|
-
* (see Benjamini, Krieger and Yekuteli)
|
|
26
|
-
*
|
|
27
|
-
* The two-step method of Benjamini, Krieger and Yekutiel that estimates the number
|
|
28
|
-
* of false hypotheses will be available (soon).
|
|
29
|
-
*
|
|
30
|
-
* Both methods exposed via this function (Benjamini/Hochberg, Benjamini/Yekutieli)
|
|
31
|
-
* are also available in the function ``multipletests``, as ``method="fdr_bh"`` and
|
|
32
|
-
* ``method="fdr_by"``, respectively.
|
|
33
|
-
*/
|
|
34
|
-
export declare function fdrcorrection(pvals: Float32Array, alpha?: number, method?: string, isSorted?: boolean): [Array<boolean>, Float32Array];
|
|
1
|
+
/**
|
|
2
|
+
* pvalue correction for false discovery rate
|
|
3
|
+
*
|
|
4
|
+
* @export
|
|
5
|
+
* @param {Float32Array} pvals Set of p-values of the individual tests.
|
|
6
|
+
* @param {number} [alpha=0.05] Family-wise error rate. Defaults to 0.05.
|
|
7
|
+
* @param {string} [method='n'] {'i', 'indep', 'p', 'poscorr', 'n', 'negcorr'}, optional
|
|
8
|
+
* Which method to use for FDR correction.
|
|
9
|
+
* ``{'i', 'indep', 'p', 'poscorr'}`` all refer to ``fdr_bh``
|
|
10
|
+
* (Benjamini/Hochberg for independent or positively
|
|
11
|
+
* correlated tests). ``{'n', 'negcorr'}`` both refer to ``fdr_by``
|
|
12
|
+
* (Benjamini/Yekutieli for general or negatively correlated tests).
|
|
13
|
+
* Defaults to 'n'.
|
|
14
|
+
* @param {boolean} [isSorted=false] If False (default), the p_values will be sorted, but the corrected
|
|
15
|
+
* pvalues are in the original order. If True, then it assumed that the
|
|
16
|
+
* pvalues are already sorted in ascending order.
|
|
17
|
+
* @return {[Array<boolean>, Float32Array]} rejected : ndarray, bool
|
|
18
|
+
* True if a hypothesis is rejected, False if not
|
|
19
|
+
* pvalue-corrected : ndarray
|
|
20
|
+
* pvalues adjusted for multiple hypothesis testing to limit FDR
|
|
21
|
+
* @see
|
|
22
|
+
* If there is prior information on the fraction of true hypothesis, then alpha
|
|
23
|
+
* should be set to ``alpha * m/m_0`` where m is the number of tests,
|
|
24
|
+
* given by the p-values, and m_0 is an estimate of the true hypothesis.
|
|
25
|
+
* (see Benjamini, Krieger and Yekuteli)
|
|
26
|
+
*
|
|
27
|
+
* The two-step method of Benjamini, Krieger and Yekutiel that estimates the number
|
|
28
|
+
* of false hypotheses will be available (soon).
|
|
29
|
+
*
|
|
30
|
+
* Both methods exposed via this function (Benjamini/Hochberg, Benjamini/Yekutieli)
|
|
31
|
+
* are also available in the function ``multipletests``, as ``method="fdr_bh"`` and
|
|
32
|
+
* ``method="fdr_by"``, respectively.
|
|
33
|
+
*/
|
|
34
|
+
export declare function fdrcorrection(pvals: Float32Array, alpha?: number, method?: string, isSorted?: boolean): [Array<boolean>, Float32Array];
|
|
35
35
|
//# sourceMappingURL=multiple-tests.d.ts.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"multiple-tests.d.ts","sourceRoot":"","sources":["multiple-tests.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"multiple-tests.d.ts","sourceRoot":"","sources":["multiple-tests.ts"],"names":[],"mappings":"AAoIA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAgCG;AACH,wBAAgB,aAAa,CAC3B,KAAK,EAAE,YAAY,EACnB,KAAK,GAAE,MAAa,EACpB,MAAM,GAAE,MAAW,EACnB,QAAQ,GAAE,OAAe,GAAG,CAAC,KAAK,CAAC,OAAO,CAAC,EAAE,YAAY,CAAC,CAyD3D"}
|