~kirkland/calc-stats/trunk

« back to all changes in this revision

Viewing changes to bin/stats

  • Committer: Dustin Kirkland
  • Date: 2017-07-26 22:42:41 UTC
  • Revision ID: kirkland@ubuntu.com-20170726224241-o97rbacgtu8mgz6f
initial commit

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
#!/bin/sh
 
2
#
 
3
# stats: perform statistical calculations on a file or standard input
 
4
#        containing numbers
 
5
#
 
6
# Copyright 2017 Dustin Kirkland <dustin.kirkland@gmail.com>
 
7
#
 
8
# Licensed under the Apache License, Version 2.0 (the "License");
 
9
# you may not use this file except in compliance with the License.
 
10
# You may obtain a copy of the License at
 
11
#
 
12
#    http://www.apache.org/licenses/LICENSE-2.0
 
13
#
 
14
# Unless required by applicable law or agreed to in writing, software
 
15
# distributed under the License is distributed on an "AS IS" BASIS,
 
16
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 
17
# See the License for the specific language governing permissions and
 
18
# limitations under the License.
 
19
 
 
20
 
 
21
set -e
 
22
 
 
23
# Work from a temp file, which we'll clean up when done
 
24
file=$(mktemp)
 
25
trap "rm -f ${file}" HUP INT QUIT ILL TRAP KILL BUS TERM
 
26
 
 
27
# Handle either a file as argument, or standard input
 
28
if [ -f "$1" ]; then
 
29
        cat "$1" >"$file"
 
30
else
 
31
        cat /dev/stdin >"$file"
 
32
fi
 
33
 
 
34
# Determine our operating mode by how we were called
 
35
case "$0" in
 
36
        *min)
 
37
                awk -v min=0 'NR == 1 || $1 < min {line = $0; min = $1} END {print min}' "$file"
 
38
        ;;
 
39
        *max)
 
40
                awk -v max=0 'NR == 1 || $1 > max {line = $0; max = $1} END {print max}' "$file"
 
41
        ;;
 
42
        *mean|*avg)
 
43
                awk '{ sum += $1; n++ } END { if (n > 0) print sum / n; }' "$file"
 
44
        ;;
 
45
        *median)
 
46
                sort -g "$file" | awk '{ a[i++]=$1; } END { print a[int(i/2)]; }'
 
47
        ;;
 
48
        *mode)
 
49
                sort -g "$file" | uniq -c | sort -r -g | head -n1 | awk '{print $2}'
 
50
        ;;
 
51
        *stdev|*stddev)
 
52
                awk '{ x+=$0; y+=$0^2 } END { print sqrt(y/NR-(x/NR)^2) }' "$file"
 
53
        ;;
 
54
        *histogram)
 
55
                sort -g "$file" | uniq -c | sort -r -g
 
56
        ;;
 
57
        *summate)
 
58
                awk '{ sum += $1 } END { print sum }' "$file"
 
59
        ;;
 
60
esac
 
61
rm -f "$file"