feat: add implementation

---
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: na
  - task: lint_package_json
    status: na
  - task: lint_repl_help
    status: na
  - task: lint_javascript_src
    status: passed
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: na
  - task: lint_javascript_tests
    status: na
  - task: lint_javascript_benchmarks
    status: na
  - 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: na
  - task: lint_typescript_tests
    status: na
  - task: lint_license_headers
    status: passed
---

Author: anandkaranubc

lib/node_modules/@stdlib/stats/base/dists/bradford/pdf/lib/factory.js

@@ -0,0 +1,81 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var constantFunction = require( '@stdlib/utils/constant-function' );
+var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var ln = require( '@stdlib/math/base/special/ln' );
+
+
+// MAIN //
+
+/**
+* Returns a function for evaluating the probability density function (PDF) for a bradford distribution with shape parameter `c`.
+*
+* @param {number} c - shape parameter
+* @returns {Function} PDF
+*
+* @example
+* var pdf = factory( 5.0 );
+* var y = pdf( 0.5 );
+* // returns ~0.797
+*
+* y = pdf( 1.0 );
+* // returns ~0.465
+*/
+function factory( c ) {
+	if (
+		isnan( c ) ||
+		c <= 0.0
+	) {
+		return constantFunction( NaN );
+	}
+	return pdf;
+
+	/**
+	* Evaluates the probability density function (PDF) for a bradford distribution.
+	*
+	* @private
+	* @param {number} x - input value
+	* @returns {number} evaluated PDF
+	*
+	* @example
+	* var y = pdf( 0.5 );
+	* // returns <number>
+	*/
+	function pdf( x ) {
+		var k;
+		if (
+			isnan( x ) ||
+			x < 0.0 ||
+			x > 1.0
+		) {
+			return NaN;
+		}
+		k = ln( 1.0 + c );
+		return c / ( ( 1.0 + ( c * x ) ) * k );
+	}
+}
+
+
+// EXPORTS //
+
+module.exports = factory;

lib/node_modules/@stdlib/stats/base/dists/bradford/pdf/lib/index.js

@@ -0,0 +1,54 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+/**
+* Bradford distribution probability density function (PDF).
+*
+* @module @stdlib/stats/base/dists/bradford/pdf
+*
+* @example
+* var pdf = require( '@stdlib/stats/base/dists/bradford/pdf' );
+*
+* var y = pdf( 0.5, 5.0 );
+* // returns ~0.797
+*
+* var myPDF = pdf.factory( 5.0 );
+* y = myPDF( 0.5 );
+* // returns ~0.797
+*
+* y = myPDF( 1.0 );
+* // returns ~0.465
+*/
+
+// MODULES //
+
+var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' );
+var main = require( './main.js' );
+var factory = require( './factory.js' );
+
+
+// MAIN //
+
+setReadOnly( main, 'factory', factory );
+
+
+// EXPORTS //
+
+module.exports = main;

lib/node_modules/@stdlib/stats/base/dists/bradford/pdf/lib/main.js

@@ -0,0 +1,86 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var ln = require( '@stdlib/math/base/special/ln' );
+
+
+// MAIN //
+
+/**
+* Returns the probability density function (PDF) for a Bradford distribution with shape parameter `c` at a value `x`.
+*
+* @param {number} x - input value
+* @param {number} c - shape parameter
+* @returns {number} evaluated PDF
+*
+* @example
+* var v = pdf( 0.1, 0.1 );
+* // returns ~1.039
+*
+* @example
+* var v = pdf( 0.5, 5.0 );
+* // returns ~0.797
+*
+* @example
+* var v = pdf( 1.0, 10.0 );
+* // returns ~0.379
+*
+* @example
+* var v = pdf( 0.5, 0.0 );
+* // returns NaN
+*
+* @example
+* var v = pdf( 2.0, 0.5 );
+* // returns NaN
+*
+* @example
+* var v = pdf( -1.0, 0.5 );
+* // returns NaN
+*
+* @example
+* var v = pdf( NaN, 1.0 );
+* // returns NaN
+*
+* @example
+* var v = pdf( 1.0, NaN );
+* // returns NaN
+*/
+function pdf( x, c ) {
+	var k;
+	if (
+		isnan( c ) ||
+		isnan( x ) ||
+		c <= 0.0 ||
+		x < 0.0 ||
+		x > 1.0
+	) {
+		return NaN;
+	}
+	k = ln( 1.0 + c );
+	return c / ( ( 1.0 + ( c * x ) ) * k );
+}
+
+
+// EXPORTS //
+
+module.exports = pdf;