X-Git-Url: https://git.piment-noir.org/?a=blobdiff_plain;f=src%2Futils%2FStatisticUtils.ts;h=9f50b8eaa12b08037ef9c49b042ac96a62463bef;hb=840ca85d7c40a6ee6f3a85a50e68dbea2f90acb8;hp=9ca4398ca76590250de3b1c5a935f42ca373cd6f;hpb=736073f0c51fc24b3d676fa2cedc6052f928fbbc;p=e-mobility-charging-stations-simulator.git diff --git a/src/utils/StatisticUtils.ts b/src/utils/StatisticUtils.ts index 9ca4398c..9f50b8ea 100644 --- a/src/utils/StatisticUtils.ts +++ b/src/utils/StatisticUtils.ts @@ -1,37 +1,10 @@ -/** - * Computes the average of the given data set. - * - * @param dataSet - Data set. - * @returns The average of the given data set. - * @internal - */ -export const average = (dataSet: number[]): number => { - if (Array.isArray(dataSet) && dataSet.length === 0) { - return 0 - } else if (Array.isArray(dataSet) && dataSet.length === 1) { - return dataSet[0] - } - return dataSet.reduce((accumulator, nb) => accumulator + nb, 0) / dataSet.length -} +import { mean } from 'rambda' -/** - * Computes the median of the given data set. - * - * @param dataSet - Data set. - * @returns The median of the given data set. - * @internal - */ -export const median = (dataSet: number[]): number => { - if (Array.isArray(dataSet) && dataSet.length === 0) { - return 0 - } else if (Array.isArray(dataSet) && dataSet.length === 1) { - return dataSet[0] - } - const sortedDataSet = dataSet.slice().sort((a, b) => a - b) - return ( - (sortedDataSet[(sortedDataSet.length - 1) >> 1] + sortedDataSet[sortedDataSet.length >> 1]) / 2 - ) -} +export const min = (...args: number[]): number => + args.reduce((minimum, num) => (minimum < num ? minimum : num), Infinity) + +export const max = (...args: number[]): number => + args.reduce((maximum, num) => (maximum > num ? maximum : num), -Infinity) // TODO: use order statistics tree https://en.wikipedia.org/wiki/Order_statistic_tree export const nthPercentile = (dataSet: number[], percentile: number): number => { @@ -48,17 +21,16 @@ export const nthPercentile = (dataSet: number[], percentile: number): number => if (percentile === 100) { return sortedDataSet[sortedDataSet.length - 1] } - const percentileIndexBase = (percentile / 100) * (sortedDataSet.length - 1) - const percentileIndexInteger = Math.floor(percentileIndexBase) + const base = (percentile / 100) * (sortedDataSet.length - 1) + const baseIndex = Math.floor(base) // eslint-disable-next-line @typescript-eslint/no-unnecessary-condition - if (sortedDataSet[percentileIndexInteger + 1] != null) { + if (sortedDataSet[baseIndex + 1] != null) { return ( - sortedDataSet[percentileIndexInteger] + - (percentileIndexBase - percentileIndexInteger) * - (sortedDataSet[percentileIndexInteger + 1] - sortedDataSet[percentileIndexInteger]) + sortedDataSet[baseIndex] + + (base - baseIndex) * (sortedDataSet[baseIndex + 1] - sortedDataSet[baseIndex]) ) } - return sortedDataSet[percentileIndexInteger] + return sortedDataSet[baseIndex] } /** @@ -70,10 +42,7 @@ export const nthPercentile = (dataSet: number[], percentile: number): number => * @see https://en.wikipedia.org/wiki/Unbiased_estimation_of_standard_deviation * @internal */ -export const stdDeviation = ( - dataSet: number[], - dataSetAverage: number = average(dataSet) -): number => { +export const stdDeviation = (dataSet: number[], dataSetAverage: number = mean(dataSet)): number => { if (Array.isArray(dataSet) && (dataSet.length === 0 || dataSet.length === 1)) { return 0 }