Используя awk
, вы можете смоделировать 2D-массив, построив составной индекс из ключа (значение первого столбца) и индекса столбца:
awk '
{
c[$1]++;
for (i=2;i<=NF;i++) {
s[$1"."i]+=$i};
}
END {
for (k in c) {
printf "%s\t", k;
for(i=2;i<NF;i++) printf "%.1f\t", s[k"."i]/c[k];
printf "%.1f\n", s[k"."NF]/c[k];
}
}' file
a 16.5 17.5 18.5
b 19.5 20.5 21.5
Аналогичный подход может быть реализован в perl
более напрямую, используя хеш массивов.
В качестве альтернативы, есть пакет данных GNU
, который (по крайней мере, начиная с версии 1.1.0) очень компактно поддерживает групповые средние значения, например.
datamash --sort --whitespace groupby 1 mean 2-4 < file
a 16.5 17.5 18.5
b 19.5 20.5 21.5
FWIW вот моя попытка решения perl, включая нормализацию к глобальному максимальному среднему значению, как запрошено в комментариях. ОТКАЗ ОТ ОТВЕТСТВЕННОСТИ: Я начинающий программист на Perl, поэтому он может продемонстрировать плохие методы программирования .
#!/usr/bin/perl
use strict;
use warnings;
use List::MoreUtils qw(pairwise minmax);
use Math::Round qw(nearest);
my @hdr;
my %sums = ();
my %count = ();
my $key;
while (defined($_ = <ARGV>)) {
chomp $_;
my @F = split(' ', $_, 0);
# UGLY: hardcoded to expect exactly 1 header row
if ($. == 1) {
@hdr = @F;
next;
}
# sum column-wise, grouped by first column
$key = shift @F;
if ( exists $sums{$key} ) {
$sums{$key} = [ pairwise { $a + $b } @{ $sums{$key} }, @F];
}
else {
$sums{$key} = \@F;
}
$count{$key}++;
}
my %avgs = ();
# NB should really initialize $maxavg to a suitably large NEGATIVE value
my $maxavg = 0.0;
# find the column averages, and the global max of those averages
for $key ( keys %sums ) {
$avgs{$key} = [ map { $_ / $count{$key} } @{ $sums{$key} } ];
# NB could use List::Util=max here, but we're alresdy using List::MoreUtils
my ($kmin, $kmax) = minmax @{ $avgs{$key} };
$maxavg = $kmax > $maxavg ? $kmax : $maxavg;
}
# normalize and print the results, rounded to nearest 0.01
print join "\t", @hdr, "\n";
for $key ( sort keys %avgs ) {
print join "\t", $key, (map { nearest (0.01, $_ / $maxavg) } @{ $avgs{$key} }), "\n";
}
Сохраняется как colavgnorm.pl
и делается исполняемым, затем запускается как
$ ./colavgnorm.pl file
K C1 C2 C3
a 0.77 0.81 0.86
b 0.91 0.95 1
, где файл
- это
K C1 C2 C3
a 12 13 14
b 15 16 17
a 21 22 23
b 24 25 26