@stdlib/stats-base-dists-signrank-pdf 0.2.2 → 0.3.0

package/NOTICE

@@ -1 +1 @@
-Copyright (c) 2016-2024 The Stdlib Authors.
+Copyright (c) 2016-2026 The Stdlib Authors.

package/README.md

@@ -141,7 +141,7 @@ for ( i = 0; i < 10; i++ ) {
     x = randu() * 30.0;
     n = ceil( randu() * 30.0 );
     y = pdf( x, n );
-    console.log( 'x: %d, n: %d, F(x;n): %d', x.toFixed( 4 ), n.toFixed( 4 ), y.toFixed( 4 ) );
+    console.log( 'x: %d, n: %d, f(x;n): %d', x.toFixed( 4 ), n.toFixed( 4 ), y.toFixed( 4 ) );
 }
 ```
 
@@ -149,6 +149,103 @@ for ( i = 0; i < 10; i++ ) {
 
 <!-- /.examples -->
 
+<!-- C interface documentation. -->
+
+* * *
+
+<section class="c">
+
+## C APIs
+
+<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
+
+<section class="intro">
+
+</section>
+
+<!-- /.intro -->
+
+<!-- C usage documentation. -->
+
+<section class="usage">
+
+### Usage
+
+```c
+#include "stdlib/stats/base/dists/signrank/pdf.h"
+```
+
+#### stdlib_base_dists_signrank_pdf( x, n )
+
+Evaluates the probability density function (PDF) of the Wilcoxon signed rank test statistic with `n` observations.
+
+```c
+double out = stdlib_base_dists_signrank_pdf( 7.0, 9 );
+// returns ~0.01
+```
+
+The function accepts the following arguments:
+
+-   **x**: `[in] double` input value.
+-   **n**: `[in] int32_t` number of observations.
+
+```c
+double stdlib_base_dists_signrank_pdf( const double x, const int32_t n );
+```
+
+</section>
+
+<!-- /.usage -->
+
+<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="notes">
+
+</section>
+
+<!-- /.notes -->
+
+<!-- C API usage examples. -->
+
+<section class="examples">
+
+### Examples
+
+```c
+#include "stdlib/stats/base/dists/signrank/pdf.h"
+#include "stdlib/math/base/special/ceil.h"
+#include <stdlib.h>
+#include <stdio.h>
+#include <stdint.h>
+
+static double random_uniform( const double min, const double max ) {
+    double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
+    return min + ( v*(max-min) );
+}
+
+int main( void ) {
+    int32_t n;
+    double x;
+    double y;
+    int i;
+
+    for ( i = 0; i < 25; i++ ) {
+        x = random_uniform( 0, 30.0 );
+        n = (int32_t)stdlib_base_ceil( random_uniform( 1.0, 30.0 ) );
+        y = stdlib_base_dists_signrank_pdf( x, n );
+        printf( "x: %lf, n: %d, f(x;n): %lf\n", x, n, y );
+    }
+}
+```
+
+</section>
+
+<!-- /.examples -->
+
+</section>
+
+<!-- /.c -->
+
 <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
 
 <section class="related">
@@ -183,7 +280,7 @@ See [LICENSE][stdlib-license].
 
 ## Copyright
 
-Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
+Copyright &copy; 2016-2026. The Stdlib [Authors][stdlib-authors].
 
 </section>
 
@@ -196,8 +293,8 @@ Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 [npm-image]: http://img.shields.io/npm/v/@stdlib/stats-base-dists-signrank-pdf.svg
 [npm-url]: https://npmjs.org/package/@stdlib/stats-base-dists-signrank-pdf
 
-[test-image]: https://github.com/stdlib-js/stats-base-dists-signrank-pdf/actions/workflows/test.yml/badge.svg?branch=v0.2.2
-[test-url]: https://github.com/stdlib-js/stats-base-dists-signrank-pdf/actions/workflows/test.yml?query=branch:v0.2.2
+[test-image]: https://github.com/stdlib-js/stats-base-dists-signrank-pdf/actions/workflows/test.yml/badge.svg?branch=v0.3.0
+[test-url]: https://github.com/stdlib-js/stats-base-dists-signrank-pdf/actions/workflows/test.yml?query=branch:v0.3.0
 
 [coverage-image]: https://img.shields.io/codecov/c/github/stdlib-js/stats-base-dists-signrank-pdf/main.svg
 [coverage-url]: https://codecov.io/github/stdlib-js/stats-base-dists-signrank-pdf?branch=main
@@ -209,8 +306,8 @@ Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 
 -->
 
-[chat-image]: https://img.shields.io/gitter/room/stdlib-js/stdlib.svg
-[chat-url]: https://app.gitter.im/#/room/#stdlib-js_stdlib:gitter.im
+[chat-image]: https://img.shields.io/badge/zulip-join_chat-brightgreen.svg
+[chat-url]: https://stdlib.zulipchat.com
 
 [stdlib]: https://github.com/stdlib-js/stdlib
 

package/dist/index.js

@@ -1,9 +1,9 @@
-"use strict";var a=function(r,e){return function(){return e||r((e={exports:{}}).exports,e),e.exports}};var s=a(function(D,v){
+"use strict";var a=function(r,e){return function(){return e||r((e={exports:{}}).exports,e),e.exports}};var s=a(function(G,v){
 var m=require('@stdlib/utils-memoize/dist'),u;function p(r,e){var i;return e===0?r===0?1:0:(i=e*(e+1)/2,r<0||r>i?0:(r>i/2&&(r=i-r),u(r-e,e-1)+u(r,e-1)))}u=m(p);v.exports=u
-});var f=a(function(E,n){
-var N=require('@stdlib/math-base-assert-is-positive-integer/dist'),l=require('@stdlib/math-base-assert-is-finite/dist'),g=require('@stdlib/math-base-assert-is-nan/dist'),d=require('@stdlib/math-base-special-exp/dist'),y=require('@stdlib/math-base-special-ln/dist'),h=require('@stdlib/constants-float64-ln-two/dist'),w=s();function z(r,e){var i;return g(r)||!N(e)||!l(e)?NaN:(i=e*(e+1)/2,r<0||r>i?0:d(y(w(r,e))-e*h))}n.exports=z
-});var o=a(function(G,q){
-var I=require('@stdlib/math-base-assert-is-positive-integer/dist'),L=require('@stdlib/utils-constant-function/dist'),P=require('@stdlib/math-base-assert-is-finite/dist'),F=require('@stdlib/math-base-assert-is-nan/dist'),O=require('@stdlib/math-base-special-exp/dist'),R=require('@stdlib/math-base-special-ln/dist'),b=require('@stdlib/constants-float64-ln-two/dist'),j=s();function k(r){var e;if(!I(r)||!P(r))return L(NaN);return e=r*(r+1)/2,i;function i(t){return F(t)?NaN:t<0||t>e?0:O(R(j(t,r))-r*b)}}q.exports=k
-});var A=require('@stdlib/utils-define-nonenumerable-read-only-property/dist'),c=f(),B=o();A(c,"factory",B);module.exports=c;
+});var f=a(function(H,n){
+var N=require('@stdlib/math-base-assert-is-positive-integer/dist'),g=require('@stdlib/math-base-assert-is-integer/dist'),l=require('@stdlib/math-base-assert-is-finite/dist'),d=require('@stdlib/math-base-assert-is-nan/dist'),y=require('@stdlib/math-base-special-exp/dist'),I=require('@stdlib/math-base-special-ln/dist'),h=require('@stdlib/constants-float64-ln-two/dist'),w=s();function z(r,e){var i;return d(r)||!N(e)||!l(e)?NaN:!g(r)||(i=e*(e+1)/2,r<0||r>i)?0:y(I(w(r,e))-e*h)}n.exports=z
+});var o=a(function(J,q){
+var L=require('@stdlib/math-base-assert-is-positive-integer/dist'),P=require('@stdlib/utils-constant-function/dist'),F=require('@stdlib/math-base-assert-is-integer/dist'),O=require('@stdlib/math-base-assert-is-finite/dist'),R=require('@stdlib/math-base-assert-is-nan/dist'),b=require('@stdlib/math-base-special-exp/dist'),j=require('@stdlib/math-base-special-ln/dist'),k=require('@stdlib/constants-float64-ln-two/dist'),A=s();function B(r){var e;if(!L(r)||!O(r))return P(NaN);return e=r*(r+1)/2,i;function i(t){return R(t)?NaN:!F(t)||t<0||t>e?0:b(j(A(t,r))-r*k)}}q.exports=B
+});var C=require('@stdlib/utils-define-nonenumerable-read-only-property/dist'),c=f(),D=o();C(c,"factory",D);module.exports=c;
 /** @license Apache-2.0 */
 //# sourceMappingURL=index.js.map

package/dist/index.js.map

@@ -1,7 +1,7 @@
 {
   "version": 3,
   "sources": ["../lib/weights.js", "../lib/main.js", "../lib/factory.js", "../lib/index.js"],
-  "sourcesContent": ["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2020 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar memoize = require( '@stdlib/utils-memoize' );\n\n\n// VARIABLES //\n\nvar memoized;\n\n\n// FUNCTIONS //\n\n/**\n* Calculates the weight for the `(x,n)` pair and memoizes the result.\n*\n* @private\n* @param {number} x - input value\n* @param {NonNegativeInteger} n - number of observations\n* @returns {number} weight\n*/\nfunction weights( x, n ) {\n\tvar mlim;\n\n\tif ( n === 0 ) {\n\t\treturn ( x === 0 ) ? 1 : 0;\n\t}\n\tmlim = n * ( n + 1 ) / 2;\n\tif ( x < 0 || x > mlim ) {\n\t\treturn 0;\n\t}\n\tif ( x > mlim / 2 ) {\n\t\tx = mlim - x;\n\t}\n\treturn memoized( x - n, n - 1 ) + memoized( x, n - 1 );\n}\n\n\n// MAIN //\n\nmemoized = memoize( weights );\n\n\n// EXPORTS //\n\nmodule.exports = memoized;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2020 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar isPositiveInteger = require( '@stdlib/math-base-assert-is-positive-integer' );\nvar isfinite = require( '@stdlib/math-base-assert-is-finite' );\nvar isnan = require( '@stdlib/math-base-assert-is-nan' );\nvar exp = require( '@stdlib/math-base-special-exp' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar LN2 = require( '@stdlib/constants-float64-ln-two' );\nvar weights = require( './weights.js' );\n\n\n// MAIN //\n\n/**\n* Evaluates the probability density function (PDF) of the Wilcoxon signed rank test statistic with `n` observations.\n*\n* @param {number} x - input value\n* @param {PositiveInteger} n - number of observations\n* @returns {Probability} evaluated PDF\n*\n* @example\n* var y = pdf( 7.0, 9 );\n* // returns ~0.01\n*\n* @example\n* var y = pdf( 7.0, 6 );\n* // returns ~0.063\n*\n* @example\n* var y = pdf( -1.0, 40 );\n* // returns 0.0\n*\n* @example\n* var y = pdf( NaN, 10 );\n* // returns NaN\n*\n* @example\n* var y = pdf( 0.0, NaN );\n* // returns NaN\n*\n* @example\n* var y = pdf( 2.0, -1 );\n* // returns NaN\n*\n* @example\n* var y = pdf( 2.0, 1.8 );\n* // returns NaN\n*/\nfunction pdf( x, n ) {\n\tvar mlim;\n\tif (\n\t\tisnan( x ) ||\n\t\t!isPositiveInteger( n ) ||\n\t\t!isfinite( n )\n\t) {\n\t\treturn NaN;\n\t}\n\tmlim = ( n * ( n + 1 ) ) / 2;\n\tif ( x < 0.0 || x > mlim ) {\n\t\treturn 0.0;\n\t}\n\treturn exp( ln( weights( x, n ) ) - ( n * LN2 ) );\n}\n\n\n// EXPORTS //\n\nmodule.exports = pdf;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2020 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar isPositiveInteger = require( '@stdlib/math-base-assert-is-positive-integer' );\nvar constantFunction = require( '@stdlib/utils-constant-function' );\nvar isfinite = require( '@stdlib/math-base-assert-is-finite' );\nvar isnan = require( '@stdlib/math-base-assert-is-nan' );\nvar exp = require( '@stdlib/math-base-special-exp' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar LN2 = require( '@stdlib/constants-float64-ln-two' );\nvar weights = require( './weights.js' );\n\n\n// MAIN //\n\n/**\n* Returns a function for evaluating the probability density function (PDF) for the distribution of the Wilcoxon signed rank test statistic with `n` observations.\n*\n* @param {PositiveInteger} n - number of observations\n* @returns {Function} PDF\n*\n* @example\n* var pdf = factory( 8 );\n* var y = pdf( 4.0 );\n* // returns ~0.008\n*\n* y = pdf( 17.0 );\n* // returns ~0.051\n*/\nfunction factory( n ) {\n\tvar mlim;\n\tif ( !isPositiveInteger( n ) || !isfinite( n ) ) {\n\t\treturn constantFunction( NaN );\n\t}\n\tmlim = n * ( n + 1 ) / 2;\n\treturn pdf;\n\n\t/**\n\t* Evaluates the probability density function (PDF) for the distribution of the Wilcoxon signed rank test statistic.\n\t*\n\t* @private\n\t* @param {number} x - input value\n\t* @returns {Probability} evaluated PDF\n\t*\n\t* @example\n\t* var y = pdf( 2 );\n\t* // returns <number>\n\t*/\n\tfunction pdf( x ) {\n\t\tif ( isnan( x ) ) {\n\t\t\treturn NaN;\n\t\t}\n\t\tif ( x < 0.0 || x > mlim ) {\n\t\t\treturn 0.0;\n\t\t}\n\t\treturn exp( ln( weights( x, n ) ) - ( n * LN2 ) );\n\t}\n}\n\n\n// EXPORTS //\n\nmodule.exports = factory;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2020 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n/**\n* Wilcoxon signed rank test statistic probability density function (PDF).\n*\n* @module @stdlib/stats-base-dists-signrank-pdf\n*\n* @example\n* var pdf = require( '@stdlib/stats-base-dists-signrank-pdf' );\n*\n* var y = pdf( 7.0, 9 );\n* // returns ~0.01\n*\n* var mypdf = pdf.factory( 8 );\n* y = mypdf( 4.0 );\n* // returns ~0.008\n*/\n\n// MODULES //\n\nvar setReadOnly = require( '@stdlib/utils-define-nonenumerable-read-only-property' );\nvar main = require( './main.js' );\nvar factory = require( './factory.js' );\n\n\n// MAIN //\n\nsetReadOnly( main, 'factory', factory );\n\n\n// EXPORTS //\n\nmodule.exports = main;\n"],
-  "mappings": "uGAAA,IAAAA,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAU,QAAS,uBAAwB,EAK3CC,EAaJ,SAASC,EAASC,EAAGC,EAAI,CACxB,IAAIC,EAEJ,OAAKD,IAAM,EACDD,IAAM,EAAM,EAAI,GAE1BE,EAAOD,GAAMA,EAAI,GAAM,EAClBD,EAAI,GAAKA,EAAIE,EACV,GAEHF,EAAIE,EAAO,IACfF,EAAIE,EAAOF,GAELF,EAAUE,EAAIC,EAAGA,EAAI,CAAE,EAAIH,EAAUE,EAAGC,EAAI,CAAE,GACtD,CAKAH,EAAWD,EAASE,CAAQ,EAK5BH,EAAO,QAAUE,IChEjB,IAAAK,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAoB,QAAS,8CAA+C,EAC5EC,EAAW,QAAS,oCAAqC,EACzDC,EAAQ,QAAS,iCAAkC,EACnDC,EAAM,QAAS,+BAAgC,EAC/CC,EAAK,QAAS,8BAA+B,EAC7CC,EAAM,QAAS,kCAAmC,EAClDC,EAAU,IAwCd,SAASC,EAAKC,EAAGC,EAAI,CACpB,IAAIC,EACJ,OACCR,EAAOM,CAAE,GACT,CAACR,EAAmBS,CAAE,GACtB,CAACR,EAAUQ,CAAE,EAEN,KAERC,EAASD,GAAMA,EAAI,GAAQ,EACtBD,EAAI,GAAOA,EAAIE,EACZ,EAEDP,EAAKC,EAAIE,EAASE,EAAGC,CAAE,CAAE,EAAMA,EAAIJ,CAAM,EACjD,CAKAN,EAAO,QAAUQ,ICvFjB,IAAAI,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAoB,QAAS,8CAA+C,EAC5EC,EAAmB,QAAS,iCAAkC,EAC9DC,EAAW,QAAS,oCAAqC,EACzDC,EAAQ,QAAS,iCAAkC,EACnDC,EAAM,QAAS,+BAAgC,EAC/CC,EAAK,QAAS,8BAA+B,EAC7CC,EAAM,QAAS,kCAAmC,EAClDC,EAAU,IAmBd,SAASC,EAASC,EAAI,CACrB,IAAIC,EACJ,GAAK,CAACV,EAAmBS,CAAE,GAAK,CAACP,EAAUO,CAAE,EAC5C,OAAOR,EAAkB,GAAI,EAE9B,OAAAS,EAAOD,GAAMA,EAAI,GAAM,EAChBE,EAaP,SAASA,EAAKC,EAAI,CACjB,OAAKT,EAAOS,CAAE,EACN,IAEHA,EAAI,GAAOA,EAAIF,EACZ,EAEDN,EAAKC,EAAIE,EAASK,EAAGH,CAAE,CAAE,EAAMA,EAAIH,CAAM,CACjD,CACD,CAKAP,EAAO,QAAUS,IC3CjB,IAAIK,EAAc,QAAS,uDAAwD,EAC/EC,EAAO,IACPC,EAAU,IAKdF,EAAaC,EAAM,UAAWC,CAAQ,EAKtC,OAAO,QAAUD",
-  "names": ["require_weights", "__commonJSMin", "exports", "module", "memoize", "memoized", "weights", "x", "n", "mlim", "require_main", "__commonJSMin", "exports", "module", "isPositiveInteger", "isfinite", "isnan", "exp", "ln", "LN2", "weights", "pdf", "x", "n", "mlim", "require_factory", "__commonJSMin", "exports", "module", "isPositiveInteger", "constantFunction", "isfinite", "isnan", "exp", "ln", "LN2", "weights", "factory", "n", "mlim", "pdf", "x", "setReadOnly", "main", "factory"]
+  "sourcesContent": ["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2020 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar memoize = require( '@stdlib/utils-memoize' );\n\n\n// VARIABLES //\n\nvar memoized;\n\n\n// FUNCTIONS //\n\n/**\n* Calculates the weight for the `(x,n)` pair and memoizes the result.\n*\n* @private\n* @param {number} x - input value\n* @param {NonNegativeInteger} n - number of observations\n* @returns {number} weight\n*/\nfunction weights( x, n ) {\n\tvar mlim;\n\n\tif ( n === 0 ) {\n\t\treturn ( x === 0 ) ? 1 : 0;\n\t}\n\tmlim = n * ( n + 1 ) / 2;\n\tif ( x < 0 || x > mlim ) {\n\t\treturn 0;\n\t}\n\tif ( x > mlim / 2 ) {\n\t\tx = mlim - x;\n\t}\n\treturn memoized( x - n, n - 1 ) + memoized( x, n - 1 );\n}\n\n\n// MAIN //\n\nmemoized = memoize( weights );\n\n\n// EXPORTS //\n\nmodule.exports = memoized;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2020 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar isPositiveInteger = require( '@stdlib/math-base-assert-is-positive-integer' );\nvar isInteger = require( '@stdlib/math-base-assert-is-integer' );\nvar isfinite = require( '@stdlib/math-base-assert-is-finite' );\nvar isnan = require( '@stdlib/math-base-assert-is-nan' );\nvar exp = require( '@stdlib/math-base-special-exp' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar LN2 = require( '@stdlib/constants-float64-ln-two' );\nvar weights = require( './weights.js' );\n\n\n// MAIN //\n\n/**\n* Evaluates the probability density function (PDF) of the Wilcoxon signed rank test statistic with `n` observations.\n*\n* @param {number} x - input value\n* @param {PositiveInteger} n - number of observations\n* @returns {Probability} evaluated PDF\n*\n* @example\n* var y = pdf( 7.0, 9 );\n* // returns ~0.01\n*\n* @example\n* var y = pdf( 7.0, 6 );\n* // returns ~0.063\n*\n* @example\n* var y = pdf( -1.0, 40 );\n* // returns 0.0\n*\n* @example\n* var y = pdf( NaN, 10 );\n* // returns NaN\n*\n* @example\n* var y = pdf( 0.0, NaN );\n* // returns NaN\n*\n* @example\n* var y = pdf( 2.0, -1 );\n* // returns NaN\n*\n* @example\n* var y = pdf( 2.0, 1.8 );\n* // returns NaN\n*/\nfunction pdf( x, n ) {\n\tvar mlim;\n\tif (\n\t\tisnan( x ) ||\n\t\t!isPositiveInteger( n ) ||\n\t\t!isfinite( n )\n\t) {\n\t\treturn NaN;\n\t}\n\tif ( !isInteger( x ) ) {\n\t\treturn 0.0;\n\t}\n\tmlim = ( n * ( n + 1 ) ) / 2;\n\tif ( x < 0.0 || x > mlim ) {\n\t\treturn 0.0;\n\t}\n\treturn exp( ln( weights( x, n ) ) - ( n * LN2 ) );\n}\n\n\n// EXPORTS //\n\nmodule.exports = pdf;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2020 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar isPositiveInteger = require( '@stdlib/math-base-assert-is-positive-integer' );\nvar constantFunction = require( '@stdlib/utils-constant-function' );\nvar isInteger = require( '@stdlib/math-base-assert-is-integer' );\nvar isfinite = require( '@stdlib/math-base-assert-is-finite' );\nvar isnan = require( '@stdlib/math-base-assert-is-nan' );\nvar exp = require( '@stdlib/math-base-special-exp' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar LN2 = require( '@stdlib/constants-float64-ln-two' );\nvar weights = require( './weights.js' );\n\n\n// MAIN //\n\n/**\n* Returns a function for evaluating the probability density function (PDF) for the distribution of the Wilcoxon signed rank test statistic with `n` observations.\n*\n* @param {PositiveInteger} n - number of observations\n* @returns {Function} PDF\n*\n* @example\n* var pdf = factory( 8 );\n* var y = pdf( 4.0 );\n* // returns ~0.008\n*\n* y = pdf( 17.0 );\n* // returns ~0.051\n*/\nfunction factory( n ) {\n\tvar mlim;\n\tif ( !isPositiveInteger( n ) || !isfinite( n ) ) {\n\t\treturn constantFunction( NaN );\n\t}\n\tmlim = n * ( n + 1 ) / 2;\n\treturn pdf;\n\n\t/**\n\t* Evaluates the probability density function (PDF) for the distribution of the Wilcoxon signed rank test statistic.\n\t*\n\t* @private\n\t* @param {number} x - input value\n\t* @returns {Probability} evaluated PDF\n\t*\n\t* @example\n\t* var y = pdf( 2 );\n\t* // returns <number>\n\t*/\n\tfunction pdf( x ) {\n\t\tif ( isnan( x ) ) {\n\t\t\treturn NaN;\n\t\t}\n\t\tif ( !isInteger( x ) ) {\n\t\t\treturn 0.0;\n\t\t}\n\t\tif ( x < 0.0 || x > mlim ) {\n\t\t\treturn 0.0;\n\t\t}\n\t\treturn exp( ln( weights( x, n ) ) - ( n * LN2 ) );\n\t}\n}\n\n\n// EXPORTS //\n\nmodule.exports = factory;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2020 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n/**\n* Wilcoxon signed rank test statistic probability density function (PDF).\n*\n* @module @stdlib/stats-base-dists-signrank-pdf\n*\n* @example\n* var pdf = require( '@stdlib/stats-base-dists-signrank-pdf' );\n*\n* var y = pdf( 7.0, 9 );\n* // returns ~0.01\n*\n* var mypdf = pdf.factory( 8 );\n* y = mypdf( 4.0 );\n* // returns ~0.008\n*/\n\n// MODULES //\n\nvar setReadOnly = require( '@stdlib/utils-define-nonenumerable-read-only-property' );\nvar main = require( './main.js' );\nvar factory = require( './factory.js' );\n\n\n// MAIN //\n\nsetReadOnly( main, 'factory', factory );\n\n\n// EXPORTS //\n\nmodule.exports = main;\n"],
+  "mappings": "uGAAA,IAAAA,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAU,QAAS,uBAAwB,EAK3CC,EAaJ,SAASC,EAASC,EAAGC,EAAI,CACxB,IAAIC,EAEJ,OAAKD,IAAM,EACDD,IAAM,EAAM,EAAI,GAE1BE,EAAOD,GAAMA,EAAI,GAAM,EAClBD,EAAI,GAAKA,EAAIE,EACV,GAEHF,EAAIE,EAAO,IACfF,EAAIE,EAAOF,GAELF,EAAUE,EAAIC,EAAGA,EAAI,CAAE,EAAIH,EAAUE,EAAGC,EAAI,CAAE,GACtD,CAKAH,EAAWD,EAASE,CAAQ,EAK5BH,EAAO,QAAUE,IChEjB,IAAAK,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAoB,QAAS,8CAA+C,EAC5EC,EAAY,QAAS,qCAAsC,EAC3DC,EAAW,QAAS,oCAAqC,EACzDC,EAAQ,QAAS,iCAAkC,EACnDC,EAAM,QAAS,+BAAgC,EAC/CC,EAAK,QAAS,8BAA+B,EAC7CC,EAAM,QAAS,kCAAmC,EAClDC,EAAU,IAwCd,SAASC,EAAKC,EAAGC,EAAI,CACpB,IAAIC,EACJ,OACCR,EAAOM,CAAE,GACT,CAACT,EAAmBU,CAAE,GACtB,CAACR,EAAUQ,CAAE,EAEN,IAEH,CAACT,EAAWQ,CAAE,IAGnBE,EAASD,GAAMA,EAAI,GAAQ,EACtBD,EAAI,GAAOA,EAAIE,GACZ,EAEDP,EAAKC,EAAIE,EAASE,EAAGC,CAAE,CAAE,EAAMA,EAAIJ,CAAM,CACjD,CAKAP,EAAO,QAAUS,IC3FjB,IAAAI,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAoB,QAAS,8CAA+C,EAC5EC,EAAmB,QAAS,iCAAkC,EAC9DC,EAAY,QAAS,qCAAsC,EAC3DC,EAAW,QAAS,oCAAqC,EACzDC,EAAQ,QAAS,iCAAkC,EACnDC,EAAM,QAAS,+BAAgC,EAC/CC,EAAK,QAAS,8BAA+B,EAC7CC,EAAM,QAAS,kCAAmC,EAClDC,EAAU,IAmBd,SAASC,EAASC,EAAI,CACrB,IAAIC,EACJ,GAAK,CAACX,EAAmBU,CAAE,GAAK,CAACP,EAAUO,CAAE,EAC5C,OAAOT,EAAkB,GAAI,EAE9B,OAAAU,EAAOD,GAAMA,EAAI,GAAM,EAChBE,EAaP,SAASA,EAAKC,EAAI,CACjB,OAAKT,EAAOS,CAAE,EACN,IAEH,CAACX,EAAWW,CAAE,GAGdA,EAAI,GAAOA,EAAIF,EACZ,EAEDN,EAAKC,EAAIE,EAASK,EAAGH,CAAE,CAAE,EAAMA,EAAIH,CAAM,CACjD,CACD,CAKAR,EAAO,QAAUU,IC/CjB,IAAIK,EAAc,QAAS,uDAAwD,EAC/EC,EAAO,IACPC,EAAU,IAKdF,EAAaC,EAAM,UAAWC,CAAQ,EAKtC,OAAO,QAAUD",
+  "names": ["require_weights", "__commonJSMin", "exports", "module", "memoize", "memoized", "weights", "x", "n", "mlim", "require_main", "__commonJSMin", "exports", "module", "isPositiveInteger", "isInteger", "isfinite", "isnan", "exp", "ln", "LN2", "weights", "pdf", "x", "n", "mlim", "require_factory", "__commonJSMin", "exports", "module", "isPositiveInteger", "constantFunction", "isInteger", "isfinite", "isnan", "exp", "ln", "LN2", "weights", "factory", "n", "mlim", "pdf", "x", "setReadOnly", "main", "factory"]
 }

package/include/stdlib/stats/base/dists/signrank/pdf.h

@@ -0,0 +1,40 @@
+/**
+* @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.
+*/
+
+#ifndef STDLIB_STATS_BASE_DISTS_SIGNRANK_PDF_H
+#define STDLIB_STATS_BASE_DISTS_SIGNRANK_PDF_H
+
+#include <stdint.h>
+
+/*
+* If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
+*/
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+/**
+* Evaluates the probability density function (PDF) of the Wilcoxon signed rank test statistic with `n` observations.
+*/
+double stdlib_base_dists_signrank_pdf( const double x, const int32_t n );
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif // !STDLIB_STATS_BASE_DISTS_SIGNRANK_PDF_H

package/lib/factory.js

@@ -22,6 +22,7 @@
 
 var isPositiveInteger = require( '@stdlib/math-base-assert-is-positive-integer' );
 var constantFunction = require( '@stdlib/utils-constant-function' );
+var isInteger = require( '@stdlib/math-base-assert-is-integer' );
 var isfinite = require( '@stdlib/math-base-assert-is-finite' );
 var isnan = require( '@stdlib/math-base-assert-is-nan' );
 var exp = require( '@stdlib/math-base-special-exp' );
@@ -69,6 +70,9 @@ function factory( n ) {
 		if ( isnan( x ) ) {
 			return NaN;
 		}
+		if ( !isInteger( x ) ) {
+			return 0.0;
+		}
 		if ( x < 0.0 || x > mlim ) {
 			return 0.0;
 		}

package/lib/main.js

@@ -21,6 +21,7 @@
 // MODULES //
 
 var isPositiveInteger = require( '@stdlib/math-base-assert-is-positive-integer' );
+var isInteger = require( '@stdlib/math-base-assert-is-integer' );
 var isfinite = require( '@stdlib/math-base-assert-is-finite' );
 var isnan = require( '@stdlib/math-base-assert-is-nan' );
 var exp = require( '@stdlib/math-base-special-exp' );
@@ -75,6 +76,9 @@ function pdf( x, n ) {
 	) {
 		return NaN;
 	}
+	if ( !isInteger( x ) ) {
+		return 0.0;
+	}
 	mlim = ( n * ( n + 1 ) ) / 2;
 	if ( x < 0.0 || x > mlim ) {
 		return 0.0;

package/lib/native.js

@@ -0,0 +1,66 @@
+/**
+* @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 addon = require( './../src/addon.node' );
+
+
+// MAIN //
+
+/**
+* Evaluates the probability density function (PDF) of the Wilcoxon signed rank test statistic with `n` observations.
+*
+* @param {number} x - input value
+* @param {PositiveInteger} n - number of observations
+* @returns {Probability} evaluated PDF
+*
+* @example
+* var y = pdf( 7.0, 9 );
+* // returns ~0.01
+*
+* @example
+* var y = pdf( 7.0, 6 );
+* // returns ~0.063
+*
+* @example
+* var y = pdf( -1.0, 40 );
+* // returns 0.0
+*
+* @example
+* var y = pdf( NaN, 10 );
+* // returns NaN
+*
+* @example
+* var y = pdf( 0.0, NaN );
+* // returns NaN
+*
+* @example
+* var y = pdf( 2.0, -1 );
+* // returns NaN
+*/
+function pdf( x, n ) {
+	return addon( x, n );
+}
+
+
+// EXPORTS //
+
+module.exports = pdf;

package/manifest.json

@@ -0,0 +1,78 @@
+{
+  "options": {
+    "task": "build",
+    "wasm": false
+  },
+  "fields": [
+    {
+      "field": "src",
+      "resolve": true,
+      "relative": true
+    },
+    {
+      "field": "include",
+      "resolve": true,
+      "relative": true
+    },
+    {
+      "field": "libraries",
+      "resolve": false,
+      "relative": false
+    },
+    {
+      "field": "libpath",
+      "resolve": true,
+      "relative": false
+    }
+  ],
+  "confs": [
+    {
+      "task": "build",
+      "wasm": false,
+      "src": ["./src/main.c"],
+      "include": ["./include"],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-assert-is-nan",
+        "@stdlib/math-base-special-exp",
+        "@stdlib/constants-float64-ln-two",
+        "@stdlib/math-base-napi-binary",
+        "@stdlib/math-base-special-ln",
+        "@stdlib/math-base-assert-is-integer"
+      ]
+    },
+    {
+      "task": "benchmark",
+      "wasm": false,
+      "src": ["./src/main.c"],
+      "include": ["./include"],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-special-ceil",
+        "@stdlib/math-base-assert-is-nan",
+        "@stdlib/math-base-special-exp",
+        "@stdlib/constants-float64-ln-two",
+        "@stdlib/math-base-special-ln",
+        "@stdlib/math-base-assert-is-integer"
+      ]
+    },
+    {
+      "task": "examples",
+      "wasm": false,
+      "src": ["./src/main.c"],
+      "include": ["./include"],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-special-ceil",
+        "@stdlib/math-base-assert-is-nan",
+        "@stdlib/math-base-special-exp",
+        "@stdlib/constants-float64-ln-two",
+        "@stdlib/math-base-special-ln",
+        "@stdlib/math-base-assert-is-integer"
+      ]
+    }
+  ]
+}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@stdlib/stats-base-dists-signrank-pdf",
-  "version": "0.2.2",
+  "version": "0.3.0",
   "description": "Wilcoxon signed rank test statistic probability density function (PDF).",
   "license": "Apache-2.0",
   "author": {
@@ -14,9 +14,12 @@
     }
   ],
   "main": "./lib",
+  "gypfile": false,
   "directories": {
     "doc": "./docs",
+    "include": "./include",
     "lib": "./lib",
+    "src": "./src",
     "dist": "./dist"
   },
   "types": "./docs/types",
@@ -32,12 +35,15 @@
   "dependencies": {
     "@stdlib/constants-float64-ln-two": "^0.2.2",
     "@stdlib/math-base-assert-is-finite": "^0.2.2",
+    "@stdlib/math-base-assert-is-integer": "^0.2.6",
     "@stdlib/math-base-assert-is-nan": "^0.2.2",
-    "@stdlib/math-base-assert-is-positive-integer": "^0.3.0",
-    "@stdlib/math-base-special-exp": "^0.2.3",
+    "@stdlib/math-base-assert-is-positive-integer": "^0.3.1",
+    "@stdlib/math-base-napi-binary": "^0.3.1",
+    "@stdlib/math-base-special-exp": "^0.2.4",
     "@stdlib/math-base-special-ln": "^0.2.4",
     "@stdlib/utils-constant-function": "^0.2.2",
     "@stdlib/utils-define-nonenumerable-read-only-property": "^0.2.2",
+    "@stdlib/utils-library-manifest": "^0.2.3",
     "@stdlib/utils-memoize": "^0.2.2"
   },
   "devDependencies": {},

package/src/addon.c

@@ -0,0 +1,22 @@
+/**
+* @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.
+*/
+
+#include "stdlib/stats/base/dists/signrank/pdf.h"
+#include "stdlib/math/base/napi/binary.h"
+
+STDLIB_MATH_BASE_NAPI_MODULE_DI_D( stdlib_base_dists_signrank_pdf )

package/src/main.c

@@ -0,0 +1,74 @@
+/**
+* @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.
+*/
+
+#include "stdlib/stats/base/dists/signrank/pdf.h"
+#include "stdlib/math/base/assert/is_integer.h"
+#include "stdlib/math/base/assert/is_nan.h"
+#include "stdlib/math/base/special/exp.h"
+#include "stdlib/constants/float64/ln_two.h"
+#include "stdlib/math/base/special/ln.h"
+#include <stdint.h>
+
+/**
+* Calculates the weight for the (x,n) pair and memoizes the result.
+*
+* @param x   input value
+* @param n   number of observations (must be a non-negative)
+* @return    weight value
+*/
+static double weights( double x, const int32_t n ) {
+	double mlim;
+
+	if ( n == 0 ) {
+		return ( x == 0 ) ? 1.0 : 0.0;
+	}
+	mlim = ( (double)n * ( (double)n + 1.0 ) ) / 2.0;
+	if ( x < 0 || x > mlim ) {
+		return 0.0;
+	}
+	if ( x > mlim / 2.0 ) {
+		x = mlim - x;
+	}
+	return weights( x-(double)n, n-1 ) + weights( x, n-1 );
+}
+
+/**
+* Evaluates the probability density function (PDF) of the Wilcoxon signed rank test statistic with `n` observations.
+*
+* @param x   input value
+* @param n   number of observations (must be a non-negative)
+* @return    evaluated PDF
+*
+* @example
+* double y = stdlib_base_dists_signrank_pdf( 7.0, 9 );
+* // returns ~0.01
+*/
+double stdlib_base_dists_signrank_pdf( const double x, const int32_t n ) {
+	double mlim;
+	if ( stdlib_base_is_nan( x ) || n <= 0 ) {
+		return 0.0/0.0;
+	}
+	if ( !stdlib_base_is_integer( x ) ) {
+		return 0.0;
+	}
+	mlim = ( (double)n * ( (double)n + 1.0 ) ) / 2.0;
+	if ( x < 0.0 || x > mlim ) {
+		return 0.0;
+	}
+	return stdlib_base_exp( stdlib_base_ln( weights( x, n ) ) - ( (double)n * STDLIB_CONSTANT_FLOAT64_LN2 ) );
+}