ruby-statistics 0.5.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 1c58bc877eb893a77eacc586ddbc2cfca7e3655d
+  data.tar.gz: e61823ec3df43e872d451668ad8148b29b3749b8
+SHA512:
+  metadata.gz: 7fe62c3c65d455fdc40c3dca386c0bce4672025c56638cadea586720b7c10ea57b6d6cf91379772f5c41c080fa4a1cb3e4b7894fe328bde45f03e616dadf0c8d
+  data.tar.gz: 5de5004f6ff87991dc41d273a05db3ab2dbe1c87c001ec5ef3b6e93d5145ac5f7168b64ec55c6049de715b7aa82ac1cdd3dec0a771c2c9f31794067734ddb125

data/.gitignore

@@ -0,0 +1,65 @@
+*.gem
+*.rbc
+/.config
+/coverage/
+/InstalledFiles
+/pkg/
+/spec/reports/
+/spec/examples.txt
+/test/tmp/
+/test/version_tmp/
+/tmp/
+
+# Used by dotenv library to load environment variables.
+# .env
+
+## Specific to RubyMotion:
+.dat*
+.repl_history
+build/
+*.bridgesupport
+build-iPhoneOS/
+build-iPhoneSimulator/
+
+## Specific to RubyMotion (use of CocoaPods):
+#
+# We recommend against adding the Pods directory to your .gitignore. However
+# you should judge for yourself, the pros and cons are mentioned at:
+# https://guides.cocoapods.org/using/using-cocoapods.html#should-i-check-the-pods-directory-into-source-control
+#
+# vendor/Pods/
+
+## Documentation cache and generated files:
+/.yardoc/
+/_yardoc/
+/doc/
+/rdoc/
+
+## Environment normalization:
+/.bundle/
+/vendor/bundle
+/lib/bundler/man/
+
+# for a library or gem, you might want to ignore these files since the code is
+# intended to run in multiple environments; otherwise, check them in:
+# Gemfile.lock
+.ruby-version
+.ruby-gemset
+
+# unless supporting rvm < 1.11.0 or doing something fancy, ignore this:
+.rvmrc
+/.bundle/
+/.yardoc
+/Gemfile.lock
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/
+
+# rspec failure tracking
+.rspec_status
+
+# byebug
+.byebug_history

data/.rspec

@@ -0,0 +1,2 @@
+--format documentation
+--color

data/.travis.yml

@@ -0,0 +1,6 @@
+sudo: false
+language: ruby
+rvm:
+  - 2.2
+  - 2.3.1
+before_install: gem install bundler

data/CODE_OF_CONDUCT.md

@@ -0,0 +1,74 @@
+# Contributor Covenant Code of Conduct
+
+## Our Pledge
+
+In the interest of fostering an open and welcoming environment, we as
+contributors and maintainers pledge to making participation in our project and
+our community a harassment-free experience for everyone, regardless of age, body
+size, disability, ethnicity, gender identity and expression, level of experience,
+nationality, personal appearance, race, religion, or sexual identity and
+orientation.
+
+## Our Standards
+
+Examples of behavior that contributes to creating a positive environment
+include:
+
+* Using welcoming and inclusive language
+* Being respectful of differing viewpoints and experiences
+* Gracefully accepting constructive criticism
+* Focusing on what is best for the community
+* Showing empathy towards other community members
+
+Examples of unacceptable behavior by participants include:
+
+* The use of sexualized language or imagery and unwelcome sexual attention or
+advances
+* Trolling, insulting/derogatory comments, and personal or political attacks
+* Public or private harassment
+* Publishing others' private information, such as a physical or electronic
+  address, without explicit permission
+* Other conduct which could reasonably be considered inappropriate in a
+  professional setting
+
+## Our Responsibilities
+
+Project maintainers are responsible for clarifying the standards of acceptable
+behavior and are expected to take appropriate and fair corrective action in
+response to any instances of unacceptable behavior.
+
+Project maintainers have the right and responsibility to remove, edit, or
+reject comments, commits, code, wiki edits, issues, and other contributions
+that are not aligned to this Code of Conduct, or to ban temporarily or
+permanently any contributor for other behaviors that they deem inappropriate,
+threatening, offensive, or harmful.
+
+## Scope
+
+This Code of Conduct applies both within project spaces and in public spaces
+when an individual is representing the project or its community. Examples of
+representing a project or community include using an official project e-mail
+address, posting via an official social media account, or acting as an appointed
+representative at an online or offline event. Representation of a project may be
+further defined and clarified by project maintainers.
+
+## Enforcement
+
+Instances of abusive, harassing, or otherwise unacceptable behavior may be
+reported by contacting the project team at ezapata@altavistaed.com. All
+complaints will be reviewed and investigated and will result in a response that
+is deemed necessary and appropriate to the circumstances. The project team is
+obligated to maintain confidentiality with regard to the reporter of an incident.
+Further details of specific enforcement policies may be posted separately.
+
+Project maintainers who do not follow or enforce the Code of Conduct in good
+faith may face temporary or permanent repercussions as determined by other
+members of the project's leadership.
+
+## Attribution
+
+This Code of Conduct is adapted from the [Contributor Covenant][homepage], version 1.4,
+available at [http://contributor-covenant.org/version/1/4][version]
+
+[homepage]: http://contributor-covenant.org
+[version]: http://contributor-covenant.org/version/1/4/

data/Gemfile

@@ -0,0 +1,6 @@
+source "https://rubygems.org"
+
+git_source(:github) {|repo_name| "https://github.com/#{repo_name}" }
+
+# Specify your gem's dependencies in statistics.gemspec
+gemspec

data/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2017 Esteban Zapata Rojas
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2017 esteban zapata
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.

data/README.md

@@ -0,0 +1,59 @@
+# Statistics
+
+Welcome to your new gem! In this directory, you'll find the files you need to be able to package up your Ruby library into a gem. Put your Ruby code in the file `lib/statistics`. To experiment with that code, run `bin/console` for an interactive prompt.
+
+TODO: Delete this and the text above, and describe your gem
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'ruby-statistics'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install ruby-statistics
+
+## Basic Usage
+
+just require the `statistics` gem in order to load it. If you don't have defined the `Distribution` namespace, the gem will assign an alias, reducing the number of namespaces needed to use a class.
+
+Right know you can load:
+
+* The whole statistics gem. `require 'statistics'`
+* A namespace. `require 'statistics/distribution'`
+* A class. `require 'statistics/distribution/normal'`
+
+Feel free to use the one that is more convenient to you.
+
+### Hello-World Example
+```ruby
+require 'statistics'
+
+poisson = Distribution::Poisson.new(l) # Using Distribution alias.
+normal = Statistics::Distribution::StandardNormal.new # Using all namespaces.
+```
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake spec` to run the tests. You can also run `bin/console` for an interactive prompt that will allow you to experiment.
+
+To install this gem onto your local machine, run `bundle exec rake install`. To release a new version, update the version number in `version.rb`, and then run `bundle exec rake release`, which will create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/[USERNAME]/statistics. This project is intended to be a safe, welcoming space for collaboration, and contributors are expected to adhere to the [Contributor Covenant](http://contributor-covenant.org) code of conduct.
+
+## License
+
+The gem is available as open source under the terms of the [MIT License](http://opensource.org/licenses/MIT).
+
+## Code of Conduct
+
+Everyone interacting in the Statistics project’s codebases, issue trackers, chat rooms and mailing lists is expected to follow the [code of conduct](https://github.com/[USERNAME]/statistics/blob/master/CODE_OF_CONDUCT.md).

data/Rakefile

@@ -0,0 +1,6 @@
+require "bundler/gem_tasks"
+require "rspec/core/rake_task"
+
+RSpec::Core::RakeTask.new(:spec)
+
+task :default => :spec

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "statistics"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start(__FILE__)

data/bin/setup

@@ -0,0 +1,8 @@
+#!/usr/bin/env bash
+set -euo pipefail
+IFS=$'\n\t'
+set -vx
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/lib/math.rb

@@ -0,0 +1,120 @@
+module Math
+  def self.factorial(n)
+    return if n < 0
+
+    n = n.to_i # Only integers.
+
+    return 1 if n == 0 || n == 1
+    Math.gamma(n + 1) # Math.gamma(x) == (n - 1)! for integer values
+  end
+
+  def self.combination(n, r)
+    self.factorial(n)/(self.factorial(r) * self.factorial(n - r)).to_f # n!/(r! * [n - r]!)
+  end
+
+  def self.permutation(n, k)
+    self.factorial(n)/self.factorial(n - k).to_f
+  end
+
+  # Function adapted from the python implementation that exists in https://en.wikipedia.org/wiki/Simpson%27s_rule#Sample_implementation
+  # Finite integral in the interval [a, b] split up in n-intervals
+  def self.simpson_rule(a, b, n, &block)
+    unless n.even?
+      puts "The composite simpson's rule needs even intervals!"
+      return
+    end
+
+    h = (b - a)/n.to_f
+    resA = yield(a)
+    resB = yield(b)
+
+    sum = resA + resB
+
+    (1..n).step(2).each do |number|
+      res = yield(a + number * h)
+      sum += 4 * res
+    end
+
+    (1..(n-1)).step(2).each do |number|
+      res = yield(a + number * h)
+      sum += 2 * res
+    end
+
+    return sum * h / 3.0
+  end
+
+  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) do |t|
+      (t ** (s - 1)) * Math.exp(-t)
+    end
+  end
+
+  def self.beta_function(x, y)
+    return 1 if x == 1 && y == 1
+
+    (Math.gamma(x) * Math.gamma(y))/Math.gamma(x + y)
+  end
+
+  ### This implementation is an adaptation of the incomplete beta function made in C by
+  ### Lewis Van Winkle, which released the code under the zlib license.
+  ### The whole math behind this code is described in the following post: https://codeplea.com/incomplete-beta-function-c
+  def self.incomplete_beta_function(x, alp, bet)
+    return if x < 0.0
+    return 1.0 if x > 1.0
+
+    tiny = 1.0E-50
+
+    if x > ((alp + 1.0)/(alp + bet + 2.0))
+      return 1.0 - self.incomplete_beta_function(1.0 - x, bet, alp)
+    end
+
+    # 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
+
+    # 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)
+    # upper_left = (x ** alp) * ((1.0 - x) ** bet)
+    # front = upper_left/down_left
+
+    f, c, d = 1.0, 1.0, 0.0
+
+    returned_value = nil
+
+    # Let's do more iterations than the proposed implementation (200 iters)
+    (0..500).each do |number|
+      m = number/2
+
+      numerator = if number == 0
+                    1.0
+                  elsif number % 2 == 0
+                    (m * (bet - m) * x)/((alp + 2.0 * m - 1.0)* (alp + 2.0 * m))
+                  else
+                    top = -((alp + m) * (alp + bet + m) * x)
+                    down = ((alp + 2.0 * m) * (alp + 2.0 * m + 1.0))
+
+                    top/down
+                  end
+
+      d = 1.0 + numerator * d
+      d = tiny if d.abs < tiny
+      d = 1.0 / d
+
+      c = 1.0 + numerator / c
+      c = tiny if c.abs < tiny
+
+      cd = (c*d).freeze
+      f = f * cd
+
+
+      if (1.0 - cd).abs < 1.0E-10
+        returned_value = front * (f - 1.0)
+        break
+      end
+    end
+
+    returned_value
+  end
+end

data/lib/statistics/distribution/beta.rb

@@ -0,0 +1,35 @@
+module Statistics
+  module Distribution
+    class Beta
+      attr_accessor :alpha, :beta
+
+      def initialize(alp, bet)
+        self.alpha = alp.to_f
+        self.beta = bet.to_f
+      end
+
+      def cumulative_function(value)
+        Math.incomplete_beta_function(value, alpha, beta)
+      end
+
+      def density_function(value)
+        return 0 if value < 0 || value > 1 # Density function defined in the [0,1] interval
+
+        num = (value**(alpha - 1)) * ((1 - value)**(beta - 1))
+        den = Math.beta_function(alpha, beta)
+
+        num/den
+      end
+
+      def mode
+        return unless alpha > 1 && beta > 1
+
+        (alpha - 1)/(alpha + beta - 2)
+      end
+
+      def mean
+        alpha / (alpha + beta)
+      end
+    end
+  end
+end

data/lib/statistics/distribution/binomial.rb

@@ -0,0 +1,49 @@
+module Statistics
+  module Distribution
+    class Binomial
+      attr_accessor :number_of_trials, :probability_per_trial
+      def initialize(n, p)
+        self.number_of_trials = n.to_i
+        self.probability_per_trial = p
+      end
+
+      def probability_mass_function(k)
+        return if k < 0 || k > number_of_trials
+        k = k.to_i
+
+        Math.combination(number_of_trials, k) *
+          (probability_per_trial ** k) * ((1 - probability_per_trial) ** (number_of_trials - k))
+      end
+
+      def cumulative_function(k)
+        return if k < 0 || k > number_of_trials
+        k = k.to_i
+
+        p = 1 - probability_per_trial
+        Math.incomplete_beta_function(p, number_of_trials - k, 1 + k)
+      end
+
+      def mean
+        number_of_trials * probability_per_trial
+      end
+
+      def variance
+        mean * (1 - probability_per_trial)
+      end
+
+      def mode
+        test = (number_of_trials + 1) * probability_per_trial
+
+        returned = if test == 0 || (test % 1 != 0)
+                     test.floor
+                   elsif (test % 1 == 0)  && (test >= 1 && test <= number_of_trials)
+                     [test, test - 1]
+                   elsif test == number_of_trials + 1
+                     number_of_trials
+                   end
+
+        returned
+      end
+    end
+  end
+end

data/lib/statistics/distribution/chi_squared.rb

@@ -0,0 +1,37 @@
+module Statistics
+  module Distribution
+    class ChiSquared
+      attr_accessor :degrees_of_freedom
+
+      alias_method :mean, :degrees_of_freedom
+
+      def initialize(k)
+        self.degrees_of_freedom = k
+      end
+
+      def cumulative_function(value)
+        k = degrees_of_freedom/2.0
+        Math.lower_incomplete_gamma_function(k, value/2.0)/Math.gamma(k)
+      end
+
+      def density_function(value)
+        return 0 if value < 0
+
+        common = degrees_of_freedom/2.0
+
+        left_down = (2 ** common) * Math.gamma(common)
+        right = (value ** (common - 1)) * Math.exp(-(value/2.0))
+
+        (1.0/left_down) * right
+      end
+
+      def mode
+        [degrees_of_freedom - 2, 0].max
+      end
+
+      def variance
+        degrees_of_freedom * 2
+      end
+    end
+  end
+end

data/lib/statistics/distribution/f.rb

@@ -0,0 +1,46 @@
+module Statistics
+  module Distribution
+    class F
+      attr_accessor :d1, :d2 # Degrees of freedom #1 and #2
+
+      def initialize(k, j)
+        self.d1 = k
+        self.d2 = j
+      end
+
+      # 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)
+
+        1 - Math.incomplete_beta_function(k, d2/2.0, d1/2.0)
+      end
+
+      def density_function(value)
+        return if d1 < 0 || d2 < 0 # F-pdf is well defined for the [0, +infinity) interval.
+
+        val = value.to_f
+        upper = ((d1 * val) ** d1) * (d2**d2)
+        lower = (d1 * val + d2) ** (d1 + d2)
+        up = Math.sqrt(upper/lower.to_f)
+        down = val * Math.beta_function(d1/2.0, d2/2.0)
+
+        up/down.to_f
+      end
+
+      def mean
+        return if d2 <= 2
+
+        d2/(d2 - 2).to_f
+      end
+
+      def mode
+        return if d1 <= 2
+
+        left = (d1 - 2)/d1.to_f
+        right = d2/(d2 + 2).to_f
+
+        left * right
+      end
+    end
+  end
+end

data/lib/statistics/distribution/normal.rb

@@ -0,0 +1,43 @@
+module Statistics
+  module Distribution
+    class Normal
+      attr_accessor :mean, :standard_deviation, :variance
+      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
+      end
+
+      def cumulative_function(value)
+        (1/2.0) * (1.0 + Math.erf((value - mean)/(standard_deviation * Math.sqrt(2.0))))
+      end
+
+      def density_function(value)
+        return 0 if standard_deviation <= 0
+
+        up_right = (value - mean)**2.0
+        down_right = 2.0 * variance
+        right = Math.exp(-(up_right/down_right))
+        left_down = Math.sqrt(2.0 * Math::PI * variance)
+        left_up = 1.0
+
+        (left_up/(left_down) * right)
+      end
+    end
+
+    class StandardNormal < Normal
+      def initialize
+        super(0, 1) # Mean = 0, Std = 1
+      end
+
+      def density_function(value)
+        pow = (value**2)/2.0
+        euler = Math.exp(-pow)
+
+        euler/Math.sqrt(2 * Math::PI)
+      end
+    end
+  end
+end

data/lib/statistics/distribution/poisson.rb

@@ -0,0 +1,38 @@
+module Statistics
+  module Distribution
+    class Poisson
+      attr_accessor :expected_number_of_occurrences
+
+      alias_method :mean, :expected_number_of_occurrences
+      alias_method :variance, :expected_number_of_occurrences
+
+      def initialize(l)
+        self.expected_number_of_occurrences = l
+      end
+
+      def probability_mass_function(k)
+        return if k < 0 || expected_number_of_occurrences < 0
+
+        k = k.to_i
+
+        upper = (expected_number_of_occurrences ** k) * Math.exp(-expected_number_of_occurrences)
+        lower = Math.factorial(k)
+
+        upper/lower.to_f
+      end
+
+      def cumulative_function(k)
+        return if k < 0 || expected_number_of_occurrences < 0
+
+        k = k.to_i
+
+        upper = Math.lower_incomplete_gamma_function((k + 1).floor, expected_number_of_occurrences)
+        lower = Math.factorial(k.floor)
+
+        # 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)
+      end
+    end
+  end
+end

data/lib/statistics/distribution/t_student.rb

@@ -0,0 +1,50 @@
+module Statistics
+  module Distribution
+    class TStudent
+      attr_accessor :degrees_of_freedom
+      attr_reader :mode
+
+      def initialize(v)
+        self.degrees_of_freedom = v
+        @mode = 0
+      end
+
+      ### Extracted from https://codeplea.com/incomplete-beta-function-c
+      ### This function is shared under zlib license and the author is Lewis Van Winkle
+      def cumulative_function(value)
+        upper = (value + Math.sqrt(value * value + degrees_of_freedom))
+        lower = (2.0 * Math.sqrt(value * value + degrees_of_freedom))
+
+        x = upper/lower
+
+        alpha = degrees_of_freedom/2.0
+        beta = degrees_of_freedom/2.0
+
+        Math.incomplete_beta_function(x, alpha, beta)
+      end
+
+      def density_function(value)
+        return if degrees_of_freedom <= 0
+
+        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)
+
+        left * right
+      end
+
+      def mean
+        0 if degrees_of_freedom > 1
+      end
+
+      def variance
+        if degrees_of_freedom > 1 && degrees_of_freedom <= 2
+          Float::INFINITY
+        elsif degrees_of_freedom > 2
+          degrees_of_freedom/(degrees_of_freedom - 2.0)
+        end
+      end
+    end
+  end
+end

data/lib/statistics/distribution/uniform.rb

@@ -0,0 +1,40 @@
+module Statistics
+  module Distribution
+    class Uniform
+      attr_accessor :left, :right
+
+      def initialize(a, b)
+        self.left = a.to_f
+        self.right = b.to_f
+      end
+
+      def density_function(value)
+        if value >= left && value <= right
+          1/(right - left)
+        else
+          0
+        end
+      end
+
+      def cumulative_function(value)
+        if value < left
+          0
+        elsif value >= left && value <= right
+          (value - left)/(right - left)
+        else
+          1
+        end
+      end
+
+      def mean
+        (1/2.0) * ( left + right )
+      end
+      alias_method :median, :mean
+
+
+      def variance
+        (1/12.0) * ( right - left ) ** 2
+      end
+    end
+  end
+end

data/lib/statistics/distribution/weibull.rb

@@ -0,0 +1,46 @@
+module Statistics
+  module Distribution
+    class Weibull
+      attr_accessor :shape, :scale # k and lambda
+
+      def initialize(k, lamb)
+        self.shape = k.to_f
+        self.scale = lamb.to_f
+      end
+
+      def cumulative_function(random_value)
+        return 0 if random_value < 0
+
+        1 - Math.exp(-((random_value/scale) ** shape))
+      end
+
+      def density_function(value)
+        return if shape <= 0 || scale <= 0
+        return 0 if value < 0
+
+        left = shape/scale
+        center = (value/scale)**(shape - 1)
+        right = Math.exp(-((value/scale)**shape))
+
+        left * center * right
+      end
+
+      def mean
+        scale * Math.gamma(1 + (1/shape))
+      end
+
+      def mode
+        return 0 if shape <= 1
+
+        scale * (((shape - 1)/shape) ** (1/shape))
+      end
+
+      def variance
+        left = Math.gamma(1 + (2/shape))
+        right = Math.gamma(1 + (1/shape)) ** 2
+
+        (scale ** 2) * (left - right)
+      end
+    end
+  end
+end

data/lib/statistics/distribution.rb

@@ -0,0 +1,11 @@
+Dir['./lib/statistics/distribution/**/*.rb'].each {|file| require file }
+
+module Statistics
+  module Distribution
+  end
+end
+
+# If Distribution is not defined, setup alias.
+if defined?(Statistics) && !(defined?(Distribution))
+  Distribution = Statistics::Distribution
+end

data/lib/statistics/version.rb

@@ -0,0 +1,3 @@
+module Statistics
+  VERSION = "0.5.0"
+end

data/lib/statistics.rb

@@ -0,0 +1,8 @@
+require 'descriptive_statistics'
+
+Dir['./lib/statistics/**/*.rb'].each {|file| require file }
+require_relative './math'
+
+module Statistics
+  # Your code goes here...
+end

data/ruby-statistics.gemspec

@@ -0,0 +1,36 @@
+# coding: utf-8
+lib = File.expand_path("../lib", __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require "statistics/version"
+
+Gem::Specification.new do |spec|
+  spec.name          = "ruby-statistics"
+  spec.version       = Statistics::VERSION
+  spec.authors       = ["esteban zapata"]
+  spec.email         = ["estebanz01@outlook.com"]
+
+  spec.summary       = %q{A ruby gem for som specific statistics. Inspired by the jStat js library.}
+  spec.description   = %q{This gem is intended to accomplish the same purpose as jStat js library:
+                          to provide ruby with statistical capabilities without the need
+                          of a statistical programming language like R or Octave. Some functions
+                          and capabilities are an implementation from other authors and are
+                          referenced properly in the class/method.}
+  spec.homepage      = "https://github.com/estebanz01/ruby-statistics"
+  spec.license       = "MIT"
+
+  # Prevent pushing this gem to RubyGems.org. To allow pushes either set the 'allowed_push_host'
+  # to allow pushing to a single host or delete this section to allow pushing to any host.
+  spec.files         = `git ls-files -z`.split("\x0").reject do |f|
+    f.match(%r{^(test|spec|features)/})
+  end
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", '>= 1.15.4'
+  spec.add_development_dependency "rake", '>= 12.0.0'
+  spec.add_development_dependency "rspec", '>= 3.6.0'
+  spec.add_development_dependency "grb", '0.4.1'
+  spec.add_development_dependency 'byebug', '9.1.0'
+  spec.add_dependency 'descriptive_statistics', '2.5.1'
+end

metadata

@@ -0,0 +1,158 @@
+--- !ruby/object:Gem::Specification
+name: ruby-statistics
+version: !ruby/object:Gem::Version
+  version: 0.5.0
+platform: ruby
+authors:
+- esteban zapata
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2017-09-11 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.15.4
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.15.4
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 12.0.0
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 12.0.0
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 3.6.0
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 3.6.0
+- !ruby/object:Gem::Dependency
+  name: grb
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - '='
+      - !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::Dependency
+  name: byebug
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - '='
+      - !ruby/object:Gem::Version
+        version: 9.1.0
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - '='
+      - !ruby/object:Gem::Version
+        version: 9.1.0
+- !ruby/object:Gem::Dependency
+  name: descriptive_statistics
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - '='
+      - !ruby/object:Gem::Version
+        version: 2.5.1
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - '='
+      - !ruby/object:Gem::Version
+        version: 2.5.1
+description: |-
+  This gem is intended to accomplish the same purpose as jStat js library:
+                            to provide ruby with statistical capabilities without the need
+                            of a statistical programming language like R or Octave. Some functions
+                            and capabilities are an implementation from other authors and are
+                            referenced properly in the class/method.
+email:
+- estebanz01@outlook.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- CODE_OF_CONDUCT.md
+- Gemfile
+- LICENSE
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- lib/math.rb
+- lib/statistics.rb
+- lib/statistics/distribution.rb
+- lib/statistics/distribution/beta.rb
+- lib/statistics/distribution/binomial.rb
+- lib/statistics/distribution/chi_squared.rb
+- lib/statistics/distribution/f.rb
+- lib/statistics/distribution/normal.rb
+- lib/statistics/distribution/poisson.rb
+- lib/statistics/distribution/t_student.rb
+- lib/statistics/distribution/uniform.rb
+- lib/statistics/distribution/weibull.rb
+- lib/statistics/version.rb
+- ruby-statistics.gemspec
+homepage: https://github.com/estebanz01/ruby-statistics
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.6.13
+signing_key: 
+specification_version: 4
+summary: A ruby gem for som specific statistics. Inspired by the jStat js library.
+test_files: []