w-statistic/src/arrGeometricStd.mjs

A tool for statistical analysis.

import size from 'lodash-es/size.js'
import isarr from 'wsemi/src/isarr.mjs'
import arrFilterByPnumAndToLog from 'wsemi/src/arrFilterByPnumAndToLog.mjs'
import arrStd from './arrStd.mjs'


/**
 * 計算陣列內有效數字之幾何標準差
 *
 * Unit Test: {@link https://github.com/yuda-lyu/w-statistic/blob/master/test/arrGeometricStd.test.js Github}
 * @memberOf w-statistic
 * @param {Array} arr 輸入陣列，只提取有效數字(或為字串的數字)進行計算
 * @returns {Number} 回傳幾何標準差
 * @example
 *
 * let arr
 *
 * arr = ['abc', '-2.5', -2.5, '-1', -1, '-0.1', -0.1, '0', 0, '0.1', 0.1, '1', 1, '2.5', 2.5, 22.5, 'xyz']
 * console.log(arrGeometricStd(arr))
 * // => 6.774243857996147
 *
 * arr = ['abc', '0', 0, '0.1', 0.1, '1', 1, '2.5', 2.5, 22.5, 'xyz']
 * console.log(arrGeometricStd(arr))
 * // => 6.774243857996147
 *
 */
function arrGeometricStd(arr) {

    //check
    if (!isarr(arr)) {
        return null
    }
    if (size(arr) === 0) {
        return null
    }

    //取大於0且轉log數值
    let rs = arrFilterByPnumAndToLog(arr)

    //n
    let n = size(rs)

    //check
    if (n === 0) {
        return null
    }

    //check
    if (n <= 1) {
        // throw new Error('至少要2點才能計算')
        return null
    }

    //r
    // let r = 0
    // each(rs, (v) => {
    //     r += (Math.log(v / geoAvg)) ** 2
    // })
    // r /= (n - 1)
    // r = Math.sqrt(r)
    // r = Math.exp(r)
    let sigma = arrStd(rs) //sigma為樣本標準差(用n-1)
    let r = Math.exp(sigma)

    return r
}


export default arrGeometricStd