statsample 0.3.3 → 0.3.4

data/test/test_combination.rb

@@ -0,0 +1,42 @@
+$:.unshift(File.dirname(__FILE__)+'/../lib/')
+require 'statsample'
+require 'test/unit'
+
+class StatsampleCombinationTestCase < Test::Unit::TestCase
+	def initialize(*args)
+        super
+	end
+    def test_basic
+        k=3
+        n=5
+        expected=[[0,1,2],[0,1,3],[0,1,4],[0,2,3],[0,2,4],[0,3,4],[1,2,3],[1,2,4],[1,3,4],[2,3,4]]
+        comb=Statsample::Combination.new(k,n)
+        a=[]
+        comb.each{|y|
+            a.push(y)
+        }
+        assert_equal(expected,a)
+    end
+    def test_gsl_versus_ruby
+        if HAS_GSL
+        k=3
+        n=10
+        gsl=Statsample::Combination.new(k,n,false)
+        gsl_array=[]
+        gsl.each{|y|
+            gsl_array.push(y)
+        }
+        rb=Statsample::Combination.new(k,n,true)
+        rb_array=[]
+        rb.each{|y|
+            rb_array.push(y)
+        }
+        assert(gsl.d.is_a? Statsample::Combination::CombinationGsl)
+        assert(rb.d.is_a?  Statsample::Combination::CombinationRuby)
+        
+        assert_equal(rb_array,gsl_array)
+    else
+            puts "Not CombinationRuby vs CombinationGSL (no gsl)"
+        end
+    end
+end

data/test/test_csv.rb

@@ -9,7 +9,7 @@ class StatsampleCSVTestCase < Test::Unit::TestCase
 		super
 	end
     def test_read
-                assert_equal(6,@ds.cases)
+        assert_equal(6,@ds.cases)
         assert_equal(%w{id name age city a1},@ds.fields)
         id=[1,2,3,4,5,6].to_vector(:scale)
         name=["Alex","Claude","Peter","Franz","George","Fernand"].to_vector(:nominal)

data/test/test_dataset.rb

@@ -148,7 +148,12 @@ class StatsampleDatasetTestCase < Test::Unit::TestCase
     def test_case_as
         assert_equal({'id'=>1,'name'=>'Alex','city'=>'New York','age'=>20,'a1'=>'a,b'},@ds.case_as_hash(0))
         assert_equal([5,'George',5,'Tome','a,b,c'],@ds.case_as_array(4))
+        # Native methods
+        assert_equal({'id'=>1,'name'=>'Alex','city'=>'New York','age'=>20,'a1'=>'a,b'},@ds._case_as_hash(0))
+        assert_equal([5,'George',5,'Tome','a,b,c'],@ds._case_as_array(4))
 
+        
+        
     end
     def test_delete_vector
         @ds.delete_vector('name')

data/test/test_statistics.rb

@@ -65,8 +65,23 @@ class StatsampleStatisicsTestCase < Test::Unit::TestCase
     end
     def test_prop_pearson
 	if HAS_GSL    
-        assert_in_delta(0.42,Statsample::Bivariate.prop_pearson(Statsample::Bivariate.t_r(0.084,94),94),0.01)        
-        assert_in_delta(0.65,Statsample::Bivariate.prop_pearson(Statsample::Bivariate.t_r(0.046,95),95),0.01)
+        assert_in_delta(0.42, Statsample::Bivariate.prop_pearson(Statsample::Bivariate.t_r(0.084,94), 94),0.01)
+        assert_in_delta(0.65, Statsample::Bivariate.prop_pearson(Statsample::Bivariate.t_r(0.046,95), 95),0.01)
+        r=0.9
+        n=100
+        t=Statsample::Bivariate.t_r(r,n)
+        assert(Statsample::Bivariate.prop_pearson(t,n,:both)<0.05)
+        assert(Statsample::Bivariate.prop_pearson(t,n,:right)<0.05)
+        assert(Statsample::Bivariate.prop_pearson(t,n,:left)>0.05)
+
+        r=-0.9
+        n=100
+        t=Statsample::Bivariate.t_r(r,n)
+        assert(Statsample::Bivariate.prop_pearson(t,n,:both)<0.05)
+        assert(Statsample::Bivariate.prop_pearson(t,n,:right)>0.05)
+        assert(Statsample::Bivariate.prop_pearson(t,n,:left)<0.05)
+
+        
 	else
 		puts "Bivariate.prop_pearson not tested (no ruby-gsl)"
 	end
@@ -79,10 +94,12 @@ class StatsampleStatisicsTestCase < Test::Unit::TestCase
 			
 		end
 	end
+
 	def test_spearman
 		v1=[86,97,99,100,101,103,106,110,112,113].to_vector(:scale)
 		v2=[0,20,28,27,50,29,7,17,6,12].to_vector(:scale)
         assert_in_delta(-0.175758,Statsample::Bivariate.spearman(v1,v2),0.0001)
+        
 	end
 	def test_point_biserial
 		c=[1,3,5,6,7,100,200,300,400,300].to_vector(:scale)

data/test/test_svg_graph.rb

@@ -1,13 +1,16 @@
 $:.unshift(File.dirname(__FILE__)+'/../lib/')
 require 'statsample'
+require 'tmpdir'
 require 'tempfile'
+require 'fileutils'
 require 'test/unit'
 begin
 	require 'statsample/graph/svggraph'
 class StatsampleSvgGraphTestCase < Test::Unit::TestCase
 
 	def initialize(*args)
-		@image_path=File.dirname(__FILE__)+"/images"
+        @image_path=Dir::tmpdir+"/images"
+        FileUtils.mkdir(@image_path) if !File.exists? @image_path 
 		super
 	end
     def test_histogram
@@ -50,7 +53,7 @@ class StatsampleSvgGraphTestCase < Test::Unit::TestCase
 		file=@image_path+"/svggraph_histogram.svg"		
 		hist=vector.svggraph_histogram(5)
 		File.open(file,"wb") {|fp|
-		fp.write(hist.burn)
+            fp.write(hist.burn)
 		}
 		assert(File.exists?(file))
 		else

data/test/test_vector.rb

@@ -87,6 +87,7 @@ class StatsampleVectorTestCase < Test::Unit::TestCase
 	def test_nominal
 		assert_equal(@c[1],5)
 		assert_equal({ 1=>1,2=>1,3=>1,4=>1,5=>5,6=>2,7=>1,8=>1, 9=>1,10=>1},@c.frequencies)
+        assert_equal({ 1=>1,2=>1,3=>1,4=>1,5=>5,6=>2,7=>1,8=>1, 9=>1,10=>1},@c._frequencies)
 		assert_equal({ 1 => 1.quo(15) ,2=>1.quo(15), 3=>1.quo(15),4=>1.quo(15),5=>5.quo(15),6=>2.quo(15),7=>1.quo(15), 8=>1.quo(15), 9=>1.quo(15),10=>1.quo(15)}, @c.proportions)
         assert_equal(@c.proportion, 1.quo(15))
         assert_equal(@c.proportion(2), 1.quo(15))
@@ -266,5 +267,9 @@ class StatsampleVectorTestCase < Test::Unit::TestCase
         assert_equal(v1.labels,v3.labels)
         assert_not_same(v1.labels,v3.labels)
     end
-     
+    def test_paired_ties
+        a=[0,0,0,1,1,2,3,3,4,4,4].to_vector(:ordinal)
+        expected=[2,2,2,4.5,4.5,6,7.5,7.5,10,10,10].to_vector(:ordinal)
+        assert_equal(expected,a.ranked)
+    end
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification 
 name: statsample
 version: !ruby/object:Gem::Version 
-  version: 0.3.3
+  version: 0.3.4
 platform: ruby
 authors: 
 - Claudio Bustos
@@ -9,7 +9,7 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2009-08-11 00:00:00 -04:00
+date: 2009-08-21 00:00:00 -04:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -82,6 +82,7 @@ files:
 - lib/statsample/bivariate.rb
 - lib/statsample/chidistribution.rb
 - lib/statsample/codification.rb
+- lib/statsample/combination.rb
 - lib/statsample/converters.rb
 - lib/statsample/crosstab.rb
 - lib/statsample/dataset.rb
@@ -95,6 +96,7 @@ files:
 - lib/statsample/htmlreport.rb
 - lib/statsample/multiset.rb
 - lib/statsample/regression.rb
+- lib/statsample/regression/logit.rb
 - lib/statsample/regression/multiple.rb
 - lib/statsample/regression/multiple/alglibengine.rb
 - lib/statsample/regression/multiple/gslengine.rb
@@ -111,6 +113,7 @@ files:
 - test/_test_chart.rb
 - test/test_anova.rb
 - test/test_codification.rb
+- test/test_combination.rb
 - test/test_crosstab.rb
 - test/test_csv.csv
 - test/test_csv.rb
@@ -162,6 +165,7 @@ test_files:
 - test/test_crosstab.rb
 - test/test_svg_graph.rb
 - test/test_csv.rb
+- test/test_combination.rb
 - test/test_resample.rb
 - test/test_stratified.rb
 - test/test_vector.rb