math.statistics: adding trim-mean and signal-to-noise.
parent
bab0f663ae
commit
4c2d01c5ec
|
@ -11,6 +11,8 @@ IN: math.statistics.tests
|
||||||
[ 1.0 ] [ { 1 1 1 } geometric-mean ] unit-test
|
[ 1.0 ] [ { 1 1 1 } geometric-mean ] unit-test
|
||||||
[ 1/3 ] [ { 1 1 1 } harmonic-mean ] unit-test
|
[ 1/3 ] [ { 1 1 1 } harmonic-mean ] unit-test
|
||||||
[ 5+1/4 ] [ { 1 3 5 7 } contraharmonic-mean ] unit-test
|
[ 5+1/4 ] [ { 1 3 5 7 } contraharmonic-mean ] unit-test
|
||||||
|
[ 18 ] [ { 4 8 15 16 23 42 } 0 trim-mean ] unit-test
|
||||||
|
[ 15+1/2 ] [ { 4 8 15 16 23 42 } 0.2 trim-mean ] unit-test
|
||||||
|
|
||||||
[ 0 ] [ { 1 } range ] unit-test
|
[ 0 ] [ { 1 } range ] unit-test
|
||||||
[ 89 ] [ { 1 2 30 90 } range ] unit-test
|
[ 89 ] [ { 1 2 30 90 } range ] unit-test
|
||||||
|
|
|
@ -32,6 +32,9 @@ IN: math.statistics
|
||||||
: contraharmonic-mean ( seq -- x )
|
: contraharmonic-mean ( seq -- x )
|
||||||
[ sum-of-squares ] [ sum ] bi / ; inline
|
[ sum-of-squares ] [ sum ] bi / ; inline
|
||||||
|
|
||||||
|
: trim-mean ( seq p -- x )
|
||||||
|
swap [ length [ * >integer ] keep over - ] keep <slice> mean ;
|
||||||
|
|
||||||
<PRIVATE
|
<PRIVATE
|
||||||
|
|
||||||
:: ((kth-object)) ( seq k nth-quot exchange-quot quot: ( x y -- ? ) -- elt )
|
:: ((kth-object)) ( seq k nth-quot exchange-quot quot: ( x y -- ? ) -- elt )
|
||||||
|
@ -255,6 +258,8 @@ ALIAS: var sample-var
|
||||||
|
|
||||||
ALIAS: std sample-std
|
ALIAS: std sample-std
|
||||||
|
|
||||||
|
: signal-to-noise ( seq -- x ) [ mean ] [ std ] bi / ;
|
||||||
|
|
||||||
: mean-dev ( seq -- x ) dup mean v-n vabs mean ;
|
: mean-dev ( seq -- x ) dup mean v-n vabs mean ;
|
||||||
|
|
||||||
: median-dev ( seq -- x ) dup median v-n vabs mean ;
|
: median-dev ( seq -- x ) dup median v-n vabs mean ;
|
||||||
|
|
Loading…
Reference in New Issue