t/Statistics.t

   1 package Math::GSL::Statistics::Test;
   2 use Math::GSL::Test qw/:all/;
   3 use base q{Test::Class};
   4 use Test::More;
   5 use Math::GSL::Statistics qw/:all/;
   6 use Math::GSL::Sort qw/:all/;
   7 use Math::GSL qw/:all/;
   8 use strict;
   9
  10 sub make_fixture : Test(setup) {
  11     my $self = shift;
  12     $self->{data} =   [
  13         17, 18, 16, 18, 12,
  14         20, 18, 20, 20, 22,
  15         20, 10, 8, 12, 16,
  16         16, 18, 20, 18, 21
  17     ];
  18     $self->{datab} = [
  19         19, 20, 22, 24, 10,
  20         25, 20, 22, 21, 23,
  21         20, 10, 12, 14, 12,
  22         20, 22, 24, 23, 17
  23     ];
  24 }
  25
  26 sub teardown : Test(teardown) {
  27 }
  28
  29 sub GSL_STATS_MEAN : Tests {
  30     ok_similar( gsl_stats_mean([2 .. 5], 1, 4), 3.5 );
  31 }
  32
  33 sub GSL_STATS_VARIANCE : Tests {
  34     ok_similar( gsl_stats_variance([2 .. 5], 1, 4), 5/3);
  35 }
  36
  37 sub GSL_STATS_SD : Tests {
  38  my $self = shift;
  39  my $x = gsl_stats_sd($self->{data}, 1, 20);
  40  ok_similar($x, 3.79750610685209);
  41 }
  42
  43 sub GSL_STATS_ABSDEV : Tests {
  44  my $self = shift;
  45  my $x = gsl_stats_absdev($self->{data},1, 20);
  46  ok_similar($x, 2.9);
  47 }
  48
  49 sub GSL_STATS_SKEW : Tests {
  50  my $self = shift;
  51  my $x = gsl_stats_skew($self->{data}, 1, 20);
  52  ok_similar($x, -0.909355923168064);
  53 }
  54
  55 sub GSL_STATS_KURTOSIS : Tests {
  56  my $self = shift;
  57  my $x = gsl_stats_kurtosis($self->{data}, 1, 20);
  58  ok_similar_relative($x, -0.233692524908094);
  59 }
  60
  61 sub GSL_STATS_COVARIANCE : Tests {
  62  my $self = shift;
  63  my $x = gsl_stats_covariance($self->{data}, 1, $self->{datab}, 1, 20);
  64  ok_similar_relative($x, 14.5263157894737);
  65 }
  66
  67 sub GSL_STATS_CORRELATION : Tests {
  68  my $self = shift;
  69  my $x = gsl_stats_correlation($self->{data}, 1, $self->{datab}, 1, 20);
  70  ok_similar_relative($x, 0.793090350710101);
  71 }
  72
  73 sub GSL_STATS_PVARIANCE : Tests {
  74  my $self = shift;
  75  my $x = gsl_stats_pvariance($self->{data}, 1, 20, $self->{datab}, 1, 20);
  76  ok_similar_relative($x, 18.8421052631579);
  77 }
  78
  79 sub GSL_STATS_TTEST : Tests {
  80  my $self = shift;
  81  my $x = gsl_stats_ttest($self->{data}, 1, 20, $self->{datab}, 1, 20);
  82  ok_similar_relative($x, -1.45701922702927);
  83 }
  84
  85 sub GSL_STATS_MAX : Tests {
  86  my $self = shift;
  87  my $x = gsl_stats_max($self->{data}, 1, 20);
  88  ok_similar_relative($x, 22);
  89 }
  90
  91 sub GSL_STATS_MIN : Tests {
  92  my $self = shift;
  93  my $x = gsl_stats_min($self->{data}, 1, 20);
  94  ok_similar_relative($x, 8);
  95 }
  96
  97 sub GSL_STATS_MINMAX : Tests {
  98  my $self = shift;
  99  my ($min, $max) = gsl_stats_minmax($self->{data}, 1, 20);
 100  ok_similar_relative($max, 22);
 101  ok_similar_relative($min, 8);
 102 }
 103
 104 sub GSL_STATS_MAX_INDEX : Tests {
 105  my $self = shift;
 106  my $x = gsl_stats_max_index($self->{data}, 1, 20);
 107  ok_similar_relative($x, 9);
 108 }
 109
 110 sub GSL_STATS_MIN_INDEX : Tests {
 111  my $self = shift;
 112  my $x = gsl_stats_min_index($self->{data}, 1, 20);
 113  ok_similar_relative($x, 12);
 114 }
 115
 116 sub GSL_STATS_MINMAX_INDEX : Tests {
 117  my $self = shift;
 118  my ($min, $max) = gsl_stats_minmax_index($self->{data}, 1, 20);
 119  ok_similar_relative($min, 12);
 120  ok_similar_relative($max, 9);
 121 }
 122
 123 sub GSL_STATS_MEDIAN_FROM_SORTED_DATA_EVEN : Tests {
 124  my $self = shift;
 125  my $sorted = gsl_sort($self->{data}, 1, 20);
 126  my $x = gsl_stats_median_from_sorted_data($sorted, 1, 20);
 127  ok_similar_relative($x, 18);
 128 }
 129
 130 sub GSL_STATS_MEDIAN_FROM_SORTED_DATA_ODD : Tests {
 131  my $self = shift;
 132  my $sorted = gsl_sort($self->{data}, 1, 20);
 133  my $x = gsl_stats_median_from_sorted_data($sorted, 1, 19);
 134  ok_similar_relative($x, 18);
 135 }
 136
 137 sub GSL_STATS_QUANTILE_FROM_SORTED_DATA_0 : Tests {
 138  my $self = shift;
 139  my $sorted = gsl_sort($self->{data}, 1, 20);
 140  my $x = gsl_stats_quantile_from_sorted_data($sorted, 1, 20, 0);
 141  ok_similar_relative($x, 8);
 142 }
 143
 144 sub GSL_STATS_QUANTILE_FROM_SORTED_DATA_100 : Tests {
 145  my $self = shift;
 146  my $sorted = gsl_sort($self->{data}, 1, 20);
 147  my $x = gsl_stats_quantile_from_sorted_data($sorted, 1, 20, 1);
 148  ok_similar_relative($x, 22);
 149 }
 150
 151 sub GSL_STATS_QUANTILE_FROM_SORTED_DATA_50_EVEN : Tests {
 152  my $self = shift;
 153  my $sorted = gsl_sort($self->{data}, 1, 20);
 154  my $x = gsl_stats_quantile_from_sorted_data($sorted, 1, 20, 0.5);
 155  ok_similar_relative($x, 18);
 156 }
 157
 158 sub GSL_STATS_QUANTILE_FROM_SORTED_DATA_50_ODD : Tests {
 159  my $self = shift;
 160  my $sorted = gsl_sort($self->{data}, 1, 20);
 161  my $x = gsl_stats_quantile_from_sorted_data($sorted, 1, 20, 0.5);
 162  ok_similar_relative($x, 18);
 163 }
 164
 165 sub GSL_STATS_LAG1_AUTOCORRELATION : Tests {
 166   my $nacc3 = 1001 ;
 167   my $i;
 168   my $numacc3;
 169   $numacc3->[0] = 1000000.2 ;
 170   for ($i = 1 ; $i < 1000  ; $i += 2)
 171     {
 172       $numacc3->[$i] = 1000000.1 ;
 173       $numacc3->[$i+1] = 1000000.3 ; }
 174   my $mean = gsl_stats_mean ($numacc3, 1, $nacc3);
 175   my $sd = gsl_stats_sd ($numacc3, 1, $nacc3);
 176   my $lag1 = gsl_stats_lag1_autocorrelation ($numacc3, 1, $nacc3);
 177
 178   my $expected_mean = 1000000.2;
 179   my $expected_sd = 0.1;
 180   my $expected_lag1 = -0.999;
 181
 182   ok_similar_relative($mean, $expected_mean);
 183   ok_similar_relative($sd, $expected_sd);
 184   ok_similar_relative($lag1, $expected_lag1);
 185 }
 186 Test::Class->runtests;