power-sums
v1.1.1
Published
sums of array items raised to incremental exponents
Downloads
27
Maintainers
Readme
power-sums
• Introduction • Limitations • Usage • Test • License •
Introduction
If a sample array and an exponent are provided, returns the sum of array values raised to one exponent
f([A],b) = ∑ai^b
If a sample array and an array of length N are provided, returns the same array filled with power sums for 1 to N.
f([A], [B]) = [ ∑ai, ∑ai^2, ∑ai^3, ...]
Limitations
While Kahan summation is used internally to reduce floating point errors, care must be taken if different sums are to be substracted for variance or skew calculations (e.g large data sets far from the origin)
Usage
single power
With a single exponent, powerSums(samples, N)
returns a single sum of values raised to exponent N.
powerSums = require('power-sums')
powerSums([0, 1, 2], 2) // return 0^2 + 1^2 + 2^2 = 5
powerSums([1/4, 1/2, 1], -1) // return 4 + 2 + 1 = 5
The internal methods for the most common cases are also exposed:
powerSums = require('power-sums')
powerSums[1]([0, 1, 2]) // direct normal sum, with floating point error correction
powerSums[2]([1/4, 1/2, 1]) // direct normal sum-square, with floating point error correction
multiple powers
The actual use case is to obtain multiple nth order sums in a single pass for calculating raw and central moments. The maximum order is the length of the parameter array and the results will be inserted in this Array before it is returned.
powerSums = require('power-sums')
powerSums([0,1,2], Array(4)) // returns [∑xi, ∑xi^2, ∑xi^3, ∑xi^4] = [3, 5, 9, 17]
The array returned is the same provided.
If the same result container is used multiple times, the content will be overwritten.
powerSums = require('power-sums')
var sums = [9,9,9]
powerSums([2,2], sums) // sums = [4, 8, 16]
edge cases
powerSums = require('power-sums')
powerSums([], Array(3)) // returns [0, 0, 0]
powerSums([0,1,2]) // returns the simple sum of 5
powerSums([0,1,2], []) // returns the empty set []
Test
In node, from the root folder type npm test
.
(test is not included with the package and must be obtained from the git repository)
License
Released under the MIT License