File Coverage

blib/lib/Math/Random/GaussianRange.pm

Criterion	Covered	Total	%
statement	39	39	100.0
branch	8	8	100.0
condition	10	13	76.9
subroutine	7	7	100.0
pod	1	1	100.0
total	65	68	95.5

use strict;
use warnings;

package Math::Random::GaussianRange;

#ABSTRACT: Given a range, generate a set of random numbers normally distributed.

use Math::Random qw(random_normal);
use Statistics::Basic qw(median);
use Carp;
use base 'Exporter';

our @EXPORT = 'generate_normal_range';


sub generate_normal_range {
    my $rh    = shift;

    my $min   = $rh->{min};
    my $max   = $rh->{max};
    my $n     = $rh->{n} || 100;
    my $round = $rh->{round} || 0;

    unless ( defined $min && defined $max ) {
        croak "Specify a range using the min and max parameters.";
    }

    if ( $min > $max ) {
        croak "The minimum cannot exceeed the maximum";
    }

    my @range  = ( $min .. $max );
    my $median = median( @range );
    my $sd     = $median/3;

    if ( $median == 0 || $sd == 0 ) {
        croak "Median and SD are both null.";
    }

    my @output =
        random_normal( $n, $median, $sd );

    if ( defined $round && $round > 0) {

        my $ra_rounded = [ ];
        foreach my $i ( @output ) {
            push( @$ra_rounded, sprintf('%.0f', $i) );
        }

        return $ra_rounded;

    }

    return \@output;

}


1;

__END__