statistics3 0.0.1

data/statistics3.gemspec

@@ -0,0 +1,117 @@
+# Generated by juwelier
+# DO NOT EDIT THIS FILE DIRECTLY
+# Instead, edit Juwelier::Tasks in Rakefile, and run 'rake gemspec'
+# -*- encoding: utf-8 -*-
+# stub: statistics3 0.0.1 ruby lib
+# stub: ext/extconf.rb
+
+Gem::Specification.new do |s|
+  s.name = "statistics3".freeze
+  s.version = "0.0.1"
+
+  s.required_rubygems_version = Gem::Requirement.new(">= 0".freeze) if s.respond_to? :required_rubygems_version=
+  s.require_paths = ["lib".freeze]
+  s.authors = ["Fred Mitchell".freeze]
+  s.date = "2016-12-12"
+  s.description = "\n  Statistics3 is a module that provides normal,\n  Chi-square, t- and F- probability distributions\n  for Ruby. It is a fork/continuation of Shin-ichiro Hara's original code.\n  It provides a native, compiled extension and a pure Ruby implementation.\n  ".freeze
+  s.email = "fred.mitchell@gmx.de".freeze
+  s.extensions = ["ext/extconf.rb".freeze]
+  s.extra_rdoc_files = [
+    "LICENSE.txt",
+    "README.org"
+  ]
+  s.files = [
+    ".document",
+    ".rspec",
+    ".ruby-version",
+    ".semver",
+    ".travis.yml",
+    "Gemfile",
+    "Gemfile.lock",
+    "History.rdoc",
+    "LICENSE.txt",
+    "README.org",
+    "Rakefile",
+    "examples/mklist.rb",
+    "examples/show.rb",
+    "ext/Makefile",
+    "ext/_statistics3.c",
+    "ext/_statistics3.o",
+    "ext/_statistics3.so",
+    "ext/extconf.rb",
+    "lib/statistics3.rb",
+    "lib/statistics3/base.rb",
+    "lib/statistics3/no_ext.rb",
+    "spec/spec_helper.rb",
+    "statistics3.gemspec",
+    "test/sample_tbl.rb",
+    "test/test_ext.rb",
+    "test/test_inv.rb"
+  ]
+  s.homepage = "http://github.com/flajann2/statistics3".freeze
+  s.licenses = ["MIT".freeze]
+  s.required_ruby_version = Gem::Requirement.new(">= 2.2.2".freeze)
+  s.rubygems_version = "2.6.7".freeze
+  s.summary = "Statistics3 library".freeze
+
+  if s.respond_to? :specification_version then
+    s.specification_version = 4
+
+    if Gem::Version.new(Gem::VERSION) >= Gem::Version.new('1.2.0') then
+      s.add_runtime_dependency(%q<semver>.freeze, ["~> 1"])
+      s.add_runtime_dependency(%q<awesome_print>.freeze, ["~> 1"])
+      s.add_runtime_dependency(%q<colorize>.freeze, ["~> 0"])
+      s.add_runtime_dependency(%q<psych>.freeze, ["~> 2"])
+      s.add_development_dependency(%q<rspec>.freeze, ["~> 3.5.0"])
+      s.add_development_dependency(%q<test-unit>.freeze, [">= 0"])
+      s.add_development_dependency(%q<yard>.freeze, ["~> 0.7"])
+      s.add_development_dependency(%q<rdoc>.freeze, ["~> 3.12"])
+      s.add_development_dependency(%q<bundler>.freeze, ["~> 1.0"])
+      s.add_development_dependency(%q<juwelier>.freeze, ["~> 2.1.0"])
+      s.add_development_dependency(%q<simplecov>.freeze, [">= 0"])
+      s.add_development_dependency(%q<pry>.freeze, ["~> 0"])
+      s.add_development_dependency(%q<pry-byebug>.freeze, ["~> 3"])
+      s.add_development_dependency(%q<pry-doc>.freeze, ["~> 0"])
+      s.add_development_dependency(%q<pry-remote>.freeze, ["~> 0"])
+      s.add_development_dependency(%q<pry-rescue>.freeze, ["~> 1"])
+      s.add_development_dependency(%q<pry-stack_explorer>.freeze, ["~> 0"])
+    else
+      s.add_dependency(%q<semver>.freeze, ["~> 1"])
+      s.add_dependency(%q<awesome_print>.freeze, ["~> 1"])
+      s.add_dependency(%q<colorize>.freeze, ["~> 0"])
+      s.add_dependency(%q<psych>.freeze, ["~> 2"])
+      s.add_dependency(%q<rspec>.freeze, ["~> 3.5.0"])
+      s.add_dependency(%q<test-unit>.freeze, [">= 0"])
+      s.add_dependency(%q<yard>.freeze, ["~> 0.7"])
+      s.add_dependency(%q<rdoc>.freeze, ["~> 3.12"])
+      s.add_dependency(%q<bundler>.freeze, ["~> 1.0"])
+      s.add_dependency(%q<juwelier>.freeze, ["~> 2.1.0"])
+      s.add_dependency(%q<simplecov>.freeze, [">= 0"])
+      s.add_dependency(%q<pry>.freeze, ["~> 0"])
+      s.add_dependency(%q<pry-byebug>.freeze, ["~> 3"])
+      s.add_dependency(%q<pry-doc>.freeze, ["~> 0"])
+      s.add_dependency(%q<pry-remote>.freeze, ["~> 0"])
+      s.add_dependency(%q<pry-rescue>.freeze, ["~> 1"])
+      s.add_dependency(%q<pry-stack_explorer>.freeze, ["~> 0"])
+    end
+  else
+    s.add_dependency(%q<semver>.freeze, ["~> 1"])
+    s.add_dependency(%q<awesome_print>.freeze, ["~> 1"])
+    s.add_dependency(%q<colorize>.freeze, ["~> 0"])
+    s.add_dependency(%q<psych>.freeze, ["~> 2"])
+    s.add_dependency(%q<rspec>.freeze, ["~> 3.5.0"])
+    s.add_dependency(%q<test-unit>.freeze, [">= 0"])
+    s.add_dependency(%q<yard>.freeze, ["~> 0.7"])
+    s.add_dependency(%q<rdoc>.freeze, ["~> 3.12"])
+    s.add_dependency(%q<bundler>.freeze, ["~> 1.0"])
+    s.add_dependency(%q<juwelier>.freeze, ["~> 2.1.0"])
+    s.add_dependency(%q<simplecov>.freeze, [">= 0"])
+    s.add_dependency(%q<pry>.freeze, ["~> 0"])
+    s.add_dependency(%q<pry-byebug>.freeze, ["~> 3"])
+    s.add_dependency(%q<pry-doc>.freeze, ["~> 0"])
+    s.add_dependency(%q<pry-remote>.freeze, ["~> 0"])
+    s.add_dependency(%q<pry-rescue>.freeze, ["~> 1"])
+    s.add_dependency(%q<pry-stack_explorer>.freeze, ["~> 0"])
+  end
+end
+

data/test/sample_tbl.rb

@@ -0,0 +1,129 @@
+#!/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
+
+  require 'statistics3'
+  show_tbl(Statistics3, *ARGV)
+end

data/test/test_ext.rb

@@ -0,0 +1,46 @@
+$test = true
+
+require 'test/unit'
+require 'sample_tbl'
+
+require 'statistics3'
+
+class TestStatistics3 < Test::Unit::TestCase
+  module Base
+    extend Statistics3::Base
+  end
+
+  module Extension
+    extend Statistics3::Extension
+  end
+
+  def test_normal
+    norm_tbl(0.0, 3.1) do |x|
+      a, b = Base.normal___x(x), Extension.normal___x(x)
+      assert_in_delta a, b, 0.000001
+    end
+  end
+
+  def test_chi
+    chi2_tbl() do |n, x|
+      a, b = Base.pchi2_x(n, x), Extension.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 = Base.ptx__x(n, x), Extension.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 = Base.pf_x(n1, n2, y), Extension.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,53 @@
+$test = true
+
+require 'test/unit'
+require 'sample_tbl'
+require 'statistics3'
+
+class TestStatistics3Inverted < Test::Unit::TestCase
+  $mod = Statistics3
+
+  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,313 @@
+--- !ruby/object:Gem::Specification
+name: statistics3
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- Fred Mitchell
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2016-12-12 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: semver
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1'
+- !ruby/object:Gem::Dependency
+  name: awesome_print
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1'
+- !ruby/object:Gem::Dependency
+  name: colorize
+  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'
+- !ruby/object:Gem::Dependency
+  name: psych
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '2'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '2'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 3.5.0
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 3.5.0
+- !ruby/object:Gem::Dependency
+  name: test-unit
+  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'
+- !ruby/object:Gem::Dependency
+  name: yard
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.7'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.7'
+- !ruby/object:Gem::Dependency
+  name: rdoc
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.12'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.12'
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.0'
+- !ruby/object:Gem::Dependency
+  name: juwelier
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 2.1.0
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 2.1.0
+- !ruby/object:Gem::Dependency
+  name: simplecov
+  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'
+- !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'
+- !ruby/object:Gem::Dependency
+  name: pry-byebug
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3'
+- !ruby/object:Gem::Dependency
+  name: pry-doc
+  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'
+- !ruby/object:Gem::Dependency
+  name: pry-remote
+  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'
+- !ruby/object:Gem::Dependency
+  name: pry-rescue
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1'
+- !ruby/object:Gem::Dependency
+  name: pry-stack_explorer
+  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: "\n  Statistics3 is a module that provides normal,\n  Chi-square, t-
+  and F- probability distributions\n  for Ruby. It is a fork/continuation of Shin-ichiro
+  Hara's original code.\n  It provides a native, compiled extension and a pure Ruby
+  implementation.\n  "
+email: fred.mitchell@gmx.de
+executables: []
+extensions:
+- ext/extconf.rb
+extra_rdoc_files:
+- LICENSE.txt
+- README.org
+files:
+- ".document"
+- ".rspec"
+- ".ruby-version"
+- ".semver"
+- ".travis.yml"
+- Gemfile
+- Gemfile.lock
+- History.rdoc
+- LICENSE.txt
+- README.org
+- Rakefile
+- examples/mklist.rb
+- examples/show.rb
+- ext/Makefile
+- ext/_statistics3.c
+- ext/_statistics3.o
+- ext/_statistics3.so
+- ext/extconf.rb
+- lib/statistics3.rb
+- lib/statistics3/base.rb
+- lib/statistics3/no_ext.rb
+- spec/spec_helper.rb
+- statistics3.gemspec
+- test/sample_tbl.rb
+- test/test_ext.rb
+- test/test_inv.rb
+homepage: http://github.com/flajann2/statistics3
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: 2.2.2
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.6.7
+signing_key: 
+specification_version: 4
+summary: Statistics3 library
+test_files: []