rubystats 0.2.3 → 0.4.1

data/test/tc_poisson.rb

@@ -0,0 +1,35 @@
+$:.unshift File.join(File.dirname(__FILE__), "..", "lib")
+require 'minitest/autorun'
+require 'rubystats/poisson_distribution'
+
+class TestPoisson < MiniTest::Unit::TestCase
+  def test_simple
+    rate = 3.5
+    k = 6
+    
+    pois = Rubystats::PoissonDistribution.new(rate)
+    cdf = pois.cdf(k)
+    pmf = pois.pmf(k)
+    pdf = pois.pdf(k)
+    mean = pois.mean
+    inv_cdf = pois.icdf(cdf)
+
+    assert_equal pmf, pdf
+    assert_in_delta(0.0770984, pdf, 0.0000001 )
+    assert_in_delta(0.9347119, cdf, 0.0000001)
+    assert_equal(k,inv_cdf)
+    assert_equal(3.5,mean)
+  end
+  
+  def test_rng
+    rate = 3.5
+    pois = Rubystats::PoissonDistribution.new(rate)
+  
+    rngs = []
+    n = 1000
+    n.times { rngs << pois.rng.to_f }
+    avg = rngs.reduce(:+) / rngs.size    
+    
+    assert_in_delta(avg, pois.mean, 0.5)
+  end
+end

data/test/tc_require_all.rb

@@ -4,17 +4,17 @@ $:.unshift File.join(File.dirname(__FILE__), "..", "lib")
 # test that we can use old api that wasn't namespaced
 #
 
-require 'test/unit'
+require 'minitest/autorun'
 require 'rubystats'
 
-class TestNormal < Test::Unit::TestCase
+class TestNormal < MiniTest::Unit::TestCase
   def test_simple
     norm = NormalDistribution.new(10,2)
     cdf = norm.cdf(11)
 
-    assert_equal("0.691462461274013",cdf.to_s)
-    assert_not_nil(norm.rng)
+    assert_equal("0.6914624612740131",cdf.to_s)
+    refute_nil(norm.rng)
 
-    expd = ExponentialDistribution.new(2)
+    ExponentialDistribution.new(2)
   end
 end

data/test/tc_studentt.rb

@@ -0,0 +1,43 @@
+$:.unshift File.join(File.dirname(__FILE__), '..', 'lib')
+require 'minitest/autorun'
+require 'rubystats/student_t_distribution'
+
+class TestStudentT < MiniTest::Unit::TestCase
+  def test_pdf
+    t = Rubystats::StudentTDistribution.new(8)
+        
+    pdf = t.pdf(0.0)
+    assert_in_epsilon pdf, 0.386699, 0.000001, "pdf"
+  
+	pdf = t.pdf(1.0)
+    assert_in_epsilon pdf, 0.2276076, 0.000001, "pdf"  
+  end
+
+  def test_properties
+    t = Rubystats::StudentTDistribution.new(8)
+    assert_in_epsilon t.mean, 0.0, 0.01, "mean"
+    assert_in_epsilon t.variance, 1.333, 0.01, "variance"
+  end
+  
+  def test_rng
+    t = Rubystats::StudentTDistribution.new(8)
+    
+    total = 100000
+    values = Array.new(total).map{t.rng}
+
+    mean = values.inject(0.0) {|sum, v| sum + v} / values.size 
+    assert_in_delta mean, 0.0, 0.1, "rng mean"    
+	
+	var = values.inject(0.0) {|sum, v| sum + (v - mean)**2} / values.size 
+    assert_in_delta var, 8.0/(8.0-2.0), 0.1, "rng mean"    
+  end
+  
+  def test_integer_input
+    dfi, dff = 8, 8.0
+	xi,xf = 1, 1.0
+    dti = Rubystats::StudentTDistribution.new(dfi).pdf(xi)
+    dtf = Rubystats::StudentTDistribution.new(dff).pdf(xf)
+    assert_in_delta dti, dtf, 0.000001   
+  end
+  
+end

data/test/tc_unif.rb

@@ -0,0 +1,48 @@
+$:.unshift File.join(File.dirname(__FILE__), '..', 'lib')
+require 'minitest/autorun'
+require 'rubystats/uniform_distribution'
+
+class TestUnif < MiniTest::Unit::TestCase
+  def test_cdf
+
+    unif = Rubystats::UniformDistribution.new(1,6)
+    cdf = unif.cdf(2.5)
+
+    assert_equal 0.3, cdf
+  end
+
+  def test_pdf
+    unif = Rubystats::UniformDistribution.new(1.0, 6.0)
+    
+    pdf = unif.pdf(2.5)
+    assert_equal 0.2, pdf
+
+    pdf_out_of_lower_bounds = unif.pdf(0.0)
+    assert_equal 0.0, pdf_out_of_lower_bounds
+ 
+    pdf_out_of_upper_bounds = unif.pdf(8.0)
+    assert_equal 0.0, pdf_out_of_upper_bounds
+  end
+  
+  def test_rng
+    lower, upper = 1,6
+    unif = Rubystats::UniformDistribution.new(lower,upper)
+    
+    total = 10000
+    values = Array.new(total).map{unif.rng}
+    
+    assert_equal values.min >= lower, true
+    assert_equal values.max <= upper, true
+  end
+  
+  def test_integer_input
+    loweri, upperi = 1, 6
+    lowerf, upperf = 1.0, 6.0
+    xi = 3
+    xf = 3.0
+    dunifi = Rubystats::NormalDistribution.new(loweri,upperi).pdf(xi)
+    duniff = Rubystats::NormalDistribution.new(lowerf,upperf).pdf(xf)
+    assert_in_delta dunifi, duniff, 0.00000001   
+  end
+  
+end

data/test/tc_weibull.rb

@@ -0,0 +1,51 @@
+$:.unshift File.join(File.dirname(__FILE__), "..", "lib")
+require 'minitest/autorun'
+require 'rubystats/weibull_distribution'
+
+class TestWeibull < MiniTest::Unit::TestCase
+  def test_simple
+    scale = 6.0
+    shape = 4.0
+    wbd = Rubystats::WeibullDistribution.new(scale,shape)
+    assert_in_delta(5.438415, wbd.mean, 0.000001)
+    assert_in_delta(2.327813, wbd.variance, 0.000001)
+    assert_in_delta(0.0243884, wbd.pdf(2.0), 0.000001)
+    assert_in_delta(0.2381909, wbd.pdf(5.0), 0.000001)
+    assert_in_delta(0.3826092, wbd.cdf(5.0), 0.000001)
+    # make sure the mean of RNG values is close to the expected mean
+    rngs = []
+    n = 10000
+    n.times {|i| rngs << wbd.rng }
+    avg = rngs.inject(0.0) {|sum,x|sum + x} / rngs.size
+    assert_in_delta(avg, 5.438415, 0.1)
+
+    x_vals = [0.5, 1, 1.5, 2, 2.5, 3, 3.5]
+    
+    p_vals = wbd.pdf(x_vals)
+    c_vals = wbd.cdf(x_vals)
+    
+    expected_pvals = [0.00038578, 0.00308404, 0.01037606, 0.02438840, 0.04679345, 0.07828442, 0.11786168]
+    expected_cvals = [0.00004822, 0.00077131, 0.00389863, 0.01226978, 0.02969111, 0.06058694, 0.10933684]
+
+    x_vals.size.times do |i|
+      assert_in_delta expected_pvals[i], p_vals[i], 0.000001
+      assert_in_delta expected_cvals[i], c_vals[i], 0.000001
+    end
+    
+    x_vals.each do |x|
+      cdf = wbd.cdf(x)
+      inv_cdf = wbd.icdf(cdf)
+      assert_in_delta(x, inv_cdf, 0.00000001)
+    end
+  end
+  
+  def test_integer_input
+    scalei,shapei = 6, 4
+    scalef,shapef = 6.0, 4.0
+    xi = 3
+    xf = 3.0
+    dwbi = Rubystats::WeibullDistribution.new(scalei,shapei).pdf(xi)
+    dwbf = Rubystats::WeibullDistribution.new(scalef,shapef).pdf(xf)
+    assert_in_delta dwbi, dwbf, 0.000001   
+  end
+end

data/test/ts_stats.rb

@@ -1,4 +1,4 @@
-require 'test/unit'
+require 'minitest/autorun'
 require 'tc_fisher'
 require 'tc_binomial'
 require 'tc_beta'

metadata

@@ -1,41 +1,79 @@
---- !ruby/object:Gem::Specification 
+--- !ruby/object:Gem::Specification
 name: rubystats
-version: !ruby/object:Gem::Version 
-  version: 0.2.3
+version: !ruby/object:Gem::Version
+  version: 0.4.1
 platform: ruby
-authors: 
+authors:
+- Ilya Scharrenbroich
 - Bryan Donovan - http://www.bryandonovan.com
+- Phillip Baker
 autorequire: 
 bindir: bin
 cert_chain: []
-
-date: 2008-07-06 00:00:00 -07:00
-default_executable: 
-dependencies: 
-- !ruby/object:Gem::Dependency 
+date: 2023-05-21 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: minitest
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '4.2'
+    - - "<"
+      - !ruby/object:Gem::Version
+        version: '5.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '4.2'
+    - - "<"
+      - !ruby/object:Gem::Version
+        version: '5.0'
+- !ruby/object:Gem::Dependency
   name: hoe
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.7.0
   type: :development
-  version_requirement: 
-  version_requirements: !ruby/object:Gem::Requirement 
-    requirements: 
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
     - - ">="
-      - !ruby/object:Gem::Version 
+      - !ruby/object:Gem::Version
         version: 1.7.0
-    version: 
-description: Ruby Stats is a port of the statistics libraries from PHPMath. Probability distributions include binomial, beta, and normal distributions with PDF, CDF and inverse CDF as well as Fisher's Exact Test.
-email: b.dondo+rubyforge@gmail.com
+- !ruby/object:Gem::Dependency
+  name: matrix
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+description: Ruby Stats is a port of the statistics libraries from PHPMath. Probability
+  distributions include binomial, beta, and normal distributions with PDF, CDF and
+  inverse CDF as well as Fisher's Exact Test.
+email: 
 executables: []
-
 extensions: []
-
-extra_rdoc_files: 
-- History.txt
-- Manifest.txt
-- README.txt
-files: 
+extra_rdoc_files: []
+files:
+- ".github/workflows/test.yml"
+- ".gitignore"
+- Gemfile
 - History.txt
 - Manifest.txt
-- README.txt
+- README.rdoc
 - Rakefile
 - examples/beta.rb
 - examples/binomial.rb
@@ -43,47 +81,77 @@ files:
 - examples/failrate_vs_goal.rb
 - examples/fisher.rb
 - examples/norm.rb
+- examples/uniform.rb
 - lib/rubystats.rb
 - lib/rubystats/beta_distribution.rb
 - lib/rubystats/binomial_distribution.rb
+- lib/rubystats/cauchy_distribution.rb
 - lib/rubystats/exponential_distribution.rb
 - lib/rubystats/fishers_exact_test.rb
+- lib/rubystats/gamma_distribution.rb
+- lib/rubystats/lognormal_distribution.rb
 - lib/rubystats/modules.rb
+- lib/rubystats/multivariate_normal_distribution.rb
 - lib/rubystats/normal_distribution.rb
+- lib/rubystats/poisson_distribution.rb
 - lib/rubystats/probability_distribution.rb
+- lib/rubystats/student_t_distribution.rb
+- lib/rubystats/uniform_distribution.rb
+- lib/rubystats/version.rb
+- lib/rubystats/weibull_distribution.rb
+- rubystats.gemspec
 - test/tc_beta.rb
 - test/tc_binomial.rb
+- test/tc_cauchy.rb
 - test/tc_exponential.rb
 - test/tc_fisher.rb
+- test/tc_gamma.rb
+- test/tc_lnorm.rb
+- test/tc_multivariate_normal.rb
 - test/tc_norm.rb
+- test/tc_poisson.rb
 - test/tc_require_all.rb
+- test/tc_studentt.rb
+- test/tc_unif.rb
+- test/tc_weibull.rb
 - test/ts_stats.rb
-has_rdoc: true
-homepage: http://rubyforge.org/projects/rubystats/
+homepage: https://github.com/phillbaker/rubystats
+licenses:
+- MIT
+metadata: {}
 post_install_message: 
-rdoc_options: 
-- --main
-- README.txt
-require_paths: 
+rdoc_options: []
+require_paths:
 - lib
-required_ruby_version: !ruby/object:Gem::Requirement 
-  requirements: 
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
   - - ">="
-    - !ruby/object:Gem::Version 
-      version: "0"
-  version: 
-required_rubygems_version: !ruby/object:Gem::Requirement 
-  requirements: 
+    - !ruby/object:Gem::Version
+      version: 1.9.3
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
   - - ">="
-    - !ruby/object:Gem::Version 
-      version: "0"
-  version: 
+    - !ruby/object:Gem::Version
+      version: '0'
 requirements: []
-
-rubyforge_project: rubystats
-rubygems_version: 1.2.0
+rubyforge_project: 
+rubygems_version: 2.4.3
 signing_key: 
-specification_version: 2
-summary: Port of PHPMath to Ruby
-test_files: []
-
+specification_version: 4
+summary: ''
+test_files:
+- test/tc_beta.rb
+- test/tc_binomial.rb
+- test/tc_cauchy.rb
+- test/tc_exponential.rb
+- test/tc_fisher.rb
+- test/tc_gamma.rb
+- test/tc_lnorm.rb
+- test/tc_multivariate_normal.rb
+- test/tc_norm.rb
+- test/tc_poisson.rb
+- test/tc_require_all.rb
+- test/tc_studentt.rb
+- test/tc_unif.rb
+- test/tc_weibull.rb
+- test/ts_stats.rb