alglib 1.0.0 → 1.0.1

data/test/test_alglib.rb

@@ -1,16 +1,52 @@
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../lib"))
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../ext"))
+
 require "test/unit"
+begin
+    require 'statsample'
+    HAS_SS=true
+rescue
+    HAS_SS=false
+end
+
+if(HAS_GSL.nil?)
+begin
+    require 'rbgsl'
+    HAS_GSL=true
+rescue
+    HAS_GSL=false
+end
+end
+
 require "alglib"
 require 'matrix'
 class TestAlglib < Test::Unit::TestCase
     def test_linear_regression
-a=[1,3,2,4,3,5,4,6,5,7]
-b=[3,3,4,4,5,5,6,6,4,4]
-c=[11,22,30,40,50,65,78,79,99,100]
-y=[3,4,5,6,7,8,9,10,20,30]
-matrix=Matrix.columns([a,b,c,y])
-lr=Alglib::LinearRegression.build_from_matrix(matrix)
-assert_in_delta(0.695,lr.coeffs[0],0.001)
-assert_in_delta(11.027,lr.constant,0.001)
-assert_in_delta(1.785,lr.process([1,3,11]),0.001)
+        a=[1,3,2,4,3,5,4,6,5,7]
+        b=[3,3,4,4,5,5,6,6,4,4]
+        c=[11,22,30,40,50,65,78,79,99,100]
+        y=[3,4,5,6,7,8,9,10,20,30]
+        matrix=Matrix.columns([a,b,c,y])
+        lr=Alglib::LinearRegression.build_from_matrix(matrix)
+        assert_in_delta(0.695,lr.coeffs[0],0.001)
+        assert_in_delta(11.027,lr.constant,0.001)
+        assert_in_delta(1.785,lr.process([1,3,11]),0.001)
+  end
+  def test_array_to_real1darray
+      v=[1,2,3]
+      r1d=Alglib.array_to_real1darray(v)
+      assert_equal(3,r1d.size)
+      assert_equal(1,r1d[0])
+      assert_equal(2,r1d[1])
+      assert_equal(3,r1d[2])
   end
+  def test_real2darray_to_array
+      m=Matrix[[0,3,4],[1,2,5],[3,5,6]]
+      v2=m.to_alglib_matrix
+      m2=Alglib.real2darray_to_array(v2)
+      assert_equal(3,m2.size)
+      assert_equal([0,3,4],m2[0])
+      assert_equal([1,2,5],m2[1])
+      assert_equal([3,5,6],m2[2])
+      end
 end

data/test/test_correlation.rb

@@ -0,0 +1,44 @@
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../lib"))
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../ext"))
+
+require "test/unit"
+begin
+    require 'statsample'
+    HAS_SS=true
+rescue
+    HAS_SS=false
+end
+
+if(HAS_GSL.nil?)
+begin
+    require 'rbgsl'
+    HAS_GSL=true
+rescue
+    HAS_GSL=false
+end
+end
+
+require "alglib"
+require 'matrix'
+class TestAlglibCorrelation < Test::Unit::TestCase
+  def test_correlation
+      a=[1,2,3,6,7,9,0,10]
+      b=[2,5,  7,8,4,8, 0,  5]
+
+      if(HAS_GSL)
+          a_gsl=GSL::Vector.alloc(a)
+          b_gsl=GSL::Vector.alloc(b)
+          corr_expected=GSL::Stats::correlation(a_gsl, b_gsl)
+          corr_alglib=Alglib.pearson_correlation(a,b)
+          assert_equal(corr_expected,corr_alglib)
+      end
+      if(HAS_SS)
+          va=a.to_vector(:scale)
+          vb=b.to_vector(:scale)
+          corr_expected=Statsample::Bivariate::spearman(va, vb)
+          corr_alglib=Alglib.spearman_correlation(a,b)
+          assert_in_delta(corr_expected,corr_alglib,0.001)
+      end
+  end
+
+end

data/test/test_correlationtest.rb

@@ -0,0 +1,45 @@
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../lib"))
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../ext"))
+
+require "test/unit"
+begin
+    require 'statsample'
+    HAS_SS=true
+rescue
+    HAS_SS=false
+end
+
+if(HAS_GSL.nil?)
+begin
+    require 'rbgsl'
+    HAS_GSL=true
+rescue
+    HAS_GSL=false
+end
+end
+
+require "alglib"
+require 'matrix'
+class TestAlglibCorrelationTests < Test::Unit::TestCase
+    def test_pearson
+      a=[2,3,4,5,6,7,8,9 ,10,22,20]
+      b=[1,4,3,5,6,8,9,10,10,20,30]
+      r=Alglib.pearson_correlation(a,b)
+      t= Statsample::Bivariate.t_r(r, a.size)
+      out  = Alglib.pearson_correlation_significance(r, a.size)
+      assert_in_delta(out[:both], Statsample::Bivariate.prop_pearson(t,a.size), 0.0001)
+      assert_in_delta(out[:left], Statsample::Bivariate.prop_pearson(t, a.size,:left), 0.0001)
+      assert_in_delta(out[:right], Statsample::Bivariate.prop_pearson(t,a.size, :right), 0.0001)
+    end
+    def test_spearman
+      a=[2,3,4,5,6,7,8,9 ,10,22,20]
+      b=[1,4,3,5,6,8,9,10,10,20,30]
+      r=Alglib.spearman_correlation(a,b)
+      t= Statsample::Bivariate.t_r(r, a.size)
+      out  = Alglib.spearman_correlation_significance(r, a.size)
+      
+      assert_in_delta(out[:both], Statsample::Bivariate.prop_pearson(t,a.size), 0.0001)
+      assert_in_delta(out[:left], Statsample::Bivariate.prop_pearson(t, a.size,:left), 0.0001)
+      assert_in_delta(out[:right], Statsample::Bivariate.prop_pearson(t,a.size, :right), 0.0001)
+    end
+end

data/test/test_linreg.rb

@@ -0,0 +1,35 @@
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../lib"))
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../ext"))
+
+require "test/unit"
+begin
+    require 'statsample'
+    HAS_SS=true
+rescue
+    HAS_SS=false
+end
+
+if(HAS_GSL.nil?)
+begin
+    require 'rbgsl'
+    HAS_GSL=true
+rescue
+    HAS_GSL=false
+end
+end
+
+require "alglib"
+require 'matrix'
+class TestAlglibLinreg < Test::Unit::TestCase
+    def test_linear_regression
+        a=[1,3,2,4,3,5,4,6,5,7]
+        b=[3,3,4,4,5,5,6,6,4,4]
+        c=[11,22,30,40,50,65,78,79,99,100]
+        y=[3,4,5,6,7,8,9,10,20,30]
+        matrix=Matrix.columns([a,b,c,y])
+        lr=Alglib::LinearRegression.build_from_matrix(matrix)
+        assert_in_delta(0.695,lr.coeffs[0],0.001)
+        assert_in_delta(11.027,lr.constant,0.001)
+        assert_in_delta(1.785,lr.process([1,3,11]),0.001)
+  end
+end

data/test/test_logit.rb

@@ -0,0 +1,43 @@
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../lib"))
+$:.unshift(File.expand_path(File.dirname(__FILE__)+"/../ext"))
+
+require "test/unit"
+begin
+    require 'statsample'
+    HAS_SS=true
+rescue
+    HAS_SS=false
+end
+
+if(HAS_GSL.nil?)
+begin
+    require 'rbgsl'
+    HAS_GSL=true
+rescue
+    HAS_GSL=false
+end
+end
+
+require "alglib"
+require 'matrix'
+class TestAlglib < Test::Unit::TestCase
+  def test_logit
+        a=[1,2,3,4,5,6,7,8,9,10]
+        b=[0,0,1,0,0,0,1,1,1,1]
+        matrix=Matrix.columns([a,b])
+        lr=Alglib::Logit.build_from_matrix(matrix)
+  end
+  def atest_logit_1
+      ds=Statsample::CSV.read(File.dirname(__FILE__)+"/sample_logit.csv")
+      ds['female'].recode! {|c|
+          c=='female' ? 1: 0
+      }
+      constant=[1]*(ds.cases)
+      matrix=Matrix.columns([ds['female'].data,ds['hon'].data])
+      lr=Alglib::Logit.build_from_matrix(matrix)
+      lr.unpack
+      matrix=Matrix.columns([ds['math'].data, ds['female'].data,ds['read'].data, ds['hon'].data])
+      lr=Alglib::Logit.build_from_matrix(matrix)
+
+  end
+end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification 
 name: alglib
 version: !ruby/object:Gem::Version 
-  version: 1.0.0
+  version: 1.0.1
 platform: ruby
 authors: 
 - Claudio Bustos
@@ -9,7 +9,7 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2009-08-03 00:00:00 -04:00
+date: 2009-08-20 00:00:00 -04:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -20,7 +20,7 @@ dependencies:
     requirements: 
     - - ">="
       - !ruby/object:Gem::Version 
-        version: 2.3.2
+        version: 2.3.3
     version: 
 description: ALGLIB - is a multilingual collection of algorithms designed to solve problems in the field of numeric analysis and data processing. ALGLIB is distributed under a 3-clause BSD license. This license is unrestrictive and allows using the package both in open and proprietary programs and even relicensing the code (if the new license doesn't conflict with the terms of the 3-clause BSD license).
 email: 
@@ -274,9 +274,19 @@ files:
 - ext/alglib/variancetests.h
 - ext/alglib/wsr.cpp
 - ext/alglib/wsr.h
+- ext/ap.i
+- ext/correlation.i
 - ext/extconf.rb
+- ext/logit.i
 - lib/alglib.rb
+- lib/alglib/correlation.rb
+- lib/alglib/linearregression.rb
+- lib/alglib/logit.rb
 - test/test_alglib.rb
+- test/test_correlation.rb
+- test/test_correlationtest.rb
+- test/test_linreg.rb
+- test/test_logit.rb
 has_rdoc: true
 homepage: http://rubyforge.org/projects/ruby-statsample/
 licenses: []
@@ -308,4 +318,8 @@ signing_key:
 specification_version: 3
 summary: ALGLIB - is a multilingual collection of algorithms designed to solve problems in the field of numeric analysis and data processing
 test_files: 
+- test/test_correlation.rb
+- test/test_correlationtest.rb
 - test/test_alglib.rb
+- test/test_logit.rb
+- test/test_linreg.rb