statsample 0.3.1 → 0.3.2

data/test/test_anova.rb

@@ -10,10 +10,10 @@ class StatsampleAnovaTestCase < Test::Unit::TestCase
 		super
 	end
     def test_basic
-	assert_in_delta(72.933, @anova.sst,0.001)
+        assert_in_delta(72.933, @anova.sst,0.001)
         assert_in_delta(14.8,@anova.sswg,0.001)
         assert_in_delta(58.133,@anova.ssbg,0.001)
-	assert_in_delta(@anova.sst,@anova.sswg+@anova.ssbg,0.00001)
+        assert_in_delta(@anova.sst,@anova.sswg+@anova.ssbg,0.00001)
         assert_equal(14,@anova.df_total)
         assert_equal(12,@anova.df_wg)
         assert_equal(2,@anova.df_bg)

data/test/test_regression.rb

@@ -27,6 +27,33 @@ class StatsampleRegressionTestCase < Test::Unit::TestCase
         assert_in_delta(3.035,lr.f,0.001)
         
     end
+
+
+    def test_multiple_regression_gsl
+	    if HAS_GSL
+        @a=[1,3,2,4,3,5,4,6,5,7].to_vector(:scale)
+        @b=[3,3,4,4,5,5,6,6,4,4].to_vector(:scale)
+        @c=[11,22,30,40,50,65,78,79,99,100].to_vector(:scale)
+        @y=[3,4,5,6,7,8,9,10,20,30].to_vector(:scale)
+        ds={'a'=>@a,'b'=>@b,'c'=>@c,'y'=>@y}.to_dataset
+        lr=Statsample::Regression::Multiple::GslEngine.new(ds,'y')
+        model_test(lr)
+        predicted=[1.7857, 6.0989, 3.2433, 7.2908, 4.9667, 10.3428, 8.8158, 10.4717, 23.6639, 25.3198]
+        c_predicted=lr.predicted
+        predicted.each_index{|i|
+            assert_in_delta(predicted[i],c_predicted[i],0.001)
+        }
+        residuals=[1.2142, -2.0989, 1.7566, -1.29085, 2.033, -2.3428, 0.18414, -0.47177, -3.66395, 4.6801]
+        c_residuals=lr.residuals
+        residuals.each_index{|i|
+            assert_in_delta(residuals[i],c_residuals[i],0.001)
+        }
+	else
+		puts "Regression::Multiple::AlglibEngine not tested (no Alglib)"
+	end
+    end
+
+
     def test_multiple_regression_alglib
 	    if HAS_ALGIB
         @a=[1,3,2,4,3,5,4,6,5,7].to_vector(:scale)
@@ -59,7 +86,7 @@ class StatsampleRegressionTestCase < Test::Unit::TestCase
         s_coeffs={'a'=>0.151,'b'=>-0.547,'c'=>0.997}
         cs_coeefs=lr.standarized_coeffs
         s_coeffs.each_key{|k|
-            assert_in_delta(s_coeffs[k],cs_coeefs[k],0.001)
+            assert_in_delta(s_coeffs[k], cs_coeefs[k],0.001)
         }
         assert_in_delta(639.6,lr.sst,0.001)
         assert_in_delta(583.76,lr.ssr,0.001)
@@ -86,7 +113,7 @@ class StatsampleRegressionTestCase < Test::Unit::TestCase
         assert_in_delta(4.559, lr.constant_se,0.001)
         assert_in_delta(2.419, lr.constant_t,0.001)
     end
-    def test_regression_pairwise
+    def test_regression_rubyengine
         @a=[nil,1,3,2,4,3,5,4,6,5,7].to_vector(:scale)
         @b=[nil,3,3,4,4,5,5,6,6,4,4].to_vector(:scale)
         @c=[nil,11,22,30,40,50,65,78,79,99,100].to_vector(:scale)

data/test/test_statistics.rb

@@ -127,6 +127,15 @@ class StatsampleStatisicsTestCase < Test::Unit::TestCase
 		
 		end
     end
+    def test_ml
+        if(true)
+		real=[1,1,1,1].to_vector(:scale)
+		
+		pred=[0.0001,0.0001,0.0001,0.0001].to_vector(:scale)
+        # puts  Statsample::Bivariate.maximum_likehood_dichotomic(pred,real)
+
+        end
+    end
     def test_simple_linear_regression
 		a=[1,2,3,4,5,6].to_vector(:scale)
 		b=[6,2,4,10,12,8].to_vector(:scale)

data/test/test_vector.rb

@@ -11,6 +11,10 @@ class StatsampleVectorTestCase < Test::Unit::TestCase
     def test_enumerable
         val=@c.collect {|v| v}
         assert_equal(val,[5,5,5,5,5,6,6,7,8,9,10,1,2,3,4,nil,-99,-99])
+    end
+    def test_product
+        a=[1,2,3,4,5].to_vector(:scale)
+        assert_equal(120,a.product)
     end
 	def test_missing_values
 		@c.missing_values=[10]

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification 
 name: statsample
 version: !ruby/object:Gem::Version 
-  version: 0.3.1
+  version: 0.3.2
 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-04 00:00:00 -04:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -51,13 +51,14 @@ extensions: []
 
 extra_rdoc_files: 
 - History.txt
+- LICENSE.txt
 - Manifest.txt
 - README.txt
 files: 
 - History.txt
+- LICENSE.txt
 - Manifest.txt
 - README.txt
-- Rakefile
 - bin/statsample
 - demo/benchmark.rb
 - demo/chi-square.rb
@@ -95,6 +96,7 @@ files:
 - lib/statsample/regression.rb
 - lib/statsample/regression/multiple.rb
 - lib/statsample/regression/multiple/alglibengine.rb
+- lib/statsample/regression/multiple/gslengine.rb
 - lib/statsample/regression/multiple/rubyengine.rb
 - lib/statsample/regression/simple.rb
 - lib/statsample/reliability.rb
@@ -102,6 +104,7 @@ files:
 - lib/statsample/srs.rb
 - lib/statsample/test.rb
 - lib/statsample/vector.rb
+- setup.rb
 - test/_test_chart.rb
 - test/test_anova.rb
 - test/test_codification.rb

data/Rakefile

@@ -1,22 +0,0 @@
-#!/usr/bin/ruby
-# -*- ruby -*-
-
-require 'rubygems'
-require 'hoe'
-require './lib/statsample'
-
-if File.exists? './local_rakefile.rb'
-	require './local_rakefile'
-end
-	
-Hoe.spec('statsample') do |p|
-	p.version=Statsample::VERSION
-	p.rubyforge_name = "ruby-statsample"
-	p.developer('Claudio Bustos', 'clbustos@gmail.com')
-	p.extra_deps << ["spreadsheet","=0.6.4"] << "svg-graph"
-	p.clean_globs << "test/images/*" 
-	#  p.rdoc_pattern = /^(lib|bin|ext\/distributions)|txt$/
-end
-
-
-# vim: syntax=Ruby