如何计算每一列的百分比?

如何计算每一列的百分比?

我正在尝试将一些数据转换为每列总数的百分比,与此线程非常相似,只是我需要对每列执行此操作: 使用 AWK 计算并除以总计

数据会喜欢这样的东西(但更多的列和行):

ID     Sample1     Sample2      Sample3
One      10          0            5
Two      3           6            8
Three    3           4            7

所需的输出如下所示:

ID     Sample1     Sample2     Sample3
One     62.50        0.0        25.0
Two     18.75       60.0        40.0
Three   18.75       40.0        35.0   

以下内容适用于单列,但我想对除第一列之外的每一列执行此操作。

gawk -F"\t" '{a[NR]=$1;x+=(b[NR]=$2)}END{while(++i<=NR)print a[i]"\t"100*b[i]/x}' file.txt 

非常感谢您提供的任何帮助。

答案1

输出与您要求的不 100% 相同,但希望足够接近:

function percent(value, total) {
    return sprintf("%.2f", 100 * value / total);
}
{
    label[NR] = $1
    for (i = 2; i <= NF; ++i) {
        sum[i] += col[i][NR] = $i
    }
}
END {
    title = label[1]
    for (i = 2; i <= length(col) + 1; ++i) {
        title = title "\t" col[i][1];
    }
    print title
    for (j = 2; j <= NR; ++j) {
        line = label[j]
        for (i = 2; i <= length(col) + 1; ++i) {
            line = line "\t" percent(col[i][j], sum[i]);
        }
        print line
    }
}

产生输出:

ID    Sample1 Sample2 Sample3
One   62.50   0.00    25.00
Two   18.75   60.00   40.00
Three 18.75   40.00   35.00

运行它gawk -f script.awk file.txt

当然,你可以将脚本扁平化为一行,但我认为最好将其保存在这样的脚本文件中,这样更容易阅读和维护。

一个更简单、更好的版本,也可以与 BSD AWK 一起使用,而不仅仅是 GNU AWK:

function percent(value, total) {
    return sprintf("%.2f", 100 * value / total)
}
BEGIN { OFS = "\t" }
NR == 1 { gsub(/ +/, OFS); print; next }
{
    label[NR] = $1
    for (i = 2; i <= NF; ++i) {
        sum[i] += col[i, NR] = $i
    }
}
END {
    for (j = 2; j <= NR; ++j) {
        $1 = label[j]
        for (i = 2; i <= NF; ++i) {
            $i = percent(col[i, j], sum[i])
        }
        print
    }
}

相关内容