Basic contracts of functions:

• Functions do not modify their arguments e.g. change their order
• Invalid input, like empty lists to functions that need 1+ items to work, will cause functions to return null.

Basic Array Operations

.mixin()

Optionally mix in the following functions into the Array prototype. Otherwise you can use them off of the simple-statistics object itself.

.mean(x)

Mean of a single-dimensional Array of numbers. Also available as .average(x)

.sum(x)

Sum of a single-dimensional Array of numbers.

.variance(x)

Variance of a single-dimensional Array of numbers.

.standard_deviation(x)

Standard Deviation of a single-dimensional Array of numbers.

.median(x)

Median of a single-dimensional array of numbers.

.geometric_mean(x)

Geometric mean of a single-dimensional array of positive numbers.

.min(x)

Finds the minimum of a single-dimensional array of numbers. This runs in linear O(n) time.

.max(x)

Finds the maximum of a single-dimensional array of numbers. This runs in linear O(n) time.

.t_test(sample, x)

Does a student's t-test of a dataset sample, represented by a single-dimensional array of numbers. x is the known value, and the result is a measure of statistical significance.

.sample_variance(x)

Produces sample variance of a single-dimensional array of numbers.

.sample_covariance(x)

Produces sample covariance of two single-dimensional arrays of numbers.

.sample_correlation(x)

Produces sample correlation of two single-dimensional arrays of numbers.

.quantile(sample, p)

Does a quantile of a dataset sample, at p. For those familiary with the k/q syntax, p == k/q. sample must be a single-dimensional array of numbers, and p must be a number greater than zero and less than one.

.iqr(sample)

Calculates the Interquartile range of a sample - the difference between the upper and lower quartiles. Useful as a measure of dispersion.

Also available as .interquartile_range(x)

.jenks(data, number_of_classes)

Find the Jenks Natural Breaks for a single-dimensional array of numbers as input and a desired number_of_classes. The result is a single-dimensional with class breaks, including the minimum and maximum of the input array.

.r_squared(data, function)

Find the r-squared value of a particular dataset, expressed as a two-dimensional Array of numbers, against a Function.

var r_squared = ss.r_squared([[1, 1]], function(x) { return x * 2; });

Regression

.linear_regression()

Create a new linear regression solver.

.data([[1, 1], [2, 2]])

Set the data of a linear regression. The input is a two-dimensional array of numbers, which are treated as coordinates, like [[x, y], [x1, y1]].

.line()

Get the linear regression line: this returns a function that you can give x values and it will return y values. Internally, this uses the m() and b() values and the classic y = mx + b equation.

var linear_regression_line = ss.linear_regression()
    .data([[0, 1], [2, 2], [3, 3]]).line();
linear_regression_line(5);

.m()

Just get the slope of the fitted regression line, the m component of the full line equation. Returns a number.

.b()

Just get the y-intercept of the fitted regression line, the b component of the line equation. Returns a number.

Classification

.bayesian()

Create a naïve bayesian classifier.

.train(item, category)

Train the classifier to classify a certain item, given as an object with keys, to be in a certain category, given as a string.

.score(item)

Get the classifications of a certain item, given as an object of category -> score mappings.

var bayes = ss.bayesian();
bayes.train({ species: 'Cat' }, 'animal');
bayes.score({ species: 'Cat' });
// { animal: 1 }