statistical 0.1.1 → 0.1.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 292baa7cd3d35aaa61324625eaaf0558e8436145
-  data.tar.gz: 67d5997eb4a9d622e5467eed349dd448ff0fb4bd
+  metadata.gz: 4dbc8862a8de1c1bd3f982635434403c1664683e
+  data.tar.gz: a1a52cffbb67236dcb224fe72634d8eb1434a4d8
 SHA512:
-  metadata.gz: 4fae7c95b560274b9251969e57ebe2559b8eb59fba27090accb5540b869ff9c8df5a4395008011efc7878cdb70ba24afe1a4a6224aac19d616f09603e072547a
-  data.tar.gz: c52abd3464ed8511b6a76950f7c1d827453109e77f41073d2a1937d6d0b208b3f3d8f3fc40e8283ca64f9f7de038b92a98dd30553ce9c87ed4119c572b1710a8
+  metadata.gz: 70b321f214492e41ba6b0f18be837404d2ba8d6b3687b21345825865c0e6d0ceead63c9b905fca095ab68980add4e3ae15405437e2afcb8fb9bdbd81e97c4c2e
+  data.tar.gz: ea3e9138092628d5d1551ad0d2eb8a26d6c929001c3a84a1d949a5afa1750a8ac2a00041794317b9c294396cae6bd4d7e7497ef799722f7fcfc12235f33e5d66

data/.gitignore

@@ -12,4 +12,6 @@
 /coverage/**
 *.so
 *.o
-.gitmsg
+.gitmsg
+/profile/**
+spec/tmp/**

data/Gemfile

@@ -1,6 +1,10 @@
 source 'https://rubygems.org'
 
-gem 'codeclimate-test-reporter', :group => :test, :require => nil
+group :test do
+  gem 'codeclimate-test-reporter', require: nil
+  gem 'rspec-prof'
+end
+
 gem 'yard', require: false
 
 # Other dependencies go in statistical.gemspec

data/Rakefile

@@ -10,6 +10,7 @@ require 'statistical/version'
 
 RSpec::Core::RakeTask.new(:spec) do |t|
   t.verbose = false
+  t.rspec_opts = '-f p'
 end
 
 # Add a cop task for code linting
@@ -23,7 +24,7 @@ end
 
 # Build the gem
 task :build do
-  system "gem build statistical.gemspec"
+  system 'gem build statistical.gemspec'
 end
 
 # Install gem locally
@@ -34,7 +35,7 @@ end
 # Release gem to github
 task :release => :build do
   system "git tag -a v#{Statistical::VERSION} -m 'Version #{Statistical::VERSION}'"
-  system "git push --tags"
+  system 'git push --tags'
   system "gem push statistical-#{Statistical::VERSION}.gem"
 end
 

data/data/template/distribution.erb

@@ -1,4 +1,4 @@
-require 'statistical/exceptions'
+require 'statistical/helpers'
 
 module Statistical
   module Distribution
@@ -70,9 +70,7 @@ module Statistical
       # @return [Boolean] true if-and-only-if two instances are of the same
       #   class and have the same parameters.
       def eql?(other)
-        return false unless other.is_a?(self.class)
-        # Compare parameters and other stuff here
-        return true
+        return other.is_a?(self.class)
       end
 
       alias :== :eql?

data/data/template/rng.erb

@@ -1,14 +1,13 @@
-require 'statistical/exceptions'
 require 'statistical/distribution/<%= distribution%>'
 
 module Statistical
-  # Companion RNG class for the continuous uniform distribution. Requires a
+  # Companion RNG class for the continuous <%= distribution.capitalize.camelcase %> distribution. Requires a
   #   distrbution object of the corresponding distribution
   #
   # @author Vaibhav Yenamandra
   #
-  # @attr_reader [Numeric] lower The lower bound of the uniform distribution.
-  # @attr_reader [Numeric] upper The upper bound of the uniform distribution.
+  # @attr_reader [Numeric] lower The lower bound of the <%= distribution.capitalize.camelcase %> distribution.
+  # @attr_reader [Numeric] upper The upper bound of the <%= distribution.capitalize.camelcase %> distribution.
   module Rng
     class <%= distribution.capitalize.camelcase %>
       attr_reader :generator, # other params go here
@@ -18,15 +17,16 @@ module Statistical
           raise TypeError, "Expected Distribution object or nil, found #{dobj.class}"
         end
         dobj = Statistical::Distribution::<%= distribution.capitalize.camelcase %>.new if dobj.nil?
-        @generator = Rng::Uniform.new(Distribution::Uniform.new, seed)
+        @generator = Random.new(seed)
         # Map other parameters here
         @sdist = dobj
       end
 
       # Return the next random number from the sequence
       #
-      # @return next random number in the sequence
+      # @return [Float] next random number in the sequence
       def rand
+        @sdist.quantile(@generator.rand)
       end
 
       # Compare against another rng to see if they are the same
@@ -34,9 +34,8 @@ module Statistical
       # @return true if and only if, source distributions are the same and the
       #   prng has the same initial state
       def eql?(other)
-        return false unless other.is_a?(self.class)
-        # Compare parameters and other stuff here
-        @generator == other.generator
+        return other.is_a?(self.class) &&
+               @generator == other.generator
       end
 
       # Return the type of the source distribution

data/data/template/spec.erb

@@ -3,17 +3,11 @@ require 'statistical/rng/<%= distribution%>'
 require 'statistical/distribution/<%= distribution %>'
 
 describe Statistical::Rng::<%= distribution.capitalize.camelcase %> do
-  it 'passes the G-test at a 95% significance level' do
-  end
 
   describe '.new' do
-    context 'when called with no arguments' do
-      fail
-    end
+    context 'when called with no arguments'
 
-    context 'when parameters are specified' do
-      fail
-    end
+    context 'when parameters are specified'
 
     context 'when initialized with a seed' do
       it 'should be repeatable for the same arguments' do
@@ -23,17 +17,11 @@ describe Statistical::Rng::<%= distribution.capitalize.camelcase %> do
   end
 
   describe '#rand' do
-    it 'returns a number between 0 and 1 by default' do
-      fail
-    end
-
-    it 'returns a bounded value when bounds are specified' do
-      fail
-    end
+    it 'passes the G-test at a 95% significance level'
   end
 
   describe '#==' do
-    context 'when compared against another uniform distribution' do
+    context 'when compared against another <%= distribution.capitalize.camelcase %> distribution' do
       it 'should return true if the bounds and seed are the same' do
         fail
       end
@@ -42,101 +30,70 @@ describe Statistical::Rng::<%= distribution.capitalize.camelcase %> do
         fail
       end
     end
-
-    context 'when compared against any other distribution class' do
-      it 'should return false if classes are different' do
-        fail
-      end
-    end
   end
 end
 
 
 describe Statistical::Distribution::<%= distribution.capitalize.camelcase %> do
   describe '.new' do
-    context 'when called with no arguments' do
-      fail
-    end
-    
-    context 'when upper and lower bounds are specified' do
-      fail
-    end
+    context 'when called with no arguments'
+
+    context 'when upper and lower bounds are specified'
   end
-  
-  
+
+
   describe '#pdf' do
-    context 'when called with x < lower_bound' do
-      fail
-    end
-    
-    context 'when called with x > upper_bound' do
-      fail
-    end
-    
-    context 'when called with x in [lower_bound, upper_bound]' do
-      fail
-    end
+    context 'when called with x < lower_bound'
+
+    context 'when called with x > upper_bound'
+
+    context 'when called with x in [lower_bound, upper_bound]'
   end
-  
+
   describe '#cdf' do
-    context 'when called with x < lower' do
-      fail
-    end
-    
-    context 'when called with x > upper' do
-      fail
-    end
-    
-    context 'when called with x in [lower, upper]' do
-      fail
-    end
+    context 'when called with x < lower'
+
+    context 'when called with x > upper'
+
+    context 'when called with x in [lower, upper]'
   end
-  
+
   describe '#quantile' do
-    context 'when called for x > 1' do
-      fail
-    end
-    
     context 'when called for x < 0' do
-      fail
+      let(:<%= distribution[0] %>dist) {Statistical::Distribution::<%= distribution.capitalize.camelcase %>.new}
+      it {
+        expect {<%= distribution[0] %>dist.quantile(-Float::EPSILON)}.to raise_error(RangeError)
+      }
     end
-    
+
+    context 'when called for x > 1' do
+      let(:<%= distribution[0] %>dist) {Statistical::Distribution::<%= distribution.capitalize.camelcase %>.new}
+      it {  
+        expect {<%= distribution[0] %>dist.quantile(1 + Float::EPSILON)}.to raise_error(RangeError)
+      }
+    end
+
     context 'when called for x in [0, 1]' do
-      fail
+      let(:<%= distribution[0] %>dist) {Statistical::Distribution::<%= distribution.capitalize.camelcase %>.new}
     end
   end
-  
+
   describe '#p_value' do
-    it 'should be the same as #quantile' do
-      fail
-    end
-  end  
-  
+    it 'should be the same as #quantile'
+  end
+
   describe '#mean' do
-    it 'should return the correct mean' do
-      fail
-    end
+    it 'should return the correct mean'
   end
-  
+
   describe '#variance' do
-    it 'should return the correct variance' do
-      fail
-    end
+    it 'should return the correct variance'
   end
-  
+
   describe '#==' do
-    context 'when compared against another Uniform distribution' do
-      it 'returns `true` if they have the same parameters' do
-        fail
-      end
-      
-      it 'returns `false` if they have different parameters' do
-        fail
-      end
-    end
-    
-    context 'when compared against any distribution type' do
-      fail
+    context 'when compared against another <%= distribution.capitalize.camelcase %> distribution' do
+      it 'returns `true` if they have the same parameters'
+      it 'returns `false` if they have different parameters'
     end
   end
 end

data/lib/statistical.rb

@@ -1,7 +1,12 @@
 require 'statistical/version'
 require 'statistical/distribution'
 require 'statistical/rng'
-require 'core_extensions'
+require 'statistical/core_extensions'
 
 module Statistical
+  using Statistical::StringExtensions
+  using Statistical::ArrayExtensions
+
+  # Truncated Euler-Mascheroni constant
+  EULER_GAMMA = 0.5772156649015328
 end

data/lib/statistical/core_extensions.rb

@@ -0,0 +1,44 @@
+# All core class modifications go here
+module Statistical
+  # Module to contain all String refinements
+  module StringExtensions
+    refine String do
+      # Convert CamelCase to snake_case
+      def snakecase
+        gsub(/::/, '/')
+          .gsub(/([A-Z]+)([A-Z][a-z])/, '\1_\2')
+          .gsub(/([a-z\d])([A-Z])/, '\1_\2')
+          .tr('-', '_')
+          .downcase
+      end
+
+      # convert snake_case to CamelCase
+      def camelcase
+        split('_').map(&:capitalize).join
+      end
+    end
+  end
+
+  # Module to contain all Array refinements
+  module ArrayExtensions
+    refine Array do
+      def sum
+        inject(:+)
+      end
+
+      def mean
+        sum.fdiv(count)
+      end
+
+      # Population variance
+      def pvariance
+        map { |x| (x - mean)**2}.sum.fdiv(count)
+      end
+
+      # Sample variance
+      def svariance
+        map { |x| (x - mean)**2}.sum.fdiv(count - 1)
+      end
+    end
+  end
+end

data/lib/statistical/distribution.rb

@@ -5,11 +5,15 @@ require 'statistical/distribution/bernoulli'
 require 'statistical/distribution/exponential'
 require 'statistical/distribution/laplace'
 require 'statistical/distribution/weibull'
+require 'statistical/distribution/frechet'
+require 'statistical/distribution/gumbel'
 
 module Statistical
   # Factory module used to create instances of various distributions classes
   # nested under itself
   module Distribution
+    using Statistical::StringExtensions
+
     # @private
     # No need to document this
     # Dynamically add constants when called

data/lib/statistical/distribution/frechet.rb

@@ -0,0 +1,125 @@
+require 'statistical/helpers'
+
+module Statistical
+  module Distribution
+    # This class models the Frechet (or inverse Weibull) distribution which is
+    #   a type-2 extreme value distribution.
+    #
+    # @author Vaibhav Yenamandra
+    # @attr_reader [Numeric] alpha Mandatory shape parameter of the distribution
+    # @attr_reader [Numeric] location Optional location parameter
+    # @attr_reader [Numeric] scale Optional scale parameter
+    # @attr_reader [Numeric] support Support / domain of the distribution's
+    #   PDF / CDF
+    class Frechet
+      attr_reader :alpha, :location, :scale, :support
+
+      # Returns a new `Statistical::Distribution::Frechet` instance
+      #
+      # @param [Numeric] alpha Mandatory shape parameter of the distribution
+      # @param [Numeric] location Optional location parameter
+      # @param [Numeric] scale Optional scale parameter
+      # @return [Frechet] `Statistical::Distribution::Frechet` instance
+      def initialize(alpha = nil, location = 0, scale = 1)
+        raise ArgumentError if alpha.nil?
+
+        @alpha = alpha.to_f
+        @location = location.to_f
+        @scale = scale.to_f
+        @support = Domain[@location, Float::INFINITY, :full_open]
+      end
+
+      # Returns value of probability density function at a point.
+      #
+      # @param [Numeric] x A real valued point
+      # @return [Float] Probability density function evaluated at `x`
+      def pdf(x)
+        xs = (x - @location) / @scale
+
+        # Not following the usual array lookup idiom here since the PDF
+        # expression is ill-defined (Complex) for x < location. Lazily evaluated
+        # expressions would help retain the array lookup idiom, but that's for
+        # later
+        if (@support <=> x).zero?
+          return (@alpha / @scale) * xs**(-1 - @alpha) * Math.exp(-(xs**-@alpha))
+        else
+          return 0
+        end
+      end
+
+      # Returns value of cumulative density function upto the specified point
+      #
+      # @param [Numeric] x A real valued point
+      # @return [Float] Cumulative density function evaluated for F(u <= x)
+      def cdf(x)
+        xs = (x - @location) / @scale
+
+        case @support <=> x
+        when 0
+          return Math.exp(-(xs**-@alpha))
+        when 1
+          return 1
+        when -1
+          return 0
+        end
+      end
+
+      # Returns value of inverse CDF for a given probability
+      #
+      # @see #p_value
+      #
+      # @param [Numeric] p a value within [0, 1]
+      # @return [Float] Inverse CDF for valid p
+      # @raise [RangeError] if p > 1 or p < 0
+      def quantile(p)
+        raise RangeError, "`p` must be in [0, 1], found: #{p}" if p < 0 || p > 1
+        return @location + @scale * ((-Math.log(p))**(-1 / @alpha))
+      end
+
+      # Returns the mean value for the calling instance. Calculated mean, and
+      #   not inferred from simulations
+      #
+      # @return Mean of the distribution
+      def mean
+        return [
+          Float::INFINITY,
+          @location + @scale * Math.gamma(1 - 1 / @alpha),
+          Float::INFINITY
+        ][@alpha <=> 1.0]
+      end
+
+      # Returns the expected value of variance for the calling instance.
+      #
+      # @return Variance of the distribution
+      def variance
+        return [
+          Float::INFINITY,
+          @scale * @scale * (
+            Math.gamma(1 - 2 / @alpha) -
+            Math.gamma(1 - 1 / @alpha)**2
+          ),
+          Float::INFINITY
+        ][@alpha <=> 2.0]
+      end
+
+      # Compares two distribution instances and returns a boolean outcome
+      #   Available publicly as #==
+      #
+      # @private
+      #
+      # @return [Boolean] true if-and-only-if two instances are of the same
+      #   class and have the same parameters.
+      def eql?(other)
+        return other.is_a?(self.class) &&
+               @alpha == other.alpha &&
+               @location == other.location &&
+               @scale == other.scale
+      end
+
+      alias :== :eql?
+      alias :p_value :quantile
+
+      private :eql?
+    end
+  end
+end

data/lib/statistical/distribution/gumbel.rb

@@ -0,0 +1,98 @@
+require 'statistical/helpers'
+
+module Statistical
+  module Distribution
+    # Distribution class for Extreme value type-3 or Gumbel distribution
+    #
+    # @author Vaibhav Yenamandra
+    #
+    # @attr_reader location [Float] The distributions location parameter
+    # @attr_reader scale [Float] The dstribution's scale parameter
+    # @attr_reader support [Float] The region of the real line where this
+    #   distribution is defined to exist
+    class Gumbel
+      attr_reader :location, :scale, :support
+
+      # Returns a new `Statistical::Distribution::Gumbel` instance
+      #
+      # @param [Types] location Description
+      # @return `Statistical::Distribution::Gumbel` instance
+      def initialize(location = 0, scale = 1)
+        @location = location.to_f
+        @scale = scale.to_f
+        @support = Domain[-Float::INFINITY, Float::INFINITY, :full_open]
+      end
+
+      # Returns value of probability density function at a point.
+      #
+      # @param [Numeric] x A real valued point
+      # @return Probability density function evaluated at `x`
+      def pdf(x)
+        xa = (x - @location) / @scale
+        return [
+          Math.exp(-xa - Math.exp(-xa)) / @scale,
+          0.0,
+          0.0
+        ][@support <=> x]
+      end
+
+      # Returns value of cumulative density function upto a point.
+      #
+      # @param [Numeric] x A real valued point
+      # @return Cumulative density function evaluated for F(u <= x)
+      def cdf(x)
+        xa = (x - @location) / @scale
+        return [
+          Math.exp(-Math.exp(-xa)),
+          1.0,
+          0.0
+        ][@support <=> x]
+      end
+
+      # Returns value of inverse CDF for a given probability
+      #
+      # @see #p_value
+      #
+      # @param [Numeric] p a value within [0, 1]
+      # @return Inverse CDF for valid p
+      # @raise [RangeError] if p > 1 or p < 0
+      def quantile(p)
+        raise RangeError, "`p` must be in [0, 1], found: #{p}" if p < 0 || p > 1
+        return @location - @scale * Math.log(-Math.log(p))
+      end
+
+      # Returns the expected value of the mean for the calling instance.
+      #
+      # @return Mean of the distribution
+      def mean
+        return @location + @scale * Statistical::EULER_GAMMA
+      end
+
+      # Returns the expected value of variance for the calling instance.
+      #
+      # @return Variance of the distribution
+      def variance
+        return ((Math::PI * @scale)**2) / 6
+      end
+
+      # Compares two distribution instances and returns a boolean outcome
+      #   Available publicly as #==
+      #
+      # @private
+      #
+      # @param other A distribution object (preferred)
+      # @return [Boolean] true if-and-only-if two instances are of the same
+      #   class and have the same parameters.
+      def eql?(other)
+        return other.is_a?(self.class) &&
+               @location == other.location &&
+               @scale == other.scale
+      end
+
+      alias :== :eql?
+      alias :p_value :quantile
+
+      private :eql?
+    end
+  end
+end

data/lib/statistical/distribution/laplace.rb

@@ -43,8 +43,7 @@ module Statistical
       def cdf(x)
         return [0.5,
                 1.0 - 0.5 * Math.exp((@location - x).fdiv(@scale)),
-                0.5 * Math.exp((x - @location).fdiv(@scale))
-        ][x <=> @location]
+                0.5 * Math.exp((x - @location).fdiv(@scale))][x <=> @location]
       end
 
       # Returns value of inverse CDF for a given probability
@@ -59,8 +58,7 @@ module Statistical
 
         return [@location,
                 @location - @scale * Math.log(2 * (1.0 - p)),
-                @scale * Math.log(2 * p) + @location
-        ][p <=> 0.5]
+                @scale * Math.log(2 * p) + @location][p <=> 0.5]
       end
 
       # Returns the expected mean value for the calling instance.

data/lib/statistical/distribution/uniform_discrete.rb

@@ -1,3 +1,5 @@
+require 'statistical/core_extensions'
+
 module Statistical
   # Module to collect all abstractions of distributions
   module Distribution
@@ -8,6 +10,8 @@ module Statistical
     # @attr_reader [Array, Numeric] support The support set of valid values a
     # random variate from the distribution can take. Must have at least 1 value
     class UniformDiscrete
+      using Statistical::ArrayExtensions
+
       attr_reader :count, :support, :lower, :upper
       # Returns a model for the discrete uniform distribution on all elements
       # present in the given set of elemets `elems`
@@ -108,7 +112,7 @@ module Statistical
       #
       # @return [Float]  Variance of the distribution
       def variance
-        return @support.variance
+        return @support.pvariance
       end
 
       # Compares two distribution instances and returns a boolean

data/lib/statistical/distribution/weibull.rb

@@ -33,8 +33,7 @@ module Statistical
         return [(@shape / @scale) * ((x / @scale)**(@shape - 1)) *
                 Math.exp(-((x / @scale)**@shape)),
                 0.0,
-                0.0
-        ][@support <=> x]
+                0.0][@support <=> x]
       end
 
       # Returns value of cumulative density function at a point. Calculated
@@ -45,8 +44,7 @@ module Statistical
       def cdf(x)
         return [1 - Math.exp(-((x / @scale)**@shape)),
                 1.0,
-                0.0
-        ][@support <=> x]
+                0.0][@support <=> x]
       end
 
       # Returns value of inverse CDF for a given probability

data/lib/statistical/helpers.rb

@@ -83,9 +83,9 @@ module Statistical
       @exclusions.each do |e|
         case e
         when Fixnum, Bignum, Float
-          has_val = has_val || (e == val)
+          has_val = (has_val || (e == val))
         when Range
-          has_val ||= e.include?(val)
+          has_val = (has_val || e.include?(val))
         end
       end
       return has_val

data/lib/statistical/rng.rb

@@ -5,11 +5,15 @@ require 'statistical/rng/bernoulli'
 require 'statistical/rng/exponential'
 require 'statistical/rng/laplace'
 require 'statistical/rng/weibull'
+require 'statistical/rng/frechet'
+require 'statistical/rng/gumbel'
 
 module Statistical
   # Factory module to create instances of the various classes
   # nested under itself
   module Rng
+    using Statistical::StringExtensions
+
     # @private
     # No need to document this
     # Dynamically add constants when called

data/lib/statistical/rng/frechet.rb

@@ -0,0 +1,59 @@
+require 'statistical/distribution/frechet'
+
+module Statistical
+  # Companion RNG class for the continuous uniform distribution. Requires a
+  #   distrbution object of the corresponding distribution
+  #
+  # @author Vaibhav Yenamandra
+  #
+  # @attr_reader [Numeric] lower The lower bound of the uniform distribution.
+  # @attr_reader [Numeric] upper The upper bound of the uniform distribution.
+  module Rng
+    class Frechet
+      attr_reader :generator, :alpha, :location, :scale
+
+      def initialize(dobj = nil, seed = Random.new_seed)
+        if dobj.nil?
+          raise ArgumentError, 'Need a alpha-parametrized Frechet object!'
+        elsif !dobj.is_a?(Statistical::Distribution::Frechet)
+          raise TypeError, "Expected Frechet Distribution found #{dobj.class}"
+        end
+
+        @generator = Random.new(seed)
+        @alpha = dobj.alpha
+        @location = dobj.location
+        @scale = dobj.scale
+        @sdist = dobj
+      end
+
+      # Return the next random number from the sequence
+      #
+      # @return [Float] next random number in the sequence
+      def rand
+        return @sdist.quantile(@generator.rand)
+      end
+
+      # Compare against another rng to see if they are the same
+      #
+      # @return true if and only if, source distributions are the same and the
+      #   prng has the same initial state
+      def eql?(other)
+        return other.is_a?(self.class) &&
+               @generator == other.generator &&
+               @alpha == other.alpha &&
+               @location == other.location &&
+               @scale == other.scale
+      end
+
+      # Return the type of the source distribution
+      #
+      # @return source distribution's type
+      def type
+        @sdist.class
+      end
+
+      alias_method :==, :eql?
+      private :eql?
+    end
+  end
+end

data/lib/statistical/rng/gumbel.rb

@@ -0,0 +1,59 @@
+require 'statistical/distribution/gumbel'
+
+module Statistical
+  # Companion RNG class for the continuous Gumbel distribution. Requires a
+  #   distrbution object of the corresponding distribution
+  #
+  # @author Vaibhav Yenamandra
+  #
+  # @attr_reader [Float] location The location parameter of the Gumbel
+  #   distribution
+  # @attr_reader [Float] scale The scale parameter of the Gumbel distribution
+  # @attr_reader [Float] generator The underlying uniform variate source used
+  #   to power `Gumbel#rand`
+  module Rng
+    class Gumbel
+      attr_reader :generator, :location, :scale
+
+      def initialize(dobj = nil, seed = Random.new_seed)
+        unless dobj.nil? || dobj.is_a?(Statistical::Distribution::Gumbel)
+          raise TypeError,
+                "Expected Distribution object or nil, found #{dobj.class}"
+        end
+        dobj = Statistical::Distribution::Gumbel.new if dobj.nil?
+        @generator = Random.new(seed)
+        @location = dobj.location
+        @scale = dobj.scale
+        @sdist = dobj
+      end
+
+      # Return the next random number from the sequence
+      #
+      # @return [Float] next random number in the sequence
+      def rand
+        @sdist.quantile(@generator.rand)
+      end
+
+      # Compare against another rng to see if they are the same
+      #
+      # @return true if and only if, source distributions are the same and the
+      #   prng has the same initial state
+      def eql?(other)
+        return other.is_a?(self.class) &&
+               @generator == other.generator &&
+               @location == other.location &&
+               @scale == other.scale
+      end
+
+      # Return the type of the source distribution
+      #
+      # @return source distribution's type
+      def type
+        @sdist.class
+      end
+
+      alias_method :==, :eql?
+      private :eql?
+    end
+  end
+end

data/lib/statistical/version.rb

@@ -1,3 +1,3 @@
 module Statistical
-  VERSION = '0.1.1'.freeze
+  VERSION = '0.1.2'.freeze
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: statistical
 version: !ruby/object:Gem::Version
-  version: 0.1.1
+  version: 0.1.2
 platform: ruby
 authors:
 - Vaibhav Yenamandra
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2016-06-16 00:00:00.000000000 Z
+date: 2016-06-20 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -116,11 +116,13 @@ files:
 - data/template/distribution.erb
 - data/template/rng.erb
 - data/template/spec.erb
-- lib/core_extensions.rb
 - lib/statistical.rb
+- lib/statistical/core_extensions.rb
 - lib/statistical/distribution.rb
 - lib/statistical/distribution/bernoulli.rb
 - lib/statistical/distribution/exponential.rb
+- lib/statistical/distribution/frechet.rb
+- lib/statistical/distribution/gumbel.rb
 - lib/statistical/distribution/laplace.rb
 - lib/statistical/distribution/two_point.rb
 - lib/statistical/distribution/uniform.rb
@@ -130,6 +132,8 @@ files:
 - lib/statistical/rng.rb
 - lib/statistical/rng/bernoulli.rb
 - lib/statistical/rng/exponential.rb
+- lib/statistical/rng/frechet.rb
+- lib/statistical/rng/gumbel.rb
 - lib/statistical/rng/laplace.rb
 - lib/statistical/rng/two_point.rb
 - lib/statistical/rng/uniform.rb

data/lib/core_extensions.rb

@@ -1,35 +0,0 @@
-# All core class modifications go here
-String.class_eval do
-  # Convert CamelCase to snake_case
-  def snakecase
-    gsub(/::/, '/')
-      .gsub(/([A-Z]+)([A-Z][a-z])/, '\1_\2')
-      .gsub(/([a-z\d])([A-Z])/, '\1_\2')
-      .tr('-', '_')
-      .downcase
-  end
-
-  # convert snake_case to CamelCase
-  def camelcase
-    split('_').map(&:capitalize).join
-  end
-end
-
-Array.class_eval do
-  def sum
-    inject(:+)
-  end
-
-  def mean
-    sum.fdiv(count)
-  end
-
-  def variance
-    map { |x| (x - mean)**2}.sum.fdiv(count)
-  end
-
-  def comprehend(&block)
-    return self if block.nil?
-    collect(&block).compact
-  end
-end