statistics2 0.54

data/lib/statistics2/no_ext.rb

@@ -0,0 +1,8 @@
+$:.unshift File.join(File.dirname(__FILE__), '..')
+
+# Indicate that we should skip loading of the C extension.
+module Statistics2
+  NO_EXT = true
+end
+
+require 'statistics2'

data/lib/statistics2/version.rb

@@ -0,0 +1,3 @@
+module Statistics2
+  VERSION = 0.54
+end

data/statistics2.gemspec

@@ -0,0 +1,35 @@
+# -*- encoding: utf-8 -*-
+
+Gem::Specification.new do |s|
+  s.name = %q{statistics2}
+  s.version = "0.54"
+
+  s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
+  s.authors = ["Brendan Ribera", "Shin-ichiro Hara"]
+  s.date = %q{2010-01-26}
+  s.description = %q{Statistics2 is a module that provides normal, Chi-square, t- and F- probability distributions for Ruby. It is a fork/continuation of Shin-ichiro Hara's original code. It can provides a native, compiled extension and a pure Ruby implementation.}
+  s.email = ["brendan.ribera@gmail.com", "sinara@blade.nagaokaut.ac.jp"]
+  s.extensions = ["ext/extconf.rb"]
+  s.extra_rdoc_files = ["Manifest.txt"]
+  s.files = ["History.rdoc", "Manifest.txt", "README.rdoc", "Rakefile", "examples/mklist.rb", "examples/show.rb", "ext/extconf.rb", "ext/statistics2.c", "lib/statistics2.rb", "lib/statistics2/no_ext.rb", "lib/statistics2/version.rb", "statistics2.gemspec", "test/sample_tbl.rb", "test/test_ext.rb", "test/test_inv.rb"]
+  s.homepage = %q{http://github.com/abscondment/statistics2}
+  s.rdoc_options = ["--main", "README.rdoc"]
+  s.require_paths = ["lib", "ext"]
+  s.rubyforge_project = %q{statistics2}
+  s.rubygems_version = %q{1.3.5}
+  s.summary = %q{Statistical Distributions for Ruby. Based on Shin-ichiro Hara's original library, updated for Ruby 1.9}
+  s.test_files = ["test/test_ext.rb", "test/test_inv.rb"]
+
+  if s.respond_to? :specification_version then
+    current_version = Gem::Specification::CURRENT_SPECIFICATION_VERSION
+    s.specification_version = 3
+
+    if Gem::Version.new(Gem::RubyGemsVersion) >= Gem::Version.new('1.2.0') then
+      s.add_development_dependency(%q<hoe>, [">= 2.4.0"])
+    else
+      s.add_dependency(%q<hoe>, [">= 2.4.0"])
+    end
+  else
+    s.add_dependency(%q<hoe>, [">= 2.4.0"])
+  end
+end

data/test/sample_tbl.rb

@@ -0,0 +1,134 @@
+#!/usr/bin/env ruby
+$test ||= false
+if $test
+  def nop(*a); ""; end
+  alias _sprintf nop
+  alias _printf nop
+  alias _puts nop
+else
+  alias _sprintf sprintf
+  alias _printf printf
+  alias _puts puts
+end
+
+def norm_line(x, n)
+  s = sprintf("%1.1f|", x)
+  (0.00).step(0.09, 0.01) do |y|
+    s << _sprintf(" %1.#{n}f", yield(x + y))
+  end
+  s
+end
+
+def norm_tbl(s, e, ln = nil, tn = nil, &b)
+  n = 4
+  unless ln
+    _printf("    " + (" "*(n-1) + "%1.2f")*10 + "\n", * (0..9).map {|x| x*0.01})
+  end
+  i = 0
+  (s.to_f).step(e.to_f, 0.1) do |x|
+    next if ln && ln.to_i != i+1
+    _puts norm_line(x, n, &b)
+    i += 1
+  end
+end
+
+def chi2_tbl(ln = nil, tn = nil)
+  pers = [0.995, 0.99, 0.975, 0.95, 0.05, 0.025, 0.01, 0.005]
+  arbi = (1..30).to_a + [40, 60, 80, 100]
+  form = "  %7.5f"
+  unless ln
+    _printf("     "); pers.each do |a|; _printf(form, a); end;  _puts
+  end
+  arbi.each_with_index do |n, i|
+    next if ln && ln.to_i != i+1
+    _printf("%4d|", n) unless tn
+    pers.each_with_index do |a, j|
+      next if tn && tn.to_i != j+1
+      form = case n
+             when 1; a >= 0.95 ? " %.4e" : " %6.3f"
+             when 2..5; a >= 0.95 ? "  %7.5f" : "  %7.3f"
+             when 6..24; "  %7.3f"
+             else; "  %7.2f" # 26..100
+             end
+      _printf(form, yield(n, a))
+    end
+    _puts
+  end
+end
+
+def t_tbl(ln = nil, tn = nil)
+  pers = [0.5, 0.4, 0.3, 0.2, 0.1, 0.05, 0.02, 0.01, 0.001]
+  arbi = (1..30).to_a + [40, 60, 120]#, 9999]
+  form = " %7.3f"
+  unless ln
+    _printf("     "); pers.each do |a|; _printf(form, a); end;  _puts
+  end
+  arbi.each_with_index do |n, i|
+    next if ln && ln.to_i != i+1
+    _printf("%4d|", n) unless tn
+    pers.each_with_index do |a, j|
+      next if tn && tn.to_i != j+1
+      _printf(form, yield(n, a))
+    end
+    _puts
+  end
+end
+
+def f_tbl(a, k1 = nil, k2 = nil)
+  arbi1 = (1..10).to_a + [12, 15, 20, 24, 30, 40, 60, 120]#, 9999]
+  arbi2 = (1..30).to_a + [40, 60, 120]#, 9999]
+  unless k1
+    _printf("     "); arbi1.each do |n1|; _printf(" %4d", n1); end;  _puts
+  end
+  arbi2.each do |n2|
+    form = n2 == 1 ? " %4d" : n2 == 2 ? " %4.1f" : " %4.2f"
+    next if k2 && k2.to_i != n2
+    _printf("%4d|", n2) unless k2
+    arbi1.each do |n1|
+      next if k1 && k1.to_i != n1
+      _printf(form, yield(n1, n2, a.to_f))
+    end
+    _puts
+  end
+end
+
+def show_tbl(mod, dist, *opts)
+  case dist
+  when nil, "norm"
+    norm_tbl(0.0, 3.1, *opts) do |x|
+      mod.normal___x(x)
+    end
+  when "chi", "chi2"
+    chi2_tbl(*opts) do |n, x|
+      mod.pchi2_x(n, x)
+    end
+  when "t"
+    t_tbl(*opts) do |n, x|
+      mod.ptx__x(n, x)
+    end
+  when "f"
+    f_tbl(*opts) do |n1, n2, a|
+      mod.pf_x(n1, n2, a)
+    end
+  end
+end
+
+
+if $0 == __FILE__
+  if ARGV.empty?
+    puts "-- This script makes tables of distributions"
+    puts "Example:"
+    puts "      #$0 norm"
+    puts "      #$0 chi2"
+    puts "      #$0 t"
+    puts "      #$0 f 0.01"
+    exit
+  end
+  
+  $:.unshift File.dirname(__FILE__)
+  $:.unshift File.join(File.dirname(__FILE__), '..', 'lib')
+  $:.unshift File.join(File.dirname(__FILE__), '..', 'ext')
+
+  require 'statistics2'
+  show_tbl(Statistics2, *ARGV)
+end

data/test/test_ext.rb

@@ -0,0 +1,53 @@
+$:.unshift File.dirname(__FILE__)
+$:.unshift File.join(File.dirname(__FILE__), '..', 'ext')
+$:.unshift File.join(File.dirname(__FILE__), '..', 'lib')
+
+$test = true
+
+require 'test/unit' unless defined?(Hoe)
+require 'sample_tbl'
+
+require 'statistics2.so'
+$mod = Statistics2
+
+
+# Request extension to be skipped.
+module Statistics20
+  NO_EXT = true
+end
+eval(File.read('lib/statistics2.rb').gsub(/Statistics2/, 'Statistics20'))
+$mod0 = Statistics20
+
+class T_Statistics2 < Test::Unit::TestCase
+
+  def test_normal
+    norm_tbl(0.0, 3.1) do |x|
+      a, b = $mod.normal___x(x), $mod0.normal___x(x)
+      assert_in_delta a, b, 0.000001
+    end
+  end
+
+  def test_chi
+    chi2_tbl() do |n, x|
+      a, b = $mod.pchi2_x(n, x), $mod0.pchi2_x(n, x)
+      assert_in_delta a/b, 1.0, 0.01
+    end
+  end
+
+  def test_t
+    t_tbl() do |n, x|
+      a, b = $mod.ptx__x(n, x), $mod.ptx__x(n, x)
+      assert_in_delta a, b, 0.001
+    end
+  end
+
+  def test_f
+    [0.01, 0.025, 0.05].each do |opt|
+      f_tbl(opt) do |n1, n2, y|
+        a, b = $mod.pf_x(n1, n2, y), $mod0.pf_x(n1, n2, y)
+        assert_in_delta a/b, 1.0, 0.01
+      end
+    end
+  end
+  
+end

data/test/test_inv.rb

@@ -0,0 +1,57 @@
+$:.unshift File.dirname(__FILE__)
+$:.unshift File.join(File.dirname(__FILE__), '..', 'ext')
+$:.unshift File.join(File.dirname(__FILE__), '..', 'lib')
+
+$test = true
+
+require 'test/unit' unless defined?(Hoe)
+require 'sample_tbl'
+require 'lib/statistics2'
+
+class T_Statistics2 < Test::Unit::TestCase
+  $mod = Statistics2
+  
+  def inv(s); "p" + s; end
+
+  def test_inv_normal
+    delta = 1.0e-6
+    meth = "normal___x"
+    norm_tbl(0.00001, 4.0) do |x|
+      pr = $mod.send(meth, x)
+      x0 = $mod.send(inv(meth), pr)
+      assert_in_delta x0/x, 1.0, delta
+    end
+  end
+
+  def test_inv_chi
+    delta = 1.0e-4
+    meth = "chi2_x"
+    chi2_tbl() do |n, pr|
+      x = $mod.send(inv(meth), n, pr)
+      pr0 = $mod.send(meth, n, x)
+      assert_in_delta pr0/pr, 1.0, delta
+    end
+  end
+
+  def test_inv_t
+    delta = 1.0e-3
+    meth = "tx__x"
+    t_tbl() do |n, pr|
+      x = $mod.send(inv(meth), n, pr)
+      pr0 = $mod.send(meth, n, x)
+      assert_in_delta pr0/pr, 1.0, delta
+    end
+  end
+
+  def test_inv_f
+    delta = 1.0e-3
+    meth = "f_x"
+    [0.01, 0.025, 0.05].each do |opt|
+      f_tbl(opt) do |n1, n2, pr|
+        x = $mod.send(inv(meth), n1, n2, pr)
+        pr0 = $mod.send(meth, n1, n2, x)
+        assert_in_delta pr0/pr, 1.0, delta
+      end
+    end
+  end
+end

metadata

@@ -0,0 +1,84 @@
+--- !ruby/object:Gem::Specification 
+name: statistics2
+version: !ruby/object:Gem::Version 
+  version: "0.54"
+platform: ruby
+authors: 
+- Brendan Ribera
+- Shin-ichiro Hara
+autorequire: 
+bindir: bin
+cert_chain: []
+
+date: 2010-01-26 00:00:00 -08:00
+default_executable: 
+dependencies: 
+- !ruby/object:Gem::Dependency 
+  name: hoe
+  type: :development
+  version_requirement: 
+  version_requirements: !ruby/object:Gem::Requirement 
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        version: 2.4.0
+    version: 
+description: Statistics2 is a module that provides normal, Chi-square, t- and F- probability distributions for Ruby. It is a fork/continuation of Shin-ichiro Hara's original code. It provides a native, compiled extension and a pure Ruby implementation.
+email: 
+- brendan.ribera@gmail.com
+- sinara@blade.nagaokaut.ac.jp
+executables: []
+
+extensions: 
+- ext/extconf.rb
+extra_rdoc_files: 
+- Manifest.txt
+files: 
+- History.rdoc
+- Manifest.txt
+- README.rdoc
+- Rakefile
+- examples/mklist.rb
+- examples/show.rb
+- ext/extconf.rb
+- ext/statistics2.c
+- lib/statistics2.rb
+- lib/statistics2/no_ext.rb
+- lib/statistics2/version.rb
+- statistics2.gemspec
+- test/sample_tbl.rb
+- test/test_ext.rb
+- test/test_inv.rb
+has_rdoc: true
+homepage: http://github.com/abscondment/statistics2
+licenses: []
+
+post_install_message: 
+rdoc_options: 
+- --main
+- README.rdoc
+require_paths: 
+- lib
+- ext
+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: "0"
+  version: 
+requirements: []
+
+rubyforge_project: statistics2
+rubygems_version: 1.3.5
+signing_key: 
+specification_version: 3
+summary: Statistical Distributions for Ruby. Based on Shin-ichiro Hara's original library, updated for Ruby 1.9
+test_files: 
+- test/test_ext.rb
+- test/test_inv.rb