1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
|
#!/usr/bin/gawk -f
BEGIN {
erriter = -1;
}
{
iter = $1;
if (iter == erriter)
next;
s = "";
i = 2;
while (index($i,"-")) {
s = sprintf("%s %s", s, $i);
i++;
}
s = sprintf("%s %s %s %s", s, $i, $(i+1), $(i+2));
i += 3;
error = 0;
for (j = i; j <= NF; j++) {
if ($j == -1)
error = 1;
}
if (!error) {
for (j = i; j <= NF; j++) {
sum[s,j-i+1] += $j;
sum2[s,j-i+1] += ($j*$j);
}
count[s]++;
}
else
erriter = iter;
numval = NF - i + 1;
}
END {
for (s in count) {
printf "%-45s ", s;
for (j = 1; j <= numval; j++) {
mean = sum[s,j]/count[s];
var = sum2[s,j]/count[s] - mean*mean;
printf "%6.2f %6.2f ", mean, sqrt(var/count[s]);
}
printf "\n";
}
}
|
