ruby-statistics 2.1.2 → 2.1.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: f02b2824672be63fdfac39ad23d7ef4c6c60c9832ac1cdeae665e5e8ccefb5f3
-  data.tar.gz: 7dd359dc331405961c1b5739a70c19831c8c12908bd09152e8d30ef6f3a42c08
+  metadata.gz: 6612502f03d8077d0158d997a42dfbc4d1002f2ab01ce2b7bdb5fbd510187e3e
+  data.tar.gz: 14fb04073b5b788dfa9e93aa586daef050dd105c2d2f8bdd17db30ad1fbcf144
 SHA512:
-  metadata.gz: 259fc76594f462aa63770351f8851cdd694d3a045e7d31d7a58a7a578c1fd6d805870b8bec8f2a7f1f2df4995537ab7eef2efb76542aa77e0371311e99e81680
-  data.tar.gz: a06d8899bcfcbc05b8ad1da6c09695d0247c11756ea26087334b229c182f5c11831b3692ccad30c75ef3dc9fc8a2699269ff7aaba22224a2d03d1280363ff1a9
+  metadata.gz: '09590f836a59563819a1a847830e5dc2ee3554415cadc81c35b2a0f43ab1af87204f028659e8aa2f30a14b58c69c3e4f65db5e722d0a00ced5d92faa1e7dce82'
+  data.tar.gz: 2e66a26c23bf1f05cb9de40e992b302c4f0fef13aa70b4e509de479cb15b9700d4032f5d548aa45110f161ef9dac417f9b1872479a02dca0e729a051be2a4fc8

data/.github/dependabot.yml

@@ -0,0 +1,15 @@
+# To get started with Dependabot version updates, you'll need to specify which
+# package ecosystems to update and where the package manifests are located.
+# Please see the documentation for all configuration options:
+# https://help.github.com/github/administering-a-repository/configuration-options-for-dependency-updates
+
+version: 2
+updates:
+  - package-ecosystem: "bundler" # See documentation for possible values
+    directory: "/" # Location of package manifests
+    schedule:
+      interval: "weekly"
+  - package-ecosystem: "github-actions" # See documentation for possible values
+    directory: "/" # Location of package manifests
+    schedule:
+      interval: "weekly"

data/.github/workflows/ruby.yml

@@ -8,9 +8,9 @@ jobs:
     runs-on: ubuntu-latest
 
     steps:
-    - uses: actions/checkout@v1
+    - uses: actions/checkout@v2.3.4
     - name: Set up Ruby 2.6
-      uses: actions/setup-ruby@v1
+      uses: actions/setup-ruby@v1.1.2
       with:
         ruby-version: 2.6.x
     - name: Build and test with Rake
@@ -23,9 +23,9 @@ jobs:
     runs-on: ubuntu-latest
 
     steps:
-    - uses: actions/checkout@v1
+    - uses: actions/checkout@v2.3.4
     - name: Set up Ruby 2.7
-      uses: actions/setup-ruby@v1
+      uses: actions/setup-ruby@v1.1.2
       with:
         ruby-version: 2.7.x
     - name: Build and test with Rake

data/lib/math.rb

@@ -9,11 +9,11 @@ module Math
   end
 
   def self.combination(n, r)
-    self.factorial(n)/(self.factorial(r) * self.factorial(n - r)).to_f # n!/(r! * [n - r]!)
+    self.factorial(n)/(self.factorial(r) * self.factorial(n - r)).to_r # n!/(r! * [n - r]!)
   end
 
   def self.permutation(n, k)
-    self.factorial(n)/self.factorial(n - k).to_f
+    self.factorial(n)/self.factorial(n - k).to_r
   end
 
   # Function adapted from the python implementation that exists in https://en.wikipedia.org/wiki/Simpson%27s_rule#Sample_implementation
@@ -24,7 +24,8 @@ module Math
       return
     end
 
-    h = (b - a)/n.to_f
+    h = (b - a)/n.to_r
+
     resA = yield(a)
     resB = yield(b)
 
@@ -45,7 +46,7 @@ module Math
 
   def self.lower_incomplete_gamma_function(s, x)
     # The greater the iterations, the better. That's why we are iterating 10_000 * x times
-    self.simpson_rule(0, x, (10_000 * x.round).round) do |t|
+    self.simpson_rule(0, x.to_r, (10_000 * x.round).round) do |t|
       (t ** (s - 1)) * Math.exp(-t)
     end
   end
@@ -72,7 +73,7 @@ module Math
     # To avoid overflow problems, the implementation applies the logarithm properties
     # to calculate in a faster and safer way the values.
     lbet_ab = (Math.lgamma(alp)[0] + Math.lgamma(bet)[0] - Math.lgamma(alp + bet)[0]).freeze
-    front = (Math.exp(Math.log(x) * alp + Math.log(1.0 - x) * bet - lbet_ab) / alp.to_f).freeze
+    front = (Math.exp(Math.log(x) * alp + Math.log(1.0 - x) * bet - lbet_ab) / alp.to_r).freeze
 
     # This is the non-log version of the left part of the formula (before the continuous fraction)
     # down_left = alp * self.beta_function(alp, bet)

data/lib/statistics/distribution/bernoulli.rb

@@ -24,7 +24,7 @@ module Statistics
       end
 
       def self.skewness(p)
-        (1.0 - 2.0*p).to_f / Math.sqrt(p * (1.0 - p))
+        (1.0 - 2.0*p).to_r / Math.sqrt(p * (1.0 - p))
       end
 
       def self.kurtosis(p)

data/lib/statistics/distribution/beta.rb

@@ -4,8 +4,8 @@ module Statistics
       attr_accessor :alpha, :beta
 
       def initialize(alp, bet)
-        self.alpha = alp.to_f
-        self.beta = bet.to_f
+        self.alpha = alp.to_r
+        self.beta = bet.to_r
       end
 
       def cumulative_function(value)

data/lib/statistics/distribution/empirical.rb

@@ -19,7 +19,7 @@ module Statistics
           summation
         end
 
-        cumulative_sum / samples.size.to_f
+        cumulative_sum / samples.size.to_r
       end
     end
   end

data/lib/statistics/distribution/f.rb

@@ -10,7 +10,7 @@ module Statistics
 
       # Formula extracted from http://www.itl.nist.gov/div898/handbook/eda/section3/eda3665.htm#CDF
       def cumulative_function(value)
-        k = d2/(d2 + d1 * value.to_f)
+        k = d2/(d2 + d1 * value.to_r)
 
         1 - Math.incomplete_beta_function(k, d2/2.0, d1/2.0)
       end
@@ -18,28 +18,28 @@ module Statistics
       def density_function(value)
         return if d1 < 0 || d2 < 0 # F-pdf is well defined for the [0, +infinity) interval.
 
-        val = value.to_f
+        val = value.to_r
         upper = ((d1 * val) ** d1) * (d2**d2)
         lower = (d1 * val + d2) ** (d1 + d2)
-        up = Math.sqrt(upper/lower.to_f)
+        up = Math.sqrt(upper/lower.to_r)
         down = val * Math.beta_function(d1/2.0, d2/2.0)
 
-        up/down.to_f
+        up/down.to_r
       end
 
       def mean
         return if d2 <= 2
 
-        d2/(d2 - 2).to_f
+        d2/(d2 - 2).to_r
       end
 
       def mode
         return if d1 <= 2
 
-        left = (d1 - 2)/d1.to_f
-        right = d2/(d2 + 2).to_f
+        left = (d1 - 2)/d1.to_r
+        right = d2/(d2 + 2).to_r
 
-        left * right
+        (left * right).to_f
       end
     end
   end

data/lib/statistics/distribution/geometric.rb

@@ -4,7 +4,7 @@ module Statistics
       attr_accessor :probability_of_success, :always_success_allowed
 
       def initialize(p, always_success: false)
-        self.probability_of_success = p.to_f
+        self.probability_of_success = p.to_r
         self.always_success_allowed = always_success
       end
 

data/lib/statistics/distribution/logseries.rb

@@ -6,7 +6,7 @@ module Statistics
         k = k.to_i
 
         left = (-1.0 / Math.log(1.0 - p))
-        right = (p ** k).to_f
+        right = (p ** k).to_r
 
         left * right / k
       end
@@ -44,7 +44,7 @@ module Statistics
         up = p + Math.log(1.0 - p)
         down = ((1.0 - p) ** 2) * (Math.log(1.0 - p) ** 2)
 
-        (-1.0 * p) * (up / down.to_f)
+        (-1.0 * p) * (up / down.to_r)
       end
     end
   end

data/lib/statistics/distribution/negative_binomial.rb

@@ -25,21 +25,21 @@ module Statistics
       end
 
       def mean
-        (probability_per_trial * number_of_failures)/(1 - probability_per_trial).to_f
+        (probability_per_trial * number_of_failures)/(1 - probability_per_trial).to_r
       end
 
       def variance
-        (probability_per_trial * number_of_failures)/((1 - probability_per_trial) ** 2).to_f
+        (probability_per_trial * number_of_failures)/((1 - probability_per_trial) ** 2).to_r
       end
 
       def skewness
-        (1 + probability_per_trial).to_f / Math.sqrt(probability_per_trial * number_of_failures)
+        (1 + probability_per_trial).to_r / Math.sqrt(probability_per_trial * number_of_failures)
       end
 
       def mode
         if number_of_failures > 1
           up = probability_per_trial * (number_of_failures - 1)
-          down = (1 - probability_per_trial).to_f
+          down = (1 - probability_per_trial).to_r
 
           (up/down).floor
         elsif number_of_failures <= 1

data/lib/statistics/distribution/normal.rb

@@ -5,9 +5,9 @@ module Statistics
       alias_method :mode, :mean
 
       def initialize(avg, std)
-        self.mean = avg.to_f
-        self.standard_deviation = std.to_f
-        self.variance = std.to_f**2
+        self.mean = avg.to_r
+        self.standard_deviation = std.to_r
+        self.variance = std.to_r**2
       end
 
       def cumulative_function(value)

data/lib/statistics/distribution/poisson.rb

@@ -18,7 +18,7 @@ module Statistics
         upper = (expected_number_of_occurrences ** k) * Math.exp(-expected_number_of_occurrences)
         lower = Math.factorial(k)
 
-        upper/lower.to_f
+        upper/lower.to_r
       end
 
       def cumulative_function(k)
@@ -31,7 +31,7 @@ module Statistics
 
         # We need the right tail, i.e.: The upper incomplete gamma function. This can be
         # achieved by doing a substraction between 1 and the lower incomplete gamma function.
-        1 - (upper/lower.to_f)
+        1 - (upper/lower.to_r)
       end
     end
   end

data/lib/statistics/distribution/t_student.rb

@@ -29,7 +29,7 @@ module Statistics
         upper = Math.gamma((degrees_of_freedom + 1)/2.0)
         lower = Math.sqrt(degrees_of_freedom * Math::PI) * Math.gamma(degrees_of_freedom/2.0)
         left = upper/lower
-        right = (1 + ((value ** 2)/degrees_of_freedom.to_f)) ** -((degrees_of_freedom + 1)/2.0)
+        right = (1 + ((value ** 2)/degrees_of_freedom.to_r)) ** -((degrees_of_freedom + 1)/2.0)
 
         left * right
       end
@@ -64,8 +64,8 @@ module Statistics
           results << Math.simpson_rule(threshold, y, 10_000) do |t|
             up = Math.gamma((v+1)/2.0)
             down = Math.sqrt(Math::PI * v) * Math.gamma(v/2.0)
-            right = (1 + ((y ** 2)/v.to_f)) ** ((v+1)/2.0)
-            left = up/down.to_f
+            right = (1 + ((y ** 2)/v.to_r)) ** ((v+1)/2.0)
+            left = up/down.to_r
 
             left * right
           end

data/lib/statistics/distribution/uniform.rb

@@ -4,8 +4,8 @@ module Statistics
       attr_accessor :left, :right
 
       def initialize(a, b)
-        self.left = a.to_f
-        self.right = b.to_f
+        self.left = a.to_r
+        self.right = b.to_r
       end
 
       def density_function(value)

data/lib/statistics/distribution/weibull.rb

@@ -4,8 +4,8 @@ module Statistics
       attr_accessor :shape, :scale # k and lambda
 
       def initialize(k, lamb)
-        self.shape = k.to_f
-        self.scale = lamb.to_f
+        self.shape = k.to_r
+        self.scale = lamb.to_r
       end
 
       def cumulative_function(random_value)

data/lib/statistics/spearman_rank_coefficient.rb

@@ -14,7 +14,7 @@ module Statistics
         if rankings.fetch(value, false)
           rankings[value][:rank] += (temporal_ranking + rankings[value][:counter])
           rankings[value][:counter] += 1
-          rankings[value][:tie_rank] = rankings[value][:rank] / rankings[value][:counter].to_f
+          rankings[value][:tie_rank] = rankings[value][:rank] / rankings[value][:counter].to_r
         else
           rankings[value] = { counter: 1, rank: temporal_ranking, tie_rank: temporal_ranking }
         end
@@ -35,7 +35,7 @@ module Statistics
       return if set_one.size == 0 && set_two.size == 0
 
       set_one_mean, set_two_mean = set_one.mean, set_two.mean
-      have_tie_ranks = (set_one + set_two).any? { |rank| rank.is_a?(Float) }
+      have_tie_ranks = (set_one + set_two).any? { |rank| rank.is_a?(Float) || rank.is_a?(Rational) }
 
       if have_tie_ranks
         numerator = 0
@@ -54,7 +54,7 @@ module Statistics
 
         denominator = Math.sqrt(squared_differences_set_one * squared_differences_set_two)
 
-        numerator / denominator.to_f # This is rho or spearman's coefficient.
+        numerator / denominator.to_r # This is rho or spearman's coefficient.
       else
         sum_squared_differences = set_one.each_with_index.reduce(0) do |memo, (rank_one, index)|
           memo += ((rank_one - set_two[index]) ** 2)
@@ -64,7 +64,7 @@ module Statistics
         numerator = 6 * sum_squared_differences
         denominator = ((set_one.size ** 3) - set_one.size)
 
-        1.0 - (numerator / denominator.to_f) # This is rho or spearman's coefficient.
+        1.0 - (numerator / denominator.to_r) # This is rho or spearman's coefficient.
       end
     end
   end

data/lib/statistics/statistical_test/chi_squared_test.rb

@@ -8,12 +8,12 @@ module Statistics
         statistic = if expected.is_a? Numeric
                       observed.reduce(0) do |memo, observed_value|
                         up = (observed_value - expected) ** 2
-                        memo += (up/expected.to_f)
+                        memo += (up/expected.to_r)
                       end
                     else
                       expected.each_with_index.reduce(0) do |memo, (expected_value, index)|
                         up = (observed[index] - expected_value) ** 2
-                        memo += (up/expected_value.to_f)
+                        memo += (up/expected_value.to_r)
                       end
                     end
 

data/lib/statistics/statistical_test/f_test.rb

@@ -19,7 +19,7 @@ module Statistics
         if args.size == 2
           variances = [args[0].variance, args[1].variance]
 
-          f_score = variances.max/variances.min.to_f
+          f_score = variances.max/variances.min.to_r
           df1 = 1 # k-1 (k = 2)
           df2 = args.flatten.size - 2 # N-k (k = 2)
         elsif args.size > 2
@@ -37,18 +37,18 @@ module Statistics
           variance_between_groups = iterator.reduce(0) do |summation, (size, index)|
             inner_calculation = size * ((sample_means[index] - overall_mean) ** 2)
 
-            summation += (inner_calculation / (total_groups - 1).to_f)
+            summation += (inner_calculation / (total_groups - 1).to_r)
           end
 
           # Variance within groups
           variance_within_groups = (0...total_groups).reduce(0) do |outer_summation, group_index|
             outer_summation += args[group_index].reduce(0) do |inner_sumation, observation|
               inner_calculation = ((observation - sample_means[group_index]) ** 2)
-              inner_sumation += (inner_calculation / (total_elements - total_groups).to_f)
+              inner_sumation += (inner_calculation / (total_elements - total_groups).to_r)
             end
           end
 
-          f_score = variance_between_groups/variance_within_groups.to_f
+          f_score = variance_between_groups/variance_within_groups.to_r
           df1 = total_groups - 1
           df2 = total_elements - total_groups
         end

data/lib/statistics/statistical_test/kolmogorov_smirnov_test.rb

@@ -17,7 +17,7 @@ module Statistics
 
         # TODO: Validate calculation of Common alpha.
         common_alpha = Math.sqrt((-0.5 * Math.log(alpha)))
-        radicand = (group_one.size + group_two.size) / (group_one.size * group_two.size).to_f
+        radicand = (group_one.size + group_two.size) / (group_one.size * group_two.size).to_r
 
         critical_d = common_alpha * Math.sqrt(radicand)
         # critical_d = self.critical_d(alpha: alpha, n: samples.size)

data/lib/statistics/statistical_test/t_test.rb

@@ -23,7 +23,7 @@ module Statistics
                     comparison_mean = args[0]
                     degrees_of_freedom = args[1].size - 1
 
-                    (data_mean - comparison_mean)/(data_std / Math.sqrt(args[1].size).to_f).to_f
+                    (data_mean - comparison_mean)/(data_std / Math.sqrt(args[1].size).to_r).to_r
                   else
                     sample_left_mean = args[0].mean
                     sample_left_variance = args[0].variance
@@ -31,12 +31,12 @@ module Statistics
                     sample_right_mean = args[1].mean
                     degrees_of_freedom = args.flatten.size - 2
 
-                    left_root = sample_left_variance/args[0].size.to_f
-                    right_root = sample_right_variance/args[1].size.to_f
+                    left_root = sample_left_variance/args[0].size.to_r
+                    right_root = sample_right_variance/args[1].size.to_r
 
                     standard_error = Math.sqrt(left_root + right_root)
 
-                    (sample_left_mean - sample_right_mean).abs/standard_error.to_f
+                    (sample_left_mean - sample_right_mean).abs/standard_error.to_r
                   end
 
         t_distribution = Distribution::TStudent.new(degrees_of_freedom)
@@ -72,7 +72,7 @@ module Statistics
 
         down = difference_std/Math.sqrt(differences.size)
 
-        t_score = (differences.mean - 0)/down.to_f
+        t_score = (differences.mean - 0)/down.to_r
 
         probability = Distribution::TStudent.new(degrees_of_freedom).cumulative_function(t_score)
 

data/lib/statistics/statistical_test/wilcoxon_rank_sum_test.rb

@@ -73,7 +73,7 @@ module Statistics
                     memo += ((t[:counter] ** 3) - t[:counter])/12.0
                   end
 
-        left = (total_group_one * total_group_two)/(n * (n - 1)).to_f
+        left = (total_group_one * total_group_two)/(n * (n - 1)).to_r
         right = (((n ** 3) - n)/12.0) - rank_sum
 
         Math.sqrt(left * right)
@@ -82,7 +82,7 @@ module Statistics
       private def ranked_sum_for(total, group)
         # sum rankings per group
         group.reduce(0) do |memo, element|
-          rank_of_element = total[element][:rank] / total[element][:counter].to_f
+          rank_of_element = total[element][:rank] / total[element][:counter].to_r
           memo += rank_of_element
         end
       end

data/lib/statistics/version.rb

@@ -1,3 +1,3 @@
 module Statistics
-  VERSION = "2.1.2"
+  VERSION = "2.1.3"
 end

data/ruby-statistics.gemspec

@@ -27,8 +27,9 @@ Gem::Specification.new do |spec|
   spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
   spec.require_paths = ["lib"]
 
-  spec.add_development_dependency "rake", '~> 12.0', '>= 12.0.0'
+  spec.add_development_dependency "rake", '>= 12.0.0', '~> 13.0'
   spec.add_development_dependency "rspec", '>= 3.6.0'
   spec.add_development_dependency "grb", '~> 0.4.1', '>= 0.4.1'
   spec.add_development_dependency 'byebug', '>= 9.1.0'
+  spec.add_development_dependency 'pry'
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: ruby-statistics
 version: !ruby/object:Gem::Version
-  version: 2.1.2
+  version: 2.1.3
 platform: ruby
 authors:
 - esteban zapata
-autorequire: 
+autorequire:
 bindir: exe
 cert_chain: []
-date: 2020-03-01 00:00:00.000000000 Z
+date: 2021-02-04 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake
@@ -19,7 +19,7 @@ dependencies:
         version: 12.0.0
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '12.0'
+        version: '13.0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
@@ -29,7 +29,7 @@ dependencies:
         version: 12.0.0
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '12.0'
+        version: '13.0'
 - !ruby/object:Gem::Dependency
   name: rspec
   requirement: !ruby/object:Gem::Requirement
@@ -48,20 +48,20 @@ dependencies:
   name: grb
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ">="
+    - - "~>"
       - !ruby/object:Gem::Version
         version: 0.4.1
-    - - "~>"
+    - - ">="
       - !ruby/object:Gem::Version
         version: 0.4.1
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ">="
+    - - "~>"
       - !ruby/object:Gem::Version
         version: 0.4.1
-    - - "~>"
+    - - ">="
       - !ruby/object:Gem::Version
         version: 0.4.1
 - !ruby/object:Gem::Dependency
@@ -78,6 +78,20 @@ dependencies:
     - - ">="
       - !ruby/object:Gem::Version
         version: 9.1.0
+- !ruby/object:Gem::Dependency
+  name: pry
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
 description: |-
   This gem is intended to accomplish the same purpose as jStat js library:
                             to provide ruby with statistical capabilities without the need
@@ -90,6 +104,7 @@ executables: []
 extensions: []
 extra_rdoc_files: []
 files:
+- ".github/dependabot.yml"
 - ".github/workflows/ruby.yml"
 - ".gitignore"
 - ".rspec"
@@ -134,7 +149,7 @@ homepage: https://github.com/estebanz01/ruby-statistics
 licenses:
 - MIT
 metadata: {}
-post_install_message: 
+post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -149,8 +164,8 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.0.8
-signing_key: 
+rubygems_version: 3.1.4
+signing_key:
 specification_version: 4
 summary: A ruby gem for som specific statistics. Inspired by the jStat js library.
 test_files: []