Skip to content

Commit 3437c9c

Browse files
committed
feat: add protocol support to stats/base/varianceyc
--- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown status: passed - task: lint_package_json status: na - task: lint_repl_help status: passed - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: passed - task: lint_javascript_tests status: passed - task: lint_javascript_benchmarks status: passed - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: na - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: passed - task: lint_license_headers status: passed ---
1 parent b688f77 commit 3437c9c

File tree

13 files changed

+429
-139
lines changed

13 files changed

+429
-139
lines changed

lib/node_modules/@stdlib/stats/base/varianceyc/README.md

+16-24
Original file line numberDiff line numberDiff line change
@@ -98,7 +98,7 @@ The use of the term `n-1` is commonly referred to as Bessel's correction. Note,
9898
var varianceyc = require( '@stdlib/stats/base/varianceyc' );
9999
```
100100

101-
#### varianceyc( N, correction, x, stride )
101+
#### varianceyc( N, correction, x, strideX )
102102

103103
Computes the [variance][variance] of a strided array `x` using a one-pass algorithm proposed by Youngs and Cramer.
104104

@@ -114,17 +114,16 @@ The function has the following parameters:
114114
- **N**: number of indexed elements.
115115
- **correction**: degrees of freedom adjustment. Setting this parameter to a value other than `0` has the effect of adjusting the divisor during the calculation of the [variance][variance] according to `N-c` where `c` corresponds to the provided degrees of freedom adjustment. When computing the [variance][variance] of a population, setting this parameter to `0` is the standard choice (i.e., the provided array contains data constituting an entire population). When computing the unbiased sample [variance][variance], setting this parameter to `1` is the standard choice (i.e., the provided array contains data sampled from a larger population; this is commonly referred to as Bessel's correction).
116116
- **x**: input [`Array`][mdn-array] or [`typed array`][mdn-typed-array].
117-
- **stride**: index increment for `x`.
117+
- **strideX**: stride length for `x`.
118118

119-
The `N` and `stride` parameters determine which elements in `x` are accessed at runtime. For example, to compute the [variance][variance] of every other element in `x`,
119+
The `N` and `stride` parameters determine which elements in the strided array are accessed at runtime. For example, to compute the [variance][variance] of every other element in `x`,
120120

121121
```javascript
122122
var floor = require( '@stdlib/math/base/special/floor' );
123123

124124
var x = [ 1.0, 2.0, 2.0, -7.0, -2.0, 3.0, 4.0, 2.0 ];
125-
var N = floor( x.length / 2 );
126125

127-
var v = varianceyc( N, 1, x, 2 );
126+
var v = varianceyc( 4, 1, x, 2 );
128127
// returns 6.25
129128
```
130129

@@ -139,13 +138,11 @@ var floor = require( '@stdlib/math/base/special/floor' );
139138
var x0 = new Float64Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0 ] );
140139
var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element
141140

142-
var N = floor( x0.length / 2 );
143-
144-
var v = varianceyc( N, 1, x1, 2 );
141+
var v = varianceyc( 4, 1, x1, 2 );
145142
// returns 6.25
146143
```
147144

148-
#### varianceyc.ndarray( N, correction, x, stride, offset )
145+
#### varianceyc.ndarray( N, correction, x, strideX, offsetX )
149146

150147
Computes the [variance][variance] of a strided array using a one-pass algorithm proposed by Youngs and Cramer and alternative indexing semantics.
151148

@@ -158,17 +155,14 @@ var v = varianceyc.ndarray( x.length, 1, x, 1, 0 );
158155

159156
The function has the following additional parameters:
160157

161-
- **offset**: starting index for `x`.
158+
- **offsetX**: starting index for `x`.
162159

163-
While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying `buffer`, the `offset` parameter supports indexing semantics based on a starting index. For example, to calculate the [variance][variance] for every other value in `x` starting from the second value
160+
While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameter supports indexing semantics based on a starting index. For example, to calculate the [variance][variance] for every other value in `x` starting from the second value
164161

165162
```javascript
166-
var floor = require( '@stdlib/math/base/special/floor' );
167-
168163
var x = [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0 ];
169-
var N = floor( x.length / 2 );
170164

171-
var v = varianceyc.ndarray( N, 1, x, 2, 1 );
165+
var v = varianceyc.ndarray( 4, 1, x, 2, 1 );
172166
// returns 6.25
173167
```
174168

@@ -183,6 +177,7 @@ var v = varianceyc.ndarray( N, 1, x, 2, 1 );
183177
- If `N <= 0`, both functions return `NaN`.
184178
- If `N - c` is less than or equal to `0` (where `c` corresponds to the provided degrees of freedom adjustment), both functions return `NaN`.
185179
- Depending on the environment, the typed versions ([`dvarianceyc`][@stdlib/stats/base/dvarianceyc], [`svarianceyc`][@stdlib/stats/base/svarianceyc], etc.) are likely to be significantly more performant.
180+
- Both functions support array-like objects having getter and setter accessors for array element access (e.g., [`@stdlib/array/base/accessor`][@stdlib/array/base/accessor]).
186181

187182
</section>
188183

@@ -195,18 +190,13 @@ var v = varianceyc.ndarray( N, 1, x, 2, 1 );
195190
<!-- eslint no-undef: "error" -->
196191

197192
```javascript
198-
var randu = require( '@stdlib/random/base/randu' );
199-
var round = require( '@stdlib/math/base/special/round' );
193+
var uniform = require( '@stdlib/random/array/uniform' );
200194
var Float64Array = require( '@stdlib/array/float64' );
201195
var varianceyc = require( '@stdlib/stats/base/varianceyc' );
202196

203-
var x;
204-
var i;
205-
206-
x = new Float64Array( 10 );
207-
for ( i = 0; i < x.length; i++ ) {
208-
x[ i ] = round( (randu()*100.0) - 50.0 );
209-
}
197+
var x = uniform( 10, -50.0, 50.0, {
198+
'dtype': 'float64'
199+
});
210200
console.log( x );
211201

212202
var v = varianceyc( x.length, 1, x, 1 );
@@ -256,6 +246,8 @@ console.log( v );
256246

257247
[mdn-typed-array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/TypedArray
258248

249+
[@stdlib/array/base/accessor]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/array/base/accessor
250+
259251
[@stdlib/stats/base/svarianceyc]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/stats/base/svarianceyc
260252

261253
[@youngs:1971a]: https://doi.org/10.1080/00401706.1971.10488826

lib/node_modules/@stdlib/stats/base/varianceyc/benchmark/benchmark.js

+9-8
Original file line numberDiff line numberDiff line change
@@ -21,13 +21,20 @@
2121
// MODULES //
2222

2323
var bench = require( '@stdlib/bench' );
24-
var randu = require( '@stdlib/random/base/randu' );
2524
var isnan = require( '@stdlib/math/base/assert/is-nan' );
25+
var uniform = require( '@stdlib/random/array/uniform' );
2626
var pow = require( '@stdlib/math/base/special/pow' );
2727
var pkg = require( './../package.json' ).name;
2828
var varianceyc = require( './../lib/varianceyc.js' );
2929

3030

31+
// VARIABLES //
32+
33+
var options = {
34+
'dtype': 'generic'
35+
};
36+
37+
3138
// FUNCTIONS //
3239

3340
/**
@@ -38,13 +45,7 @@ var varianceyc = require( './../lib/varianceyc.js' );
3845
* @returns {Function} benchmark function
3946
*/
4047
function createBenchmark( len ) {
41-
var x;
42-
var i;
43-
44-
x = [];
45-
for ( i = 0; i < len; i++ ) {
46-
x.push( ( randu()*20.0 ) - 10.0 );
47-
}
48+
var x = uniform( len, -10.0, 10.0, options );
4849
return benchmark;
4950

5051
function benchmark( b ) {

lib/node_modules/@stdlib/stats/base/varianceyc/benchmark/benchmark.ndarray.js

+9-8
Original file line numberDiff line numberDiff line change
@@ -21,13 +21,20 @@
2121
// MODULES //
2222

2323
var bench = require( '@stdlib/bench' );
24-
var randu = require( '@stdlib/random/base/randu' );
24+
var uniform = require( '@stdlib/random/array/uniform' );
2525
var isnan = require( '@stdlib/math/base/assert/is-nan' );
2626
var pow = require( '@stdlib/math/base/special/pow' );
2727
var pkg = require( './../package.json' ).name;
2828
var varianceyc = require( './../lib/ndarray.js' );
2929

3030

31+
// VARIABLES //
32+
33+
var options = {
34+
'dtype': 'generic'
35+
};
36+
37+
3138
// FUNCTIONS //
3239

3340
/**
@@ -38,13 +45,7 @@ var varianceyc = require( './../lib/ndarray.js' );
3845
* @returns {Function} benchmark function
3946
*/
4047
function createBenchmark( len ) {
41-
var x;
42-
var i;
43-
44-
x = [];
45-
for ( i = 0; i < len; i++ ) {
46-
x.push( ( randu()*20.0 ) - 10.0 );
47-
}
48+
var x = uniform( len, -10.0, 10.0, options );
4849
return benchmark;
4950

5051
function benchmark( b ) {

lib/node_modules/@stdlib/stats/base/varianceyc/docs/repl.txt

+15-19
Original file line numberDiff line numberDiff line change
@@ -1,10 +1,10 @@
11

2-
{{alias}}( N, correction, x, stride )
2+
{{alias}}( N, correction, x, strideX )
33
Computes the variance of a strided array using a one-pass algorithm proposed
44
by Youngs and Cramer.
55

6-
The `N` and `stride` parameters determine which elements in `x` are accessed
7-
at runtime.
6+
The `N` and stride parameters determine which elements in the strided array
7+
are accessed at runtime.
88

99
Indexing is relative to the first index. To introduce an offset, use a typed
1010
array view.
@@ -31,8 +31,8 @@
3131
x: Array<number>|TypedArray
3232
Input array.
3333

34-
stride: integer
35-
Index increment.
34+
strideX: integer
35+
Stride length.
3636

3737
Returns
3838
-------
@@ -46,27 +46,24 @@
4646
> {{alias}}( x.length, 1, x, 1 )
4747
~4.3333
4848

49-
// Using `N` and `stride` parameters:
49+
// Using `N` and stride parameters:
5050
> x = [ -2.0, 1.0, 1.0, -5.0, 2.0, -1.0 ];
51-
> var N = {{alias:@stdlib/math/base/special/floor}}( x.length / 2 );
52-
> var stride = 2;
53-
> {{alias}}( N, 1, x, stride )
51+
> {{alias}}( 3, 1, x, 2 )
5452
~4.3333
5553

5654
// Using view offsets:
5755
> var x0 = new {{alias:@stdlib/array/float64}}( [ 1.0, -2.0, 3.0, 2.0, 5.0, -1.0 ] );
5856
> var x1 = new {{alias:@stdlib/array/float64}}( x0.buffer, x0.BYTES_PER_ELEMENT*1 );
59-
> N = {{alias:@stdlib/math/base/special/floor}}( x0.length / 2 );
60-
> stride = 2;
61-
> {{alias}}( N, 1, x1, stride )
57+
> {{alias}}( 3, 1, x1, 2 )
6258
~4.3333
6359

64-
{{alias}}.ndarray( N, correction, x, stride, offset )
60+
61+
{{alias}}.ndarray( N, correction, x, strideX, offsetX )
6562
Computes the variance of a strided array using a one-pass algorithm proposed
6663
by Youngs and Cramer and alternative indexing semantics.
6764

6865
While typed array views mandate a view offset based on the underlying
69-
buffer, the `offset` parameter supports indexing semantics based on a
66+
buffer, the offset parameter supports indexing semantics based on a
7067
starting index.
7168

7269
Parameters
@@ -89,10 +86,10 @@
8986
x: Array<number>|TypedArray
9087
Input array.
9188

92-
stride: integer
93-
Index increment.
89+
strideX: integer
90+
Stride length.
9491

95-
offset: integer
92+
offsetX: integer
9693
Starting index.
9794

9895
Returns
@@ -109,8 +106,7 @@
109106

110107
// Using offset parameter:
111108
> var x = [ 1.0, -2.0, 3.0, 2.0, 5.0, -1.0 ];
112-
> var N = {{alias:@stdlib/math/base/special/floor}}( x.length / 2 );
113-
> {{alias}}.ndarray( N, 1, x, 2, 1 )
109+
> {{alias}}.ndarray( 3, 1, x, 2, 1 )
114110
~4.3333
115111

116112
See Also

lib/node_modules/@stdlib/stats/base/varianceyc/docs/types/index.d.ts

+14-7
Original file line numberDiff line numberDiff line change
@@ -20,7 +20,14 @@
2020

2121
/// <reference types="@stdlib/types"/>
2222

23-
import { NumericArray } from '@stdlib/types/array';
23+
import { NumericArray, Collection, AccessorArrayLike } from '@stdlib/types/array';
24+
25+
26+
/**
27+
* Input array.
28+
*/
29+
type InputArray = NumericArray | Collection<number> | AccessorArrayLike<number>;
30+
2431

2532
/**
2633
* Interface describing `varianceyc`.
@@ -32,7 +39,7 @@ interface Routine {
3239
* @param N - number of indexed elements
3340
* @param correction - degrees of freedom adjustment
3441
* @param x - input array
35-
* @param stride - stride length
42+
* @param strideX - stride length for `x`
3643
* @returns variance
3744
*
3845
* @example
@@ -41,16 +48,16 @@ interface Routine {
4148
* var v = varianceyc( x.length, 1, x, 1 );
4249
* // returns ~4.3333
4350
*/
44-
( N: number, correction: number, x: NumericArray, stride: number ): number;
51+
( N: number, correction: number, x: InputArray, stride: number ): number;
4552

4653
/**
4754
* Computes the variance of a strided array using a one-pass algorithm proposed by Youngs and Cramer and alternative indexing semantics.
4855
*
4956
* @param N - number of indexed elements
5057
* @param correction - degrees of freedom adjustment
5158
* @param x - input array
52-
* @param stride - stride length
53-
* @param offset - starting index
59+
* @param strideX - stride length for `x`
60+
* @param offsetX - starting index
5461
* @returns variance
5562
*
5663
* @example
@@ -59,7 +66,7 @@ interface Routine {
5966
* var v = varianceyc.ndarray( x.length, 1, x, 1, 0 );
6067
* // returns ~4.3333
6168
*/
62-
ndarray( N: number, correction: number, x: NumericArray, stride: number, offset: number ): number;
69+
ndarray( N: number, correction: number, x: InputArray, stride: number, offset: number ): number;
6370
}
6471

6572
/**
@@ -68,7 +75,7 @@ interface Routine {
6875
* @param N - number of indexed elements
6976
* @param correction - degrees of freedom adjustment
7077
* @param x - input array
71-
* @param stride - stride length
78+
* @param strideX - stride length
7279
* @returns variance
7380
*
7481
* @example

lib/node_modules/@stdlib/stats/base/varianceyc/docs/types/test.ts

+3
Original file line numberDiff line numberDiff line change
@@ -16,6 +16,7 @@
1616
* limitations under the License.
1717
*/
1818

19+
import AccessorArray = require( '@stdlib/array/base/accessor' );
1920
import varianceyc = require( './index' );
2021

2122

@@ -26,6 +27,7 @@ import varianceyc = require( './index' );
2627
const x = new Float64Array( 10 );
2728

2829
varianceyc( x.length, 1, x, 1 ); // $ExpectType number
30+
varianceyc( x.length, 1, new AccessorArray( x ), 1 ); // $ExpectType number
2931
}
3032

3133
// The compiler throws an error if the function is provided a first argument which is not a number...
@@ -101,6 +103,7 @@ import varianceyc = require( './index' );
101103
const x = new Float64Array( 10 );
102104

103105
varianceyc.ndarray( x.length, 1, x, 1, 0 ); // $ExpectType number
106+
varianceyc.ndarray( x.length, 1, new AccessorArray( x ), 1, 0 ); // $ExpectType number
104107
}
105108

106109
// The compiler throws an error if the `ndarray` method is provided a first argument which is not a number...

lib/node_modules/@stdlib/stats/base/varianceyc/examples/index.js

+4-10
Original file line numberDiff line numberDiff line change
@@ -18,18 +18,12 @@
1818

1919
'use strict';
2020

21-
var randu = require( '@stdlib/random/base/randu' );
22-
var round = require( '@stdlib/math/base/special/round' );
23-
var Float64Array = require( '@stdlib/array/float64' );
21+
var uniform = require( '@stdlib/random/array/uniform' );
2422
var varianceyc = require( './../lib' );
2523

26-
var x;
27-
var i;
28-
29-
x = new Float64Array( 10 );
30-
for ( i = 0; i < x.length; i++ ) {
31-
x[ i ] = round( (randu()*100.0) - 50.0 );
32-
}
24+
var x = uniform( 10, -50.0, 50.0, {
25+
'dtype': 'float64'
26+
});
3327
console.log( x );
3428

3529
var v = varianceyc( x.length, 1, x, 1 );

0 commit comments

Comments
 (0)