perl - tendencia - modo estadistico calculadora casio fx-350ms

¿Usando perl para encontrar la mediana, el modo, la desviación estándar?

Question

Tengo una matriz de números. ¿Cuál es la forma más fácil de calcular la mediana, el modo y el desarrollo estándar para el conjunto de datos?

Answer 1

#!/usr/bin/perl
#
# stdev - figure N, min, max, median, mode, mean, & std deviation
#
# pull out all the real numbers in the input
# stream and run standard calculations on them.
# they may be intermixed with other test, need
# not be on the same or different lines, and
# can be in scientific notion (avagadro=6.02e23).
# they also admit a leading + or -.
#
# Tom Christiansen
# [email protected]
use strict;
use warnings;
use List::Util qw< min max >;
sub by_number {
if ($a < $b){ -1 } elsif ($a > $b) { 1 } else { 0 }
}
#
my $number_rx = qr{
# leading sign, positive or negative
(?: [+-] ? )
# mantissa
(?= [0123456789.] )
(?:
# "N" or "N." or "N.N"
(?:
(?: [0123456789] + )
(?:
(?: [.] )
(?: [0123456789] * )
) ?
|
# ".N", no leading digits
(?:
(?: [.] )
(?: [0123456789] + )
)
)
)
# abscissa
(?:
(?: [Ee] )
(?:
(?: [+-] ? )
(?: [0123456789] + )
)
|
)
}x;
my $n = 0;
my $sum = 0;
my @values = ();
my %seen = ();
while (<>) {
while (/($number_rx)/g) {
$n++;
my $num = 0 + $1; # 0+ is so numbers in alternate form count as same
$sum += $num;
push @values, $num;
$seen{$num}++;
}
}
die "no values" if $n == 0;
my $mean = $sum / $n;
my $sqsum = 0;
for (@values) {
$sqsum += ( $_ ** 2 );
}
$sqsum /= $n;
$sqsum -= ( $mean ** 2 );
my $stdev = sqrt($sqsum);
my $max_seen_count = max values %seen;
my @modes = grep { $seen{$_} == $max_seen_count } keys %seen;
my $mode = @modes == 1
? $modes[0]
: "(" . join(", ", @modes) . ")";
$mode .= '' @ '' . $max_seen_count;
my $median;
my $mid = int @values/2;
my @sorted_values = sort by_number @values;
if (@values % 2) {
$median = $sorted_values[ $mid ];
} else {
$median = ($sorted_values[$mid-1] + $sorted_values[$mid])/2;
}
my $min = min @values;
my $max = max @values;
printf "n is %d, min is %g, max is %g/n", $n, $min, $max;
printf "mode is %s, median is %g, mean is %g, stdev is %g/n",
$mode, $median, $mean, $stdev;

Answer 2

• Statistics::Basic::Mean
• Statistics::Basic::Median
• Statistics::Basic::Mode
• Statistics::Basic::StdDev

Answer 3

Dependiendo de qué tan profundo necesita ir, la respuesta de erickb puede funcionar. Sin embargo, para la funcionalidad numérica en Perl hay PDL . pdl crear un piddle (el objeto que contiene sus datos) utilizando la función pdl . Desde allí puede usar las operaciones en esta página para hacer las estadísticas que necesita.

Editar: mirando a su alrededor encontré dos llamadas a funciones que hacen EXACTAMENTE lo que necesita. statsover proporciona estadísticas sobre una dimensión de un piddle, mientras que las stats hacen lo mismo sobre todo el piddle.

my $piddle = pdl @data;
my ($mean,$prms,$median,$min,$max,$adev,$rms) = statsover $piddle;