statsample 0.13.1 → 0.14.0

data/test/test_rserve_extension.rb

@@ -12,7 +12,7 @@ class StatsampleRserveExtensionTestCase < MiniTest::Unit::TestCase
     should "return a valid rexp for numeric vector" do
       a=100.times.map {|i| rand()>0.9 ? nil : i+rand() }.to_scale
       rexp=a.to_REXP
-      assert_instance_of(rexp, Rserve::REXP::Double)
+      assert(rexp.is_a? Rserve::REXP::Double)
       assert_equal(rexp.to_ruby,a.data_with_nils)
       @r.assign 'a',rexp
       assert_equal(a.data_with_nils, @r.eval('a').to_ruby)
@@ -23,7 +23,7 @@ class StatsampleRserveExtensionTestCase < MiniTest::Unit::TestCase
       c=100.times.map {|i| rand()>0.9 ? nil : i+rand() }.to_scale
       ds={'a'=>a,'b'=>b,'c'=>c}.to_dataset
       rexp=ds.to_REXP
-      assert_instance_of(rexp, Rserve::REXP::GenericVector)
+      assert(rexp.is_a? Rserve::REXP::GenericVector)
       ret=rexp.to_ruby
       assert_equal(a.data_with_nils, ret['a']) 
       @r.assign 'df', rexp

data/test/test_stest.rb

@@ -2,15 +2,29 @@ require(File.dirname(__FILE__)+'/helpers_tests.rb')
 
 
 class StatsampleTestTestCase < MiniTest::Unit::TestCase
-  def test_chi_square
+  def test_chi_square_matrix_with_expected
     real=Matrix[[95,95],[45,155]]
     expected=Matrix[[68,122],[72,128]]
     assert_nothing_raised do
       chi=Statsample::Test.chi_square(real,expected)
     end
-    chi=Statsample::Test.chi_square(real,expected)
+    chi=Statsample::Test.chi_square(real,expected).chi_square
     assert_in_delta(32.53,chi,0.1)
+    
   end
+  def test_chi_square_matrix_only_observed
+    observed=Matrix[[20,30,40],[30,40,50],[60,70,80],[10,20,40]]
+    assert_nothing_raised do
+      chi=Statsample::Test.chi_square(observed)
+    end
+    chi=Statsample::Test.chi_square(observed)
+    assert_in_delta(9.5602, chi.chi_square, 0.0001)
+    assert_in_delta(0.1444, chi.probability, 0.0001)
+
+    assert_equal(6, chi.df)
+    
+  end
+  
   def test_u_mannwhitney
     a=[1,2,3,4,5,6].to_scale
     b=[0,5,7,9,10,11].to_scale

data/test/test_vector.rb

@@ -241,6 +241,11 @@ class StatsampleTestVector < MiniTest::Unit::TestCase
     v2=[1,2,3].to_vector()
     assert_equal(v1,v2)
   end
+  def test_vector_percentil
+    a=[1,2,2,3,4,5,5,5,6,10].to_scale
+    expected=[10,25,25,40,50,70,70,70,90,100].to_scale
+    assert_equal(expected, a.vector_percentil)
+  end
   def test_ordinal
     @c.type=:ordinal
     assert_equal(5,@c.median)
@@ -249,7 +254,11 @@ class StatsampleTestVector < MiniTest::Unit::TestCase
 
     v=[200000, 200000, 210000, 220000, 230000, 250000, 250000, 250000, 270000, 300000, 450000, 130000, 140000, 140000, 140000, 145000, 148000, 165000, 170000, 180000, 180000, 180000, 180000, 180000, 180000 ].to_scale
     assert_equal(180000,v.median)
-
+    a=[7.0, 7.0, 7.0, 7.0, 7.0, 8.0, 8.0, 8.0, 9.0, 9.0, 10.0, 10.0, 10.0, 10.0, 10.0, 12.0, 12.0, 13.0, 14.0, 14.0, 2.0, 3.0, 3.0, 3.0, 3.0, 4.0, 4.0, 4.0, 4.0, 4.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 6.0, 6.0, 6.0].to_scale
+    assert_equal(4.5, a.percentil(25))
+    assert_equal(6.5, a.percentil(50))
+    assert_equal(9.5, a.percentil(75))
+    assert_equal(3.0, a.percentil(10))
   end
   def test_ranked
     v1=[0.8,1.2,1.2,2.3,18].to_vector(:ordinal)

metadata

@@ -4,9 +4,9 @@ version: !ruby/object:Gem::Version
   prerelease: false
   segments: 
   - 0
-  - 13
-  - 1
-  version: 0.13.1
+  - 14
+  - 0
+  version: 0.14.0
 platform: ruby
 authors: 
 - Claudio Bustos
@@ -35,7 +35,7 @@ cert_chain:
   rpP0jjs0
   -----END CERTIFICATE-----
 
-date: 2010-07-03 00:00:00 -04:00
+date: 2010-08-16 00:00:00 -04:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -126,9 +126,9 @@ dependencies:
       - !ruby/object:Gem::Version 
         segments: 
         - 0
-        - 2
-        - 0
-        version: 0.2.0
+        - 3
+        - 1
+        version: 0.3.1
   type: :runtime
   version_requirements: *id007
 - !ruby/object:Gem::Dependency 
@@ -191,14 +191,15 @@ description: |-
   Include:
   * Descriptive statistics: frequencies, median, mean, standard error, skew, kurtosis (and many others).
   * Imports and exports datasets from and to Excel, CSV and plain text files.
-  * Correlations: Pearson's r, Spearman's rank correlation (rho), point biserial, tau a, tau b, gamma,  Tetrachoric and Polychoric.
+  * Correlations: Pearson's r, Spearman's rank correlation (rho), point biserial, tau a, tau b and  gamma.  Tetrachoric and Polychoric correlation provides by +statsample-bivariate-extension+ gem.
   * Anova: generic and vector-based One-way ANOVA and Two-way ANOVA
   * Tests: F, T, Levene, U-Mannwhitney.
   * Regression: Simple, Multiple (OLS), Probit  and Logit
-  * Factorial Analysis: Extraction (PCA and Principal Axis), Rotation (Varimax, Equimax, Quartimax) and Parallel Analysis, for estimation of number of factors.
+  * Factorial Analysis: Extraction (PCA and Principal Axis), Rotation (Varimax, Equimax, Quartimax) and Parallel Analysis and Velicer's MAP test, for estimation of number of factors.
   * Reliability analysis for simple scale and a DSL to easily analyze multiple scales using factor analysis and correlations, if you want it.
   * Dominance Analysis, with multivariate dependent and bootstrap (Azen & Budescu)
   * Sample calculation related formulas
+  * Structural Equation Modeling (SEM), using R libraries +sem+ and +OpenMx+
   * Creates reports on text, html and rtf, using ReportBuilder gem
 email: 
 - clbustos@gmail.com
@@ -219,6 +220,7 @@ files:
 - Rakefile
 - bin/statsample
 - data/crime.txt
+- data/hartman_23.matrix
 - data/locale/es/LC_MESSAGES/statsample.mo
 - data/repeated_fields.csv
 - data/test_binomial.csv
@@ -241,6 +243,7 @@ files:
 - examples/tetrachoric.rb
 - examples/u_test.rb
 - examples/vector.rb
+- examples/velicer_map_test.rb
 - lib/distribution.rb
 - lib/distribution/chisquare.rb
 - lib/distribution/f.rb
@@ -265,6 +268,7 @@ files:
 - lib/statsample/dominanceanalysis.rb
 - lib/statsample/dominanceanalysis/bootstrap.rb
 - lib/statsample/factor.rb
+- lib/statsample/factor/map.rb
 - lib/statsample/factor/parallelanalysis.rb
 - lib/statsample/factor/pca.rb
 - lib/statsample/factor/principalaxis.rb
@@ -301,6 +305,7 @@ files:
 - lib/statsample/rserve_extension.rb
 - lib/statsample/srs.rb
 - lib/statsample/test.rb
+- lib/statsample/test/chisquare.rb
 - lib/statsample/test/f.rb
 - lib/statsample/test/levene.rb
 - lib/statsample/test/t.rb