如何将 150 1536 孔板中的值压缩为 perl 中每个坐标的一个值

Question

我有 150 个 32 x 48 的值矩阵，代表酵母平板图像的颜色强度。我想以每个板的坐标（x，y）为例生成一个平均值，然后从中创建一个新板。现在，我的文件中的所有值均采用以下格式：

• EA_D01_5-8 30,22 -0.397914165526517

• EA_D01_5-8 30,23 -0.326759277147352

• EA_D01_5-8 31,22 -0.172657520010773

• EA_D01_5-8 31,23 -0.103405885199075

第一组代表车牌名称，第二组代表 x,y 坐标，第三组代表 z 分数值。

while($line = <DATAFILE>){                 
     chomp($line); 
     my @temp = split(/\t/, $line);
     $long_name = $temp[0];
     $coords = $temp[1];
     $zscores = $temp[2];
     $stats{$long_name}{$coords}[0] = $zscores;

}

我创建了一个哈希 $stats{$long_name}{$coords}[0] 来存储值 然后我创建了一个循环来尝试访问这些值：

foreach $long_name ( sort keys %stats ) {

        foreach $coords( keys %{$stats{$long_name}} ){

                my $zscoreV = $stats{$long_name}{$coords}[0];
                $totalV = $totalV + $zscoreV;

        }

}

但是我做错了一些事情，因为当我尝试将总数除以 150（即 $long_name 中唯一值的数量）时，我没有剩下 1536 个值。有谁知道我做错了什么？谢谢

Answer 1

也许我误解了你想要的。我认为你想要一个“平均”板，其中 0,0 是所有 150 个板的 0,0 值的平均值，0,1 是 0,1 值的平均值，等等。

如果是这样，你就不需要想要一个长名称的散列，你只需要一个坐标的散列：（

my %coords_total;
my %coords_count;
while(my $line = <DATAFILE>) {
    chomp($line); 
    my ($long_name, $coords, $zscores) = split(/\t/, $line);
    $coords_total{$coords} += $zscores;
    $coords_count{$coords}++;
}

my %coords_average = map $coords_total{$_} / $coords_count{$_}, keys %coords_total;

假设任意数量的车牌；如果你总是想假设 150 个车牌，你可以这样做 $coords_average{$coords} += $zscores / 150; 在初始循环中。）