diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/README.md b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/README.md new file mode 100644 index 000000000000..da0f72d3487d --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/README.md @@ -0,0 +1,267 @@ + + +# Quantile Function + +> [Log-logistic][log-logistic-distribution] distribution [quantile function][quantile-function]. + +
+ +The [quantile function][quantile-function] for a [log-logistic][log-logistic-distribution] random variable is + + + +```math +Q(p;\alpha,\beta) = \alpha\left(\dfrac{p}{1-p}\right)^{1/\beta} +``` + + + + + +for `0 <= p <= 1`, where `α > 0` is the scale parameter and `β > 0` is the shape parameter. + +
+ + + +
+ +## Usage + +```javascript +var quantile = require( '@stdlib/stats/base/dists/log-logistic/quantile' ); +``` + +#### quantile( p, alpha, beta ) + +Evaluates the [quantile function][quantile-function] for a [log-logistic][log-logistic-distribution] distribution with parameters `alpha` (scale parameter) and `beta` (shape parameter). + +```javascript +var y = quantile( 0.5, 1.0, 1.0 ); +// returns 1.0 + +y = quantile( 0.5, 4.0, 2.0 ); +// returns 4.0 + +y = quantile( 0.8, 1.0, 1.0 ); +// returns 4.0 +``` + +If provided a probability `p` outside the interval `[0,1]`, the function returns `NaN`. + +```javascript +var y = quantile( 1.9, 1.0, 1.0 ); +// returns NaN + +y = quantile( -0.1, 1.0, 1.0 ); +// returns NaN +``` + +If provided `NaN` as any argument, the function returns `NaN`. + +```javascript +var y = quantile( NaN, 1.0, 1.0 ); +// returns NaN + +y = quantile( 0.5, NaN, 1.0 ); +// returns NaN + +y = quantile( 0.5, 1.0, NaN ); +// returns NaN +``` + +If provided `alpha <= 0`, the function returns `NaN`. + +```javascript +var y = quantile( 0.5, -1.0, 1.0 ); +// returns NaN +``` + +If provided `beta <= 0`, the function returns `NaN`. + +```javascript +var y = quantile( 0.5, 1.0, -1.0 ); +// returns NaN +``` + +#### quantile.factory( alpha, beta ) + +Returns a function for evaluating the [quantile function][quantile-function] of a [log-logistic][log-logistic-distribution] distribution with parameters `alpha` (scale parameter) and `beta` (shape parameter). + +```javascript +var myQuantile = quantile.factory( 1.0, 1.0 ); + +var y = myQuantile( 0.5 ); +// returns 1.0 + +y = myQuantile( 0.8 ); +// returns 4.0 +``` + +
+ + + +
+ +## Examples + + + +```javascript +var uniform = require( '@stdlib/random/array/uniform' ); +var logEachMap = require( '@stdlib/console/log-each-map' ); +var quantile = require( '@stdlib/stats/base/dists/log-logistic/quantile' ); + +var opts = { + 'dtype': 'float64' +}; +var p = uniform( 10, 0.0, 1.0, opts ); +var alpha = uniform( 10, 0.1, 10.0, opts ); +var beta = uniform( 10, 0.1, 10.0, opts ); + +logEachMap( 'p: %0.4f, alpha: %0.4f, beta: %0.4f, Q(p;alpha,beta): %0.4f', p, alpha, beta, quantile ); +``` + +
+ + + + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/stats/base/dists/log-logistic/quantile.h" +``` + +#### stdlib_base_dists_log_logistic_quantile( p, alpha, beta ) + +Evaluates the [quantile function][quantile-function] for a [log-logistic][log-logistic-distribution] distribution with parameters `alpha` (scale parameter) and `beta` (shape parameter). + +```c +double out = stdlib_base_dists_log_logistic_quantile( 0.5, 1.0, 1.0 ); +// returns 1.0 +``` + +The function accepts the following arguments: + +- **p**: `[in] double` input probability. +- **alpha**: `[in] double` scale parameter. +- **beta**: `[in] double` shape parameter. + +```c +double stdlib_base_dists_log_logistic_quantile( const double p, const double alpha, const double beta ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/stats/base/dists/log-logistic/quantile.h" +#include +#include + +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 ) { + double alpha; + double beta; + double p; + double y; + int i; + + for ( i = 0; i < 25; i++ ) { + p = random_uniform( 0.0, 1.0 ); + alpha = random_uniform( 0.1, 10.0 ); + beta = random_uniform( 0.1, 10.0 ); + y = stdlib_base_dists_log_logistic_quantile( p, alpha, beta ); + printf( "p: %lf, alpha: %lf, beta: %lf, Q(p;alpha,beta): %lf\n", p, alpha, beta, y ); + } +} +``` + +
+ + + +
+ + + + + +
+ +
+ + + + + +
+ +[log-logistic-distribution]: https://en.wikipedia.org/wiki/Log-logistic_distribution + +[quantile-function]: https://en.wikipedia.org/wiki/Quantile_function + +
+ + diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/benchmark.js new file mode 100644 index 000000000000..36957309c7dc --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/benchmark.js @@ -0,0 +1,95 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var EPS = require( '@stdlib/constants/float64/eps' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var quantile = require( './../lib' ); + + +// MAIN // + +bench( pkg, function benchmark( b ) { + var alpha; + var beta; + var opts; + var p; + var y; + var i; + + opts = { + 'dtype': 'float64' + }; + p = uniform( 100, EPS, 1.0 - EPS, opts ); + alpha = uniform( 100, EPS, 5.0, opts ); + beta = uniform( 100, EPS, 5.0, opts ); + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = quantile( p[ i % p.length ], alpha[ i % alpha.length ], beta[ i % beta.length ] ); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); + +bench( format( '%s::factory', pkg ), function benchmark( b ) { + var myquantile; + var alpha; + var beta; + var opts; + var p; + var y; + var i; + + alpha = 2.0; + beta = 3.0; + myquantile = quantile.factory( alpha, beta ); + + opts = { + 'dtype': 'float64' + }; + p = uniform( 100, EPS, 1.0 - EPS, opts ); + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = myquantile( p[ i % p.length ] ); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/benchmark.native.js new file mode 100644 index 000000000000..347acf698341 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/benchmark.native.js @@ -0,0 +1,71 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var EPS = require( '@stdlib/constants/float64/eps' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var quantile = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( quantile instanceof Error ) +}; + + +// MAIN // + +bench( format( '%s::native', pkg ), opts, function benchmark( b ) { + var alpha; + var beta; + var opts; + var p; + var y; + var i; + + opts = { + 'dtype': 'float64' + }; + p = uniform( 100, EPS, 1.0 - EPS, opts ); + alpha = uniform( 100, EPS, 5.0, opts ); + beta = uniform( 100, EPS, 5.0, opts ); + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = quantile( p[ i % p.length ], alpha[ i % alpha.length ], beta[ i % beta.length ] ); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/c/Makefile b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/c/Makefile new file mode 100644 index 000000000000..979768abbcec --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 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. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := benchmark.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled benchmarks. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/c/benchmark.c new file mode 100644 index 000000000000..f4f5c4abe6b2 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/benchmark/c/benchmark.c @@ -0,0 +1,141 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/log-logistic/quantile.h" +#include "stdlib/constants/float64/eps.h" +#include +#include +#include +#include +#include + +#define NAME "log-logistic-quantile" +#define ITERATIONS 1000000 +#define REPEATS 3 + +/** +* Prints the TAP version. +*/ +static void print_version( void ) { + printf( "TAP version 13\n" ); +} + +/** +* Prints the TAP summary. +* +* @param total total number of tests +* @param passing total number of passing tests +*/ +static void print_summary( int total, int passing ) { + printf( "#\n" ); + printf( "1..%d\n", total ); + printf( "# total %d\n", total ); + printf( "# pass %d\n", passing ); + printf( "#\n" ); + printf( "# ok\n" ); +} + +/** +* Prints benchmark results. +* +* @param elapsed elapsed time in seconds +*/ +static void print_results( double elapsed ) { + double rate = (double)ITERATIONS / elapsed; + printf( " ---\n" ); + printf( " iterations: %d\n", ITERATIONS ); + printf( " elapsed: %0.9f\n", elapsed ); + printf( " rate: %0.9f\n", rate ); + printf( " ...\n" ); +} + +/** +* Returns a clock time. +* +* @return clock time +*/ +static double tic( void ) { + struct timeval now; + gettimeofday( &now, NULL ); + return (double)now.tv_sec + (double)now.tv_usec/1.0e6; +} + +/** +* Generates a random number on the interval [min,max). +* +* @param min minimum value (inclusive) +* @param max maximum value (exclusive) +* @return random number +*/ +static double random_uniform( const double min, const double max ) { + double v = (double)rand() / ( (double)RAND_MAX + 1.0 ); + return min + ( v*(max-min) ); +} + +/** +* Runs a benchmark. +* +* @return elapsed time in seconds +*/ +static double benchmark( void ) { + double elapsed; + double alpha[ 100 ]; + double beta[ 100 ]; + double p[ 100 ]; + double y; + double t; + int i; + + for ( i = 0; i < 100; i++ ) { + p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 1.0 - STDLIB_CONSTANT_FLOAT64_EPS ); + alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 ); + beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 ); + } + t = tic(); + for ( i = 0; i < ITERATIONS; i++ ) { + y = stdlib_base_dists_log_logistic_quantile( p[ i%100 ], alpha[ i%100 ], beta[ i%100 ] ); + if ( y != y ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( y != y ) { + printf( "should not return NaN\n" ); + } + return elapsed; +} + +/** +* Main execution sequence. +*/ +int main( void ) { + double elapsed; + int i; + + srand( time( NULL ) ); + + print_version(); + for ( i = 0; i < REPEATS; i++ ) { + printf( "# c::%s\n", NAME ); + elapsed = benchmark(); + print_results( elapsed ); + printf( "ok %d benchmark finished\n", i+1 ); + } + print_summary( REPEATS, REPEATS ); +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/binding.gyp b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/binding.gyp new file mode 100644 index 000000000000..0d6508a12e99 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/binding.gyp @@ -0,0 +1,170 @@ +# @license Apache-2.0 +# +# Copyright (c) 2026 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. + +# A `.gyp` file for building a Node.js native add-on. +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # List of files to include in this file: + 'includes': [ + './include.gypi', + ], + + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Target name should match the add-on export name: + 'addon_target_name%': 'addon', + + # Set variables based on the host OS: + 'conditions': [ + [ + 'OS=="win"', + { + # Define the object file suffix: + 'obj': 'obj', + }, + { + # Define the object file suffix: + 'obj': 'o', + } + ], # end condition (OS=="win") + ], # end conditions + }, # end variables + + # Define compile targets: + 'targets': [ + + # Target to generate an add-on: + { + # The target name should match the add-on export name: + 'target_name': '<(addon_target_name)', + + # Define dependencies: + 'dependencies': [], + + # Define directories which contain relevant include headers: + 'include_dirs': [ + # Local include directory: + '<@(include_dirs)', + ], + + # List of source files: + 'sources': [ + '<@(src_files)', + ], + + # Settings which should be applied when a target's object files are used as linker input: + 'link_settings': { + # Define libraries: + 'libraries': [ + '<@(libraries)', + ], + + # Define library directories: + 'library_dirs': [ + '<@(library_dirs)', + ], + }, + + # C/C++ compiler flags: + 'cflags': [ + # Enable commonly used warning options: + '-Wall', + + # Aggressive optimization: + '-O3', + ], + + # C specific compiler flags: + 'cflags_c': [ + # Specify the C standard to which a program is expected to conform: + '-std=c99', + ], + + # C++ specific compiler flags: + 'cflags_cpp': [ + # Specify the C++ standard to which a program is expected to conform: + '-std=c++11', + ], + + # Linker flags: + 'ldflags': [], + + # Apply conditions based on the host OS: + 'conditions': [ + [ + 'OS=="mac"', + { + # Linker flags: + 'ldflags': [ + '-undefined dynamic_lookup', + '-Wl,-no-pie', + '-Wl,-search_paths_first', + ], + }, + ], # end condition (OS=="mac") + [ + 'OS!="win"', + { + # C/C++ flags: + 'cflags': [ + # Generate platform-independent code: + '-fPIC', + ], + }, + ], # end condition (OS!="win") + ], # end conditions + }, # end target <(addon_target_name) + + # Target to copy a generated add-on to a standard location: + { + 'target_name': 'copy_addon', + + # Declare that the output of this target is not linked: + 'type': 'none', + + # Define dependencies: + 'dependencies': [ + # Require that the add-on be generated before building this target: + '<(addon_target_name)', + ], + + # Define a list of actions: + 'actions': [ + { + 'action_name': 'copy_addon', + 'message': 'Copying addon...', + + # Explicitly list the inputs in the command-line invocation below: + 'inputs': [], + + # Declare the expected outputs: + 'outputs': [ + '<(addon_output_dir)/<(addon_target_name).node', + ], + + # Define the command-line invocation: + 'action': [ + 'cp', + '<(PRODUCT_DIR)/<(addon_target_name).node', + '<(addon_output_dir)/<(addon_target_name).node', + ], + }, + ], # end actions + }, # end target copy_addon + ], # end targets +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/img/equation_log_logistic_quantile_function.svg b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/img/equation_log_logistic_quantile_function.svg new file mode 100644 index 000000000000..491bc2a3458c --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/img/equation_log_logistic_quantile_function.svg @@ -0,0 +1,2 @@ + + \ No newline at end of file diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/repl.txt b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/repl.txt new file mode 100644 index 000000000000..5d6f535bc272 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/repl.txt @@ -0,0 +1,72 @@ + +{{alias}}( p, alpha, beta ) + Evaluates the quantile function for a log-logistic distribution with + scale parameter `alpha` and shape parameter `beta` at a probability `p`. + + If provided `alpha <= 0` or `beta <= 0`, the function returns `NaN`. + + If `p < 0` or `p > 1`, the function returns `NaN`. + + Parameters + ---------- + p: number + Input probability. + + alpha: number + Scale parameter. + + beta: number + Shape parameter. + + Returns + ------- + out: number + Evaluated quantile function. + + Examples + -------- + > var y = {{alias}}( 0.5, 1.0, 1.0 ) + 1.0 + > y = {{alias}}( 0.5, 4.0, 2.0 ) + 4.0 + > y = {{alias}}( 0.8, 1.0, 1.0 ) + 4.0 + > y = {{alias}}( NaN, 1.0, 1.0 ) + NaN + > y = {{alias}}( 0.5, NaN, 1.0 ) + NaN + > y = {{alias}}( 0.5, 1.0, NaN ) + NaN + > y = {{alias}}( 0.5, -1.0, 1.0 ) + NaN + > y = {{alias}}( 0.5, 1.0, -1.0 ) + NaN + + +{{alias}}.factory( alpha, beta ) + Returns a function for evaluating the quantile function of a + log-logistic distribution with scale parameter `alpha` and + shape parameter `beta`. + + Parameters + ---------- + alpha: number + Scale parameter. + + beta: number + Shape parameter. + + Returns + ------- + quantile: Function + Quantile function. + + Examples + -------- + > var myQuantile = {{alias}}.factory( 1.0, 1.0 ); + > var y = myQuantile( 0.5 ) + 1.0 + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/types/index.d.ts b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/types/index.d.ts new file mode 100644 index 000000000000..960fe1056124 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/types/index.d.ts @@ -0,0 +1,128 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 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. +*/ + +// TypeScript Version: 4.1 + +/** +* Evaluates the quantile function for a log-logistic distribution. +* +* @param p - input probability +* @returns evaluated quantile function +*/ +type Unary = ( p: number ) => number; + +/** +* Interface for the quantile function of a log-logistic distribution. +*/ +interface Quantile { + /** + * Evaluates the quantile function for a log-logistic distribution with scale parameter `alpha` and shape parameter `beta` at a probability `p`. + * + * ## Notes + * + * - If provided `alpha <= 0` or `beta <= 0`, the function returns `NaN`. + * - If `p < 0` or `p > 1`, the function returns `NaN`. + * + * @param p - input probability + * @param alpha - scale parameter + * @param beta - shape parameter + * @returns evaluated quantile function + * + * @example + * var y = quantile( 0.5, 1.0, 1.0 ); + * // returns 1.0 + * + * @example + * var y = quantile( 0.5, 4.0, 2.0 ); + * // returns 4.0 + * + * @example + * var y = quantile( 0.8, 1.0, 1.0 ); + * // returns 4.0 + * + * @example + * var y = quantile( 0.0, 1.0, 1.0 ); + * // returns 0.0 + * + * @example + * var y = quantile( 1.0, 1.0, 1.0 ); + * // returns Infinity + * + * @example + * var y = quantile( NaN, 1.0, 1.0 ); + * // returns NaN + * + * @example + * var y = quantile( 0.5, NaN, 1.0 ); + * // returns NaN + * + * @example + * var y = quantile( 0.5, 1.0, NaN ); + * // returns NaN + * + * @example + * var y = quantile( 0.5, -1.0, 1.0 ); + * // returns NaN + * + * @example + * var y = quantile( 0.5, 1.0, -1.0 ); + * // returns NaN + */ + ( p: number, alpha: number, beta: number ): number; + + /** + * Returns a function for evaluating the quantile function for a log-logistic distribution. + * + * @param alpha - scale parameter + * @param beta - shape parameter + * @returns quantile function + * + * @example + * var myquantile = quantile.factory( 1.0, 1.0 ); + * var y = myquantile( 0.5 ); + * // returns 1.0 + * + * y = myquantile( 0.8 ); + * // returns 4.0 + */ + factory( alpha: number, beta: number ): Unary; +} + +/** +* Log-logistic distribution quantile function. +* +* @param p - input probability +* @param alpha - scale parameter +* @param beta - shape parameter +* @returns evaluated quantile function +* +* @example +* var y = quantile( 0.5, 1.0, 1.0 ); +* // returns 1.0 +* +* @example +* var myquantile = quantile.factory( 1.0, 1.0 ); +* var y = myquantile( 0.5 ); +* // returns 1.0 +*/ +declare var quantile: Quantile; + + +// EXPORTS // + +export = quantile; diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/types/test.ts b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/types/test.ts new file mode 100644 index 000000000000..fe24843415d0 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/docs/types/test.ts @@ -0,0 +1,119 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 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. +*/ + +import quantile = require( './index' ); + + +// TESTS // + +// The function returns a number... +{ + quantile( 0.5, 2, 4 ); // $ExpectType number + quantile( 0.8, 2, 8 ); // $ExpectType number +} + +// The compiler throws an error if the function is provided values other than three numbers... +{ + quantile( true, 3, 6 ); // $ExpectError + quantile( false, 2, 4 ); // $ExpectError + quantile( '5', 1, 2 ); // $ExpectError + quantile( [], 1, 2 ); // $ExpectError + quantile( {}, 2, 4 ); // $ExpectError + quantile( ( x: number ): number => x, 2, 4 ); // $ExpectError + + quantile( 0.5, true, 12 ); // $ExpectError + quantile( 0.5, false, 12 ); // $ExpectError + quantile( 0.5, '5', 10 ); // $ExpectError + quantile( 0.5, [], 16 ); // $ExpectError + quantile( 0.5, {}, 18 ); // $ExpectError + quantile( 0.5, ( x: number ): number => x, 16 ); // $ExpectError + + quantile( 0.5, 5, true ); // $ExpectError + quantile( 0.5, 5, false ); // $ExpectError + quantile( 0.5, 2, '5' ); // $ExpectError + quantile( 0.5, 4, [] ); // $ExpectError + quantile( 0.5, 4, {} ); // $ExpectError + quantile( 0.5, 5, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided an unsupported number of arguments... +{ + quantile(); // $ExpectError + quantile( 0.5 ); // $ExpectError + quantile( 0.5, 0 ); // $ExpectError + quantile( 0.5, 0, 4, 1 ); // $ExpectError +} + +// Attached to main export is a `factory` method which returns a function... +{ + quantile.factory( 3, 4 ); // $ExpectType Unary +} + +// The `factory` method returns a function which returns a number... +{ + const fcn = quantile.factory( 3, 4 ); + fcn( 0.5 ); // $ExpectType number +} + +// The compiler throws an error if the function returned by the `factory` method is provided invalid arguments... +{ + const fcn = quantile.factory( 3, 4 ); + fcn( true ); // $ExpectError + fcn( false ); // $ExpectError + fcn( '5' ); // $ExpectError + fcn( [] ); // $ExpectError + fcn( {} ); // $ExpectError + fcn( ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function returned by the `factory` method is provided an unsupported number of arguments... +{ + const fcn = quantile.factory( 3, 4 ); + fcn(); // $ExpectError + fcn( 0.5, 0 ); // $ExpectError + fcn( 0.5, 0, 1 ); // $ExpectError +} + +// The compiler throws an error if the `factory` method is provided values other than two numbers... +{ + quantile.factory( true, 3 ); // $ExpectError + quantile.factory( false, 2 ); // $ExpectError + quantile.factory( '5', 1 ); // $ExpectError + quantile.factory( [], 1 ); // $ExpectError + quantile.factory( {}, 2 ); // $ExpectError + quantile.factory( ( x: number ): number => x, 2 ); // $ExpectError + + quantile.factory( 9, true ); // $ExpectError + quantile.factory( 9, false ); // $ExpectError + quantile.factory( 5, '5' ); // $ExpectError + quantile.factory( 8, [] ); // $ExpectError + quantile.factory( 9, {} ); // $ExpectError + quantile.factory( 8, ( x: number ): number => x ); // $ExpectError + + quantile.factory( [], true ); // $ExpectError + quantile.factory( {}, false ); // $ExpectError + quantile.factory( false, '5' ); // $ExpectError + quantile.factory( {}, [] ); // $ExpectError + quantile.factory( '5', ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the `factory` method is provided an unsupported number of arguments... +{ + quantile.factory( 0 ); // $ExpectError + quantile.factory( 0, 4, 8 ); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/c/Makefile b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/c/Makefile new file mode 100644 index 000000000000..c8f8e9a1517b --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 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. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := example.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled examples. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/c/example.c new file mode 100644 index 000000000000..4a63a2a5bcb5 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/c/example.c @@ -0,0 +1,42 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/log-logistic/quantile.h" +#include +#include + +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 ) { + double alpha; + double beta; + double p; + double y; + int i; + + for ( i = 0; i < 25; i++ ) { + p = random_uniform( 0.0, 1.0 ); + alpha = random_uniform( 0.1, 10.0 ); + beta = random_uniform( 0.1, 10.0 ); + y = stdlib_base_dists_log_logistic_quantile( p, alpha, beta ); + printf( "p: %lf, alpha: %lf, beta: %lf, Q(p;alpha,beta): %lf\n", p, alpha, beta, y ); + } +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/index.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/index.js new file mode 100644 index 000000000000..67a5fe3a6ac6 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/examples/index.js @@ -0,0 +1,32 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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'; + +var uniform = require( '@stdlib/random/array/uniform' ); +var logEachMap = require( '@stdlib/console/log-each-map' ); +var quantile = require( './../lib' ); + +var opts = { + 'dtype': 'float64' +}; +var p = uniform( 10, 0.0, 1.0, opts ); +var alpha = uniform( 10, 0.1, 10.0, opts ); +var beta = uniform( 10, 0.1, 10.0, opts ); + +logEachMap( 'p: %0.4f, alpha: %0.4f, beta: %0.4f, Q(p;alpha,beta): %0.4f', p, alpha, beta, quantile ); diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/include.gypi b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/include.gypi new file mode 100644 index 000000000000..bee8d41a2caf --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/include.gypi @@ -0,0 +1,53 @@ +# @license Apache-2.0 +# +# Copyright (c) 2026 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. + +# A GYP include file for building a Node.js native add-on. +# +# Main documentation: +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Source directory: + 'src_dir': './src', + + # Include directories: + 'include_dirs': [ + ' 1.0 ) { + return NaN; + } + if ( p === 0.0 ) { + return 0.0; + } + if ( p === 1.0 ) { + return Infinity; + } + return alpha * pow( p / ( 1.0 - p ), 1.0 / beta ); + } +} + + +// EXPORTS // + +module.exports = factory; diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/index.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/index.js new file mode 100644 index 000000000000..975f1dc888f2 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/index.js @@ -0,0 +1,51 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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'; + +/** +* Log-logistic distribution quantile function. +* +* @module @stdlib/stats/base/dists/log-logistic/quantile +* +* @example +* var quantile = require( '@stdlib/stats/base/dists/log-logistic/quantile' ); +* +* var y = quantile( 0.5, 1.0, 1.0 ); +* // returns 1.0 +* +* var myQuantile = quantile.factory( 1.0, 1.0 ); +* var y = myQuantile( 0.5 ); +* // returns 1.0 +*/ + +// 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; diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/main.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/main.js new file mode 100644 index 000000000000..3c01ddf6bfc3 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/main.js @@ -0,0 +1,101 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 pow = require( '@stdlib/math/base/special/pow' ); + + +// MAIN // + +/** +* Evaluates the quantile function for a log-logistic distribution with scale parameter `alpha` and shape parameter `beta` at a probability `p`. +* +* @param {Probability} p - input probability +* @param {PositiveNumber} alpha - scale parameter +* @param {PositiveNumber} beta - shape parameter +* @returns {number} evaluated quantile function +* +* @example +* var y = quantile( 0.5, 1.0, 1.0 ); +* // returns 1.0 +* +* @example +* var y = quantile( 0.5, 4.0, 2.0 ); +* // returns 4.0 +* +* @example +* var y = quantile( 0.8, 1.0, 1.0 ); +* // returns 4.0 +* +* @example +* var y = quantile( -0.5, 1.0, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( 1.5, 1.0, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( NaN, 1.0, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( 0.5, NaN, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( 0.5, 1.0, NaN ); +* // returns NaN +* +* @example +* var y = quantile( 0.5, -1.0, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( 0.5, 1.0, -1.0 ); +* // returns NaN +*/ +function quantile( p, alpha, beta ) { + if ( + isnan( p ) || + isnan( alpha ) || + isnan( beta ) || + alpha <= 0.0 || + beta <= 0.0 || + p < 0.0 || + p > 1.0 + ) { + return NaN; + } + if ( p === 0.0 ) { + return 0.0; + } + if ( p === 1.0 ) { + return Infinity; + } + return alpha * pow( p / ( 1.0 - p ), 1.0 / beta ); +} + + +// EXPORTS // + +module.exports = quantile; diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/native.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/native.js new file mode 100644 index 000000000000..67342fd981e9 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/lib/native.js @@ -0,0 +1,83 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 quantile function for a log-logistic distribution with scale parameter `alpha` and shape parameter `beta` at a probability `p`. +* +* @param {Probability} p - input probability +* @param {PositiveNumber} alpha - scale parameter +* @param {PositiveNumber} beta - shape parameter +* @returns {number} evaluated quantile function +* +* @example +* var y = quantile( 0.5, 1.0, 1.0 ); +* // returns 1.0 +* +* @example +* var y = quantile( 0.5, 4.0, 2.0 ); +* // returns 4.0 +* +* @example +* var y = quantile( 0.8, 1.0, 1.0 ); +* // returns 4.0 +* +* @example +* var y = quantile( -0.5, 1.0, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( 1.5, 1.0, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( NaN, 1.0, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( 0.5, NaN, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( 0.5, 1.0, NaN ); +* // returns NaN +* +* @example +* var y = quantile( 0.5, -1.0, 1.0 ); +* // returns NaN +* +* @example +* var y = quantile( 0.5, 1.0, -1.0 ); +* // returns NaN +*/ +function quantile( p, alpha, beta ) { + return addon( p, alpha, beta ); +} + + +// EXPORTS // + +module.exports = quantile; diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/manifest.json b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/manifest.json new file mode 100644 index 000000000000..444a3a08bea9 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/manifest.json @@ -0,0 +1,80 @@ +{ + "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/napi/ternary", + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/pow" + ] + }, + { + "task": "benchmark", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/pow", + "@stdlib/constants/float64/eps" + ] + }, + { + "task": "examples", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/pow" + ] + } + ] +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/package.json b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/package.json new file mode 100644 index 000000000000..63cd03805d3f --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/package.json @@ -0,0 +1,70 @@ +{ + "name": "@stdlib/stats/base/dists/log-logistic/quantile", + "version": "0.0.0", + "description": "Log-logistic distribution quantile function.", + "license": "Apache-2.0", + "author": { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + }, + "contributors": [ + { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + } + ], + "main": "./lib", + "gypfile": true, + "directories": { + "benchmark": "./benchmark", + "doc": "./docs", + "example": "./examples", + "include": "./include", + "lib": "./lib", + "src": "./src", + "test": "./test" + }, + "types": "./docs/types", + "scripts": {}, + "homepage": "https://github.com/stdlib-js/stdlib", + "repository": { + "type": "git", + "url": "git://github.com/stdlib-js/stdlib.git" + }, + "bugs": { + "url": "https://github.com/stdlib-js/stdlib/issues" + }, + "dependencies": {}, + "devDependencies": {}, + "engines": { + "node": ">=0.10.0", + "npm": ">2.7.0" + }, + "os": [ + "aix", + "darwin", + "freebsd", + "linux", + "macos", + "openbsd", + "sunos", + "win32", + "windows" + ], + "keywords": [ + "stdlib", + "stdmath", + "statistics", + "stats", + "distribution", + "dist", + "probability", + "quantile", + "inverse", + "location-scale", + "log-logistic function", + "log-logistic", + "univariate", + "continuous" + ] +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/Makefile b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/Makefile new file mode 100644 index 000000000000..2caf905cedbe --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 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. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + + +# RULES # + +#/ +# Removes generated files for building an add-on. +# +# @example +# make clean-addon +#/ +clean-addon: + $(QUIET) -rm -f *.o *.node + +.PHONY: clean-addon + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: clean-addon + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/addon.c b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/addon.c new file mode 100644 index 000000000000..9c66348ea5ad --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/addon.c @@ -0,0 +1,23 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/log-logistic/quantile.h" +#include "stdlib/math/base/napi/ternary.h" + +// cppcheck-suppress shadowFunction +STDLIB_MATH_BASE_NAPI_MODULE_DDD_D( stdlib_base_dists_log_logistic_quantile ) diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/main.c new file mode 100644 index 000000000000..44fe1989984b --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/src/main.c @@ -0,0 +1,54 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/log-logistic/quantile.h" +#include "stdlib/math/base/assert/is_nan.h" +#include "stdlib/math/base/special/pow.h" + +/** +* Evaluates the quantile function for a log-logistic distribution with scale parameter `alpha` and shape parameter `beta` at a probability `p`. +* +* @param p input probability +* @param alpha scale parameter +* @param beta shape parameter +* @return evaluated quantile function +* +* @example +* double y = stdlib_base_dists_log_logistic_quantile( 0.5, 1.0, 1.0 ); +* // returns 1.0 +*/ +double stdlib_base_dists_log_logistic_quantile( const double p, const double alpha, const double beta ) { + if ( + stdlib_base_is_nan( p ) || + stdlib_base_is_nan( alpha ) || + stdlib_base_is_nan( beta ) || + alpha <= 0.0 || + beta <= 0.0 || + p < 0.0 || + p > 1.0 + ) { + return 0.0 / 0.0; + } + if ( p == 0.0 ) { + return 0.0; + } + if ( p == 1.0 ) { + return 1.0 / 0.0; + } + return alpha * stdlib_base_pow( p / ( 1.0 - p ), 1.0 / beta ); +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/REQUIRE b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/REQUIRE new file mode 100644 index 000000000000..98be20b58ed3 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/REQUIRE @@ -0,0 +1,3 @@ +Distributions 0.23.8 +julia 1.5 +JSON 0.21 diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/data.json b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/data.json new file mode 100644 index 000000000000..a007196a2ec9 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/data.json @@ -0,0 +1 @@ +{"p":[0.776474186285538,0.02405624520990446,0.9416507715357758,0.18777713513344874,0.27598017810109643,0.5254118369008591,0.09524951515394808,0.11108201860899758,0.5042970800205191,0.28066120794282146,0.17282353853578125,0.6460121639656222,0.9829171239381662,0.7598036063194165,0.32960098109230973,0.9373893109626752,0.7438147095726735,0.1434010812806501,0.779956503517022,0.5946059800457781,0.46101226359300496,0.971151958514181,0.00855142467780512,0.5666529441902006,0.11881846808487762,0.5382033617883939,0.16331514640864664,0.3010785491708008,0.35812893484054076,0.5703898094220099,0.41679079975062017,0.4945928483102975,0.6677120620157184,0.43660089174220973,0.13689624058677397,0.5011386631595889,0.059040030632449225,0.12076575390055866,0.8425971000756862,0.7524141118640424,0.738359530542634,0.7474696351000754,0.7670834986504242,0.542179339466659,0.4210151607773531,0.3337590342756467,0.7433246974760934,0.5014079117237105,0.7053930212323387,0.8776699101545384,0.4851218648446063,0.024680567350874538,0.5699608230666762,0.9168118820778224,0.7635899447600549,0.23598815018911123,0.9746630222378116,0.5182316607106998,0.18806949865436184,0.7365925349850517,0.15521820403251962,0.21897383180444474,0.7094127732249729,0.5222437419638117,0.26225510881631653,0.8624222516733602,0.13161012452825283,0.6493733354358049,0.7799113095468373,0.15136090254729773,0.19016108742248972,0.8145074576941029,0.7971282211310773,0.2977040189531197,0.2062906869286647,0.7236505456825613,0.74515058998833,0.3171310828413043,0.2396745470878765,0.7427414442949256,0.6545230322254969,0.15652706554234497,0.4473558650812781,0.8265236062013759,0.058890662159235996,0.11664745274275057,0.9425368176415565,0.1294027332741312,0.7327688243172665,0.9276977143233638,0.5215422899505371,0.7221984640338568,0.501334159266249,0.08234526851828838,0.1589412381653046,0.1592344145306286,0.10253609451324985,0.9993804456519582,0.44297499421752407,0.13680239659986637,0.36457253210981344,0.9145968651626029,0.25847529869600794,0.8883871191163184,0.27656577878845323,0.4245989066885445,0.5173142766891496,0.8043680364341041,0.21606689637275944,0.5951959355252763,0.46460614818863744,0.200703485780803,0.6624201443144035,0.9872713565283531,0.12228206283958198,0.6489008990722536,0.9452947716585245,0.5684685438891701,0.07259874926982857,0.9962711912104545,0.8112989823262684,0.9702943306761487,0.6346481543371933,0.9073994990139977,0.7976897386949444,0.6152047780345602,0.9890828935891394,0.6634491428803502,0.9943192993801058,0.798907673918325,0.7250315600925963,0.6370010180632484,0.1264387643160747,0.4956920136062275,0.016403478819103645,0.3506878777702438,0.29623223381421004,0.6501617079551927,0.14081225249023643,0.729520108894915,0.04179842234330167,0.5203559646703086,0.4203063119716106,0.29364283639625643,0.5995852620315512,0.04516504419325784,0.4357944714755906,0.10366445713979022,0.4757542669882153,0.6475926370082505,0.24837778675720634,0.9171388673510104,0.41341695028297054,0.7402333311684006,0.13683688253064874,0.536789428011762,0.8907609350936261,0.8485427742594944,0.8900412167046374,0.4881928293499205,0.5782041601614101,0.30285056190210136,0.5485884351411068,0.60350172013683,0.05594731511842377,0.7678621191286097,0.39646323758453117,0.5162717279858053,0.12111496847851111,0.5098909259022126,0.16214760581279286,0.4580635867483054,0.7969251353101863,0.5840928763058053,0.3394709897182473,0.6570568584503405,0.06278752477710947,0.47469541506831514,0.06295424399016603,0.526230446614778,0.11009115470444386,0.18666299680000198,0.5141325818729205,0.620798611423953,0.47485999990564176,0.6541848480776447,0.23090170601236681,0.61616921511084,0.6201519036604763,0.14279188585534275,0.6335078591392511,0.6749140221599625,0.4408564122456369,0.04570688403670298,0.6428263038269153,0.6880057002458846,0.7381065675397211,0.29909268034130077,0.8494136524797087,0.8835167250583067,0.10172273281130662,0.32148202248299573,0.9471619996460583,0.9464786436759356,0.6118460430673791,0.6614210343673879,0.28636818972225764,0.1261428298150442,0.5650694340836111,0.24899096012632937,0.41128287015498355,0.3189597163936495,0.8876422475188144,0.8560988203149555,0.007373240759798394,0.47737391339476076,0.914927192931926,0.45487581402894905,0.9386141812490012,0.7395091324206309,0.30416917837159496,0.14048490834984406,0.5692759750327188,0.4852702386855494,0.24936650056172804,0.8096243798199427,0.3745082240971953,0.6750268024962365,0.10659073681874043,0.781349479922643,0.3861155144825108,0.45980225476386716,0.7396722632525732,0.9357823316854235,0.5776390736942416,0.5472216533848966,0.1161183129578498,0.7284136022613081,0.7255296849120727,0.27176085409723183,0.8109433891801019,0.8795431159125757,0.8097841599925748,0.1402303043315245,0.5183895119070834,0.06874617505312242,0.8277225824015177,0.18574225061517202,0.999795557748324,0.21568726342496847,0.06126079308785504,0.3376617948098528,0.8162620706305626,0.20472785222004464,0.339026772572661,0.48459745882038496,0.5023662940814801,0.6437574672889841,0.3750775178548287,0.2017912824533925,0.8506376195851548,0.2552061878353744,0.24195142341813858,0.9685688667553429,0.08793543117312796,0.995908871069987,0.3185399444676942,0.444970353771988,0.8094522875509627,0.2456578613253002,0.02046608407925854,0.5605497008038924,0.8732735315351428,0.8349890139926044,0.582491352836197,0.26831329075502053,0.21986066266009718,0.19744089918757757,0.07489932791061493,0.17866445151414023,0.8923204964947675,0.31270022404818776,0.19736367784973408,0.09092264892745727,0.5599765379926418,0.15896531226022814,0.3669846030416879,0.28165033336203993,0.30442158506428885,0.5535644927345742,0.9012288036437601,0.0016999744981673892,0.5157667858597903,0.504175307912661,0.9025533059999704,0.3234052904128544,0.8372208343704102,0.2851953271827432,0.8573168833677614,0.9274005289512877,0.6691557901773455,0.4438957234466483,0.025850994717384168,0.4695297549736517,0.5585108922525637,0.15043637915197716,0.8280435183597952,0.7155009816983646,0.4329191492911528,0.0864666476156366,0.4290261367094945,0.4816943410085792,0.7351760467000625,0.3019956993736759,0.9902733624787667,0.8841516164148113,0.8714598122801285,0.8723459104406746,0.5382570668857998,0.34950924658985816,0.42883310825340404,0.9065065982293459,0.6724283890741145,0.5970017746470606,0.7198706328494258,0.4731176316313852,0.11785696423057335,0.6798894150286614,0.3451657185696302,0.8609018522358209,0.8489504295508046,0.12792362531644152,0.35348673331678016,0.11430864385858319,0.7698864350862427,0.5821209313901312,0.943171218679971,0.010391283255508643,0.4096941478699514,0.9332709930065148,0.9756971000686432,0.3859536004606574,0.39246964428864173,0.654625848460832,0.3093479335591416,0.7991280724368033,0.23746892782270224,0.13366649413163634,0.5598076068321698,0.38881740669902953,0.22413835937791338,0.0200388877597244,0.08952786525190515,0.33488253058906636,0.7913375542351022,0.6662303459044687,0.018522190233664748,0.13138503561259984,0.11962750344599105,0.16300166702421026,0.8644122976976836,0.39970547828675196,0.30394660999839496,0.023791930425938412,0.25471413644229046,0.9339617082077225,0.12996739866721496,0.2453365254108244,0.19451906451987644,0.985487097544548,0.7754006437102086,0.03261175126211491,0.8373386974019834,0.7706627972814695,0.7094587863228052,0.771631563716499,0.7928577339884229,0.40271716211534603,0.0005318735535131403,0.6956196824006733,0.3733203325669485,0.34220608795190666,0.27860156416375137,0.1799602801208371,0.40079044501349026,0.7272257661609023,0.9485768295794812,0.8030881310870854,0.36920337461698693,0.9563177115092472,0.8290148370127062,0.8293615527117424,0.5346019904553998,0.8810668871866758,0.4527536754934174,0.26315973010075133,0.7717338749962734,0.6298623271011174,0.9333911586372408,0.5546364496414142,0.9428995337510888,0.023703155362768236,0.32495732803483524,0.02453056881227833,0.06094936790316208,0.5021169825185716,0.7357403278763315,0.6005459602142229,0.7546427123384816,0.5539938012944343,0.3559835067777619,0.5399744942707962,0.5734050589685233,0.8182077154186511,0.7039897402358481,0.09631329677615974,0.9295069283892576,0.7042972152694233,0.876791444515586,0.5725600452475728,0.8732634398186643,0.8870573854281729,0.6930626099721409,0.8471745898040854,0.6384552470355758,0.6150989415591586,0.1187693031413859,0.6992439255433451,0.8526833612398712,0.6390449753099967,0.007452067092083903,0.32719563771995097,0.08909550537220379,0.4538513650811178,0.2913435885932629,0.8952666351302325,0.7309601402469932,0.18398123282769308,0.006780813225261806,0.9718292757986347,0.1367921654329911,0.32193822118259785,0.9140637644712313,0.357738958668093,0.017962778477789598,0.5479247361819726,0.7465998195049288,0.02088251255456597,0.5069112809091776,0.24236901937827038,0.3276885655581008,0.09992092311742828,0.7458014435178333,0.48234472982971366,0.38535841144151406,0.9017218346909457,0.23084506365023416,0.3505651375477994,0.955901454238856,0.5146874407247335,0.8285166786207274,0.301886451698699,0.24561293781508553,0.9466148028103847,0.5585272620935957,0.23506727838615005,0.23588818996637761,0.3962808126759503,0.14259332305176797,0.5280583830838401,0.8056709086195635,0.31532777559154024,0.74083687521866,0.9794616809237551,0.287782795126061,0.9674101195892826,0.4155859744198125,0.2783425871946037,0.6512410098604244,0.5058712216386474,0.3833290660235553,0.9944988625176714,0.40796422435728075,0.7436053942575065,0.593892968328064,0.4122656955726096,0.838416594607001,0.42528090810660624,0.14300210722804207,0.23901706443007487,0.46188741157992164,0.7252829012622775,0.9524605554878507,0.5788243852251123,0.2984919750895918,0.32804519484947137,0.8934500194919892,0.8618289734648645,0.700982039522174,0.7607321354666776,0.7282974476644528,0.6261502320516654,0.9580782575035498,0.717365966488821,0.3924036205669881,0.8179895890470339,0.5283180708881088,0.6753878568961411,0.6555479224917156,0.44245725611687914,0.7105702624878096,0.24043976943012324,0.756967290940189,0.5791432369260963,0.998228271495301,0.5277801662850015,0.7551031978095435,0.10695521403182107,0.27164049430085324,0.8309505629919438,0.9042361333547035,0.7978621719039569,0.41744247166413273,0.4519191780957037,0.8916008459973772,0.2844358530701929,0.1252012908850959,0.7422990607347189,0.3614903402059333,0.8387187962577589,0.8586627626463068,0.37563340985537197,0.5042645013599547,0.9704134582198956,0.2078118679517782,0.45463630242239694,0.20634644735239993,0.6448207970875899,0.10328072853955184,0.7020270884280637,0.5237421020196222,0.19630754362667857,0.7949334978226867,0.58756010964111,0.45736374379965117,0.4897185108308979,0.040573121123083755,0.019708686607673442,0.5296413344298418,0.6446433251733965,0.11661341753593213,0.8356999667249347,0.5885919772692618,0.027643831588785694,0.5379002206277679,0.6633435671407959,0.2267489556381741,0.412726876538531,0.23022179738420934,0.1323206958337153,0.15410926764103183,0.19875445417325222,0.8333954516770867,0.8900566121488256,0.9339078485921366,0.8807768736528787,0.6120817153121596,0.9266669233070466,0.2737241050228403,0.9355795703207992,0.2317363707182607,0.936837462253221,0.19411228839950123,0.08818945705974857,0.5364276983045966,0.05534143090029997,0.9592286950695004,0.0443271616107328,0.2130448558662128,0.37174682457737096,0.09319215942300763,0.6402318380271734,0.4226590208971903,0.45460364022775135,0.2716865976668308,0.589131194743362,0.8998285778931039,0.35598089344386374,0.6364928630861659,0.30950614387090347,0.7128847072400721,0.3332869615812799,0.4686817623844579,0.4780179676431511,0.49230986259663956,0.2136404751331611,0.7892906082415139,0.11071036798269618,0.4544990817141088,0.013540317993561368,0.5988192476240914,0.04089612526002684,0.4133812165365254,0.8836299784576347,0.7741536264227422,0.32152050283233735,0.03808756062350516,0.7014014609159727,0.6351255294154883,0.4478489625010944,0.4850508776099762,0.38762181144801455,0.143213061497522,0.5226988505401586,0.11370095447246409,0.6029979002383397,0.6851996370477214,0.1162308645261354,0.6657716151533704,0.9966259340653139,0.6883662815272229,0.3071268804012346,0.8817889835255008,0.11570032813674613,0.29316314644283725,0.36010873812853994,0.9551535713113792,0.9985747742869868,0.2672341239214713,0.24315053501900352,0.5082654917132121,0.6912572443631064,0.4802824594521562,0.4281558578914886,0.6797331242899176,0.28469497499604857,0.14140460026166934,0.1872272804553492,0.8602576069040799,0.7604394189468633,0.41648412202789986,0.2070347102707193,0.7870866762825908,0.5008934122821231,0.3338741195565078,0.2734720200540104,0.13081588690632495,0.08962672276277361,0.0060296310254213115,0.8427209559711932,0.9492962533601637,0.8331291610347356,0.562611068725714,0.6139329487525893,0.9073549385168294,0.6907411542316557,0.8118664730761156,0.6199340679192091,0.4222417164085357,0.5570910555362192,0.13910457543106514,0.5817576916590177,0.9271043958808949,0.7474628420916478,0.2561731429809,0.8440764131231799,0.2932047632326602,0.9333165278504689,0.4866475879500607,0.3545745353841878,0.6331707855113256,0.8557213711372682,0.3755478976020312,0.8866380871214492,0.9746599771793791,0.2025023888251778,0.33963914059598266,0.6040682790364831,0.9491960102027739,0.6272860508457825,0.26340342360540214,0.9506807404761195,0.10727461207666122,0.9540261160927013,0.7764505341335068,0.5565054102503548,0.7454923346720415,0.8597470432138095,0.3192631245983577,0.7558949348190955,0.5440319904528854,0.17791085285045316,0.016971425426737685,0.8590770895821187,0.48269887461728445,0.4116271478677269,0.23426522735557787,0.8686078250277063,0.5035392517499657,0.565911389140362,0.3602409716707764,0.11301074178745374,0.6434804072733058,0.8142347820767719,0.45279890854081195,0.19084475507296306,0.58771953977244,0.13910597805052638,0.3503585962352638,0.6922514192305318,0.09714745746289388,0.8443637423945942,0.9466959474017724,0.5390786679700508,0.700114225679888,0.16646660289224435,0.6338464300582636,0.9121995425591707,0.09209291466215896,0.24793300099276605,0.5436923074386186,0.5879372520516208,0.4398394139520516,0.1393800702431065,0.44414629159384234,0.7725535146232967,0.5330239813804034,0.35614861111204377,0.49605461100471393,0.6728982373951777,0.13103591800172043,0.5743451261670325,0.7980374688335008,0.05825273025754951,0.025793047737739272,0.6723086407575234,0.7133010756743343,0.3264357191943257,0.7358048378199167,0.4095786246883568,0.9090690893339939,0.13235791097497596,0.28959053916340305,0.7107620081895052,0.40079338536632947,0.24660383738068226,0.7843651256826198,0.19132849904388172,0.4846016784293261,0.40035297467251985,0.5359614581240872,0.5967382749275607,0.8169719878679342,0.8276146176366168,0.7124913579869667,0.8027513810565663,0.9255460012143488,0.17768681610592463,0.40830512765882543,0.8868763032082603,0.27465150401114946,0.5203616770229624,0.8369032243622911,0.3301922118088627,0.8476710746627549,0.6178124042620567,0.22819524875956065,0.3486541495708695,0.7394267277260225,0.9240999872887086,0.7566583829151698,0.22724593758414868,0.38495094751999637,0.6099708995458204,0.38882477447864217,0.6527703221481984,0.6303333518082264,0.591992088496641,0.7781797581840184,0.5109976317809741,0.3834521955863781,0.8752833609100004,0.7788313950191101,0.06661811028264286,0.8279057767310138,0.5387891880567943,0.7977345688992485,0.4976544653475874,0.7482653915941948,0.916069169725179,0.39165642173317794,0.22562261997195865,0.824633903151869,0.23838172474588148,0.14299390210163754,0.4060843579509452,0.6074069552437282,0.41857964900979994,0.30011316971515933,0.040162107769432076,0.6875135854878032,0.16658444000843242,0.7308718123231572,0.6347488405909014,0.9598419577172548,0.08898642705968651,0.37915498131746883,0.7187006176307444,0.9140425022093668,0.7074286840367496,0.6254527554724378,0.8923811844810912,0.17951553016999888,0.378058989211586,0.5635680991565175,0.08373318886600067,0.7262544914881672,0.43602094198853325,0.0059169344448224415,0.7920553259768923,0.1816133136048188,0.08926212153805002,0.7505890096096091,0.9389611662496766,0.10550429461562705,0.47525424648415937,0.8688902541902142,0.8457340594750451,0.9235470194835291,0.38580159385386836,0.4292628631412465,0.5854876762119674,0.14353375467162777,0.7977416425327544,0.7161336391338284,0.13634054062146306,0.8947827148432326,0.5340043675454142,0.4202364711969757,0.3105700819904218,0.2425929429293452,0.9388602384170401,0.17844294194537658,0.061513045003192435,0.7069600331205539,0.2913168602824845,0.588192273388495,0.040700773903685095,0.48755870766844644,0.29587390880197484,0.3121607864870761,0.6654786681964275,0.5955355752966911,0.7485355797868289,0.9319601996137528,0.29677571114302803,0.452685575260934,0.7083756389370888,0.4025896348381828,0.2330843107247822,0.7632154970657294,0.4660094452654535,0.5729220092856364,0.15621337392538237,0.8731427690082663,0.4496569090535607,0.6978463075295711,0.31424827211647943,0.48226995469831957,0.010833878889683346,0.79411240805853,0.7600852302719532,0.9628039596797486,0.3065553282300564,0.00310012362625689,0.5590921349681632,0.5537430058063757,0.362706585890723,0.43530449197028065,0.6978310947585913,0.26337952852494695,0.8107031873726007,0.35362801998043847,0.9257468552979708,0.8236768714391798,0.19786108123454227,0.09066855502356241,0.22168029719809823,0.730870077816441,0.16654904627821665,0.4065997571126114,0.3354026188186323,0.49259212790732687,0.26248793720887065,0.5700609296462853,0.8517336886556763,0.11078582616029187,0.48246455984499903,0.8623808301535341,0.606915382713259,0.7243144507351296,0.26173941846933557,0.8723976785422018,0.24143584147326047,0.6746955281767044,0.8393531179353158,0.32953395612170255,0.6044079178131319,0.7481663755054331,0.11453178172138267,0.18763203307542287,0.4729530011722385,0.9888912086535728,0.8032789573034722,0.24373937346153385,0.2304506046797523,0.24286678187489286,0.06837786237065058,0.17845356557978342,0.30310600330688553,0.8779630069428201,0.5529258404106883,0.6662138438594545,0.8912031367172082,0.6061814825653526,0.6996279713088716,0.8970088723984844,0.12016109026858102,0.0989205698462674,0.9175468596884091,0.3928484428457524,0.20890089639430176,0.3697854351509784,0.9062143024882057,0.8972571328723913,0.6570432980023633,0.5322744575139249,0.40942708059623323,0.4120648549550221,0.22920766240071955,0.4256117637933341,0.6055792334348711,0.38766295686630103,0.7407089796137869,0.06602058478925166,0.7420710402105579,0.07547733108511334,0.3990573624667888,0.4760279714358142,0.8108811156736222,0.5481068166445684,0.8063182298494882,0.26414597111095084,0.7047977715704654,0.26428088002490413,0.4919890808974825,0.8305956027268209,0.04608185711418211,0.27938730684478275,0.016304645871548562,0.22676661284652289,0.12224247884063288,0.42558005904159457,0.48017323451679217,0.4622250051654171,0.02882326200814944,0.28575184572336965,0.9704318894835084,0.23802620808078667,0.0593522543506968,0.3191917400843751,0.1240468329171911,0.677909599610478,0.36849319009203585,0.41318101906046545,0.31223934539113274,0.04912736021536945,0.6522545691249264,0.8277342885549674,0.9385365868794688,0.29362862210590457,0.9149326464713201,0.6533712178083901,0.39069976354687985,0.5653282795074919,0.4577382214257819,0.18167040268337242,0.3748012801469347,0.519803876237211,0.5120087088731253,0.0470241354980232,0.48039066320121115,0.4318026418213965,0.16958203869910093,0.6122981129967044,0.9991089441244481,0.37435319865408334,0.925015723735686,0.9434601378717091,0.5067479149174383],"alpha":[2.5759815593814213,3.18776856927723,2.4891563198997697,1.4150610924990032,3.4991405144377348,4.974531603768562,1.4212013048512024,5.057778917743153,3.4267873327669225,1.8436920761745674,0.7856271703272717,2.280815165242684,1.2424911570122705,1.1307657747478457,0.7850932088101876,1.7932636888974647,4.747364750708571,3.0839633472499326,1.5967227190605433,4.4903470504614855,3.2028801991762497,3.923322386620348,4.027769189623935,2.8217179929470397,2.1848241709726586,1.5672698944996517,4.528579403283651,0.8157495912218622,5.172539348998365,1.7369672960283218,3.6243057912096805,1.4032807094342998,3.0531656747727647,1.3006158575347677,3.096706260727263,3.6559819299591183,0.7462537074300475,4.231313985835204,0.5671307393268328,3.984321269216842,1.6680445326034583,2.5322924145655485,3.8823698734084506,0.6782112108727707,1.5563335245341443,1.2697010608076837,3.4158255556054034,0.643737023402486,2.403906435282839,4.603470864285331,2.0491628558318715,3.013558338340122,3.835880276690192,4.545959052988271,1.5397495721484296,2.3115428994790195,2.9628288110083627,4.0981909066566375,1.3027815428128793,3.029093065214678,2.5356999346167157,2.785672790037201,0.8893554555885759,1.4208984901470085,2.8679424960958926,2.9834836489464376,2.335723137261861,0.648742881785886,4.580226147860374,2.496705348852348,0.6125614043885426,3.5536633721261097,2.928030906696951,1.4006969958657751,2.487841854661989,2.365155744929231,0.772896418473727,0.5862624284906721,3.279351789211013,3.246270910688672,2.6889504482298126,3.1040801262991753,2.9667075533663585,3.1729437424103453,4.110353652117766,4.478615989253371,1.0113174981753665,0.9907612705719966,4.209958835513399,0.8270265364061209,1.9324038571414277,2.717961047788518,5.4081675791765065,2.807069890343152,3.7929856465031726,0.6298414732752706,0.9826532084043411,2.4070360057233824,4.744799987068868,4.1210632360671,3.1918625316604983,1.39833266402364,1.7734907239774962,5.162857359366601,1.7059203213771519,2.1510024842438265,3.519884796014738,4.769487905688219,0.5854809684776625,2.3402382749020365,2.957720133368645,4.262442998836852,4.606692930705686,0.5282091357950687,2.3465599737266727,2.0987360806354465,1.7402690593602317,2.464172149987294,0.727943250748952,3.6655626305371487,4.829094561089772,5.048575634460679,2.429605211034724,4.419871604542778,5.384864279042755,5.20225200618347,3.324259759523544,5.2069237023776935,0.6171561461076043,5.49215068899863,2.808870880661998,1.0232146834026308,3.6899170942587265,3.2954561792628208,3.3186325638857768,1.581929300179526,1.7046263241099642,1.0503821906101125,2.622566879663994,4.124793163251478,1.9105508406850653,2.4487479757699075,0.7221185728712427,2.177492404077914,1.300395287285468,0.7064814966758266,5.037023180484376,4.592654705332289,4.902583149867516,4.670549450014615,0.8234256162124338,3.6578066308560357,4.954297320202706,2.264330960477385,4.235542895964868,2.1255281636915555,4.438624350855556,0.6732430288097742,2.6007100748181706,4.664029784461812,0.9565791270265862,5.2014601223432155,3.905131868797129,2.34414804108164,4.380230911515204,5.141775762225208,5.45166411277973,1.8411678851492237,5.222212126745299,2.9984477705579566,2.3493079033409874,4.92768857805337,2.657562864980622,0.5461063128849748,4.948731876742247,4.061209048631978,1.3008500853963219,1.4991733275597403,1.661861497280049,2.4057437260056496,3.196124205583329,4.804372202554026,3.1571785621208988,0.7812155197508814,1.7980257070456136,3.482282537577153,3.658849169266331,4.033168765562596,4.409192004874626,3.949204443842696,1.2783188656895457,3.6887061614880308,3.3472811481079887,1.3488697417734792,2.14513941514864,3.132982119706627,1.031520717222267,4.785193250098802,1.1230973445708343,2.9276231061489284,1.4594084757850974,5.468755662431478,3.667975203345982,3.464460727697263,3.5746600977065963,4.400177360297127,0.5951947970665812,4.73175639419926,0.9137324014615396,1.0862401880513612,4.849458714969644,2.209571236762005,1.8499222979180403,5.470200017300364,2.874517179485319,4.692925880851314,3.960548123927465,4.752074264168919,5.33525194004387,5.281595382274305,2.1389109582112726,1.6797169488114057,1.5244174953212204,2.1624631663927127,4.600196772516674,1.979255443669062,1.9018218908246804,2.715646796484909,1.7500479601079528,3.126797983081268,5.146104675742634,4.336825899621969,4.251293368436568,0.8843790846881001,1.8922444474049567,2.7042644867710646,1.4490096732505549,4.933922551771131,3.4677213830336124,2.4337237439255874,4.780507079913484,1.6754885409733384,4.412943115736969,0.7538198058596441,5.257309684592014,4.980782211823073,3.5753687293049135,0.759626780946216,1.0165776850222747,0.7119070987242944,4.585941169399482,2.8409603074992136,1.5114989176682554,4.911778836372066,1.4506314306741548,2.916369717045208,1.7495735028781463,5.134247219761276,5.16779126845426,0.7618676769267414,1.5951226100267097,0.6947809674375867,3.4529420445772003,2.016731873576369,4.897288699952915,4.499254927761108,1.2270217148289637,2.2042569786308572,4.4749024668060144,2.516482343897211,1.0189476678434934,1.56053727652632,3.1160363321617117,2.6473985125930106,4.079259140541789,5.2816331274389565,0.7594943690461637,0.7124905294633062,3.6360814888571253,3.452811264834585,4.256158099058352,1.3527781763002573,1.28485653697555,5.28002724989015,2.2039844665632495,5.473637678527522,5.2180559300689815,1.8811715168570238,3.160254982947089,2.0360006693660653,2.996680452648381,4.457369318075377,3.7417792470416695,0.8529068768424071,4.73275615871688,3.9262427077384503,5.158798336423078,2.1517371439315234,5.1079190494025335,4.905319911752333,2.1964747680827816,4.749754209953216,5.115538947063172,3.8689424994722095,3.7169633465757386,2.2619301820921174,1.2581525117134735,4.43203062403181,0.8938013214201008,0.9886089369795477,3.4730263843328597,1.462280445899476,1.8425331882927027,1.2606240060116298,3.507558574189688,5.00451573746947,5.109595880318976,0.958962390810056,3.5196829929511315,1.6568949904154113,5.124571980451218,4.662561304142194,2.1996863085527494,0.8216657470490558,1.223436704030096,0.7651066912294597,1.2555055949200407,2.6962817649120323,4.309839835866631,4.117686653402928,2.16875731048472,1.779563712396734,2.9787035074685706,0.8998137523358181,3.1074975334667294,2.8316692714559863,4.226715622905696,4.925063233932106,0.6348068442178063,2.80150435499091,3.7026742789694285,4.577337400360708,2.259567952679728,4.228085281818468,1.3514580843604247,5.156529106770179,0.9195301619669567,2.227261923945189,3.445610200457401,2.875609705039312,2.298912704506609,1.1646639288812672,5.187918159398359,2.6401660134269673,4.999609759003826,3.8057159993241934,0.7503992709151792,5.216017443927165,4.054554750098628,0.6490177287870442,1.5566255549773085,2.7275652077718124,1.0947699267530437,1.4779473357497426,1.0779958934555052,0.5676960799554838,3.085151555801434,4.992857086332339,4.713442033929039,1.3303064740930721,5.045750072055763,4.656851821827596,3.3247610777369427,4.53990718219548,0.6143463600376982,3.2743215975982114,2.9588179360321085,2.042813026045946,1.5321176232715956,4.146975012054764,4.071846616697587,4.605132144100248,3.8555440954283,3.2079891862113605,2.300747209574569,2.119787524636439,5.062419595355793,1.826768381040917,1.1566450277039593,2.9699026998611338,2.1855067239446275,3.0733755062577557,2.7057960307789735,4.3118286058710265,1.61896070316568,3.887763702962263,3.9833861291492862,1.6816973324844005,4.224655621720353,4.5641458755148,4.224371679331782,5.219558132360406,2.4228316410494783,1.2419590423118183,2.534928313804749,3.3994589153339403,4.973113476474277,4.058547364706063,4.902181291293173,1.0129019668305146,3.457239403182732,3.0331080253286267,0.9948068027201047,1.2932798202965707,2.1644784178116137,2.2940510215683902,2.3475081409811356,1.487780224226007,0.6068248820659101,2.970210961934777,2.532933875667026,5.222332301672349,1.7659536885601785,1.229714977253898,1.1195196374531031,2.373345610355357,1.6244740936757203,3.446362459646,2.9357365537231312,3.2958013985558816,1.4364908104380454,2.7250889642883,4.772440384684089,2.3525600189843363,3.7830656609614053,4.021774566315318,4.32058374549613,2.4592917422171316,5.127842482292939,2.839922784442874,0.6716003286732992,4.096575127756935,0.8091655653420677,1.7594336074353008,4.799391851665904,1.6823449302953848,0.8741860098779826,3.528563085666182,4.949179423785964,4.448736011364549,0.6731980620993734,5.109317279079193,0.5130181842158972,1.346910327758995,1.458989325295132,3.1567361184007057,1.0229125261599765,3.1834595847262643,1.4725693669559998,4.446807651634919,3.9018536423222234,5.080834244175285,1.586504727172909,3.405823760517945,1.2279244208201485,4.14870431674885,1.3749548796298112,1.475779685919037,5.49070518358783,0.5384681288757371,0.6233553704397594,4.002762143030182,2.0611643971143243,2.3285721264332566,2.6592709351988146,1.0200559871710486,3.5635106027303536,4.21137870592047,1.66825944059481,4.794898729304396,1.1574152447698767,3.526943431945981,3.3160644590114217,2.745644610796985,0.8643589849836881,4.028280754039363,5.353531938299435,1.0450775283453013,0.8464263370970233,3.1760336834521814,2.4929259985475545,2.4255391080239006,1.1327510741025921,5.2337381185011385,4.979769597298868,5.25482278606547,5.105303997868669,4.217587382852198,1.587036957520786,3.6402195459675095,1.706581989916783,2.4002091798364673,5.275939378993059,1.5463354745260456,2.1078880878126145,1.3708800639448455,3.6369482687342076,2.4316160281233454,0.9463176909763755,1.012916702984636,0.6151979190668284,0.5057590751064085,5.2466707012164395,1.4835283483678758,3.7499607820614074,4.087589346731378,3.459515765072541,3.0961206219150834,5.018550606923463,2.227721234273758,3.83604281275641,2.526658406176125,4.470162883125623,4.17102986073247,2.5793980983442273,2.238884766859734,2.564149591313872,5.163024465771468,4.060420389653091,4.243680062116832,2.9925085680641104,5.14800845317681,4.34749551909419,5.45072175862823,0.7476759002612838,2.0329961470842446,1.5692722742033334,2.321467980582232,5.274138702949276,3.0742639109466694,3.8214514648365823,3.1214848085856515,0.9152484085724473,5.308296797181496,3.288131492389777,4.1509977327813825,5.384482506269982,3.9092312642940747,3.164608277748272,3.147421139814712,2.328504396829641,2.734548126526322,5.204433603145533,2.9368894063975133,4.574648834123412,0.5910762905766671,1.040824419641922,0.5763697106004269,5.190273349370754,3.21610347754506,2.5736760299127157,2.8859307800397067,3.394617030838802,1.8252727263024966,2.811072130078562,4.314430811377504,2.9537097121193874,5.2096963594908745,3.98031475836951,4.0458671828338115,3.8047763596648574,3.9232984203656125,2.144805861354027,5.053130311957157,1.3695325666984044,3.4789938710547217,4.251711519629001,5.32869069945634,2.487854435391111,4.570484197193526,3.245282186963276,4.607661579415084,4.054394792854577,4.148791232576674,1.7431801495755823,3.3233433281417764,2.966229687977032,3.7184012467331646,2.6402255983533958,5.385996699938429,1.6371467577318728,3.4988380372679133,4.121303501676256,4.714864352675008,1.2730157887278521,4.677800504824029,2.9269005986838286,0.9263965702746075,3.3777044889609495,4.9986977447791485,3.8361704827195835,2.7004640492113845,4.725581889621406,2.8453002071291627,5.454400877248841,2.623608535665019,1.7952743919049092,1.0216152695274967,5.367284674708862,2.868914249320535,3.751199810872979,4.165740128054011,4.190861242661292,4.448419572552366,1.2258590030836531,4.815772773611812,1.435507319480482,2.0636131984932797,3.9133105379598536,4.752083393448127,2.3734688512431386,5.053741476206951,0.8470281675427811,1.062032396306128,5.019256303037473,3.823918770869103,4.480145767611452,5.272209277595332,4.835017135428389,4.121028998752363,2.655043507816805,5.341914849425953,5.436179753952618,2.8049321709526494,3.420941783134526,0.8858822218881652,5.16570461327685,0.9814916232587163,0.9391401266803274,1.2616574506738938,4.935311452419216,3.894913916054204,2.4664997849779695,3.696135412095804,4.078503989150682,4.476473731492472,3.576491880095026,4.420650851511752,0.742188184917711,4.463296446467694,2.7375450749155124,4.8341674784902136,1.0452261598570118,2.9657420633248894,2.3394811357834797,1.2304752000245807,3.5516403342178258,1.1374291190398818,1.3052859214493524,4.926347564957603,1.9631757173654603,1.0251423703916025,1.772872726284584,1.8920891949840286,4.798218516775822,2.4011277021951756,5.089497055246739,1.1626061588911392,2.382101696200392,4.063980441031772,0.8563188111992599,2.778778443393885,5.1211323375618605,0.985773386690822,1.2949166674223396,4.962427394664315,4.626496388457584,0.725867174297969,5.094031693124104,2.98413778745996,1.9536052160295643,0.7885301308691013,2.592035523873923,3.870547167243859,2.252394806156129,3.7155616989795464,5.217716603986983,2.1386431613602817,2.5798288474078186,1.577587042732635,1.9753604467243164,5.2775301528763015,1.7088820406689937,4.695337369349112,2.0011868811053053,4.719638859163442,2.5878577190365064,1.6127004262236184,1.5897924886522903,1.6936105465486035,1.5505540647962053,2.6389498685911223,3.2194433367374113,3.3489753403258424,1.698676917093675,3.0896478632604927,2.5338313698612422,4.726776484571241,2.0404117074536616,1.3754257727323156,5.057199802584739,2.446787625009201,3.2138990164811085,3.9922309861791034,1.296535054601173,5.291696330803008,1.7346308740508427,5.227685582892331,3.824654938194053,3.491929503856988,0.7717943466025685,1.6750409959911245,2.2272212581528756,1.3322723555695843,2.6576864191846212,1.3738901660726872,1.2742216563631843,4.462030265890707,4.369459909550132,2.9093482831247224,5.150170997055118,3.060259751852927,0.7156361333002677,3.9590215818561845,2.474454445368967,0.7572601129040342,1.6988290151667058,1.4587097199678563,2.889342463892292,1.0207755088700572,4.604289828973766,1.9455613552985613,1.277779072397532,3.613003938226181,2.6767215983950647,2.9124202401805404,5.088156934180572,1.8392495130798843,4.338627587409308,5.350551766003324,0.5774939338513303,0.6492822078061318,5.453216895900082,3.997227600753651,3.378853599409145,5.381269246113665,1.2362647709315662,1.0896268841177676,4.743760636943858,2.53401785187163,0.5254194070192892,2.294401525319632,4.246021327909174,1.8338951030129487,3.2514996923895896,4.372520362332059,4.596968535800394,4.667977429175641,1.5762561572481633,0.867927515564801,1.5586768853734174,1.2316523827449017,3.6679658252915353,3.312219282448715,1.2514262472110156,1.800747941970396,1.9620245579413893,2.358374679906911,2.8489322833391935,4.398613976762139,4.352036881446004,0.9070121572684282,1.6839973494458866,3.5064169728351233,3.2704120629219786,0.5650770016465119,1.7441412761092019,2.271123051381948,3.603448739211272,4.157549915512197,2.6872355285189604,4.247771632279927,2.2875159107366323,1.551038233928,4.913311815909581,4.894769196273609,1.632282747532587,5.06808340059709,0.8780327813125561,4.252946394077724,0.6629679536930617,3.9056165212035405,1.2985088575919008,1.6678568959835178,2.024521925378979,5.099511396917749,1.757531580216834,1.408723784619041,4.402703917609319,5.476459200518693,0.7707816399184404,2.616386069227436,3.956973921218852,2.8603389516743665,2.5315650853634404,4.2434063466647105,3.7469233534702013,2.916441073140007,4.763104218443024,3.6683766929679535,0.5225338185074586,3.3842551649997397,5.040676760808502,3.476535492791559,4.4716242143675515,4.662729098745929,2.381496007957165,4.837176930612028,1.7483003774558277,5.1217191530034265,0.6578342061554561,4.053198765375829,0.6245693759858812,0.8683928430881633,3.775759040883096,2.158993888595698,5.3362690322875,3.8074015617207775,2.728011593341713,0.6019637996444189,3.5844425660947516,3.1173151257348373,5.219169685730545,2.1235942503979683,3.7022631198490648,4.672719804197635,5.235660264624228,2.6418167894691464,0.8795484003840044,5.340857705837677,2.876248010234623,5.384822639669362,5.424304186535934,1.1907255327038575,3.8436891167425364,0.9168442572058062,1.3368138103284672,4.362999560948603,0.7274953683533043,4.254659328421805,1.0919479342118272,4.019963869579369,4.211074042289733,3.9216479329293925,1.5680312840769657,2.580424994001935,1.7469097319739353,1.1257596364063789,5.259244873301684,2.779095442693715,3.082365281805572,1.5418753462676955,0.6797412157996181,4.439642882599199,4.03110161679929,4.2425236470998176,2.8116961053074925,4.7118389121681234,2.2946320635064685,2.7775810790608384,3.9222015745149514,2.3457076073870837,5.171635698824622,3.9700775148965524,5.47950174009404,3.521755716606175,2.162491078852587,3.3406055989796175,0.9573445464857928,1.0549386759509334,1.6182821967188317,2.559861599401931,4.578210892578723,3.0532491389028777,3.4381446635937847,1.1814761231646516,3.8696652018658426,4.540553123948589,4.2699852217652445,0.7217406145042855,4.034895476508545,2.8187534723450867,1.1417190344613588,3.1616912291867156,4.033784813354265,2.8531969003907505,3.2063973388155347,1.4885786484630477,2.064188543946188,4.57007139853087,1.6066864264723997,0.8786975044334728,4.311255639673757,3.980094900436901,4.118092407010357,3.922960425524896,2.573995164375233,0.6249317896580693,4.020255592677688,4.923616715105949,3.8093907603028727,3.9402211725185374,2.196499771343344,3.211922778164093,3.6566223446497546,4.420856052942852,0.7583792156258125,2.856749617463858,3.1339256467654515,3.6563355771848203,1.1309479639116533,2.717530571784454,3.897529712386632,1.8051710247126347,2.1854606415736195,4.876706454339682,5.060393316499514,2.1205509649014758,4.213721646065473,3.5530455348711207,4.612079704621793,4.88270073091136,5.313187579877879,1.820585211711931,2.93274199191036,5.295411700808615,0.7514568107257875,4.867165655718369,0.7918069756387853,3.4190297520626904,4.332478746170068,4.351993429613555,1.5619075395463284,4.597513041169191,3.05617642182575,2.2902046640122347,5.4385488908058806,3.2979893428302036,1.1007667941698414,3.521728422515876,1.7245063361991848,3.240273569458979,4.4473900033764595,5.214648554123599,3.5714968736927983,3.2495854368864703,1.5869074755804364,1.8394308089254914,0.8887220491495558,0.6558675256998088,2.7406395679488833,1.9789904630302664,4.203576635147924,2.3275592748681326,4.623479842843244,1.4727071208212557,3.6090322090035905,1.6372456440000167,4.6781285055861055,2.780896071113898,4.032753637440871,4.333535919967771,3.3667520787099043,1.1798550214369303,3.4246001168375146,4.68910767745527,2.4330442725801635,2.3591135266994407,1.0754930540574732,0.5542701795317491,1.9848550122955728,1.5617270893771158,2.230198941187495,1.5212369880530558,3.7310067343028592,5.248099675043681,5.284272691499362,4.774004616837386,1.5895032293780782,2.0227780512509623],"beta":[5.4589601341826555,2.7559689068337,1.116981130895926,4.408437276359957,3.0589444221872237,1.550585803624624,5.22697274013284,2.9779550569673905,0.671987125835873,3.636652611517862,3.7607990844312242,4.460499643983624,1.0477037658974127,1.577096392182847,2.91146046748327,2.771850835563396,1.827529728576513,1.959301405401796,0.944262698150869,3.2156358002878065,4.338923781228795,1.0701037187075386,4.560151629061326,4.133587008471331,1.6049463168628193,5.304904732395694,3.669746056319193,1.4878119240965666,3.3457856211484667,4.45319426723424,4.091740681354267,2.387641616989834,0.9718952239753715,3.1333081684423156,3.175756863369083,3.088386858500511,2.447174332555405,4.277554720562364,4.4053430528965905,3.128571689201741,0.8040595700470927,5.209155243181012,3.36910261128532,3.2506768311619876,5.463551685377822,1.9412225672068202,3.9782076021467514,1.9039270492323208,1.1934386046713847,4.28141045307828,0.9875999571036582,5.066367977827314,4.779477694370732,4.848704769372301,3.207364892424539,5.270520703120259,0.6368763900244818,3.6915159899434076,2.7062366189494,4.383531739157354,3.9864074099732054,2.79629874641454,2.0222916260574424,4.519804572913403,2.636189298468558,3.2015365168227925,2.9048105538420494,3.8714171450668373,2.934055607612291,0.5446086057732598,1.3014293570826392,0.5203330309372851,4.604102240268085,3.002094104864363,2.9647397828427198,4.161978406203853,4.079784897220034,3.22837908829229,1.266818587975262,5.0769806360886545,4.475499883803906,3.9572640454265793,3.75186244436125,4.530978589136066,2.4467896495266546,2.9616419239400846,4.379312787277583,1.7116381680120503,1.8212536163547928,1.01255011831836,4.677127807103332,3.6401851380998655,2.0684760486245306,3.7833906090415055,4.636993390295652,2.8365396550540805,4.93586962163518,2.193712328662522,5.054460264283955,4.313158564120166,4.7127692646534785,0.8885197812467378,1.1899946346534613,2.5439018869570136,1.537921008702667,2.838850596552746,2.294906913982067,2.3728890729630168,5.395288941060202,1.4732661130264353,2.6974204386531704,1.3727034506625198,0.9190403767183981,4.530096010621512,4.401730247085736,2.5744373734265853,4.86882343693977,1.74683650748728,2.6746507640491632,3.3241722288092963,4.075530221843613,1.9206817355167134,1.6683802411274513,5.203616210573426,0.9628401172874608,1.3216218963291966,4.524614066991784,3.501801664204561,1.4823974360374221,5.359041169306315,1.0799495224158846,3.618306684353671,0.9281029662773983,3.6104106200290293,5.109145152710592,1.8390744854725123,0.7990428363836728,4.7048481771231385,1.7643686512676833,1.1983940028673363,2.304664877959239,2.5636004723189307,2.3606331445126925,4.642042926916657,1.2869241458979481,1.1574492853714664,5.453497621440373,4.7883688389681565,2.438874882869838,3.7582402928488383,3.9581920238999233,1.502749289713497,2.6406656125547254,4.791119290226031,5.492165141081677,1.8314678520511136,0.8264925885527525,3.57734536258054,0.6543554916808676,1.3627967309882523,2.1610661119925183,1.0994524983378058,3.7223003712871967,3.9097922334237363,2.1997828274259175,2.0280090078172437,3.9958026831797904,2.1388634945844824,1.8194993943969213,2.7967019070462746,1.9080002858764251,3.8643982227697844,2.053782028740277,0.701762266447798,4.633825652533618,1.0619768919856127,3.4317972392719884,3.688485531475921,1.6544823979026686,1.3901159227999955,4.652937387548822,4.773086870050662,4.277861481107626,1.5859502391833464,2.0237937703226034,2.7521454347771153,2.038028987021364,2.651090101859009,0.5740356076125548,4.327488231756436,1.5382538779706196,4.681358890588653,4.675789484160376,1.5050119102915844,2.8647097490427305,2.7871829557416854,3.301782198909363,2.2332271552428278,5.110650678792743,1.03984563578517,2.452624666250239,4.313514549064479,2.391583652177491,4.999663277461524,3.931037164006349,4.315524900713108,0.5983902299070338,5.248619919284653,2.3808356448836863,1.68044893982449,3.059674901877779,2.9243258018217864,2.137952134358219,2.739883759342743,1.7517090919219065,1.7953647619951814,1.4881443922058764,2.980765021815008,3.0902657421390254,2.6032831449392564,4.806051850646545,4.420896350098687,2.0779312049377925,3.2742088364895223,1.0236511027001856,5.331509226007746,1.4649919647946348,1.187169370006857,4.170663632716328,4.428508949441553,2.705173378354198,3.814177432145459,4.765720851689872,4.895949586419789,0.5164402919687778,5.429835269411368,3.9799976986171077,2.5148761044773513,5.3731185042287715,5.415347334496603,4.22214212765341,2.5441743046581538,3.0322653903024084,5.3983514769164636,4.637927090504061,3.4512552871680486,2.3787611391552206,1.185832398510835,2.7692592232814968,3.3438149547644462,0.514619764499344,1.6737404326944845,2.087420957815413,1.1812079617371447,1.4032519850847054,4.362415424362872,1.5228645047845697,0.7216070768854583,4.2965243994365565,4.117993854029757,0.8750392245957201,1.65648959116186,2.294868381405447,1.8199979888184987,3.9984878716900116,2.9247822460572044,4.3918447691883475,2.881126737569198,1.0159185192380122,1.133395472376731,2.829544008843113,0.8655596208100709,2.087687257494199,1.6039551987116158,0.6821079119437985,5.320524321627506,5.144775433234848,1.4397305059486272,1.6479716912061844,3.1656249226781616,4.340404427873926,1.8119948281985805,3.268336558228881,2.751037205142982,2.4098732601972013,2.1349156229716613,2.2309465512877953,2.297331178522991,0.8735453442825343,4.246314609691657,4.631421743446545,1.6360277756738535,4.26671264288237,3.933048327469276,4.53969206333905,4.82310505835027,3.5688116294669077,3.7697760479817948,5.157553927986427,2.0215742703852744,4.82891815481776,1.6864821170033124,1.1539194295325013,1.6229498049469722,0.75973765554824,4.872200997873263,4.151560788011183,2.033494693403985,4.577349390615512,0.7667102681054756,3.2621104197041104,0.8606895048481418,4.2189091030823205,2.5044677209761663,2.1202131152073465,5.233778699392673,2.7275929587266785,1.3175555838582682,5.080781276738664,3.9420618121771884,1.546186348269279,2.4226800833340416,4.225548333893981,1.9916500834831736,4.176805679279021,3.5754819079581086,1.9298679866244923,2.0807774771803618,0.9849612737456654,1.0424441862682872,2.631477254880705,3.9158759377983694,3.3111024484034965,3.7695296556567577,3.2364001529469095,2.88767821635724,5.423257255497083,3.1605362027995305,4.7944614211805066,3.433336957281052,4.946146894289572,3.4126804337506904,4.335234462390268,2.8413104717845235,3.7634038214328154,1.3186167318466797,4.847638078968715,0.8523671009535101,3.0587107603698813,3.6603879871965135,4.088523020636764,5.2673582341130185,3.1421623046435103,0.6318099860372668,5.159704825418176,1.3595185079100376,1.357361502637195,0.7845394831637438,4.7034661450056845,2.5588235203391902,3.154805317757143,1.294417723451274,2.311963290711691,2.1072697335368957,1.8721175696953625,3.071682686992385,1.5316925896942497,4.028975853329949,2.052406500441421,2.5781215991476287,1.8818796814067535,0.7480304988069626,1.798327464712502,4.116034369767887,4.410160184444468,3.345148967430864,3.1487325830926856,1.5297090619681981,0.5166903156738485,5.279158087504902,3.030067837522054,3.2958178383240595,1.2272476972447428,4.482730735909796,3.8854807994796685,3.336894373211324,3.106749965623224,1.1624342259800264,0.7727774749889289,4.8737977462101725,3.257616860832827,4.687131236841674,0.8107782190676467,4.106143223977498,4.516104030494214,1.2603808209010843,3.2403202076316955,0.6331637344506564,0.7541884957417229,3.340928199350846,1.776935257687185,0.8559963456053716,5.108293009613891,1.8636056988138034,1.3924278700004618,4.796110094492005,4.745114436172669,0.5276464311738189,3.738746816944226,4.655285024251484,4.453541953752826,1.770073278982916,1.3981088702348616,4.975161754201428,1.3601486049986473,5.306398896589536,0.7144109106161496,5.382463824680401,5.435489333467658,2.591763276405001,1.1485913475457883,3.2222560376423086,2.6139060110649552,3.839485706833085,0.5993488818347232,1.6092618646799135,1.2859510692741756,0.6544523777439843,3.4928414312438356,3.7914764954142717,3.858928920424118,4.858738577784303,2.4694414987082522,5.185093694775249,0.795341918635702,2.9389783202714286,3.8388315724524054,2.719145092060262,3.8535250669161942,1.8224837780214083,5.166036383455526,1.344514682956281,4.6681003959442275,1.0589836155907313,0.8976171409643081,4.874895596954136,3.9235122169574064,4.7990390506240415,3.7464699946006315,2.7735184763152394,2.1228994757019244,3.9725347690617383,4.060465771218579,4.98271267584707,4.555831796438395,0.9675199039469513,4.845917717489748,1.59209715977162,0.5701087224125485,4.907932660133359,5.203826948382196,1.25068069409999,3.905271453963646,2.0701328784285487,2.558816308460485,1.933455013701166,1.003298178746322,2.6920634075260814,4.049501812696839,1.2003974317516837,3.198791689016724,3.0264675372009817,0.9480475874863921,1.4123961467134756,2.5467307058848316,4.610949980160607,4.867348191245034,3.463045002738911,3.83406203101383,0.9675383375529674,1.1653527154832528,1.0295194289320826,1.0321831307366987,0.5675350505375161,1.6621976732155348,0.6797146764104807,3.960016556979272,1.996549222023399,4.052782704062121,3.572346086166786,0.6057286135895352,3.1232499592441174,4.783181634768787,3.018449346408058,5.463269023319005,2.2888177313673004,4.913380240231037,2.205121238229487,2.8686634376327853,4.471306023590352,2.9001209646006663,3.5447591336492073,0.9616522571097381,5.208566957441085,1.634717382213487,5.0311237428078694,1.2049161054758046,3.760472056053082,3.5038971342241307,2.553301534634282,0.6331656943939341,3.924812407846546,3.4480405139476664,2.581400945830568,1.0443584354238318,5.401537028445121,3.307372469895478,5.069607742008602,1.1926794472126274,4.194734255934053,2.70324989075979,1.00597538242515,1.5688879331024859,4.970985410078083,4.271383807253134,3.443613946597811,4.257382984751402,0.7216049226897232,2.8146991378874104,1.097290031867812,0.7113839303401988,2.88420929370644,3.354350141463784,3.608433762249084,3.335454112342244,5.466278635656886,2.0070314012892916,4.7092260088022035,3.035693113389067,3.646539906251892,3.7907274568212763,1.745715745332846,4.268358433835957,1.929478353888456,2.1943108278722194,2.8394527479148697,2.1804191270921764,4.5116381195138135,5.198559893263137,0.6252807574898795,1.4500877707468303,1.9833089644976527,3.0982207020718984,3.9412696928886755,2.6447255298818035,4.670931597358373,2.0844249198326343,4.110127112520227,3.1357939448861503,2.549800359731763,4.337699079106413,2.4744789119969806,3.0714732025709877,3.9472527139639997,0.5113983155686725,4.394140402659769,4.432107721108741,3.6093994178823783,4.371001937229615,1.9726305436149745,1.4786604328550208,2.374033898475518,2.308953996426712,1.1870568549528067,3.4774826200581197,4.461610462382157,2.034618204944367,0.8912088653393575,1.4608551386108566,5.34114636412829,0.9230104835719428,5.20346216728036,4.408434615928632,1.7501260470537845,4.292146021052436,4.2846644314399125,1.3118779456641612,4.4487535564173974,3.9983064621437645,3.042080271738727,1.5058724943767283,1.9709204316783386,4.484243094909661,3.0906354283720723,3.670086697386311,2.290511303711726,1.657787895401876,4.804092369760769,1.463426731272069,2.1747353800944307,2.738479308878303,3.7162603818789326,3.3967978115597717,4.972443955036612,4.194927574867991,4.0333758155646375,4.441147774602421,3.4236494949249856,4.362024307653068,3.356747364415544,1.2219084795519681,4.29026320932215,1.8183990993408667,3.5659044474858375,4.068686493226573,4.687220922146654,3.008314179076787,4.616823386342789,2.3923653342726174,3.9751540510322476,4.473006667552034,1.4823461081808138,4.084021512953677,0.7708605120077512,0.8430023738952115,4.934595498079869,0.8191530025437939,1.2470380268020638,1.499069822280253,2.2080648299762524,2.1496850180738964,2.463719370687035,1.1553459573639808,1.113056631739639,2.025836849922122,5.471441440188884,3.257671618963886,4.31348055082312,4.527093351747732,1.6927975353124687,4.18441183027231,3.759355517431297,4.700810076013193,1.127923798344615,2.9216005504293636,1.4999384938463363,1.4900793305046318,1.6303235412577306,4.7155694140092415,2.7971347616508018,0.9988895257709469,4.2560472924951345,4.8471408061212795,3.0931785118232833,2.282491209639907,4.685320757820213,5.336270498554935,1.3968979405461452,3.19880075091299,0.8290377350547424,1.8579969840315809,4.021599697571093,3.7696374288794097,1.8836177205561195,5.202613990439023,1.542717405234751,4.8605054051241074,3.3126171094562435,2.9685521077517025,1.2648377336846635,5.42111124737669,3.765019347801983,1.731597298322924,1.4129898031412598,2.4310441320605634,4.700012449677237,1.58882099953732,2.5218762922769926,1.046614108088569,2.5180202163016348,5.1233132856333645,1.5270784782410525,1.9734939190533458,5.268167463535727,1.6693093286721512,1.4984576629500475,1.3960134655408614,1.9495237655681736,2.999905764781132,1.434864584403885,4.9050858581048296,3.7493844070365308,1.1403108851787465,4.928798193908873,4.799188326328682,0.596315280980187,1.8220355549173135,0.5536067920689631,1.0364825447886379,3.855107125839031,1.557947591372278,2.1586511047966823,1.5493775258934068,4.069292928818733,2.1318594304204845,5.027104503901445,4.939701640702831,3.9648398822173365,4.122842978904157,4.9553428423216435,4.388310553636021,1.1157870510550056,3.2635939929810736,0.8378412213866444,3.329810410771346,0.8451708331952397,1.6886553410060996,3.741275866783581,1.6881051351796967,4.420518349877908,0.9218994871751447,1.1433922001775738,5.361115695780324,2.3231329445079094,2.9480252548268986,2.327365452203308,0.7848214989013587,0.9052956585196066,4.209332141438775,3.3616614539233263,1.621896556854285,3.359531775144001,1.8513749422458032,4.070905554917687,4.719108922124041,4.092865376599209,2.9186468593445944,3.9598162700894406,2.557018661003722,2.45500119689068,2.668086996033347,4.945649330801799,3.2539110221537078,3.941318791994128,1.7355375346710329,1.769048935588144,0.611362622194684,4.160201426071843,3.269305825030278,1.0209894012981198,1.62047355210396,2.1605374006586295,3.6132930843678217,1.567749887353372,1.8704794270487208,2.079476670291904,2.59588080115772,2.189597026128087,4.686326813886332,2.165444871359475,3.5924726214265466,0.6737158092550473,5.338001357908586,1.588890304240282,3.2071030692423625,5.097612224163581,4.40088298888087,1.3105647688750994,0.8180074655553133,3.158392008926158,5.023432552359466,1.6921789184229947,3.296970209235201,5.038837285631185,2.3178512022656426,1.918831343373752,1.7235773612673406,5.25834937214883,1.5902142404375348,2.7849075136569628,1.2560917120556288,0.9777010023823913,2.786561308002094,5.180254763862012,2.9033492099080145,3.0542200133601023,1.6398421311593643,2.765525337755977,4.857465260478575,2.4898443134072394,5.250135116771875,4.614108841036201,1.797465199592982,3.300188734151803,4.006555852307213,3.0528256557661413,3.8603576505423405,2.4944782756518586,2.088746670845734,2.5853086502992344,1.8013461074805355,2.7248937784727776,5.147593420812308,5.4942775901081395,3.8180020647858677,0.9172445904399846,1.1891689448463922,2.745174806794956,0.7842458319300711,1.9443702752557774,1.962337406419798,3.2466842894304087,2.7042564829980704,2.535907852182585,1.9071454501855574,4.596465077123468,1.427552140931357,4.071017660843868,4.428594180752192,4.711670552193667,3.2791384253702356,5.172993846838756,3.898815279503232,2.3643630253937093,0.9095994809965742,3.722231826425098,2.785387378665422,2.857760465953471,3.0513413899563986,4.792094786273677,2.16322591781614,5.442809613378637,1.18695339161443,5.2905708371488736,4.615955840881691,4.535738653869609,3.6061962318570435,2.4225354208136114,4.999486152138416,4.292021437749859,4.141183744373731,5.073036177107447,1.8736254491436795,0.8672783047157091,0.8445406379820626,1.1646049800573195,5.491093157293084,2.9448368712980217,5.455805239981622,1.9929680826087233,3.7289559129168923,3.8516819838535215,0.547196775295858,4.463596455562546,3.978271054511607,0.9300561151209592,0.7382383684767375,1.4748936720675783,5.068117755076396,4.979283237302941,4.328304877313049,5.474866434659147,3.6054435025315286,1.8606811451945588,2.004862658566408,2.87431790962197,3.1895628803376637,3.2820423206570166,3.5023012274220475,3.0996140084922503,2.7165897065422047,2.2705736653064905,3.2467433537095753,1.0518631876141677,3.216971700945713,1.9052193254674539,0.8146695274841804,4.532434582730561,1.5163092109824257,3.7308586665648917,2.3906090028362055,4.612848586128037,3.607963906364529,3.377912444329268,1.9594502904277618,2.298315956135598,0.6548165485161546,3.625678887509946,5.034955806639777,4.6116801526805515,2.7101153703669607,1.7834762830826085,3.983214682869474,0.5208317929402202,4.029783505946088,4.499710365622131,3.342905470153565,1.7539083451390665,1.816950066978464,5.342349860583442,1.5814413518421762,0.6063631608980553,1.1723910948633616,2.206714295988575,2.0510661769535306,2.842405579751711,5.035795340196291,4.13856085953461,4.969404049793779,1.452343601644775,3.9320345420556198,4.263597234580125,3.051568107826161,5.3747257784265665,2.9386871045274816,4.993035309142321,5.187882314483974,3.306449890473676,0.9667445465767905,5.051435376905429,1.7795989202157187,4.2376624900705435,3.919178318746201,4.077662500290293,1.7629141484985018,1.8229338551164904,2.932920633604806,1.0108965304975865,4.745613067982978,3.4716869337281815,4.334510159412838,2.347354776403154,2.0834622360544834,3.7930723743409738,1.5744349997059166,1.8392620418161494,4.9350901630290025,1.132992107229824,1.3977502178254415,1.074519077676292,0.792832996661182,1.53516647261379,3.0667755133074865,5.168490498857615,1.2112222148054836,4.054369929621286,5.324490189950622,2.8948933007171136,4.126484992423684,1.4412008221593837,3.156154070022414,1.848514526083443,2.666639371093445,1.2574338026682625,1.3290913010791312,5.124185488709119,3.473487597332282,1.0374613019433372,0.7930394111368395,2.457363748697012,0.7701296817753973,2.5591527236560427,5.329802261520074,5.380646112790172,1.741179223059734,3.0955875277165696,3.3779185754107006,1.4748960434054288,3.0822276774014696,2.6086888717464602,3.4260769348376305,0.9276805644839377,1.2482800783557209,1.4970462337338137,3.8511601361532897,3.1852506815747725,0.8980345803689422,1.04017157924526,5.060504496824197,2.313566493928419,4.203538794606311,1.0863308889146246,2.8860078158356246,1.3476539312202236,4.299288020905487,1.4192304474767687,1.45607987930179,1.7861897755042206,1.420556406563734,3.4571300596235695,4.883953419049859,2.3381256821020964,1.478731244011003,4.936650609503813,2.3239108015249075,3.358502967437012,2.1011623264240793,5.0870972645602475,3.440213190590079,5.226029964338533],"expected":[3.23600448477042,0.831676553308887,30.019857044091022,1.0150770463207444,2.5528640450064026,5.311859198382411,0.9238747487031688,2.5156976560621374,3.515571713353154,1.4232813910012392,0.5180898860462665,2.6101189039557013,59.44500648486634,2.346910335941724,0.6151966124468988,4.760434492703767,8.506504171202046,1.23861059009712,6.098605165599269,5.058387719007925,3.0895796457831195,104.90064125800782,1.4203560751480417,3.0108757759135067,0.6269447023331094,1.6131640627580082,2.9014473912226033,0.463155721045574,4.3447443405668595,1.8511179316722322,3.3386132388325875,1.3906259102901801,6.260213949215299,1.198971634918966,1.7341833126060462,3.6613776404343805,0.2407322463300831,2.6602390299686545,0.8299938982008261,5.683972018357419,6.061328488986405,3.1187842128450267,5.530200823949757,0.7144297208943677,1.4681713856900789,0.8893189450742254,4.462484539682976,0.6456439642002442,4.996259029602964,7.294110039988699,1.929293333880426,1.458494376506902,4.068753239776553,7.457159376281399,2.2192585998122656,1.8496838307300776,913.2129600574741,4.179993109321653,0.7588471723625235,3.829955800979615,1.657754303613058,1.767781323366129,1.3827709371235037,1.4491656714750474,1.937223354461703,5.2932240518627465,1.2199128353590596,0.7606894092869472,7.0494049039750895,0.10534196490416209,0.20119725401482733,61.03715160483853,3.941464487943175,1.0524117840633853,1.5792186299333684,2.9806407626452867,1.0053870821069497,0.4622892594518397,1.3182902626891042,4.00020216853453,3.1016182953331723,2.02809180915312,2.8041994970241557,4.4781867876196495,1.3242558749455127,2.2607972383786668,1.9156115737241435,0.32530865868126796,7.325130863050013,10.28105239834984,1.9683579802867606,3.5336779684234934,5.422138618054966,1.4842483702427218,2.648097785232984,0.35032477304376874,0.633175408962442,69.77261331811806,4.534539824010378,2.688602057684687,2.8369172708107273,20.16355656625984,0.7314755458761124,11.668805479427704,0.9128994494619904,1.9326183764721319,3.627772939597587,8.654348318276684,0.46107952978265077,3.0401528261004644,2.806239511818972,1.557591056065261,9.592585284249116,1.3801914174789618,1.499562062322142,2.6642328579697043,3.124585922116417,2.8853032807578076,0.28084003636120886,19.68744226410875,6.906892739572919,31.007306735226074,3.3828360409289893,6.853141943093153,22.38647950300117,7.419733536207257,9.000065890983914,6.32053954015382,20.117230293638148,7.104524807563305,6.893398572574879,1.195271302086156,0.45980847333228086,3.2797645521273626,1.4892891824570924,1.1316592062119126,0.577191816510451,1.1982712956245358,0.9409239235897988,9.439926914614192,0.4908199945683819,2.527816668378224,0.6301704045567221,1.8023393605411697,1.779604228401042,0.05060966517426851,4.804056814738615,2.9269449293431147,4.711308051793279,5.4913860998329564,0.6224906479602608,18.113545196518658,4.339533815614755,2.8174903973521013,3.0287757752593207,2.3036973912665517,56.228987366819545,1.0898650594341301,63.531960195636344,4.5051342251266,1.1068933593213168,2.436586344893353,4.115129020012653,2.6100368952561244,1.212316549236754,9.274575900299794,4.907450895896313,1.8980773369562978,1.7570878856295018,3.0411725488433543,0.9933803971242798,4.717881428626222,5.171182448431235,0.8860163706360974,4.286543094849593,7.491287279803677,0.5917560865578579,1.4585585730648665,0.3249162695161609,2.59453019575183,2.0396965184338507,3.529527256325427,3.199187632337342,1.0659987558334563,1.7107954536002263,4.38996038524954,2.027414310767995,4.82152455768166,10.356696360633089,2.609998210613581,1.8245544425502789,4.311643958054227,3.181379307528015,0.1791017209052759,2.633568010568804,4.160840186673719,1.4117776721514248,3.2680220346065427,1.5755446072405315,20.546965691433698,0.6004415333718824,4.599196258726264,12.261591344436184,6.154165751488408,4.013387434715503,5.138767614268564,0.12940928607374305,3.2724287899833473,1.0199238763392002,0.5631295409129126,4.313044305985491,1.7047237804424609,4.864193452387466,10.487323548907483,0.1750280407003503,4.462065550459644,19.54180863402103,4.472118219532502,12.895230495446278,7.885589826218663,1.8005881635493197,1.1150722165462583,1.7433897067487687,2.123886835623365,1.5676323929102156,2.5966787507890774,1.3400291119445527,5.026825663046705,1.0511456182276862,4.168635880071714,4.335507920009639,4.1574229830507425,5.292779250995788,1.5285795037234295,3.469491403613429,2.800283759451928,0.8701425341083036,7.304097119202111,4.155402132131313,2.028712738936839,6.749248820758713,3.660286069336771,7.115483722130639,0.5387437183771515,5.34139416877008,2.3407403597301943,6.916402731126223,0.2184408703460125,21.84684601315025,0.4838950017334491,0.02280622464525952,1.8995398299721984,3.0878904685510737,1.557076170438842,0.9014258197911451,2.875458608919333,1.760481706612965,11.657185963453607,4.588850684492608,0.5455744971596431,11.646294745900702,0.3639538804086588,2.099270780278996,13.263398442793655,2.7283133092398177,29.448049455502957,1.031929486597411,2.04148924884033,18.583519271267615,0.9351942727646488,0.2596701538883136,2.067267664454123,7.854873060747301,7.27510059846443,6.646705083876212,4.374022416247195,0.5937658271941105,0.2690020255382095,0.7909994601678731,2.132548928688583,6.928004041543427,0.8759400596917787,0.8364583947804355,2.2864557250064923,2.4358599045776446,2.508310922079349,4.086772763128608,1.2514954131180263,1.22716019593618,2.1417841331678447,4.830192134856137,0.09050558902675751,3.797516640743391,0.8565364317139055,7.72786194756275,3.3690641589990435,8.162855729172772,1.6863023864476037,7.23168946182544,17.295410400806052,2.541408671825441,4.1556213401324955,0.22028406909032316,3.5886985360123282,5.065147686869012,1.5855005020934512,1.837222334322593,6.975395300201127,0.842613304674443,0.04566782444887097,3.181671000574575,1.3429715498613413,2.3470462988361165,0.9021988541477853,31.04428107997178,7.379117308625278,10.306879457499935,4.123711929320683,3.6275190345822907,1.4153312410155803,4.2574862486976865,11.908453891289716,2.6078207377410685,1.0008965585094505,1.5336197774834526,0.7424177878309809,0.44242811305152363,3.872438331519493,2.249630523667177,23.66210992449676,4.179572323485541,1.0900227458128557,2.4822061356503173,0.5227093194764488,4.513035660140383,3.1761179892331963,7.095192462934724,1.164976895991711,0.5882448717662608,6.040709777025095,7.811581721012458,3.995027375125998,2.0429295957479012,5.295189389083361,1.0917351954439858,14.694339843716058,0.7228546983813894,0.24861424354823947,3.7273223200607033,2.54136778651247,1.696775794883328,0.5565227018600403,2.479790625116565,0.8911712422488247,6.473450568476686,6.327517310913624,0.04027565642074357,0.4696574324578695,2.6524427240872335,0.3424342319786016,2.800220102268768,1.992160772515188,0.7650291418741464,0.25360629866582385,0.6075149670608064,1.344860145427511,0.8916593045279768,3.7777111385143933,2.3586858643078195,6.831349429036249,9.747037326071094,0.05011355041974141,8.269340252252452,6.09445120973477,0.7521913470278486,4.711868795696705,4.531570871502294,1.5787924939236768,7.061294040081794e-7,4.849826676591386,3.4320020681510126,3.776872517274855,1.7758412241156556,2.2871809798622307,2.074515028654163,2.8439094493591544,12.936897934489172,6.121629181888728,0.5783319924630023,5.5942609554254075,3.548247759882904,4.30641731096086,3.210348415309961,7.022109022819481,1.5524156563885112,1.717603925923962,5.801174525554857,3.894005422933523,139.96210814054493,4.873967722212562,20.469814936826822,0.06779621507536518,2.099754130368849,0.17211614672055914,0.35561566678527967,3.405466279524826,6.170842753114404,8.789679981551332,6.620624953558959,1.0611939116362596,3.0263398109319892,3.320499268132158,1.2291640802286554,1.7498564688560967,4.09246032927395,1.5044075419705236,86.79225527036209,1.7480851326870757,0.8706788430744328,3.324817567292634,13.596435565138677,9.900322531381402,2.411585810242355,1.9209781531037575,2.8912981553745722,3.1759675288095566,0.3418778872465343,12.509506667451042,4.853231681879661,3.8317006082418574,0.4043598562802314,2.349329548442857,1.8616422420129988,2.270049163095772,1.237314149338705,8.346215425907088,5.605543623931096,1.421983652589313,1.4057669764259968,19.81963455925373,0.47015690359564166,2.354052235822415,1.3427600821503811,1.0124719797737116,0.055618881358777,1.7500308331659409,1.1513602129939708,1.5826771468744996,4.9858367302625455,2.94964642952033,0.47986890356605694,2.9380338815412648,0.6687365410089686,1.3279470906542374,1.316882896298408,31.20105503115844,0.7979504171219209,2.1612964131973653,324.68997211835847,4.500372947009708,5.281143324042166,2.599127570204681,1.1902831987719291,13.659836371222774,1.346135324588207,2.253609245229092,0.42610414187573564,1.2621366028641117,3.525633250519837,0.5913000676673906,0.9723325365594506,3.098141170407527,6.241065141538017,35.92864361690873,1.8630823925088977,2.1280339876375653,3.322454784191478,3.1985129133078742,1.9633703587176206,4.912713940592748,0.7696713268993519,549.3269834758223,2.3117435386655862,17.924830210161755,1.0864282160609369,2.3908285129705877,8.113554077652658,0.8987639051825284,0.544141426391731,2.296658305408585,1.937284048776896,3.3098082893481124,2.119800249716533,5.815117890516572,4.258757274584757,3.8415148133592,7.870116982532798,9.673649467699244,2.135847885160202,4.714945942068055,2.3976256056083236,2.7761040236263175,136.60062256036292,1.8491296717405592,1.6132125305474685,1.8480864294841486,3.9958410862201057,2.9546699260372757,1.1370991357141338,0.9252293752705081,2.5413720877553416,0.37727863198746797,7.2942923546432334,1.6788335693039051,1614.4314732132257,4.172638825463579,4.862659199900005,2.037110557211282,2.194966643734856,3.2563046744507234,8.802176752980525,9.892019687228903,3.6146439501530874,4.012256843915968,4.224433864482914,1.7127020346380182,1.62416145305121,22.368034207025364,3.3173572881943025,19.067291126782795,37.80093601955066,4.316459395080641,4.36966085363667,14.339790690655704,0.5005838587630127,1.9664380295300068,0.8020586311123764,2.634870752581821,2.58792240518913,3.888698069730152,3.918472810187143,1.392198301936746,1.2571845198410407,6.376901381708499,3.041676391570436,4.091300619360704,1.2621000680811592,1.6444478105161822,3.2377006393595416,8.15862657962295,0.5762745545548534,6.209611574283717,5.842207679166758,1.190064730477296,4.845071851353014,0.683445845976812,0.5778023330862139,0.5289722238661997,3.531974193610429,1.53821870513329,1.7381014795908305,1.6428905625156711,5.733533540337912,3.1004457244964154,498.78085674951444,6.801024132115413,3.27383975926559,10.520154629500958,3.1839251773046167,15.707306410864478,1.6916498280137309,12.217738785663224,1.1578129157811154,0.7062133769001953,1.428243104554202,1.841916228068834,20.07609451950539,0.16989552619931347,1.017110864934279,4.142817702598488,0.2758638941121097,5.147376751551234,3.7774837372090384,3.73883801645184,1.3853738575397352,3.614954950313208,15.81107386158101,3.25447657962481,3.0372948079509543,4.13722772749292,2.9948110633873237,2.461172839847427,4.007616291035329,4.582533176474025,1.262389859221526,2.6482716605839003,6.49206557426488,0.6004074568836076,2.9816688409609187,0.6957506392929039,4.440384954096149,1.1554172914230496,4.262905722868923,4.277497701257619,7.316196771624775,2.1801579928179193,0.8676954741530069,1.3110399510732842,6.094514597480653,2.6954416864919843,3.571994725135746,3.744541304987373,1.5669845851274333,4.563219891004939,0.7400317102750165,5.264937239867499,1.8590377046225026,1.3298439723098847,5.219670349923316,19.876452116082767,2.8335341952163873,2.919109684363881,1.3854342924011713,0.07591385255809888,1.767047520452343,3.4033946236636003,187.45649114452107,1008.7815006711855,2.466993006499504,2.464179996568128,2.6961972520916415,7.409280147579027,5.077280257048109,2.1627775786176695,4.959921137883569,0.7485981245341887,2.9694367604481005,0.6983464567541012,1.4030745746629014,2.4962310998810593,4.5531777938203035,2.7249714725964593,3.257414080202042,3.707864619478517,3.2197822538528635,2.333627189742376,1.0034734791609903,1.0664761102375717,0.25139219597939033,8.133625943193723,51.42064139417776,7.053450279956029,1.100951021266181,3.445586926829071,6.35732671396686,1.4606944431664606,4.671210038291077,1.6144979684940652,1.1834008389767718,6.4964857152540345,0.7360445998803046,1.1128081413096027,3.4806253126435966,3.366143827090984,3.909292484085817,7.175555164984693,4.2467529313067995,2.5786114006938265,2.3396165939872535,2.5309414693974843,0.9470304956180506,4.458610949581111,3.817984443071777,4.226374079726363,5.810908332990899,3.7071067175187054,3.0444240145511388,0.8582362467792182,83.53943151500495,3.669506488716319,1.5983247976158768,5.474018358943827,0.8858306837483996,6.8829360017065655,4.748702778063687,4.323274709804871,11.26730095717071,5.420781862443811,2.0043671627300017,3.4682109556385523,2.0477699194241996,3.5086810719621977,0.04861563174602143,6.775558348958254,1.97252538467975,2.592566659272044,1.3509365765512082,48.888716560073505,1.6116563328408093,1.8142127585100163,1.0724808344461105,1.0160466575163905,4.713029014353238,4.815330966799729,1.5542947860889138,2.3180091950739308,2.722382227389799,2.9847578978741973,1.7566089350580005,1.6198890828443429,3.042859063378352,11.137897834173378,7.760275124657169,4.812900308552831,1.6724958348375403,0.78675250184975,2.4006982389598406,9.7731168755275,0.9859975904875219,2.7167275178449803,0.933347913845102,2.2857772393742897,2.1289940849346545,0.6085161452447143,2.462935167231985,2.323422485920443,1.508165863534636,2.319940373957666,4.353108337574666,3.6056512173052746,1.6041579527993906,3.3457114925051226,1.5032167562904695,1.9984632475003756,1.1462377953435445,0.9026110221513136,2.3215417476389755,1.2148578098224485,4.312870095509429,0.8795000291852579,10.91239443411368,1.330246086857289,0.9698056558066307,4.538779115075314,2.1230883566591165,1.54910945072583,42.058571242463564,1.3006755471996048,4.257628212862269,3.602076104460946,0.6311977803550011,0.7784106334956357,8.250080194987707,10.872975924027333,5.488824727961615,10.568651470446941,3.2639602997369597,0.541245719158569,4.3827166663230015,6.558443345117297,0.4009625012494629,2.5894117294099894,5.7681137604715875,1.1750124646176245,5.552893393492044,4.8045093622947705,3.485160170948179,2.8975412405037595,5.641167897288066,1.9149778617715154,1.953588758344042,0.5975432935762259,3.1820032727053382,3.619612844662986,1.029595201744354,2.502205617455114,2.6740465759237066,2.5313598044469283,6.272674735699477,4.468657732906528,2.9818677193102436,6.65484145274868,2.6456975717202407,2.1064383703608236,5.617945531048344,0.5945858237157103,4.027102727375932,2.2634311877245135,4.5093863502339895,10.857597007328913,2.471037313075313,3.2515385304071316,5.412519666472381,1.0908455046094694,3.1425040673741935,4.321638031251242,1.8276529515502433,4.4425522405829305,0.5853985929840844,1.2460360055635722,1.0270709508022249,2.1631292477594823,1.5766425433884583,1.8443463214566609,4.648963413335335,0.4038177294409866,1.1609287579272118,1.9825499332558332,89.71253985366867,8.624038948902589,1.0009498214209678,5.019446092502198,2.2558759042675183,2.3505317794529343,2.894722592851537,2.5214148302576054,7.421742444756916,2.7377997807474803,1.4976133973446395,4.872411327741998,0.3301637940377969,2.16005381195624,6.686748749970111,11.04592635677144,0.42647865604732965,4.540271415686098,4.6959217758178875,8.773510427163764,3.95580946435505,4.64809024518444,0.5766700705278327,4.318705403509945,0.1386792432409755,1.1255405052872525,4.613893086435092,1.4371289681067032,9.661058522944527,4.027639600622145,2.557949160764498,0.49989486331376676,2.722846208565226,5.340973173173273,2.3102988427317133,0.09172814835453735,10.503603788133539,2.1779736811462174,5.58685097151958,0.9034064541174462,0.8715604306284469,3.4568046029497737,2.3270963222330385,6.4376093434494095,11.000525408633544,1.5203546085420476,7.4209394056229065,0.3626216571063751,1.0337115113436368,7.963777875122144,0.6729905494410807,3.3495577908821406,1.4309906541654898,3.9212115096029803,4.568572983219121,1.5841190298433194,4.104155082783144,2.405256935897165,2.271150524382454,0.8875414035470159,5.153774544519396,0.647754008353031,5.066271144546993,2.5621975001977844,1.85166154256396,2.0432627637824674,0.6699381236306344,4.805708549208346,3.664484101566349,4.160861492151186,1.9327515415373537,3.4761422921720313,2.5508783616143242,3.215299878739553,4.375508881169582,8.379031718667408,12.033495001166916,1.9154283516592558,0.06395528246849452,2.3625971374725547,1.1674585127143664,0.7440160223194183,1.407581617089268,1.7445851238658086,4.544277066311236,0.42008088630833,3.6874519161026225,1.742450913003301,2.0753614571650543,4.36248015630187,11.724238241977837,0.7828739493770709,7.431987061132116,0.5097602347464121,5.88377454995608,1.8819642419924907,5.756753835027621,5.104581654100198,2.7845855028320057,1.6491176222158643,2.569854615052296,1.117705484609397,1.1068006922503228,0.8566599876107821,18.769368168220733,5.171023877288888,2.8013081566592875,3.0813150762558017,2.067395534973979,0.2836406748653321,0.8285835100285762,4.175477326392046,11.545589743514755,4.142843076634321,2.620082112895413,5.379677606499003,4.6701233635560895,7.029818535417523,1.5862988573656194,0.39861431103517087,1.9674645459769509,7.319169056976123,1.0228747343742148,1.541047513921192,3.017572868255801,3.282669503196198,8.656107915828134,6.944513444876854,5.194704787126141,1.5347115519211951,3.267596944157391,1.149265348369166,3.1599905657816243,6.455905733116663,4.577384260459624,2.230544385319871,0.32906104441396905,6.872226059894311,0.46940340811024756,4.225308534983218,0.773605904934406,9.388156516709216,4.60571501829914,9.414075789432355,1.0636487474329825,9.18525863064836,1.414585887258405,2.275926524075159,8.595368503466984,0.17774049382522455,0.33328180477548514,0.6640164885538107,0.3506890656907691,1.4998106012474348,4.20404630996479,5.138312112222066,3.274086398916585,1.0431958564548962,1.2099543664872328,19.617715843550535,0.6092859128205418,0.22741486158391022,2.197007575693711,0.2406417348649485,7.630103294966888,1.6241543702539591,4.220907936874101,1.1493395361014083,0.13319410323878408,2.997231085086213,6.379409434306156,9.03394037740535,3.272713133333114,38.59091751925124,4.193733437920293,0.8484486477848734,3.640477002974659,4.161375266556199,0.865452557866078,1.7714970742191736,1.1372034285834374,0.5620266981757687,1.071930522860589,1.5101787135166493,1.8523564016042937,1.1026630797010204,4.541792973457086,42.46760978365481,4.138384084840919,7.8231589817136715,3.6022940818801024,2.033253059637609]} diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/runner.jl b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/runner.jl new file mode 100644 index 000000000000..b3b660fe32ef --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/fixtures/julia/runner.jl @@ -0,0 +1,77 @@ +#!/usr/bin/env julia +# +# @license Apache-2.0 +# +# Copyright (c) 2026 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. + +import Distributions: quantile, LogLogistic +import JSON + +""" + gen( p, alpha, beta, name ) + +Generate fixture data and write to file. + +# Arguments + +* `p`: input probability +* `alpha`: scale parameter +* `beta`: shape parameter +* `name::AbstractString`: output filename + +# Examples + +``` julia +julia> p = rand( 1000 ); +julia> alpha = rand( 1000 ) .* 5.0 .+ 0.5; +julia> beta = rand( 1000 ) .* 5.0 .+ 0.5; +julia> gen( p, alpha, beta, "data.json" ); +``` +""" +function gen( p, alpha, beta, name ) + z = Array{Float64}( undef, length(p) ); + for i in eachindex(p) + z[ i ] = quantile( LogLogistic( alpha[i], beta[i] ), p[i] ); + end + + # Store data to be written to file as a collection: + data = Dict([ + ("p", p), + ("alpha", alpha), + ("beta", beta), + ("expected", z) + ]); + + # Based on the script directory, create an output filepath: + filepath = joinpath( dir, name ); + + # Write the data to the output filepath as JSON: + outfile = open( filepath, "w" ); + write( outfile, JSON.json(data) ); + write( outfile, "\n" ); + close( outfile ); +end + +# Get the filename: +file = @__FILE__; + +# Extract the directory in which this file resides: +dir = dirname( file ); + +# Generate fixture data: +p = rand( 1000 ); +alpha = rand( 1000 ) .* 5.0 .+ 0.5; +beta = rand( 1000 ) .* 5.0 .+ 0.5; +gen( p, alpha, beta, "data.json" ); diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.factory.js new file mode 100644 index 000000000000..a454e9698a55 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.factory.js @@ -0,0 +1,182 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var abs = require( '@stdlib/math/base/special/abs' ); +var PINF = require( '@stdlib/constants/float64/pinf' ); +var NINF = require( '@stdlib/constants/float64/ninf' ); +var EPS = require( '@stdlib/constants/float64/eps' ); +var factory = require( './../lib/factory.js' ); + + +// FIXTURES // + +var data = require( './fixtures/julia/data.json' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof factory, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function returns a function', function test( t ) { + var quantile = factory( 1.0, 1.0 ); + t.strictEqual( typeof quantile, 'function', 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `NaN` for any parameter, the returned function always returns `NaN`', function test( t ) { + var quantile; + + quantile = factory( NaN, 1.0 ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + quantile = factory( 1.0, NaN ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + quantile = factory( NaN, NaN ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided a non-positive `alpha`, the returned function always returns `NaN`', function test( t ) { + var quantile; + + quantile = factory( 0.0, 1.0 ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + quantile = factory( -1.0, 1.0 ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + quantile = factory( NINF, 1.0 ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided a non-positive `beta`, the returned function always returns `NaN`', function test( t ) { + var quantile; + + quantile = factory( 1.0, 0.0 ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + quantile = factory( 1.0, -1.0 ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + quantile = factory( 1.0, NINF ); + t.strictEqual( isnan( quantile( 0.5 ) ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the returned function returns `NaN` if provided `NaN` for `p`', function test( t ) { + var quantile; + var y; + + quantile = factory( 1.0, 1.0 ); + + y = quantile( NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the returned function returns `NaN` if provided `p < 0` or `p > 1`', function test( t ) { + var quantile; + var y; + + quantile = factory( 1.0, 1.0 ); + + y = quantile( -0.1 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 1.1 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( NINF ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( PINF ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the returned function returns `0` if provided `p = 0`', function test( t ) { + var quantile; + var y; + + quantile = factory( 1.0, 1.0 ); + + y = quantile( 0.0 ); + t.strictEqual( y, 0.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'the returned function returns `Infinity` if provided `p = 1`', function test( t ) { + var quantile; + var y; + + quantile = factory( 1.0, 1.0 ); + + y = quantile( 1.0 ); + t.strictEqual( y, PINF, 'returns expected value' ); + + t.end(); +}); + +tape( 'the returned function evaluates the log-logistic quantile for `p`, `alpha`, and `beta`', function test( t ) { + var expected; + var quantile; + var alpha; + var delta; + var beta; + var tol; + var p; + var y; + var i; + + expected = data.expected; + p = data.p; + alpha = data.alpha; + beta = data.beta; + for ( i = 0; i < p.length; i++ ) { + quantile = factory( alpha[ i ], beta[ i ] ); + y = quantile( p[ i ] ); + if ( expected[ i ] !== null && !( expected[ i ] === 0.0 && y < EPS ) ) { + if ( y === expected[ i ] ) { + t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha: '+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] ); + } else { + delta = abs( y - expected[ i ] ); + tol = 3.0 * EPS * abs( expected[ i ] ); + t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' ); + } + } + } + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.js new file mode 100644 index 000000000000..ea2a910239b6 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.js @@ -0,0 +1,38 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var quantile = require( './../lib' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof quantile, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'attached to the main export is a factory method for generating `quantile` functions', function test( t ) { + t.strictEqual( typeof quantile.factory, 'function', 'exports a factory method' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.native.js new file mode 100644 index 000000000000..f3a0dfb10a3e --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.native.js @@ -0,0 +1,165 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 resolve = require( 'path' ).resolve; +var tape = require( 'tape' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var abs = require( '@stdlib/math/base/special/abs' ); +var PINF = require( '@stdlib/constants/float64/pinf' ); +var NINF = require( '@stdlib/constants/float64/ninf' ); +var EPS = require( '@stdlib/constants/float64/eps' ); + + +// VARIABLES // + +var quantile = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( quantile instanceof Error ) +}; + + +// FIXTURES // + +var data = require( './fixtures/julia/data.json' ); + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof quantile, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'if provided `NaN` for any parameter, the function returns `NaN`', opts, function test( t ) { + var y; + + y = quantile( NaN, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, NaN, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, 1.0, NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( NaN, NaN, NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided a non-positive `alpha`, the function returns `NaN`', opts, function test( t ) { + var y; + + y = quantile( 0.5, 0.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, -1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, NINF, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided a non-positive `beta`, the function returns `NaN`', opts, function test( t ) { + var y; + + y = quantile( 0.5, 1.0, 0.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, 1.0, -1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, 1.0, NINF ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided `p < 0` or `p > 1`, the function returns `NaN`', opts, function test( t ) { + var y; + + y = quantile( -0.1, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 1.1, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( NINF, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( PINF, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided `p = 0`, the function returns `0`', opts, function test( t ) { + var y; + + y = quantile( 0.0, 1.0, 1.0 ); + t.strictEqual( y, 0.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided `p = 1`, the function returns `Infinity`', opts, function test( t ) { + var y; + + y = quantile( 1.0, 1.0, 1.0 ); + t.strictEqual( y, PINF, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function evaluates the log-logistic quantile for `p`, `alpha`, and `beta`', opts, function test( t ) { + var expected; + var alpha; + var delta; + var beta; + var tol; + var p; + var y; + var i; + + expected = data.expected; + p = data.p; + alpha = data.alpha; + beta = data.beta; + for ( i = 0; i < p.length; i++ ) { + y = quantile( p[ i ], alpha[ i ], beta[ i ] ); + if ( expected[ i ] !== null && !( expected[ i ] === 0.0 && y < EPS ) ) { + if ( y === expected[ i ] ) { + t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha: '+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] ); + } else { + delta = abs( y - expected[ i ] ); + tol = 3.0 * EPS * abs( expected[ i ] ); + t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' ); + } + } + } + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.quantile.js b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.quantile.js new file mode 100644 index 000000000000..293b760b16f8 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/log-logistic/quantile/test/test.quantile.js @@ -0,0 +1,156 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var abs = require( '@stdlib/math/base/special/abs' ); +var PINF = require( '@stdlib/constants/float64/pinf' ); +var NINF = require( '@stdlib/constants/float64/ninf' ); +var EPS = require( '@stdlib/constants/float64/eps' ); +var quantile = require( './../lib' ); + + +// FIXTURES // + +var data = require( './fixtures/julia/data.json' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof quantile, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'if provided `NaN` for any parameter, the function returns `NaN`', function test( t ) { + var y; + + y = quantile( NaN, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, NaN, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, 1.0, NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( NaN, NaN, NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided a non-positive `alpha`, the function returns `NaN`', function test( t ) { + var y; + + y = quantile( 0.5, 0.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, -1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, NINF, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided a non-positive `beta`, the function returns `NaN`', function test( t ) { + var y; + + y = quantile( 0.5, 1.0, 0.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, 1.0, -1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 0.5, 1.0, NINF ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided `p < 0` or `p > 1`, the function returns `NaN`', function test( t ) { + var y; + + y = quantile( -0.1, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( 1.1, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( NINF, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = quantile( PINF, 1.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided `p = 0`, the function returns `0`', function test( t ) { + var y; + + y = quantile( 0.0, 1.0, 1.0 ); + t.strictEqual( y, 0.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided `p = 1`, the function returns `Infinity`', function test( t ) { + var y; + + y = quantile( 1.0, 1.0, 1.0 ); + t.strictEqual( y, PINF, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function evaluates the log-logistic quantile for `p`, `alpha`, and `beta`', function test( t ) { + var expected; + var alpha; + var delta; + var beta; + var tol; + var p; + var y; + var i; + + expected = data.expected; + p = data.p; + alpha = data.alpha; + beta = data.beta; + for ( i = 0; i < p.length; i++ ) { + y = quantile( p[ i ], alpha[ i ], beta[ i ] ); + if ( expected[ i ] !== null && !( expected[ i ] === 0.0 && y < EPS ) ) { + if ( y === expected[ i ] ) { + t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha: '+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] ); + } else { + delta = abs( y - expected[ i ] ); + tol = 3.0 * EPS * abs( expected[ i ] ); + t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' ); + } + } + } + t.end(); +});