UNPKG

@stdlib/stats

Version:

Standard library statistical functions.

88 lines (78 loc) 1.9 kB
/** * @license Apache-2.0 * * Copyright (c) 2018 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 exp = require( '@stdlib/math/base/special/exp' ); var NINF = require( '@stdlib/constants/float64/ninf' ); // MAIN // /** * Evaluates the probability density function (PDF) for a Gumbel distribution with location parameter `mu` and scale parameter `beta` at a value `x`. * * @param {number} x - input value * @param {number} mu - location parameter * @param {PositiveNumber} beta - scale parameter * @returns {number} evaluated PDF * * @example * var y = pdf( 0.0, 0.0, 2.0 ); * // returns ~0.184 * * @example * var y = pdf( 0.0, 0.0, 1.0 ); * // returns ~0.368 * * @example * var y = pdf( 1.0, 3.0, 2.0 ); * // returns ~0.09 * * @example * var y = pdf( NaN, 0.0, 1.0 ); * // returns NaN * * @example * var y = pdf( 0.0, NaN, 1.0 ); * // returns NaN * * @example * var y = pdf( 0.0, 0.0, NaN ); * // returns NaN * * @example * // Negative scale parameter: * var y = pdf( 2.0, 0.0, -1.0 ); * // returns NaN */ function pdf( x, mu, beta ) { var z; if ( isnan( x ) || isnan( mu ) || isnan( beta ) || beta <= 0.0 ) { return NaN; } if ( x === NINF ) { return 0.0; } z = ( x - mu ) / beta; return ( 1.0 / beta ) * exp( -z - exp( -z ) ); } // EXPORTS // module.exports = pdf;