Perl 脚本中的统计

Question

我有以下问题： 我想创建一个 perl 脚本，从文本文件（包含几列数字的文件）读取并计算一些统计数据（平均值、中位数、标准差、方差）。我已经构建了一个脚本，但是由于我还没有爱上 perl，所以我无法解决它的语法问题...... 这是我的 Perl 脚本..

#!/usr/bin/perl -w

use strict;
open(FILEHANDLE, data.txt);
while (<FILEHANDLE>) {
shift @ARGV;
my @array = split(\t,$_);
}
close(FILEHANDLE);
###### mean, sum and size
$N = $sum = 0;
$array[$x-1]; 
$N++;
$sum += $array[$x-1];
###### minimum and the maximum
($min = 0, $max = 0);
$max = $array[$x-1] if ($max < $array[$x-1]), (my@sorted = sort { $a <=> $b } @samples) {
print join(" ",@sorted);
}
##### median
if ($N % 2==1) {
print "$median = $sorted[int($N/2)]\n"; ## check this out
};
else ($median = ($sorted[$N/2] + $sorted[($N/2)-1]) / 2)) { 
print "$median\n"; # check this out
};
##### quantiles 1º and 3º
if $qt1 = $sorted[$r25-1] {
print "\n"; # check this out
};
else $qt1 = $fr*($sorted[$ir] - $sorted[$ir-1]) + $sorted[$ir-1] {
print "\n"; # check this out
};
##### variance
for (my $i=0; 
$i<scalar(@samples); 
$i++) 
{
$Var += ($samples[$i]-$mean)**2;
$Var = $Var/($N-1);
};
###### standard error
($Std = sqrt($Var)/ sqrt($N));
############################################################
print "$min\n";
print "$max\n";
print "$mean\n";
print "$median\n";
print "$qt1\n";
print "$var\n";
print "$std\n";

exit(0);

我想让它工作。请帮忙。提前致谢！

Answer 1

您的代码中的错误：

open(FILEHANDLE, data.txt);

需要引用 data.txt。您没有检查打开的返回值，例如 ... 或 die $!。您应该使用词法文件句柄和三个参数open，例如open my $fh、'<'、"data.txt"或die $!。

shift @ARGV;

除了从参数列表中删除第一个值（然后立即将其丢弃）之外，这不会执行任何操作。

my @array = split(\t,$_);

您使用 \t 作为裸字，它应该是正则表达式，/\t/。您的 @array 是在 while 循环的词法范围内声明的，并且在此块之外将是未定义的。

$N = $sum = 0;

这两个变量都没有声明，这将导致脚本在使用 strict 时死掉（这是一个非常好的主意）。使用 my $N 来解决这个问题。另外，$N 也不是一个很好的变量名。

$array[$x-1];

这将无济于事。 $x 未声明（见上文），也未定义。整个语句什么也不做，就像有一行3;。我相信您会得到一个错误，例如在 void context 中无用地使用变量。

$N++;

这将 $N 增加到 1，这是一件无用的事情，因为您只用上面的几行将其初始化为 0。

嗯......这样的例子还在继续。我建议你从小处开始，使用严格和警告，因为它们是非常好的工具，并一一解决错误。一个非常好的主意是为您的计算创建子例程，例如：

sub sum {
   # code here
   return $sum;
}

转到 perldoc.perl.org 并阅读文档。特别有用的是与语法相关的和 perlfunc。

另外，您应该知道此功能可以在模块中找到，您可以在 CPAN 中找到这些模块。

Answer 2

您的主要问题是您尚未声明变量，例如 $N、$max 等。

您需要使用 my 引入所有新变量当你第一次引用它们时。就像您对 $array 和 $i 所做的那样。例如

$N = $sum = 0;

应该变成

my( $N, $sum ) = ( 0, 0 );