#!/usr/bin/perl -w use strict; my $sum; my $square; my $n; my $do_geom = 1; my $geom; while (<>) { chomp; s#_##g; my $x; foreach $x (split) { $sum += $x; $square += $x*$x; ++$n; if ($x == 0) { $do_geom = 0; } else { $geom += log($x); } } } my $mean = $sum / $n; my $mean_squares = $square / $n; my $sd = sqrt($mean_squares - $mean * $mean); printf "total: %15.4f\n", $sum; printf "mean: %15.4f\n", $mean; printf "sd: %15.4f %5.2f%%\n", $sd, 100 * ($sd / $mean); if ($do_geom && $geom != 0) { printf "geom: %15.4f\n", exp($geom/$n); }