diff --git a/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/lib/main.js b/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/lib/main.js
index 0b3a10605b21..298ff5188ba2 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/lib/main.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/lib/main.js
@@ -20,8 +20,8 @@
 
 // MODULES //
 
-var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var sqrt = require( '@stdlib/math/base/special/sqrt' );
+var variance = require( '@stdlib/stats/base/dists/bernoulli/variance' );
 
 
 // MAIN //
@@ -49,14 +49,7 @@ var sqrt = require( '@stdlib/math/base/special/sqrt' );
 * // returns NaN
 */
 function stdev( p ) {
-	if (
-		isnan( p ) ||
-		p < 0.0 ||
-		p > 1.0
-	) {
-		return NaN;
-	}
-	return sqrt( ( 1.0-p ) * p );
+	return sqrt( variance( p ) );
 }
 
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/manifest.json b/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/manifest.json
index 546b58881775..0b79e0701b2b 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/manifest.json
+++ b/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/manifest.json
@@ -39,8 +39,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/napi/unary",
-        "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/sqrt"
+        "@stdlib/math/base/special/sqrt",
+        "@stdlib/stats/base/dists/bernoulli/variance"
       ]
     },
     {
@@ -55,8 +55,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/sqrt"
+        "@stdlib/math/base/special/sqrt",
+        "@stdlib/stats/base/dists/bernoulli/variance"
       ]
     },
     {
@@ -71,8 +71,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/sqrt"
+        "@stdlib/math/base/special/sqrt",
+        "@stdlib/stats/base/dists/bernoulli/variance"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/src/main.c
index 0db3e3c4cf20..c3fd826e4188 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/src/main.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/bernoulli/stdev/src/main.c
@@ -17,8 +17,8 @@
 */
 
 #include "stdlib/stats/base/dists/bernoulli/stdev.h"
-#include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/sqrt.h"
+#include "stdlib/stats/base/dists/bernoulli/variance.h"
 
 /**
 * Returns the standard deviation of a Bernoulli distribution.
@@ -31,12 +31,5 @@
 * // returns ~0.3
 */
 double stdlib_base_dists_bernoulli_stdev( const double p ) {
-	if (
-		stdlib_base_is_nan( p ) ||
-		p < 0.0 ||
-		p > 1.0
-	) {
-		return 0.0/0.0; // NaN
-	}
-	return stdlib_base_sqrt( ( 1.0-p ) * p );
+	return stdlib_base_sqrt( stdlib_base_dists_bernoulli_variance( p ) );
 }