distributions-lognormal-pdf
Version:
Lognormal distribution probability density function (PDF)
399 lines (298 loc) • 9.88 kB
Markdown
Probability Density Function
===
[![NPM version][npm-image]][npm-url] [![Build Status][travis-image]][travis-url] [![Coverage Status][codecov-image]][codecov-url] [![Dependencies][dependencies-image]][dependencies-url]
> [Lognormal](https://en.wikipedia.org/wiki/Lognormal_distribution) distribution probability density function (PDF).
The [probability density function](https://en.wikipedia.org/wiki/Probability_density_function) (PDF) for a [lognormal](https://en.wikipedia.org/wiki/Lognormal_distribution) random variable is
<div class="equation" align="center" data-raw-text="f(x;\mu,\sigma) = \frac{1}{x\sqrt{2\pi\sigma^2}} e^{-\frac{\left(\ln x-\mu\right)^2}{2\sigma^2}}" data-equation="eq:pdf_function">
<img src="https://cdn.rawgit.com/distributions-io/lognormal-pdf/c1d82cb66e4000ee374d7d4aa9f9c41e36d58d48/docs/img/eqn.svg" alt="Probability density function (PDF) for a lognormal distribution.">
<br>
</div>
where `mu` is the location parameter and `sigma > 0` is the scale parameter. According to the definition, the *natural logarithm* of a random variable from a
[](https://en.wikipedia.org/wiki/Lognormal_distribution) follows a [normal distribution](https://en.wikipedia.org/wiki/Normal_distribution).
``` bash
$ npm install distributions-lognormal-pdf
```
For use in the browser, use [browserify](https://github.com/substack/node-browserify).
``` javascript
var pdf = require( 'distributions-lognormal-pdf' );
```
Evaluates the [probability density function](https://en.wikipedia.org/wiki/Probability_density_function) (PDF) for the [lognormal](https://en.wikipedia.org/wiki/Lognormal_distribution) distribution. `x` may be either a [`number`](https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Number), an [`array`](https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Array), a [`typed array`](https://developer.mozilla.org/en-US/docs/Web/JavaScript/Typed_arrays), or a [`matrix`](https://github.com/dstructs/matrix).
``` javascript
var matrix = require( 'dstructs-matrix' ),
mat,
out,
x,
i;
out = pdf( 1 );
// returns 0.399
out = pdf( -1 );
// returns 0
x = [ 0, 0.5, 1, 1.5, 2, 2.5 ];
out = pdf( x );
// returns [ 0, ~0.627, ~0.399, ~0.245, ~0.157, ~0.105 ]
x = new Int8Array( x );
out = pdf( x );
// returns Float64Array( [0,0,~0.399,~0.399,~0.157,~0.157] )
x = new Float32Array( 6 );
for ( i = 0; i < 6; i++ ) {
x[ i ] = i * 0.5;
}
mat = matrix( x, [3,2], 'float32' );
/*
[ 0 0.5
1 1.5
2 2.5 ]
*/
out = pdf( mat );
/*
[ 0 ~0.627
~0.399 ~0.245
~0.157 ~0.105 ]
*/
```
The function accepts the following `options`:
* __mu__: location parameter. Default: `0`.
* __sigma__: scale parameter. Default: `1`.
* __accessor__: accessor `function` for accessing `array` values.
* __dtype__: output [`typed array`](https://developer.mozilla.org/en-US/docs/Web/JavaScript/Typed_arrays) or [`matrix`](https://github.com/dstructs/matrix) data type. Default: `float64`.
* __copy__: `boolean` indicating if the `function` should return a new data structure. Default: `true`.
* __path__: [deepget](https://github.com/kgryte/utils-deep-get)/[deepset](https://github.com/kgryte/utils-deep-set) key path.
* __sep__: [deepget](https://github.com/kgryte/utils-deep-get)/[deepset](https://github.com/kgryte/utils-deep-set) key path separator. Default: `'.'`.
A [Lognormal](https://en.wikipedia.org/wiki/Lognormal_distribution) distribution is a function of two parameters: `mu`(location parameter) and `sigma > 0`(scale parameter). By default, `mu` is equal to `0` and `sigma` is equal to `1`. To adjust either parameter, set the corresponding option.
``` javascript
var x = [ 0, 0.5, 1, 1.5, 2, 2.5 ];
var out = pdf( x, {
'mu': 8,
'sigma': 2,
});
// returns [ 0, 0, 0, 0, 0, 0 ]
```
For non-numeric `arrays`, provide an accessor `function` for accessing `array` values.
``` javascript
var data = [
[],
[],
[],
[],
[],
[]
];
function getValue( d, i ) {
return d[ 1 ];
}
var out = pdf( data, {
'accessor': getValue
});
// returns [ 0, ~0.627, ~0.399, ~0.245, ~0.157, ~0.105 ]
```
To [deepset](https://github.com/kgryte/utils-deep-set) an object `array`, provide a key path and, optionally, a key path separator.
``` javascript
var data = [
{'x':[0,0]},
{'x':[1,0.5]},
{'x':[2,1]},
{'x':[3,1.5]},
{'x':[4,2]},
{'x':[5,2.5]}
];
var out = pdf( data, {
'path': 'x/1',
'sep': '/'
});
/*
[
{'x':[0,0]},
{'x':[1,~0.627]},
{'x':[2,~0.399]},
{'x':[3,~0.245]},
{'x':[4,~0.157]},
{'x':[5,~0.105]}
]
*/
var bool = ( data === out );
// returns true
```
By default, when provided a [`typed array`](https://developer.mozilla.org/en-US/docs/Web/JavaScript/Typed_arrays) or [`matrix`](https://github.com/dstructs/matrix), the output data structure is `float64` in order to preserve precision. To specify a different data type, set the `dtype` option (see [`matrix`](https://github.com/dstructs/matrix) for a list of acceptable data types).
``` javascript
var x, out;
x = new Int8Array( [0,1,2,3,4] );
out = pdf( x, {
'dtype': 'int32'
});
// returns Int32Array( [0,0,0,0,0] )
// Works for plain arrays, as well...
out = pdf( [0,0.5,1,1.5,2], {
'dtype': 'uint8'
});
// returns Uint8Array( [0,0,0,0,0] )
```
By default, the function returns a new data structure. To mutate the input data structure (e.g., when input values can be discarded or when optimizing memory usage), set the `copy` option to `false`.
``` javascript
var bool,
mat,
out,
x,
i;
x = [ 0, 0.5, 1, 1.5, 2 ];
out = pdf( x, {
'copy': false
});
// returns [ 0, ~0.627, ~0.399, ~0.245, ~0.157 ]
bool = ( x === out );
// returns true
x = new Float32Array( 6 );
for ( i = 0; i < 6; i++ ) {
x[ i ] = i * 0.5;
}
mat = matrix( x, [3,2], 'float32' );
/*
[ 0 0.5
1 1.5
2 2.5 ]
*/
out = pdf( mat, {
'copy': false
});
/*
[ 0 ~0.627
~0.399 ~0.245
~0.157 ~0.105 ]
*/
bool = ( mat === out );
// returns true
```
* If an element is __not__ a numeric value, the evaluated [PDF](https://en.wikipedia.org/wiki/Lognormal_distribution) is `NaN`.
``` javascript
var data, out;
out = pdf( null );
// returns NaN
out = pdf( true );
// returns NaN
out = pdf( {'a':'b'} );
// returns NaN
out = pdf( [ true, null, [] ] );
// returns [ NaN, NaN, NaN ]
function getValue( d, i ) {
return d.x;
}
data = [
{'x':true},
{'x':[]},
{'x':{}},
{'x':null}
];
out = pdf( data, {
'accessor': getValue
});
// returns [ NaN, NaN, NaN, NaN ]
out = pdf( data, {
'path': 'x'
});
/*
[
{'x':NaN},
{'x':NaN},
{'x':NaN,
{'x':NaN}
]
*/
```
* Be careful when providing a data structure which contains non-numeric elements and specifying an `integer` output data type, as `NaN` values are cast to `0`.
``` javascript
var out = pdf( [ true, null, [] ], {
'dtype': 'int8'
});
// returns Int8Array( [0,0,0] );
```
``` javascript
var pdf = require( 'distributions-lognormal-pdf' ),
matrix = require( 'dstructs-matrix' );
var data,
mat,
out,
tmp,
i;
// Plain arrays...
data = new Array( 10 );
for ( i = 0; i < data.length; i++ ) {
data[ i ] = i * 0.5;
}
out = pdf( data );
// Object arrays (accessors)...
function getValue( d ) {
return d.x;
}
for ( i = 0; i < data.length; i++ ) {
data[ i ] = {
'x': data[ i ]
};
}
out = pdf( data, {
'accessor': getValue
});
// Deep set arrays...
for ( i = 0; i < data.length; i++ ) {
data[ i ] = {
'x': [ i, data[ i ].x ]
};
}
out = pdf( data, {
'path': 'x/1',
'sep': '/'
});
// Typed arrays...
data = new Float32Array( 10 );
for ( i = 0; i < data.length; i++ ) {
data[ i ] = i * 0.5;
}
out = pdf( data );
// Matrices...
mat = matrix( data, [5,2], 'float32' );
out = pdf( mat );
// Matrices (custom output data type)...
out = pdf( mat, {
'dtype': 'uint8'
});
```
To run the example code from the top-level application directory,
``` bash
$ node ./examples/index.js
```
Unit tests use the [Mocha](http://mochajs.org/) test framework with [Chai](http://chaijs.com) assertions. To run the tests, execute the following command in the top-level application directory:
``` bash
$ make test
```
All new feature development should have corresponding unit tests to validate correct functionality.
This repository uses [Istanbul](https://github.com/gotwarlost/istanbul) as its code coverage tool. To generate a test coverage report, execute the following command in the top-level application directory:
``` bash
$ make test-cov
```
Istanbul creates a `./reports/coverage` directory. To access an HTML version of the report,
``` bash
$ make view-cov
```
---
[](http://opensource.org/licenses/MIT).
Copyright © 2015. The [Compute.io](https://github.com/compute-io) Authors.
[]: http://img.shields.io/npm/v/distributions-lognormal-pdf.svg
[]: https://npmjs.org/package/distributions-lognormal-pdf
[]: http://img.shields.io/travis/distributions-io/lognormal-pdf/master.svg
[]: https://travis-ci.org/distributions-io/lognormal-pdf
[]: https://img.shields.io/codecov/c/github/distributions-io/lognormal-pdf/master.svg
[]: https://codecov.io/github/distributions-io/lognormal-pdf?branch=master
[]: http://img.shields.io/david/distributions-io/lognormal-pdf.svg
[]: https://david-dm.org/distributions-io/lognormal-pdf
[]: http://img.shields.io/david/dev/distributions-io/lognormal-pdf.svg
[]: https://david-dm.org/dev/distributions-io/lognormal-pdf
[]: http://img.shields.io/github/issues/distributions-io/lognormal-pdf.svg
[]: https://github.com/distributions-io/lognormal-pdf/issues