求助计算平均值的Perl程序

liujie126

我有一个文件，其内容是：1 1 1 -1 1 1 1 1 -1 1 1 -1 -1 1 -1 1 1 -1 -1 1 1 1 1 1 1 -1 1 1 1 1 1 -1 -1 1 1 -1 1 1 1 -1 1 1 1 1 -1 -1 1 -1 -1 -1 1 -1 1 -1 -1 1 1 1 -1 -1 1 1 1 1 1 1 1 -1 1 -1 1 -1 1 1 1 1 1 1 -1 -1 1 1 1 1 -1 1 -1 -1 -1 1 1 1 1 1 1 -1 1 1 1 1 1 -1 1 1 1 1 1 1 1 -1 -1 1 1 1 -1 -1 1 1 -1 -1 -1 1 1 -1 1 1 1 1 -1 1 -1 1 -1 1 1 -1 -1 1 1 1 1 1 1 1 每个数字之间有空格隔开，现要逐级的计算每4个数的平均值，即第1--4位这4个数平均值，返回数值；第2--5位的4个数返回一个平均值；第3--6位这4个数返回一个平均值。我的数据较大，由于本人Perl编程水平一般，希望高手赐教，衷心谢谢！！

xiumu2280

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. while (<DATA>) {
   
6.         chomp;
   
7.         my @data = split /\s+/;
   
8.         for (my $m=0;$m<=@data-4;$m++) {
   
9.                 my $seq = $data[$m]+$data[$m+1]+$data[$m+2]+$data[$m+3];
   
10.                 my $aver = $seq/4;
    
11.                 print "$m\t$seq\t$aver\n";
    
12.         }
    
13. }
    
14. 
    
15. 
    
16. 
    
17. __DATA__
    
18. 1 1 1 1 -1 1 1 1 1 -1 1

复制代码

这和kmer这个概念很像···
程序很简陋，不过能用。

yestreenstars

1. #!/usr/bin/perl
   
2. use 5.010;
   
3. 
   
4. my @data = split /\s+/, <DATA>;
   
5. my $n = 0;
   
6. while ( $#data >= 4 ) {
   
7.         say ++$n, '--', $n+3, ' ', ($data[0]+$data[1]+$data[2]+$data[3])/4;
   
8.         shift @data;
   
9. }
   
10. 
    
11. __DATA__
    
12. 1 1 1 -1 1 1 1 1 -1 1 1 -1 -1 1 -1 1 1 -1 -1 1 1 1 1 1 1 -1 1 1 1 1 1 -1 -1 1 1 -1 1 1 1 -1 1 1 1 1 -1 -1 1 -1 -1 -1 1 -1 1 -1 -1 1 1 1 -1 -1 1 1 1 1 1 1 1 -1 1 -1 1 -1 1 1 1 1 1 1 -1 -1 1 1 1 1 -1 1 -1 -1 -1 1 1 1 1 1 1 -1 1 1 1 1 1 -1 1 1 1 1 1 1 1 -1 -1 1 1 1 -1 -1 1 1 -1 -1 -1 1 1 -1 1 1 1 1 -1 1 -1 1 -1 1 1 -1 -1 1 1 1 1 1 1 1

复制代码

sars318

这个更像perl程序.!回复 3# yestreenstars


   

pitonas

小伙伴们, 这个高端大气上档次。{:2_179:}

1. #!/usr/bin/perl -w
   
2. my $s = <DATA>;
   
3. print +( 4 - tr/-// * 2 ) / 4 . ' ' for $s =~ /(?:\S+\s*){4}/g;
   
4. 
   
5. __DATA__
   
6. 1 1 1 -1 1 1 1 1 -1 1 1 -1 -1 1 -1 1 1 -1 -1 1 1 1 1 1 1 -1 1 1 1 1 1 -1 -1 1 1 -1 1 1 1 -1 1 1 1 1 -1 -1 1 -1 -1 -1 1 -1 1 -1 -1 1 1 1 -1 -1 1 1 1 1 1 1 1 -1 1 -1 1 -1 1 1 1 1 1 1 -1 -1 1 1 1 1 -1 1 -1 -1 -1 1 1 1 1 1 1 -1 1 1 1 1 1 -1 1 1 1 1 1 1 1 -1 -1 1 1 1 -1 -1 1 1 -1 -1 -1 1 1 -1 1 1 1 1 -1 1 -1 1 -1 1 1 -1 -1 1 1 1 1 1 1 1
   

复制代码

yestreenstars

回复 5# pitonas
小伙伴你算出来的好像有问题~

   

pitonas

{:2_177:} 大神
0.5 1 0 0 0 1 0.5 0.5 0 0.5 1 -0.5 -0.5 0 0 1 0.5 0 1 0 1 -0.5 0.5 0.5 1 0.5 1 0 0 0 0 1 0 0 0.5 1

1 1 1 -1   0.5
1 1 1 1    1
-1 1 1 -1  0
-1 1 -1 1  0
1 -1 -1 1  0
1 1 1 1    1
.....
回复 6# yestreenstars


   

pitonas

没注意到  1--4, 2--5 {:2_172:} 是有问题

回复 6# yestreenstars


   

pitonas

小伙伴们, 这个正确了 {:2_172:}

1. #!/usr/bin/perl -w
   
2. my @s = map split, <DATA>;
   
3. print +( 4 - "@s[$_..$_+3]" =~ tr/-// * 2 ) / 4, " " for 0 .. @s - 4;
   
4. 
   
5. __DATA__
   
6. 1 1 1 -1 1 1 1 1 -1 1 1 -1
   

复制代码

jason680

回复 1# liujie126

The other way to assign average number by call

use strict;
use warnings;

sub avg{
  my($sAvg, @aData) = @_;
  my @aRet;
  my $sTotal = 0;

  return if(@aData < $sAvg);   # check it

  $sTotal += $aData[$_]  foreach(0 .. $sAvg-2);   # precalculate the total without last one [0,1,2]
  foreach(0 .. @aData-$sAvg){
    $sTotal += $aData[$_ + $sAvg-1];                  # summary with last one [0,1,2,3]
    push @aRet, $sTotal/$sAvg;                            # to get average
    $sTotal -= $aData[$_];                                   # remove the first one [1,2,3]
  }
  return @aRet;
}
while(<DATA>){
  chomp;
  my @aData = split;
  my @aAvg = avg(4,@aData);
  print "@aAvg\n";
}

__DATA__
1 1 1 -1 1 1 1 1 -1 1 1 -1 -1 1 -1 1 1 -1 -1 1 1 1 1 1 1 -1 1 1 1 1 1 -1 -1 1 1 -1 1 1 1 -1 1 1 1 1 -1 -1 1 -1 -1 -1 1 -1 1 -1 -1 1 1 1 -1 -1 1 1 1 1 1 1 1 -1 1 -1 1 -1 1 1 1 1 1 1 -1 -1 1 1 1 1 -1 1 -1 -1 -1 1 1 1 1 1 1 -1 1 1 1 1 1 -1 1 1 1 1 1 1 1 -1 -1 1 1 1 -1 -1 1 1 -1 -1 -1 1 1 -1 1 1 1 1 -1 1 -1 1 -1 1 1 -1 -1 1 1 1 1 1 1 1
1 1 1 1 -1 -1 -1 -1