diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/cdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/cdf/benchmark/benchmark.js
index 94cf026fe34d..1152bd27c6b2 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/cdf/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/cdf/benchmark/benchmark.js
@@ -21,8 +21,9 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var ceil = require( '@stdlib/math/base/special/ceil' );
-var randu = require( '@stdlib/random/base/randu' );
+var Float64Array = require( '@stdlib/array/float64' );
+var uniform = require( '@stdlib/random/base/uniform' );
+var discreteUniform = require( '@stdlib/random/base/discrete-uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
 var pkg = require( './../package.json' ).name;
@@ -32,18 +33,26 @@ var cdf = require( './../lib' );
 // MAIN //
 
 bench( pkg, function benchmark( b ) {
+	var len;
 	var d1;
 	var d2;
 	var x;
 	var y;
 	var i;
 
+	len = 100;
+	x = new Float64Array( len );
+	d1 = new Float64Array( len );
+	d2 = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( EPS, 100.0 );
+		d1[ i ] = discreteUniform( 1, 100 );
+		d2[ i ] = discreteUniform( 1, 100 );
+	}
+
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = ( randu()*100.0 ) + EPS;
-		d1 = ceil( randu()*100.0 );
-		d2 = ceil( randu()*100.0 );
-		y = cdf( x, d1, d2 );
+		y = cdf( x[ i % len ], d1[ i % len ], d2[ i % len ]);
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
@@ -58,6 +67,7 @@ bench( pkg, function benchmark( b ) {
 
 bench( pkg+':factory', function benchmark( b ) {
 	var mycdf;
+	var len;
 	var d1;
 	var d2;
 	var x;
@@ -67,11 +77,15 @@ bench( pkg+':factory', function benchmark( b ) {
 	d1 = 1.5;
 	d2 = 1.5;
 	mycdf = cdf.factory( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( EPS, 20.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = ( randu()*20.0 ) + EPS;
-		y = mycdf( x );
+		y = mycdf( x[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/ctor/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/ctor/benchmark/benchmark.js
index 00fe2679a7b8..92b979030ad4 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/ctor/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/ctor/benchmark/benchmark.js
@@ -21,7 +21,8 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
+var Float64Array = require( '@stdlib/array/float64' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
 var pkg = require( './../package.json' ).name;
@@ -32,15 +33,22 @@ var F = require( './../lib' );
 
 bench( pkg+'::instantiation', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
 	var i;
 
+	len = 100;
+	d1 = new Float64Array( len );
+	d2 = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( EPS, 10.0 );
+	}
+
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		d1 = ( randu() * 10.0 ) + EPS;
-		d2 = ( randu() * 10.0 ) + EPS;
-		dist = new F( d1, d2 );
+		dist = new F( d1[ i % len ], d2[ i % len ] );
 		if ( !( dist instanceof F ) ) {
 			b.fail( 'should return a distribution instance' );
 		}
@@ -81,6 +89,7 @@ bench( pkg+'::get:d1', function benchmark( b ) {
 
 bench( pkg+'::set:d1', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
 	var y;
@@ -89,12 +98,16 @@ bench( pkg+'::set:d1', function benchmark( b ) {
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	y = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		y[ i ] = uniform( EPS, 10.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		y = ( 10.0*randu() ) + EPS;
-		dist.d1 = y;
-		if ( dist.d1 !== y ) {
+		dist.d1 = y[ i % len ];
+		if ( dist.d1 !== y[ i % len ] ) {
 			b.fail( 'should return set value' );
 		}
 	}
@@ -134,6 +147,7 @@ bench( pkg+'::get:d2', function benchmark( b ) {
 
 bench( pkg+'::set:d2', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
 	var y;
@@ -142,12 +156,16 @@ bench( pkg+'::set:d2', function benchmark( b ) {
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	y = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		y[ i ] = uniform( EPS, 10.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		y = ( 10.0*randu() ) + EPS;
-		dist.d2 = y;
-		if ( dist.d2 !== y ) {
+		dist.d2 = y[ i % len ];
+		if ( dist.d2 !== y[ i % len ] ) {
 			b.fail( 'should return set value' );
 		}
 	}
@@ -161,18 +179,25 @@ bench( pkg+'::set:d2', function benchmark( b ) {
 
 bench( pkg+':entropy', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
+	var x;
 	var y;
 	var i;
 
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( EPS, 10.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		dist.d1 = ( 10.0*randu() ) + EPS;
+		dist.d1 = x[ i % len ];
 		y = dist.entropy;
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
@@ -188,18 +213,25 @@ bench( pkg+':entropy', function benchmark( b ) {
 
 bench( pkg+':kurtosis', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
+	var x;
 	var y;
 	var i;
 
 	d1 = 10.56;
 	d2 = 8.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( 8.0 + EPS, 20.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		dist.d1 = ( 10.0*randu() ) + 8.0 + EPS;
+		dist.d1 = x[ i % len ];
 		y = dist.kurtosis;
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
@@ -215,18 +247,25 @@ bench( pkg+':kurtosis', function benchmark( b ) {
 
 bench( pkg+':mean', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
+	var x;
 	var y;
 	var i;
 
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( EPS, 10.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		dist.d1 = ( 10.0*randu() ) + EPS;
+		dist.d1 = x[ i % len ];
 		y = dist.mean;
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
@@ -242,18 +281,25 @@ bench( pkg+':mean', function benchmark( b ) {
 
 bench( pkg+':mode', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
+	var x;
 	var y;
 	var i;
 
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( 2.0 + EPS, 10.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		dist.d1 = ( 10.0*randu() ) + 2.0 + EPS;
+		dist.d1 = x[ i % len ];
 		y = dist.mode;
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
@@ -269,18 +315,25 @@ bench( pkg+':mode', function benchmark( b ) {
 
 bench( pkg+':skewness', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
+	var x;
 	var y;
 	var i;
 
 	d1 = 10.56;
 	d2 = 6.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( 6.0 + EPS, 20.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		dist.d1 = ( 10.0*randu() ) + 6.0 + EPS;
+		dist.d1 = x[ i % len ];
 		y = dist.skewness;
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
@@ -296,18 +349,25 @@ bench( pkg+':skewness', function benchmark( b ) {
 
 bench( pkg+':stdev', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
+	var x;
 	var y;
 	var i;
 
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( 1.0 + EPS, 10.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		dist.d1 = ( 10.0*randu() ) + 1.0 + EPS;
+		dist.d1 = x[ i % len ];
 		y = dist.stdev;
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
@@ -323,18 +383,25 @@ bench( pkg+':stdev', function benchmark( b ) {
 
 bench( pkg+':variance', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
+	var x;
 	var y;
 	var i;
 
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( 1.0 + EPS, 10.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		dist.d1 = ( 10.0*randu() ) + 1.0 + EPS;
+		dist.d1 = x[ i % len ];
 		y = dist.variance;
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
@@ -350,6 +417,7 @@ bench( pkg+':variance', function benchmark( b ) {
 
 bench( pkg+':cdf', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
 	var x;
@@ -359,11 +427,15 @@ bench( pkg+':cdf', function benchmark( b ) {
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( 0.0, 1.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = randu();
-		y = dist.cdf( x );
+		y = dist.cdf( x[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
@@ -378,6 +450,7 @@ bench( pkg+':cdf', function benchmark( b ) {
 
 bench( pkg+':pdf', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
 	var x;
@@ -387,11 +460,15 @@ bench( pkg+':pdf', function benchmark( b ) {
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( 0.0, 1.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = randu();
-		y = dist.pdf( x );
+		y = dist.pdf( x[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
@@ -406,6 +483,7 @@ bench( pkg+':pdf', function benchmark( b ) {
 
 bench( pkg+':quantile', function benchmark( b ) {
 	var dist;
+	var len;
 	var d1;
 	var d2;
 	var x;
@@ -415,11 +493,15 @@ bench( pkg+':quantile', function benchmark( b ) {
 	d1 = 10.56;
 	d2 = 5.54;
 	dist = new F( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( 0.0, 1.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = randu();
-		y = dist.quantile( x );
+		y = dist.quantile( x[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/entropy/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/entropy/benchmark/benchmark.js
index 06c413ab13c3..0d130c420f0a 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/entropy/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/entropy/benchmark/benchmark.js
@@ -21,7 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var Float64Array = require( '@stdlib/array/float64' );
 var EPS = require( '@stdlib/constants/float64/eps' );
@@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( EPS, 10.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/entropy/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/f/entropy/benchmark/benchmark.native.js
index 1bcda255a4a6..f7b36d56ef6f 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/entropy/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/entropy/benchmark/benchmark.native.js
@@ -23,7 +23,7 @@
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var Float64Array = require( '@stdlib/array/float64' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var EPS = require( '@stdlib/constants/float64/eps' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var tryRequire = require( '@stdlib/utils/try-require' );
@@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( EPS, 10.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/kurtosis/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/kurtosis/benchmark/benchmark.js
index 8c0386a78b51..95eb51d4aef8 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/kurtosis/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/kurtosis/benchmark/benchmark.js
@@ -21,7 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var Float64Array = require( '@stdlib/array/float64' );
 var EPS = require( '@stdlib/constants/float64/eps' );
@@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + 8.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 8.0 + EPS, 20.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/kurtosis/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/f/kurtosis/benchmark/benchmark.native.js
index c2fe90957168..98159644d06f 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/kurtosis/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/kurtosis/benchmark/benchmark.native.js
@@ -24,7 +24,7 @@ var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var Float64Array = require( '@stdlib/array/float64' );
 var EPS = require( '@stdlib/constants/float64/eps' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 var pkg = require( './../package.json' ).name;
@@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + 8.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 8.0 + EPS, 20.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/mean/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/mean/benchmark/benchmark.js
index 18e9a887c4c0..ab484d9766d0 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/mean/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/mean/benchmark/benchmark.js
@@ -21,7 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var Float64Array = require( '@stdlib/array/float64' );
 var EPS = require( '@stdlib/constants/float64/eps' );
@@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + 2.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 2.0 + EPS, 10.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/mean/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/f/mean/benchmark/benchmark.native.js
index db11d88dd8c8..40774e05a5ff 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/mean/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/mean/benchmark/benchmark.native.js
@@ -24,7 +24,7 @@ var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var EPS = require( '@stdlib/constants/float64/eps' );
 var Float64Array = require( '@stdlib/array/float64' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 var pkg = require( './../package.json' ).name;
@@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + 2.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 2.0 + EPS, 10.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/mode/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/mode/benchmark/benchmark.js
index 0b07733176df..2b8de9bea2fd 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/mode/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/mode/benchmark/benchmark.js
@@ -21,7 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var Float64Array = require( '@stdlib/array/float64' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
@@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + 2.0 + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + EPS;
+		d1[ i ] = uniform( 2.0 + EPS, 10.0 );
+		d2[ i ] = uniform( EPS, 10.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/mode/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/f/mode/benchmark/benchmark.native.js
index a7b9d7604b1c..6eebdd5b13be 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/mode/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/mode/benchmark/benchmark.native.js
@@ -24,7 +24,7 @@ var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var Float64Array = require( '@stdlib/array/float64' );
 var EPS = require( '@stdlib/constants/float64/eps' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 var pkg = require( './../package.json' ).name;
@@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + 2.0 + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + EPS;
+		d1[ i ] = uniform( 2.0 + EPS, 10.0 );
+		d2[ i ] = uniform( EPS, 10.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/pdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/pdf/benchmark/benchmark.js
index a0c9300a552e..4988b19e4e4c 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/pdf/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/pdf/benchmark/benchmark.js
@@ -21,8 +21,9 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var ceil = require( '@stdlib/math/base/special/ceil' );
-var randu = require( '@stdlib/random/base/randu' );
+var Float64Array = require( '@stdlib/array/float64' );
+var uniform = require( '@stdlib/random/base/uniform' );
+var discreteUniform = require( '@stdlib/random/base/discrete-uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
 var pkg = require( './../package.json' ).name;
@@ -32,18 +33,26 @@ var pdf = require( './../lib' );
 // MAIN //
 
 bench( pkg, function benchmark( b ) {
+	var len;
 	var d1;
 	var d2;
 	var x;
 	var y;
 	var i;
 
+	len = 100;
+	x = new Float64Array( len );
+	d1 = new Float64Array( len );
+	d2 = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( EPS, 100.0 );
+		d1[ i ] = discreteUniform( 1, 100 );
+		d2[ i ] = discreteUniform( 1, 100 );
+	}
+
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = ( randu()*100.0 ) + EPS;
-		d1 = ceil( randu()*100.0 );
-		d2 = ceil( randu()*100.0 );
-		y = pdf( x, d1, d2 );
+		y = pdf( x[ i % len ], d1[ i % len ], d2[ i % len ]);
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
@@ -58,6 +67,7 @@ bench( pkg, function benchmark( b ) {
 
 bench( pkg+':factory', function benchmark( b ) {
 	var mypdf;
+	var len;
 	var d1;
 	var d2;
 	var x;
@@ -67,11 +77,15 @@ bench( pkg+':factory', function benchmark( b ) {
 	d1 = 1.5;
 	d2 = 1.5;
 	mypdf = pdf.factory( d1, d2 );
+	len = 100;
+	x = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = uniform( EPS, 20.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = ( randu()*20.0 ) + EPS;
-		y = mypdf( x );
+		y = mypdf( x[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/quantile/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/quantile/benchmark/benchmark.js
index 9afb66183639..e910bf9f075d 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/quantile/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/quantile/benchmark/benchmark.js
@@ -21,8 +21,9 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var ceil = require( '@stdlib/math/base/special/ceil' );
-var randu = require( '@stdlib/random/base/randu' );
+var Float64Array = require( '@stdlib/array/float64' );
+var uniform = require( '@stdlib/random/base/uniform' );
+var discreteUniform = require( '@stdlib/random/base/discrete-uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var pkg = require( './../package.json' ).name;
 var quantile = require( './../lib' );
@@ -31,18 +32,26 @@ var quantile = require( './../lib' );
 // MAIN //
 
 bench( pkg, function benchmark( b ) {
+	var len;
 	var d1;
 	var d2;
 	var p;
 	var y;
 	var i;
 
+	len = 100;
+	p = new Float64Array( len );
+	d1 = new Float64Array( len );
+	d2 = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		p[ i ] = uniform( 0.0, 1.0 );
+		d1[ i ] = discreteUniform( 1, 100 );
+		d2[ i ] = discreteUniform( 1, 100 );
+	}
+
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		p = randu();
-		d1 = ceil( randu()*100.0 );
-		d2 = ceil( randu()*100.0 );
-		y = quantile( p, d1, d2 );
+		y = quantile( p[ i % len ], d1[ i % len ], d2[ i % len ]);
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
@@ -57,6 +66,7 @@ bench( pkg, function benchmark( b ) {
 
 bench( pkg+':factory', function benchmark( b ) {
 	var myquantile;
+	var len;
 	var d1;
 	var d2;
 	var p;
@@ -66,11 +76,15 @@ bench( pkg+':factory', function benchmark( b ) {
 	d1 = 1.5;
 	d2 = 1.5;
 	myquantile = quantile.factory( d1, d2 );
+	len = 100;
+	p = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		p[ i ] = uniform( 0.0, 1.0 );
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		p = randu();
-		y = myquantile( p );
+		y = myquantile( p[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/benchmark.js
index 565ebf622e1a..7bafcf30a977 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/benchmark.js
@@ -21,7 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var Float64Array = require( '@stdlib/array/float64' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
@@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + 6.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 6.0 + EPS, 20.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/benchmark.native.js
index 29bd4c53ff8a..3a679c820374 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/benchmark.native.js
@@ -23,7 +23,7 @@
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var Float64Array = require( '@stdlib/array/float64' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
 var tryRequire = require( '@stdlib/utils/try-require' );
@@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + 6.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 6.0 + EPS, 20.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/benchmark.js
index 94cc83ab9897..210a8a7b3747 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/benchmark.js
@@ -21,7 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var Float64Array = require( '@stdlib/array/float64' );
 var EPS = require( '@stdlib/constants/float64/eps' );
@@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu()*10.0 ) + EPS;
-		d2[ i ] = ( randu()*10.0 ) + 4.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 4.0 + EPS, 20.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/benchmark.native.js
index e6d817b6995f..005cda94d367 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/benchmark.native.js
@@ -23,7 +23,7 @@
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var Float64Array = require( '@stdlib/array/float64' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var EPS = require( '@stdlib/constants/float64/eps' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var tryRequire = require( '@stdlib/utils/try-require' );
@@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu()*10.0 ) + EPS;
-		d2[ i ] = ( randu()*10.0 ) + 4.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 4.0 + EPS, 20.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/benchmark.js
index e0456d03116f..362ba56bab5e 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/benchmark.js
@@ -21,7 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var Float64Array = require( '@stdlib/array/float64' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
@@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + 4.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 4.0 + EPS, 20.0 );
 	}
 
 	b.tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/benchmark.native.js
index c47080cc040f..0c895d4e75f6 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/benchmark.native.js
@@ -23,7 +23,7 @@
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var Float64Array = require( '@stdlib/array/float64' );
-var randu = require( '@stdlib/random/base/randu' );
+var uniform = require( '@stdlib/random/base/uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
 var tryRequire = require( '@stdlib/utils/try-require' );
@@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	d1 = new Float64Array( len );
 	d2 = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		d1[ i ] = ( randu() * 10.0 ) + EPS;
-		d2[ i ] = ( randu() * 10.0 ) + 4.0 + EPS;
+		d1[ i ] = uniform( EPS, 10.0 );
+		d2[ i ] = uniform( 4.0 + EPS, 20.0 );
 	}
 
 	b.tic();