statsample 0.7.0 → 0.8.0

data/test/test_permutation.rb

@@ -1,37 +1,36 @@
-$:.unshift(File.dirname(__FILE__)+'/../lib/')
-require 'statsample'
-require 'test/unit'
+require(File.dirname(__FILE__)+'/test_helpers.rb')
 
-class StatsamplePermutationTestCase < Test::Unit::TestCase
-    def initialize(*args)
-          super
-    end
-    def test_number_of_permutations
-      per1=Statsample::Permutation.new(4)
-      assert_equal(24,per1.permutation_number)
-      per2=Statsample::Permutation.new([1,1,1,0,0,0])
-      assert_equal(20,per2.permutation_number)
-    end
-    def test_permutation_with_number
-      per1=Statsample::Permutation.new(2)
-      exp1=[[0,1],[1,0]]
-      assert_equal(exp1,per1.permutations)
-      per2=Statsample::Permutation.new(3)
-      exp2=[[0,1,2],[0,2,1],[1,0,2],[1,2,0],[2,0,1],[2,1,0]]
-      assert_equal(exp2,per2.permutations)
-      
-    end
-    def test_permutation_with_array_simple
-      per1=Statsample::Permutation.new(%w{a b})
-      exp1=[['a','b'],['b','a']]
-      assert_equal(exp1,per1.permutations)
-      per2=Statsample::Permutation.new(%w{a b c})
-      exp2=[%w{a b c},%w{a c b},%w{b a c} ,%w{b c a},%w{c a b},%w{c b a}]
-      assert_equal(exp2,per2.permutations)
-    end
-    def test_permutation_with_array_repeated
-      per1=Statsample::Permutation.new([0,0,1,1])
-      exp1=[[0,0,1,1],[0,1,0,1],[0,1,1,0],[1,0,0,1],[1,0,1,0],[1,1,0,0]]
-      assert_equal(exp1,per1.permutations)
-    end
-end
+
+class StatsamplePermutationTestCase < MiniTest::Unit::TestCase
+  def initialize(*args)
+    super
+  end
+  def test_number_of_permutations
+    per1=Statsample::Permutation.new(4)
+    assert_equal(24,per1.permutation_number)
+    per2=Statsample::Permutation.new([1,1,1,0,0,0])
+    assert_equal(20,per2.permutation_number)
+  end
+  def test_permutation_with_number
+    per1=Statsample::Permutation.new(2)
+    exp1=[[0,1],[1,0]]
+    assert_equal(exp1,per1.permutations)
+    per2=Statsample::Permutation.new(3)
+    exp2=[[0,1,2],[0,2,1],[1,0,2],[1,2,0],[2,0,1],[2,1,0]]
+    assert_equal(exp2,per2.permutations)
+
+  end
+  def test_permutation_with_array_simple
+    per1=Statsample::Permutation.new(%w{a b})
+    exp1=[['a','b'],['b','a']]
+    assert_equal(exp1,per1.permutations)
+    per2=Statsample::Permutation.new(%w{a b c})
+    exp2=[%w{a b c},%w{a c b},%w{b a c} ,%w{b c a},%w{c a b},%w{c b a}]
+    assert_equal(exp2,per2.permutations)
+  end
+  def test_permutation_with_array_repeated
+    per1=Statsample::Permutation.new([0,0,1,1])
+    exp1=[[0,0,1,1],[0,1,0,1],[0,1,1,0],[1,0,0,1],[1,0,1,0],[1,1,0,0]]
+    assert_equal(exp1,per1.permutations)
+  end
+end

data/test/test_promise_after.rb

@@ -0,0 +1,39 @@
+require(File.dirname(__FILE__)+'/test_helpers.rb')
+
+
+class StatsamplePromiseAfterTestCase < MiniTest::Unit::TestCase
+  class ExpensiveClass
+    extend Statsample::PromiseAfter
+    attr_reader :a, :dirty
+    def initialize
+      @a=nil
+      @b=nil
+      @dirty=false
+    end
+    def compute
+      @a="After"
+      @b="After"
+      @dirty=true
+    end
+    def a
+      @a.nil? ? nil : "@a=#{@a}"
+    end
+    def b
+      "@b=#{@b}"
+    end
+    promise_after :compute, :a
+  end
+  def setup
+    @ec=ExpensiveClass.new
+  end
+  def test_promise_after_before
+    assert_equal(nil, @ec.instance_variable_get("@a")) 
+    assert_equal(nil, @ec.instance_variable_get("@b")) 
+    assert_equal("@b=",@ec.b)
+    refute(@ec.dirty) 
+    # Calling method a active compute
+    assert_equal("@a=After", @ec.a) 
+    assert(@ec.dirty)
+    assert_equal("@b=After", @ec.b) 
+  end
+end

data/test/test_regression.rb

@@ -1,47 +1,45 @@
-$:.unshift(File.dirname(__FILE__)+'/../lib/')
-require 'statsample'
-require 'minitest/unit'
-MiniTest::Unit.autorun
+require(File.dirname(__FILE__)+'/test_helpers.rb')
+
 class StatsampleRegressionTestCase < MiniTest::Unit::TestCase
-	def test_parameters
+  def test_parameters
     @x=[13,20,10,33,15].to_vector(:scale)
-		@y=[23,18,35,10,27	].to_vector(:scale)
-		reg=Statsample::Regression::Simple.new_from_vectors(@x,@y)
+    @y=[23,18,35,10,27	].to_vector(:scale)
+    reg=Statsample::Regression::Simple.new_from_vectors(@x,@y)
     _test_simple_regression(reg)
     ds={'x'=>@x,'y'=>@y}.to_dataset
-		reg=Statsample::Regression::Simple.new_from_dataset(ds,'x','y')
+    reg=Statsample::Regression::Simple.new_from_dataset(ds,'x','y')
     _test_simple_regression(reg)
     reg=Statsample::Regression.simple(@x,@y)
     _test_simple_regression(reg)
-    
-	end
+
+  end
   def _test_simple_regression(reg)
     assert_in_delta(40.009, reg.a,0.001)
-		assert_in_delta(-0.957, reg.b,0.001)
-		assert_in_delta(4.248,reg.standard_error,0.002)
+    assert_in_delta(-0.957, reg.b,0.001)
+    assert_in_delta(4.248,reg.standard_error,0.002)
   end
-   
+
   def test_multiple_dependent
     complete=Matrix[
-    [1,0.53,0.62,0.19,-0.09,0.08,0.02,-0.12,0.08],
-    [0.53,1,0.61,0.23,0.1,0.18,0.02,-0.1,0.15],
-    [0.62,0.61,1,0.03,0.1,0.12,0.03,-0.06,0.12],
-    [0.19,0.23,0.03,1,-0.02,0.02,0,-0.02,-0.02],
-    [-0.09,0.1,0.1,-0.02,1,0.05,0.06,0.18,0.02],
-    [0.08,0.18,0.12,0.02,0.05,1,0.22,-0.07,0.36],
-    [0.02,0.02,0.03,0,0.06,0.22,1,-0.01,-0.05],
-    [-0.12,-0.1,-0.06,-0.02,0.18,-0.07,-0.01,1,-0.03],
+      [1,0.53,0.62,0.19,-0.09,0.08,0.02,-0.12,0.08],
+      [0.53,1,0.61,0.23,0.1,0.18,0.02,-0.1,0.15],
+      [0.62,0.61,1,0.03,0.1,0.12,0.03,-0.06,0.12],
+      [0.19,0.23,0.03,1,-0.02,0.02,0,-0.02,-0.02],
+      [-0.09,0.1,0.1,-0.02,1,0.05,0.06,0.18,0.02],
+      [0.08,0.18,0.12,0.02,0.05,1,0.22,-0.07,0.36],
+      [0.02,0.02,0.03,0,0.06,0.22,1,-0.01,-0.05],
+      [-0.12,-0.1,-0.06,-0.02,0.18,-0.07,-0.01,1,-0.03],
     [0.08,0.15,0.12,-0.02,0.02,0.36,-0.05,-0.03,1]]
     complete.extend Statsample::CovariateMatrix
     complete.fields=%w{adhd cd odd sex age monly mwork mage poverty}
-    
+
     lr=Statsample::Regression::Multiple::MultipleDependent.new(complete, %w{adhd cd odd})
-    
-    
+
+
     assert_in_delta(0.197, lr.r2yx,0.001)
     assert_in_delta(0.197, lr.r2yx_covariance,0.001)
     assert_in_delta(0.07, lr.p2yx,0.001)
-    
+
 
   end
   def test_multiple_regression_pairwise_2
@@ -57,7 +55,7 @@ class StatsampleRegressionTestCase < MiniTest::Unit::TestCase
     assert_in_delta(1361.130,lr.ssr,0.001)
     assert_in_delta(1046.306,lr.sse,0.001)
     assert_in_delta(3.035,lr.f,0.001)
-      
+
   end
 
 
@@ -73,12 +71,12 @@ class StatsampleRegressionTestCase < MiniTest::Unit::TestCase
       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)
+        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)
+        assert_in_delta(residuals[i],c_residuals[i],0.001)
       }
     else
       puts "Regression::Multiple::GslEngine not tested (no Gsl)"
@@ -86,54 +84,54 @@ class StatsampleRegressionTestCase < MiniTest::Unit::TestCase
   end
 
 
-  
+
   def model_test_matrix(lr,name='undefined')
-    
+
     stan_coeffs={'a'=>0.151,'b'=>-0.547,'c'=>0.997}
     unstan_coeffs={'a'=>0.695, 'b'=>-4.286, 'c'=>0.266}
-     
+
     unstan_coeffs.each_key{|k|
       assert_in_delta(unstan_coeffs[k], lr.coeffs[k],0.001,"b coeffs - #{name}")
     }
 
-   stan_coeffs.each_key{|k|
-     assert_in_delta(stan_coeffs[k], lr.standarized_coeffs[k],0.001, "beta coeffs - #{name}")
+    stan_coeffs.each_key{|k|
+      assert_in_delta(stan_coeffs[k], lr.standarized_coeffs[k],0.001, "beta coeffs - #{name}")
     }
-    
+
     assert_in_delta(11.027,lr.constant,0.001)
-    
+
     assert_in_delta(0.955,lr.r,0.001)
     assert_in_delta(0.913,lr.r2,0.001)
 
     assert_in_delta(20.908, lr.f,0.001)
     assert_in_delta(0.001, lr.significance, 0.001)
     assert_in_delta(0.226,lr.tolerance("a"),0.001)
-    
+
     coeffs_se={"a"=>1.171,"b"=>1.129,"c"=>0.072}
 
-    
-    
+
+
     ccoeffs_se=lr.coeffs_se
     coeffs_se.each_key{|k|
-        assert_in_delta(coeffs_se[k],ccoeffs_se[k],0.001)
+      assert_in_delta(coeffs_se[k],ccoeffs_se[k],0.001)
     }
     coeffs_t={"a"=>0.594,"b"=>-3.796,"c"=>3.703}
     ccoeffs_t=lr.coeffs_t
     coeffs_t.each_key{|k|
-        assert_in_delta(coeffs_t[k], ccoeffs_t[k],0.001)
+      assert_in_delta(coeffs_t[k], ccoeffs_t[k],0.001)
     }
-   
+
     assert_in_delta(639.6,lr.sst,0.001)
     assert_in_delta(583.76,lr.ssr,0.001)
     assert_in_delta(55.840,lr.sse,0.001)
-    
+
   end
   def model_test(lr,name='undefined')
     model_test_matrix(lr,name)
     assert_in_delta(4.559, lr.constant_se,0.001)
     assert_in_delta(2.419, lr.constant_t,0.001)
-    
-    assert_in_delta(1.785,lr.process([1,3,11]),0.001) 
+
+    assert_in_delta(1.785,lr.process([1,3,11]),0.001)
   end
   def test_regression_matrix
     @a=[1,3,2,4,3,5,4,6,5,7].to_vector(:scale)
@@ -141,16 +139,16 @@ class StatsampleRegressionTestCase < MiniTest::Unit::TestCase
     @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
-    
+
     cor=Statsample::Bivariate.correlation_matrix(ds)
     lr=Statsample::Regression::Multiple::MatrixEngine.new(cor,'y', :y_mean=>@y.mean, :x_mean=>{'a'=>ds['a'].mean, 'b'=>ds['b'].mean, 'c'=>ds['c'].mean}, :cases=>@a.size, :y_sd=>@y.sd , :x_sd=>{'a' => @a.sd, 'b' => @b.sd, 'c' => @c.sd})
     assert_nil(lr.constant_se)
     assert_nil(lr.constant_t)
     model_test_matrix(lr, "correlation matrix")
-    
+
     covariance=Statsample::Bivariate.covariance_matrix(ds)
     lr=Statsample::Regression::Multiple::MatrixEngine.new(covariance,'y', :y_mean=>@y.mean, :x_mean=>{'a'=>ds['a'].mean, 'b'=>ds['b'].mean, 'c'=>ds['c'].mean}, :cases=>@a.size)
-    
+
     model_test(lr , "covariance matrix")
   end
   def test_regression_rubyengine
@@ -161,10 +159,10 @@ class StatsampleRegressionTestCase < MiniTest::Unit::TestCase
     ds={'a'=>@a,'b'=>@b,'c'=>@c,'y'=>@y}.to_dataset
     lr=Statsample::Regression::Multiple::RubyEngine.new(ds,'y')
     model_test(lr, 'rubyengine with missing data')
-    
+
     predicted=[nil,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 do |i|
       if c_predicted[i].nil?
         assert(predicted[i].nil?)
@@ -181,6 +179,6 @@ class StatsampleRegressionTestCase < MiniTest::Unit::TestCase
         assert_in_delta(residuals[i],c_residuals[i],0.001)
       end
     end
-    
+
   end
-end
+end

data/test/test_reliability.rb

@@ -1,33 +1,32 @@
-$:.unshift(File.dirname(__FILE__)+'/../lib/')
-require 'statsample'
-require 'test/unit'
+require(File.dirname(__FILE__)+'/test_helpers.rb')
 
-class StatsampleReliabilityTestCase < Test::Unit::TestCase
 
-	def initialize(*args)
-		super
-        @x1=[1,1,1,1,2,2,2,2,3,3,3,30].to_vector(:scale)
-        @x2=[1,1,1,2,2,3,3,3,3,4,4,50].to_vector(:scale)
-        @x3=[2,2,1,1,1,2,2,2,3,4,5,40].to_vector(:scale)
-        @x4=[1,2,3,4,4,4,4,3,4,4,5,30].to_vector(:scale)
-        @ds={'x1'=>@x1,'x2'=>@x2,'x3'=>@x3,'x4'=>@x4}.to_dataset
-	end
-   
-    def test_general
-        ia=Statsample::Reliability::ItemAnalysis.new(@ds)
-        assert_in_delta(0.980,ia.alpha,0.001)
-        assert_in_delta(0.999,ia.alpha_standarized,0.001)
-        assert_in_delta(0.999,ia.item_total_correlation()['x1'],0.001)
-        assert_in_delta(1050.455,ia.stats_if_deleted()['x1'][:variance_sample],0.001)
-    end	
-    def test_icc
-        #p @x1.factors
-        icc=Statsample::Reliability::ItemCharacteristicCurve.new(@ds)
-        # Need to create the test!!!!
-        #p icc.curve_field('x1',1).sort
-        #p icc.curve_field('x1',2).sort
-        #p icc.curve_field('x1',3).sort
-        #p icc.curve_field('x1',30).sort
-        
-    end
+class StatsampleReliabilityTestCase < MiniTest::Unit::TestCase
+
+  def initialize(*args)
+    super
+    @x1=[1,1,1,1,2,2,2,2,3,3,3,30].to_vector(:scale)
+    @x2=[1,1,1,2,2,3,3,3,3,4,4,50].to_vector(:scale)
+    @x3=[2,2,1,1,1,2,2,2,3,4,5,40].to_vector(:scale)
+    @x4=[1,2,3,4,4,4,4,3,4,4,5,30].to_vector(:scale)
+    @ds={'x1'=>@x1,'x2'=>@x2,'x3'=>@x3,'x4'=>@x4}.to_dataset
+  end
+
+  def test_general
+    ia=Statsample::Reliability::ItemAnalysis.new(@ds)
+    assert_in_delta(0.980,ia.alpha,0.001)
+    assert_in_delta(0.999,ia.alpha_standarized,0.001)
+    assert_in_delta(0.999,ia.item_total_correlation()['x1'],0.001)
+    assert_in_delta(1050.455,ia.stats_if_deleted()['x1'][:variance_sample],0.001)
+  end
+  def test_icc
+    #p @x1.factors
+    icc=Statsample::Reliability::ItemCharacteristicCurve.new(@ds)
+    # Need to create the test!!!!
+    #p icc.curve_field('x1',1).sort
+    #p icc.curve_field('x1',2).sort
+    #p icc.curve_field('x1',3).sort
+    #p icc.curve_field('x1',30).sort
+
+  end
 end

data/test/test_resample.rb

@@ -1,24 +1,23 @@
-$:.unshift(File.dirname(__FILE__)+'/../lib/')
-require 'statsample'
-require 'test/unit'
+require(File.dirname(__FILE__)+'/test_helpers.rb')
 
-class StatsampleResampleTestCase < Test::Unit::TestCase
-	def initialize(*args)
-		super
-	end
-    def test_basic
-        r=Statsample::Resample.generate(20,1,10)
-        assert_equal(20,r.size)
-        assert(r.min>=1)
-        assert(r.max<=10)
-    end
-    def test_repeat_and_save
-        r=Statsample::Resample.repeat_and_save(400) {
-            Statsample::Resample.generate(20,1,10).count(1)
-        }
-        assert_equal(400,r.size)
-        v=Statsample::Vector.new(r,:scale)
-        a=v.count {|x|  x > 3}
-        assert(a>=30 && a<=70)
-    end
-end
+
+class StatsampleResampleTestCase < MiniTest::Unit::TestCase
+  def initialize(*args)
+    super
+  end
+  def test_basic
+    r=Statsample::Resample.generate(20,1,10)
+    assert_equal(20,r.size)
+    assert(r.min>=1)
+    assert(r.max<=10)
+  end
+  def test_repeat_and_save
+    r=Statsample::Resample.repeat_and_save(400) {
+      Statsample::Resample.generate(20,1,10).count(1)
+    }
+    assert_equal(400,r.size)
+    v=Statsample::Vector.new(r,:scale)
+    a=v.count {|x|  x > 3}
+    assert(a>=30 && a<=70)
+  end
+end

data/test/test_srs.rb

@@ -1,11 +1,10 @@
-$:.unshift(File.dirname(__FILE__)+'/../lib/')
-require 'statsample'
-require 'test/unit'
+require(File.dirname(__FILE__)+'/test_helpers.rb')
 
-class StatsampleSrsTestCase < Test::Unit::TestCase
-    def test_std_error
-        assert_equal(384,Statsample::SRS.estimation_n0(0.05,0.5,0.95).to_i)
-        assert_equal(108,Statsample::SRS.estimation_n(0.05,0.5,150,0.95).to_i)
-        assert_in_delta(0.0289,Statsample::SRS.proportion_sd_kp_wor(0.5,100,150),0.001)
-    end
+
+class StatsampleSrsTestCase < MiniTest::Unit::TestCase
+  def test_std_error
+    assert_equal(384,Statsample::SRS.estimation_n0(0.05,0.5,0.95).to_i)
+    assert_equal(108,Statsample::SRS.estimation_n(0.05,0.5,150,0.95).to_i)
+    assert_in_delta(0.0289,Statsample::SRS.proportion_sd_kp_wor(0.5,100,150),0.001)
+  end
 end