diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/pmf/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/pmf/test/test.factory.js
index 7285a136d716..2367c11b273a 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/pmf/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/pmf/test/test.factory.js
@@ -21,11 +21,10 @@
 // MODULES //
 
 var tape = require( 'tape' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
 var factory = require( './../lib/factory.js' );
 
 
@@ -174,9 +173,7 @@ tape( 'if provided a `r` which is not a positive number, the created function al
 
 tape( 'the created function evaluates the pmf for `x` given large `r` and `p`', function test( t ) {
 	var expected;
-	var delta;
 	var pmf;
-	var tol;
 	var x;
 	var r;
 	var p;
@@ -190,22 +187,14 @@ tape( 'the created function evaluates the pmf for `x` given large `r` and `p`',
 	for ( i = 0; i < x.length; i++ ) {
 		pmf = factory( r[i], p[i] );
 		y = pmf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', r: '+r[i]+', p: '+p[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[i] );
-			tol = 650.0 * EPS * abs( expected[i] );
-			t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. r: '+r[i]+'. p: '+p[i]+'. y: '+y+'. Expected: '+expected[i]+'. tol: '+tol+'. delta: '+delta+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 655 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the created function evaluates the pmf for `x` given a large parameter `r` and small `p`', function test( t ) {
 	var expected;
-	var delta;
 	var pmf;
-	var tol;
 	var i;
 	var r;
 	var p;
@@ -219,22 +208,14 @@ tape( 'the created function evaluates the pmf for `x` given a large parameter `r
 	for ( i = 0; i < x.length; i++ ) {
 		pmf = factory( r[i], p[i] );
 		y = pmf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', r: '+r[i]+', p: '+p[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[i] );
-			tol = 450.0 * EPS * abs( expected[i] );
-			t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. r: '+r[i]+'. p: '+p[i]+'. y: '+y+'. Expected: '+expected[i]+'. tol: '+tol+'. delta: '+delta+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 688 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the created function evaluates the pmf for `x` given small `r` and large parameter `p`', function test( t ) {
 	var expected;
-	var delta;
 	var pmf;
-	var tol;
 	var i;
 	var r;
 	var p;
@@ -248,22 +229,14 @@ tape( 'the created function evaluates the pmf for `x` given small `r` and large
 	for ( i = 0; i < x.length; i++ ) {
 		pmf = factory( r[i], p[i] );
 		y = pmf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', r: '+r[i]+', p: '+p[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[i] );
-			tol = 80.0 * EPS * abs( expected[i] );
-			t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. r: '+r[i]+'. p: '+p[i]+'. y: '+y+'. Expected: '+expected[i]+'. tol: '+tol+'. delta: '+delta+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 143 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the created function evaluates the pmf for `x` given small `r` and `p`', function test( t ) {
 	var expected;
-	var delta;
 	var pmf;
-	var tol;
 	var i;
 	var r;
 	var p;
@@ -277,13 +250,7 @@ tape( 'the created function evaluates the pmf for `x` given small `r` and `p`',
 	for ( i = 0; i < x.length; i++ ) {
 		pmf = factory( r[i], p[i] );
 		y = pmf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', r: '+r[i]+', p: '+p[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[i] );
-			tol = 80.0 * EPS * abs( expected[i] );
-			t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. r: '+r[i]+'. p: '+p[i]+'. y: '+y+'. Expected: '+expected[i]+'. tol: '+tol+'. delta: '+delta+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 110 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/pmf/test/test.pmf.js b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/pmf/test/test.pmf.js
index 96dbdef76e71..d62ffa9bf7df 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/pmf/test/test.pmf.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/pmf/test/test.pmf.js
@@ -21,11 +21,10 @@
 // MODULES //
 
 var tape = require( 'tape' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
 var pmf = require( './../lib' );
 
 
@@ -122,8 +121,6 @@ tape( 'if provided a success probability `p` outside of `[0,1]`, the function re
 
 tape( 'the function evaluates the pmf for `x` given large `r` and `p`', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var i;
 	var r;
 	var p;
@@ -136,21 +133,13 @@ tape( 'the function evaluates the pmf for `x` given large `r` and `p`', function
 	p = highHigh.p;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[i], r[i], p[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', r: '+r[i]+', p: '+p[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[i] );
-			tol = 650.0 * EPS * abs( expected[i] );
-			t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. r: '+r[i]+'. p: '+p[i]+'. y: '+y+'. Expected: '+expected[i]+'. tol: '+tol+'. delta: '+delta+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 655 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the pmf for `x` given large parameter `r` and small `p`', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var i;
 	var r;
 	var p;
@@ -163,21 +152,13 @@ tape( 'the function evaluates the pmf for `x` given large parameter `r` and smal
 	p = highSmall.p;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[i], r[i], p[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', r: '+r[i]+', p: '+p[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[i] );
-			tol = 450.0 * EPS * abs( expected[i] );
-			t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. r: '+r[i]+'. p: '+p[i]+'. y: '+y+'. Expected: '+expected[i]+'. tol: '+tol+'. delta: '+delta+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 688 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the pmf for `x` given small `r` and large `p`', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var i;
 	var r;
 	var p;
@@ -190,21 +171,13 @@ tape( 'the function evaluates the pmf for `x` given small `r` and large `p`', fu
 	p = smallHigh.p;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[i], r[i], p[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', r: '+r[i]+', p: '+p[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[i] );
-			tol = 80.0 * EPS * abs( expected[i] );
-			t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. r: '+r[i]+'. p: '+p[i]+'. y: '+y+'. Expected: '+expected[i]+'. tol: '+tol+'. delta: '+delta+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 143 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the pmf for `x` given small `r` and `p`', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var i;
 	var r;
 	var p;
@@ -217,13 +190,7 @@ tape( 'the function evaluates the pmf for `x` given small `r` and `p`', function
 	p = smallSmall.p;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[i], r[i], p[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', r: '+r[i]+', p: '+p[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[i] );
-			tol = 80.0 * EPS * abs( expected[i] );
-			t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. r: '+r[i]+'. p: '+p[i]+'. y: '+y+'. Expected: '+expected[i]+'. tol: '+tol+'. delta: '+delta+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 110 ), true, 'returns expected value' );
 	}
 	t.end();
 });