statsample 1.4.1 → 1.4.2

data/test/test_anova_contrast.rb

@@ -1,36 +1,36 @@
-require(File.expand_path(File.dirname(__FILE__)+'/helpers_tests.rb'))
-class StatsampleAnovaContrastTestCase < MiniTest::Unit::TestCase
+require(File.expand_path(File.dirname(__FILE__) + '/helpers_tests.rb'))
+class StatsampleAnovaContrastTestCase < Minitest::Test
   context(Statsample::Anova::Contrast) do
     setup do
-      constant=[12,13,11,12,12].to_scale
-      frequent=[9,10,9,13,14].to_scale
-      infrequent=[15,16,17,16,16].to_scale
-      never=[17,18,12,18,20].to_scale
-      @vectors=[constant, frequent, infrequent, never]
-      @c=Statsample::Anova::Contrast.new(:vectors=>@vectors)
+      constant = [12, 13, 11, 12, 12].to_scale
+      frequent = [9, 10, 9, 13, 14].to_scale
+      infrequent = [15, 16, 17, 16, 16].to_scale
+      never = [17, 18, 12, 18, 20].to_scale
+      @vectors = [constant, frequent, infrequent, never]
+      @c = Statsample::Anova::Contrast.new(vectors: @vectors)
     end
-    should "return correct value using c" do
-      @c.c([1,-1.quo(3),-1.quo(3),-1.quo(3)])
-      #@c.c([1,-0.333,-0.333,-0.333])
+    should 'return correct value using c' do
+      @c.c([1, -1.quo(3), -1.quo(3), -1.quo(3)])
+      # @c.c([1,-0.333,-0.333,-0.333])
       assert_in_delta(-2.6667, @c.psi, 0.0001)
       assert_in_delta(1.0165, @c.se, 0.0001)
       assert_in_delta(-2.623, @c.t, 0.001)
-      assert_in_delta(-4.82, @c.confidence_interval[0],0.01)
-      assert_in_delta(-0.51, @c.confidence_interval[1],0.01)
-      assert(@c.summary.size>0)
+      assert_in_delta(-4.82, @c.confidence_interval[0], 0.01)
+      assert_in_delta(-0.51, @c.confidence_interval[1], 0.01)
+      assert(@c.summary.size > 0)
     end
-    should "return correct values using c_by_index" do
-      @c.c_by_index([0],[1,2,3])
+    should 'return correct values using c_by_index' do
+      @c.c_by_index([0], [1, 2, 3])
       assert_in_delta(-2.6667, @c.psi, 0.0001)
       assert_in_delta(1.0165, @c.se, 0.0001)
       assert_in_delta(-2.623, @c.t, 0.001)
     end
-    should "return correct values using incomplete c_by_index" do
-      c1=Statsample::Anova::Contrast.new(:vectors=>@vectors, :c=>[0.5,0.5,-1,0])
-      c2=Statsample::Anova::Contrast.new(:vectors=>@vectors, :c1=>[0,1],:c2=>[2])
-      assert_equal(c1.psi,c2.psi)
-      assert_equal(c1.se,c2.se)
-      assert_equal(c1.t,c2.t)
+    should 'return correct values using incomplete c_by_index' do
+      c1 = Statsample::Anova::Contrast.new(vectors: @vectors, c: [0.5, 0.5, -1, 0])
+      c2 = Statsample::Anova::Contrast.new(vectors: @vectors, c1: [0, 1], c2: [2])
+      assert_equal(c1.psi, c2.psi)
+      assert_equal(c1.se, c2.se)
+      assert_equal(c1.t, c2.t)
     end
   end
 end

data/test/test_anovaoneway.rb

@@ -1,26 +1,26 @@
-require(File.expand_path(File.dirname(__FILE__)+'/helpers_tests.rb'))
-class StatsampleAnovaOneWayTestCase < MiniTest::Unit::TestCase
+require(File.expand_path(File.dirname(__FILE__) + '/helpers_tests.rb'))
+class StatsampleAnovaOneWayTestCase < Minitest::Test
   context(Statsample::Anova::OneWay) do
     setup do
-      @ss_num=30.08
-      @ss_den=87.88
-      @df_num=2
-      @df_den=21
-      @anova=Statsample::Anova::OneWay.new(:ss_num=>@ss_num, :ss_den=>@ss_den, :df_num=>@df_num, :df_den=>@df_den)
+      @ss_num = 30.08
+      @ss_den = 87.88
+      @df_num = 2
+      @df_den = 21
+      @anova = Statsample::Anova::OneWay.new(ss_num: @ss_num, ss_den: @ss_den, df_num: @df_num, df_den: @df_den)
     end
-    should "Statsample::Anova.oneway respond to #oneway" do
+    should 'Statsample::Anova.oneway respond to #oneway' do
       assert(Statsample::Anova.respond_to? :oneway)
     end
-    should "return correct value for ms_num and ms_den" do
+    should 'return correct value for ms_num and ms_den' do
       assert_in_delta(15.04, @anova.ms_num, 0.01)
       assert_in_delta(4.18, @anova.ms_den, 0.01)
     end
-    should "return correct value for f" do
+    should 'return correct value for f' do
       assert_in_delta(3.59, @anova.f, 0.01)
     end
-    should "respond to summary" do
+    should 'respond to summary' do
       assert(@anova.respond_to? :summary)
-      assert(@anova.summary.size>0)
+      assert(@anova.summary.size > 0)
     end
   end
 end

data/test/test_anovatwoway.rb

@@ -1,38 +1,37 @@
-require(File.expand_path(File.dirname(__FILE__)+'/helpers_tests.rb'))
-class StatsampleAnovaTwoWayTestCase < MiniTest::Unit::TestCase
+require(File.expand_path(File.dirname(__FILE__) + '/helpers_tests.rb'))
+class StatsampleAnovaTwoWayTestCase < Minitest::Test
   context(Statsample::Anova::TwoWay) do
     setup do
-      @ss_a=192.2
-      @ss_b=57.8
-      @ss_axb=168.2
-      @ss_within=75.6
-      @df_a=@df_b=1
-      @df_within=16
-      @anova=Statsample::Anova::TwoWay.new(:ss_a=>@ss_a, :ss_b=>@ss_b, :ss_axb=>@ss_axb, :ss_within=>@ss_within , :df_a=>@df_a, :df_b=>@df_b, :df_within=>@df_within)
+      @ss_a = 192.2
+      @ss_b = 57.8
+      @ss_axb = 168.2
+      @ss_within = 75.6
+      @df_a = @df_b = 1
+      @df_within = 16
+      @anova = Statsample::Anova::TwoWay.new(ss_a: @ss_a, ss_b: @ss_b, ss_axb: @ss_axb, ss_within: @ss_within, df_a: @df_a, df_b: @df_b, df_within: @df_within)
     end
-    should "Statsample::Anova.twoway respond to #twoway" do
-    assert(Statsample::Anova.respond_to? :twoway)
+    should 'Statsample::Anova.twoway respond to #twoway' do
+      assert(Statsample::Anova.respond_to? :twoway)
     end
-    should "return correct value for ms_a, ms_b and ms_axb" do
+    should 'return correct value for ms_a, ms_b and ms_axb' do
       assert_in_delta(192.2, @anova.ms_a, 0.01)
       assert_in_delta(57.8, @anova.ms_b, 0.01)
       assert_in_delta(168.2, @anova.ms_axb, 0.01)
-      
     end
-    should "return correct value for f " do
+    should 'return correct value for f ' do
       assert_in_delta(40.68, @anova.f_a, 0.01)
       assert_in_delta(12.23, @anova.f_b, 0.01)
       assert_in_delta(35.60, @anova.f_axb, 0.01)
     end
-    should "return correct value for probability for f " do
+    should 'return correct value for probability for f ' do
       assert(@anova.f_a_probability < 0.05)
       assert(@anova.f_b_probability < 0.05)
       assert(@anova.f_axb_probability < 0.05)
     end
 
-    should "respond to summary" do
+    should 'respond to summary' do
       assert(@anova.respond_to? :summary)
-      assert(@anova.summary.size>0)
+      assert(@anova.summary.size > 0)
     end
   end
 end

data/test/test_anovatwowaywithdataset.rb

@@ -1,49 +1,47 @@
-require(File.expand_path(File.dirname(__FILE__)+'/helpers_tests.rb'))
+require(File.expand_path(File.dirname(__FILE__) + '/helpers_tests.rb'))
 # Reference:
 # * http://www.uwsp.edu/psych/Stat/13/anova-2w.htm#III
-class StatsampleAnovaTwoWayWithVectorsTestCase < MiniTest::Unit::TestCase
+class StatsampleAnovaTwoWayWithVectorsTestCase < Minitest::Test
   context(Statsample::Anova::TwoWayWithVectors) do
     setup do
-      @pa=[5,4,3,4,2,18,19,14,12,15,6,7,5,8,4,6,9,5,9,3].to_scale
-      @pa.name="Passive Avoidance"
-      @a=[0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1].to_vector
-      @a.labels={0=>'0%',1=>'35%'}
-      @a.name='Diet'
-      @b=[0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1].to_vector
-      @b.labels={0=>'Young',1=>'Older'}
-      @b.name="Age"
-      @anova=Statsample::Anova::TwoWayWithVectors.new(:a=>@a,:b=>@b, :dependent=>@pa)
+      @pa = [5, 4, 3, 4, 2, 18, 19, 14, 12, 15, 6, 7, 5, 8, 4, 6, 9, 5, 9, 3].to_scale
+      @pa.name = 'Passive Avoidance'
+      @a = [0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1].to_vector
+      @a.labels = { 0 => '0%', 1 => '35%' }
+      @a.name = 'Diet'
+      @b = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1].to_vector
+      @b.labels = { 0 => 'Young', 1 => 'Older' }
+      @b.name = 'Age'
+      @anova = Statsample::Anova::TwoWayWithVectors.new(a: @a, b: @b, dependent: @pa)
     end
-    should "Statsample::Anova respond to #twoway_with_vectors" do
-    assert(Statsample::Anova.respond_to? :twoway_with_vectors)
+    should 'Statsample::Anova respond to #twoway_with_vectors' do
+      assert(Statsample::Anova.respond_to? :twoway_with_vectors)
     end
-    should "#new returns the same as Statsample::Anova.twoway_with_vectors" do
-      @anova2=Statsample::Anova.twoway_with_vectors(:a=>@a,:b=>@b, :dependent=>@pa)
+    should '#new returns the same as Statsample::Anova.twoway_with_vectors' do
+      @anova2 = Statsample::Anova.twoway_with_vectors(a: @a, b: @b, dependent: @pa)
       assert_equal(@anova.summary, @anova2.summary)
     end
-    should "return correct value for ms_a, ms_b and ms_axb" do
+    should 'return correct value for ms_a, ms_b and ms_axb' do
       assert_in_delta(192.2, @anova.ms_a, 0.01)
       assert_in_delta(57.8, @anova.ms_b, 0.01)
       assert_in_delta(168.2, @anova.ms_axb, 0.01)
-      
     end
-    should "return correct value for f " do
+    should 'return correct value for f ' do
       assert_in_delta(40.68, @anova.f_a, 0.01)
       assert_in_delta(12.23, @anova.f_b, 0.01)
       assert_in_delta(35.60, @anova.f_axb, 0.01)
     end
-    should "return correct value for probability for f " do
+    should 'return correct value for probability for f ' do
       assert(@anova.f_a_probability < 0.05)
       assert(@anova.f_b_probability < 0.05)
       assert(@anova.f_axb_probability < 0.05)
     end
 
-    should "respond to summary" do
-      
-      @anova.summary_descriptives=true
-      @anova.summary_levene=true
+    should 'respond to summary' do
+      @anova.summary_descriptives = true
+      @anova.summary_levene = true
       assert(@anova.respond_to? :summary)
-      assert(@anova.summary.size>0)
+      assert(@anova.summary.size > 0)
     end
   end
 end

data/test/test_anovawithvectors.rb

@@ -1,102 +1,100 @@
-require(File.expand_path(File.dirname(__FILE__)+'/helpers_tests.rb'))
-class StatsampleAnovaOneWayWithVectorsTestCase < MiniTest::Unit::TestCase
+require(File.expand_path(File.dirname(__FILE__) + '/helpers_tests.rb'))
+class StatsampleAnovaOneWayWithVectorsTestCase < Minitest::Test
   context(Statsample::Anova::OneWayWithVectors) do
-    
-    context("when initializing") do
+    context('when initializing') do
       setup do
-        @v1=10.times.map {rand(100)}.to_scale
-        @v2=10.times.map {rand(100)}.to_scale
-        @v3=10.times.map {rand(100)}.to_scale
+        @v1 = 10.times.map { rand(100) }.to_scale
+        @v2 = 10.times.map { rand(100) }.to_scale
+        @v3 = 10.times.map { rand(100) }.to_scale
       end
-      should "be the same using [] or args*" do
-        a1=Statsample::Anova::OneWayWithVectors.new(@v1,@v2,@v3)
-        a2=Statsample::Anova::OneWayWithVectors.new([@v1,@v2,@v3])
-        assert_equal(a1.f,a2.f)
+      should 'be the same using [] or args*' do
+        a1 = Statsample::Anova::OneWayWithVectors.new(@v1, @v2, @v3)
+        a2 = Statsample::Anova::OneWayWithVectors.new([@v1, @v2, @v3])
+        assert_equal(a1.f, a2.f)
       end
-      should "be the same using module method or object instantiation" do
-        a1=Statsample::Anova::OneWayWithVectors.new(@v1,@v2,@v3)
-        a2=Statsample::Anova.oneway_with_vectors(@v1,@v2,@v3)
-        assert_equal(a1.f,a2.f)
+      should 'be the same using module method or object instantiation' do
+        a1 = Statsample::Anova::OneWayWithVectors.new(@v1, @v2, @v3)
+        a2 = Statsample::Anova.oneway_with_vectors(@v1, @v2, @v3)
+        assert_equal(a1.f, a2.f)
       end
-      should "detect optional hash" do
-        a1=Statsample::Anova::OneWayWithVectors.new(@v1,@v2,@v3, {:name=>'aaa'})
+      should 'detect optional hash' do
+        a1 = Statsample::Anova::OneWayWithVectors.new(@v1, @v2, @v3, name: 'aaa')
         assert_equal('aaa', a1.name)
       end
-      should "omit incorrect arguments" do
-        a1=Statsample::Anova::OneWayWithVectors.new(@v1,@v2,@v3, {:name=>'aaa'})
-        a2=Statsample::Anova::OneWayWithVectors.new(@v1,nil,nil,@v2,@v3, {:name=>'aaa'})
-        assert_equal(a1.f,a2.f)
+      should 'omit incorrect arguments' do
+        a1 = Statsample::Anova::OneWayWithVectors.new(@v1, @v2, @v3, name: 'aaa')
+        a2 = Statsample::Anova::OneWayWithVectors.new(@v1, nil, nil, @v2, @v3, name: 'aaa')
+        assert_equal(a1.f, a2.f)
       end
     end
     setup do
-      @v1=[3,3,2,3,6].to_vector(:scale)
-      @v2=[7,6,5,6,7].to_vector(:scale)
-      @v3=[9,8,9,7,8].to_vector(:scale)
-      @name="Anova testing"
-      @anova=Statsample::Anova::OneWayWithVectors.new(@v1,@v2,@v3, :name=>@name)
-    end
-    should "store correctly contrasts" do
-      c1=Statsample::Anova::Contrast.new(:vectors=>[@v1,@v2,@v3], :c=>[1,-0.5, -0.5])
-      
-      c2=@anova.contrast(:c=>[1,-0.5,-0.5])
-      assert_equal(c1.t,c2.t)
-      
-    end
-    should "respond to #summary" do
+      @v1 = [3, 3, 2, 3, 6].to_vector(:scale)
+      @v2 = [7, 6, 5, 6, 7].to_vector(:scale)
+      @v3 = [9, 8, 9, 7, 8].to_vector(:scale)
+      @name = 'Anova testing'
+      @anova = Statsample::Anova::OneWayWithVectors.new(@v1, @v2, @v3, name: @name)
+    end
+    should 'store correctly contrasts' do
+      c1 = Statsample::Anova::Contrast.new(vectors: [@v1, @v2, @v3], c: [1, -0.5, -0.5])
+
+      c2 = @anova.contrast(c: [1, -0.5, -0.5])
+      assert_equal(c1.t, c2.t)
+    end
+    should 'respond to #summary' do
       assert(@anova.respond_to? :summary)
     end
-    should "have correct name of analysis on #summary" do
+    should 'have correct name of analysis on #summary' do
       assert_match(/#{@name}/, @anova.summary)
     end
-    should "returns same levene values as direct Levene creation" do
-      assert_equal(@anova.levene.f, Statsample::Test.levene([@v1,@v2,@v3]).f)
+    should 'returns same levene values as direct Levene creation' do
+      assert_equal(@anova.levene.f, Statsample::Test.levene([@v1, @v2, @v3]).f)
     end
-    should "have correct value for levene" do
-      assert_in_delta(0.604,@anova.levene.f, 0.001)
-      assert_in_delta(0.562,@anova.levene.probability, 0.001)
+    should 'have correct value for levene' do
+      assert_in_delta(0.604, @anova.levene.f, 0.001)
+      assert_in_delta(0.562, @anova.levene.probability, 0.001)
     end
-    should "have correct value for sst" do
-     assert_in_delta(72.933, @anova.sst,0.001)
+    should 'have correct value for sst' do
+      assert_in_delta(72.933, @anova.sst, 0.001)
     end
-    should "have correct value for sswg" do
-      assert_in_delta(14.8,@anova.sswg,0.001)
+    should 'have correct value for sswg' do
+      assert_in_delta(14.8, @anova.sswg, 0.001)
     end
-    should "have correct value for ssb" do
-      assert_in_delta(58.133,@anova.ssbg,0.001)
+    should 'have correct value for ssb' do
+      assert_in_delta(58.133, @anova.ssbg, 0.001)
     end
-    should "sst=sswg+ssbg" do
-      assert_in_delta(@anova.sst,@anova.sswg+@anova.ssbg,0.00001)
+    should 'sst=sswg+ssbg' do
+      assert_in_delta(@anova.sst, @anova.sswg + @anova.ssbg, 0.00001)
     end
-    should "df total equal to number of n-1" do
-      assert_equal(@v1.n+@v2.n+@v3.n-1,@anova.df_total)
+    should 'df total equal to number of n-1' do
+      assert_equal(@v1.n + @v2.n + @v3.n - 1, @anova.df_total)
     end
-    should "df wg equal to number of n-k" do
-      assert_equal(@v1.n+@v2.n+@v3.n-3,@anova.df_wg)
+    should 'df wg equal to number of n-k' do
+      assert_equal(@v1.n + @v2.n + @v3.n - 3, @anova.df_wg)
     end
-    should "df bg equal to number of k-1" do
-      assert_equal(2,@anova.df_bg)
+    should 'df bg equal to number of k-1' do
+      assert_equal(2, @anova.df_bg)
     end
-    should "f=(ssbg/df_bg)/(sswt/df_wt)" do
-      assert_in_delta((@anova.ssbg.quo(@anova.df_bg)).quo( @anova.sswg.quo(@anova.df_wg)), @anova.f, 0.001)
+    should 'f=(ssbg/df_bg)/(sswt/df_wt)' do
+      assert_in_delta((@anova.ssbg.quo(@anova.df_bg)).quo(@anova.sswg.quo(@anova.df_wg)), @anova.f, 0.001)
     end
-    should "p be correct" do
-      assert(@anova.probability<0.01)
+    should 'p be correct' do
+      assert(@anova.probability < 0.01)
     end
-    should "be correct using different test values" do
-      anova2=Statsample::Anova::OneWayWithVectors.new([@v1,@v1,@v1,@v1,@v2])
-      assert_in_delta(3.960, anova2.f,0.001)
-      assert_in_delta(0.016, anova2.probability,0.001)
+    should 'be correct using different test values' do
+      anova2 = Statsample::Anova::OneWayWithVectors.new([@v1, @v1, @v1, @v1, @v2])
+      assert_in_delta(3.960, anova2.f, 0.001)
+      assert_in_delta(0.016, anova2.probability, 0.001)
     end
-    context "with extra information on summary" do
+    context 'with extra information on summary' do
       setup do
-        @anova.summary_descriptives=true
-        @anova.summary_levene=true
-        @summary=@anova.summary
+        @anova.summary_descriptives = true
+        @anova.summary_levene = true
+        @summary = @anova.summary
       end
-      should "have section with levene statistics" do
+      should 'have section with levene statistics' do
         assert_match(/Levene/, @summary)
       end
-      should "have section with descriptives" do
+      should 'have section with descriptives' do
         assert_match(/Min/, @summary)
       end
     end

data/test/test_awesome_print_bug.rb

@@ -1,14 +1,14 @@
-require(File.expand_path(File.dirname(__FILE__)+'/helpers_tests.rb'))
-class StatsampleAwesomePrintBug < MiniTest::Test
-  context("Awesome Print integration") do
+require(File.expand_path(File.dirname(__FILE__) + '/helpers_tests.rb'))
+class StatsampleAwesomePrintBug < Minitest::Test
+  context('Awesome Print integration') do
     setup do
-      require "awesome_print"
+      require 'awesome_print'
     end
-    should "should be flawless" do
-      a=[1,2,3].to_scale
-      
-      assert(a!=[1,2,3])
-      assert_nothing_raised do 
+    should 'should be flawless' do
+      a = [1, 2, 3].to_scale
+
+      assert(a != [1, 2, 3])
+      assert_nothing_raised do
         ap a
       end
     end

data/test/test_bartlettsphericity.rb

@@ -1,25 +1,25 @@
-require(File.expand_path(File.dirname(__FILE__)+'/helpers_tests.rb'))
+require(File.expand_path(File.dirname(__FILE__) + '/helpers_tests.rb'))
 
-class StatsampleBartlettSphericityTestCase < MiniTest::Test
+class StatsampleBartlettSphericityTestCase < Minitest::Test
   include Statsample::Test
   context Statsample::Test::BartlettSphericity do
     setup do
-      @v1=[1 ,2 ,3 ,4 ,7 ,8 ,9 ,10,14,15,20,50,60,70].to_scale
-      @v2=[5 ,6 ,11,12,13,16,17,18,19,20,30,0,0,0].to_scale
-      @v3=[10,3 ,20,30,40,50,80,10,20,30,40,2,3,4].to_scale
+      @v1 = [1, 2, 3, 4, 7, 8, 9, 10, 14, 15, 20, 50, 60, 70].to_scale
+      @v2 = [5, 6, 11, 12, 13, 16, 17, 18, 19, 20, 30, 0, 0, 0].to_scale
+      @v3 = [10, 3, 20, 30, 40, 50, 80, 10, 20, 30, 40, 2, 3, 4].to_scale
       # KMO: 0.490
-      ds={'v1'=>@v1,'v2'=>@v2,'v3'=>@v3}.to_dataset
-      cor=Statsample::Bivariate.correlation_matrix(ds)
-      @bs=Statsample::Test::BartlettSphericity.new(cor, 14)
+      ds = { 'v1' => @v1, 'v2' => @v2, 'v3' => @v3 }.to_dataset
+      cor = Statsample::Bivariate.correlation_matrix(ds)
+      @bs = Statsample::Test::BartlettSphericity.new(cor, 14)
     end
-    should "have correct value for chi" do
-      assert_in_delta(9.477, @bs.value,0.001)
+    should 'have correct value for chi' do
+      assert_in_delta(9.477, @bs.value, 0.001)
     end
-    should "have correct value for df" do
+    should 'have correct value for df' do
       assert_equal(3, @bs.df)
     end
-    should "have correct value for probability" do
-      assert_in_delta(0.024,@bs.probability,0.001)
+    should 'have correct value for probability' do
+      assert_in_delta(0.024, @bs.probability, 0.001)
     end
   end
 end