Perl/펄] 표준 편차 구하기/계산 함수; Standard Deviation Function

#!/usr/bin/perl
use strict; use warnings;

  my @foo = ( 1345, 1301, 1368, 1322, 1310, 1370, 1318, 1350, 1303, 1299 );


  # 평균 구하기
  print mean(@foo), "\n";
  # 출력 결과: 1328.6


  # 표본의 표준편차 구하기 (STDEV)
  print standardDeviation(@foo, 1), "\n";
  # 출력 결과: 27.4639157198435


  # 모집단 전체의 표준편차 구하기 (STDEVP)
  print standardDeviation(@foo, 0), "\n";
  # 출력 결과: 26.0545581424825




sub mean {  # 산술 평균 구하기
  return "NaN" if ($#_ < 0); # 빈 배열일 경우의 에러 처리

  my @array = @_;
  my $sum = 0.0;

  for (my $i = 0; $i <= $#array; $i++) {
    $sum += $array[$i];
  }

  # 소괄호로 둘러싸 주지 않으면
  # 연산자 우선 순위 때문에 버그 생김
  return $sum / ($#array + 1);
}


sub standardDeviation {
  my $option = pop;
  my @array = @_;

  return "NaN" if ($#array + 1 < 2); # 배열 요소 개수가 최소한 2가 되어야 함

  my $sum = 0.0;
  my $sd = 0.0;
  my $diff;
  my $meanValue = mean(@array);

  for (my $i = 0; $i <= $#array; $i++) {
    $diff = $array[$i] - $meanValue;
    $sum += $diff * $diff;
  }
  $sd = sqrt($sum / ($#array + 1 - $option));

  return $sd;
}