perl - Perl:从多列文件的每一列计算折刀误差
问题描述
我正在尝试计算多列文件中每一列的折刀平均值和误差。我的示例数据文件如下所示:
$ cat data.HW2
1.1 2.1 3.1 4.1
1.2 2.2 3.2 4.2
1.3 2.3 3.3 4.3
1.4 2.4 3.4 4.4
我尝试的解决方案是定义最终大小与列数相同的数组(在本例中为 4)并逐行迭代它们:
cat jackkinfe.pl
#! /usr/bin/perl
use warnings; use strict;
my @n=0;
my @x;
my $j;
my $i;
my $dg;
my @x_jack;
my @x_tot=0;
my $cols;
my $col_start=0;
# read in the data
while(<>)
{
my @column = split();
$cols=@column;
foreach my $j ($col_start .. $#column) {
$x[$n[$j]][$j] = $column[$j];
$x_tot[$j] += $x[$n[$j]][$j];
$n[$j]++;
}
}
# Do the jackknife estimates
for ($j=$col_start; $j<$cols; $j++)
{
for ($i = 0; $i < $n[$j]; $i++)
{
$x_jack[$i][$j] = ($x_tot[$j] - $x[$i][$j]) / ($n[$j] - 1);
}
# Do the final jackknife estimate
my @g_jack_av=0;
my @g_jack_err=0;
for ($i = 0; $i < $n[$j]; $i++)
{
$dg = $x_jack[$i][$j];
$g_jack_av[$j] += $dg;
$g_jack_err[$j] += $dg**2;
}
$g_jack_av[$j] /= $n[$j];
$g_jack_err[$j] /= $n[$j];
$g_jack_err[$j] = sqrt(($n[$j] - 1) * abs($g_jack_err[$j] - $g_jack_av[$j]**2));
printf "%e %e ", $g_jack_av[$j], $g_jack_err[$j];
}
printf "\n";
它给了我以下两个警告:
$cat data.HW2 | perl jackknife.pl
Use of uninitialized value within @n in array element at cols_jacknife.pl line 19, <> line 1.
Use of uninitialized value within @n in array element at cols_jacknife.pl line 20, <> line 1.
它抱怨以下两行:
$x[$n[$j]][$j] = $column[$j];
$x_tot[$j] += $x[$n[$j]][$j];
但我想@n
根据数据文件的大小动态设置大小。
如何删除此警告?
由于我正在努力学习最佳实践,因此也欢迎并非常感谢有关我的 Perl 使用的任何其他建议。
解决方案
这部分代码
my @n=0;
....
foreach my $j ($col_start .. $#column) {
$x[$n[$j]][$j] = $column[$j];
$x_tot[$j] += $x[$n[$j]][$j];
$n[$j]++;
}
$j
每个大于 0的值都会触发一次警告,因为只@n
定义了第一个元素:$n[0] = 0
. 只有在循环迭代结束时,数组值才最终定义,当它被1
增量运算符设置为时$n[$j]++
。
从技术上讲,代码仍将按预期工作,因为undef
将被强制转换为0
. 所以....忽略警告应该是安全的。你可以在你的循环中做这样的事情来避免它:
$n[$j] //= 0; # $n[$j] is defined, or set to 0
这相当于
if (not defined($n[$j])) {
$n[$j] = 0;
}
推荐阅读
- excel - 检查一个范围与其他范围的精确匹配
- javascript - AssertionError [ERR_ASSERTION]: PORT 是必需的(将节点 js 作为 Windows 服务启动时)
- python - 如果当前列包含字符串,则检查其他列的值
- r - R中每第n列更改列名
- android - 任务的特定调度程序线程每个任务都应该有它自己的调度程序线程
- sql - 比较 Oracle SQL 中属性的所有条目
- javascript - 有没有办法使用 Moment.js 计算给定年份的 ISO 周数?
- c# - 更新/删除记录时如何处理多对一关系
- flutter - Flutter - 在构建器内部的开关内调用导航器
- bash - 在file2中找到离file1最近的点,shell skript