From: martinahansen Date: Wed, 12 Jan 2011 14:02:41 +0000 (+0000) Subject: posished Calc::standard_deviation X-Git-Url: https://git.donarmstrong.com/?a=commitdiff_plain;ds=sidebyside;h=dac27dcade24185a98c6b56f33f308589546a86a;hp=65ff7b75db80def2dd99f7d1c74e42d172d6d543;p=biopieces.git posished Calc::standard_deviation git-svn-id: http://biopieces.googlecode.com/svn/trunk@1201 74ccb610-7750-0410-82ae-013aeee3265d --- diff --git a/code_perl/Maasha/Calc.pm b/code_perl/Maasha/Calc.pm index 5945d4d..6f83fd0 100644 --- a/code_perl/Maasha/Calc.pm +++ b/code_perl/Maasha/Calc.pm @@ -218,24 +218,24 @@ sub standard_deviation # Returns a float. - my ( $mean_num, $num, $div, $div_sum, $mean_div, $std_div ); + my ( $mean_num, $num, $dev, $dev_sum, $mean_dev, $std_dev ); $mean_num = mean( $numbers ); - $div_sum = 0; + $dev_sum = 0; foreach $num ( @{ $numbers } ) { - $div = ( $num - $mean_num ) ** 2; + $dev = ( $num - $mean_num ) ** 2; - $div_sum += $div; + $dev_sum += $dev; } - $mean_div = $div_sum / scalar @{ $numbers }; + $mean_dev = $dev_sum / scalar @{ $numbers }; - $std_div = sqrt( abs( $mean_div ) ); + $std_dev = sqrt( $mean_dev ); - return $std_div; + return $std_dev; }