test_statistics.rb in statsample-0.3.4

- old
+ new

@@ -63,12 +63,27 @@
 	end
 #assert_equal(expected,obt)
     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
     end
 	def test_covariance
@@ -77,13 +92,15 @@
 			v2=[2,3,7,8,6,4,3,2].to_vector(:scale)
 			assert_in_delta(Statsample::Bivariate.covariance(v1,v2), Statsample::Bivariate.covariance_slow(v1,v2), 0.001)
 			
 		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)
 		d=[1,1,1,1,1,0,0,0,0,0].to_vector(:scale)
 		assert_raise TypeError do