statsample 0.11.2 → 0.12.0

data/test/test_bivariate_polychoric.rb

@@ -49,7 +49,7 @@ class StatsampleBivariatePolychoricTestCase < MiniTest::Unit::TestCase
       assert_in_delta(1.137, @poly.threshold_x[1],0.001)
       end
     end
-    if Statsample.has_gsl?
+    if true and Statsample.has_gsl?
     context "compute joint" do
       setup do
         @poly.method=:joint
@@ -57,11 +57,11 @@ class StatsampleBivariatePolychoricTestCase < MiniTest::Unit::TestCase
       end
       should "have correct values" do 
         assert_equal(:joint, @poly.method)
-      assert_in_delta(0.4192, @poly.r, 0.0001)
-      assert_in_delta(-0.2421, @poly.threshold_y[0],0.0001)
-      assert_in_delta(-0.0297, @poly.threshold_x[0],0.0001)
-      assert_in_delta(1.5938, @poly.threshold_y[1],0.0001)
-      assert_in_delta(1.1331, @poly.threshold_x[1],0.0001)
+        assert_in_delta(0.4192, @poly.r, 0.0001)
+        assert_in_delta(-0.2421, @poly.threshold_y[0],0.0001)
+        assert_in_delta(-0.0297, @poly.threshold_x[0],0.0001)
+        assert_in_delta(1.5938, @poly.threshold_y[1],0.0001)
+        assert_in_delta(1.1331, @poly.threshold_x[1],0.0001)
     end
     end
     end

data/test/test_distribution.rb

@@ -39,7 +39,7 @@ class DistributionTestCase < MiniTest::Unit::TestCase
   def test_normal_bivariate
     if Distribution.has_gsl?
       [0.2,0.4,0.6,0.8,0.9, 0.99,0.999,0.999999].each {|rho|
-        assert_equal(GSL::Ran::bivariate_gaussian_pdf(0, 0, 1,1,rho), Distribution::NormalBivariate.pdf(0,0, rho , 1,1))
+        assert_in_delta(GSL::Ran::bivariate_gaussian_pdf(0, 0, 1,1,rho), Distribution::NormalBivariate.pdf(0,0, rho , 1,1),1e-8)
 
       }
     end

data/test/test_reliability.rb

@@ -1,6 +1,4 @@
 require(File.dirname(__FILE__)+'/helpers_tests.rb')
-
-
 class StatsampleReliabilityTestCase < MiniTest::Unit::TestCase
   context Statsample::Reliability do
     context "Cronbach's alpha" do 
@@ -30,6 +28,19 @@ class StatsampleReliabilityTestCase < MiniTest::Unit::TestCase
         expected = @k.quo(@k-1) * (1-(ind_var.quo(total_sum)))
         assert_in_delta(expected, @a,1e-10)
       end
+      should "method cronbach_alpha_from_n_s2_cov return correct values" do
+        sa=Statsample::Reliability::ScaleAnalysis.new(@ds)
+        vm, cm = sa.variances_mean, sa.covariances_mean
+        assert_in_delta(sa.alpha, Statsample::Reliability.cronbach_alpha_from_n_s2_cov(@n_variables, vm,cm), 1e-10 )
+        
+      end
+      should "return correct n for desired alpha, covariance and variance" do
+        sa=Statsample::Reliability::ScaleAnalysis.new(@ds)
+        vm, cm = sa.variances_mean, sa.covariances_mean
+        n_obtained=Statsample::Reliability.n_for_desired_alpha(@a, vm,cm)
+        #p n_obtained
+        assert_in_delta(Statsample::Reliability.cronbach_alpha_from_n_s2_cov(n_obtained, vm,cm) ,@a,0.001) 
+      end
       should "standarized alpha will be equal to sum of matrix covariance less the individual variances on standarized values" do
         total_sum=@cme.total_sum
         ind_var=@dse.fields.inject(0) {|ac,v| ac+@dse[v].variance}
@@ -87,14 +98,82 @@ class StatsampleReliabilityTestCase < MiniTest::Unit::TestCase
       end
       
     end
-    context Statsample::Reliability::ItemAnalysis do
+    
+    context Statsample::Reliability::MultiScaleAnalysis do
+      
+      setup do
+        
+        size=100
+        @scales=4
+        @items_per_scale=10
+        h={}
+        @scales.times {|s|
+          @items_per_scale.times {|i|
+            h["#{s}_#{i}"] = (size.times.map {(s*2)+rand}).to_scale
+          }
+        }
+        @ds=h.to_dataset
+        @msa=Statsample::Reliability::MultiScaleAnalysis.new(:name=>'Multiple Analysis') do |m|
+          m.scale "complete", @ds
+          @scales.times {|s|
+            m.scale "scale_#{s}", @ds.clone(@items_per_scale.times.map {|i| "#{s}_#{i}"}), {:name=>"Scale #{s}"}
+          }
+        end
+      end
+        should "Retrieve correct ScaleAnalysis for whole scale" do
+          sa=Statsample::Reliability::ScaleAnalysis.new(@ds, :name=>"Scale complete") 
+          assert_equal(sa.variances_mean, @msa.scale("complete").variances_mean)
+        end
+        should "Retrieve correct ScaleAnalysis for each scale" do
+          @scales.times {|s|
+            sa=Statsample::Reliability::ScaleAnalysis.new(@ds.dup(@items_per_scale.times.map {|i| "#{s}_#{i}"}), :name=>"Scale #{s}")
+            assert_equal(sa.variances_mean,@msa.scale("scale_#{s}").variances_mean)
+          }
+        end
+        should "retrieve correct correlation matrix for each scale" do
+          vectors={'complete'=>@ds.vector_sum}
+          
+          @scales.times {|s|
+           vectors["scale_#{s}"]=@ds.dup(@items_per_scale.times.map {|i| "#{s}_#{i}"}).vector_sum 
+          }
+          ds2=vectors.to_dataset
+          assert_equal(Statsample::Bivariate.correlation_matrix(ds2), @msa.correlation_matrix)
+        end
+        should "delete scale using delete_scale" do
+          @msa.delete_scale("complete")
+          assert_equal(@msa.scales.keys.sort, @scales.times.map {|s| "scale_#{s}"})
+        end
+        should "retrieve pca for scales" do
+          @msa.delete_scale("complete")
+          vectors=Hash.new
+          @scales.times {|s|
+           vectors["scale_#{s}"]=@ds.dup(@items_per_scale.times.map {|i| "#{s}_#{i}"}).vector_sum 
+          }
+          ds2=vectors.to_dataset
+          cor_matrix=Statsample::Bivariate.correlation_matrix(ds2)
+          m=3
+          pca=Statsample::Factor::PCA.new(cor_matrix, :m=>m)
+          assert_equal(pca.component_matrix, @msa.pca(:m=>m).component_matrix)
+        end
+        should "retrieve acceptable summary" do
+          @msa.delete_scale("scale_0")
+          @msa.delete_scale("scale_1")
+          @msa.delete_scale("scale_2")
+          
+          
+          #@msa.summary_correlation_matrix=true
+          #@msa.summary_pca=true
+          assert(@msa.summary.size>0)
+        end
+    end
+    context Statsample::Reliability::ScaleAnalysis do
       setup do 
-        @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)
+        @x1=[1,1,1,1,2,2,2,2,3,3,3,30].to_scale
+        @x2=[1,1,1,2,2,3,3,3,3,4,4,50].to_scale
+        @x3=[2,2,1,1,1,2,2,2,3,4,5,40].to_scale
+        @x4=[1,2,3,4,4,4,4,3,4,4,5,30].to_scale
         @ds={'x1'=>@x1,'x2'=>@x2,'x3'=>@x3,'x4'=>@x4}.to_dataset
-        @ia=Statsample::Reliability::ItemAnalysis.new(@ds)
+        @ia=Statsample::Reliability::ScaleAnalysis.new(@ds)
         @cov_matrix=Statsample::Bivariate.covariance_matrix(@ds)
       end     
       should "return correct values for item analysis" do 

data/test/test_vector.rb

@@ -200,14 +200,6 @@ class StatsampleTestVector < MiniTest::Unit::TestCase
       
       assert_match(/#{@c.name}/, @c.summary)
     end
-    should "have output dependent of #type" do
-      @c.type=:nominal
-      assert_match(/Distribution/, @c.summary())
-      @c.type=:ordinal
-      assert_match(/median/, @c.summary())
-      @c.type=:scale
-      assert_match(/mean/, @c.summary())
-    end
     
     
   end

metadata

@@ -4,9 +4,9 @@ version: !ruby/object:Gem::Version
   prerelease: false
   segments: 
   - 0
-  - 11
-  - 2
-  version: 0.11.2
+  - 12
+  - 0
+  version: 0.12.0
 platform: ruby
 authors: 
 - Claudio Bustos
@@ -35,7 +35,7 @@ cert_chain:
   rpP0jjs0
   -----END CERTIFICATE-----
 
-date: 2010-05-05 00:00:00 -04:00
+date: 2010-06-14 00:00:00 -04:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -118,7 +118,7 @@ dependencies:
   type: :runtime
   version_requirements: *id006
 - !ruby/object:Gem::Dependency 
-  name: statistics2
+  name: extendmatrix
   prerelease: false
   requirement: &id007 !ruby/object:Gem::Requirement 
     requirements: 
@@ -126,66 +126,51 @@ dependencies:
       - !ruby/object:Gem::Version 
         segments: 
         - 0
-        - 54
-        version: "0.54"
+        - 2
+        - 0
+        version: 0.2.0
   type: :runtime
   version_requirements: *id007
 - !ruby/object:Gem::Dependency 
-  name: extendmatrix
+  name: rubyforge
   prerelease: false
   requirement: &id008 !ruby/object:Gem::Requirement 
     requirements: 
-    - - ~>
+    - - ">="
       - !ruby/object:Gem::Version 
         segments: 
-        - 0
         - 2
         - 0
-        version: 0.2.0
-  type: :runtime
+        - 4
+        version: 2.0.4
+  type: :development
   version_requirements: *id008
 - !ruby/object:Gem::Dependency 
-  name: gsl
+  name: shoulda
   prerelease: false
   requirement: &id009 !ruby/object:Gem::Requirement 
-    requirements: 
-    - - ~>
-      - !ruby/object:Gem::Version 
-        segments: 
-        - 1
-        - 12
-        - 109
-        version: 1.12.109
-  type: :runtime
-  version_requirements: *id009
-- !ruby/object:Gem::Dependency 
-  name: rubyforge
-  prerelease: false
-  requirement: &id010 !ruby/object:Gem::Requirement 
     requirements: 
     - - ">="
       - !ruby/object:Gem::Version 
         segments: 
-        - 2
         - 0
-        - 4
-        version: 2.0.4
+        version: "0"
   type: :development
-  version_requirements: *id010
+  version_requirements: *id009
 - !ruby/object:Gem::Dependency 
   name: hoe
   prerelease: false
-  requirement: &id011 !ruby/object:Gem::Requirement 
+  requirement: &id010 !ruby/object:Gem::Requirement 
     requirements: 
     - - ">="
       - !ruby/object:Gem::Version 
         segments: 
         - 2
         - 6
-        - 0
-        version: 2.6.0
+        - 1
+        version: 2.6.1
   type: :development
-  version_requirements: *id011
+  version_requirements: *id010
 description: |-
   A suite for basic and advanced statistics on Ruby. Tested on Ruby 1.8.7, 1.9.1, 1.9.2 (April, 2010) and JRuby 1.4 (Ruby 1.8.7 compatible).
   
@@ -197,6 +182,7 @@ description: |-
   * 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.
+  * Reliability analysis for simple scale and helpers to analyze multiple scales using factor analysis and correlations
   * Dominance Analysis, with multivariate dependent and bootstrap (Azen & Budescu)
   * Sample calculation related formulas
   * Creates reports on text, html and rtf, using ReportBuilder gem
@@ -224,6 +210,7 @@ files:
 - data/test_binomial.csv
 - data/tetmat_matrix.txt
 - data/tetmat_test.txt
+- doc_latex/manual/equations.tex
 - examples/correlation_matrix.rb
 - examples/dataset.rb
 - examples/dominance_analysis.rb
@@ -244,6 +231,7 @@ files:
 - lib/distribution/f.rb
 - lib/distribution/normal.rb
 - lib/distribution/normalbivariate.rb
+- lib/distribution/normalmultivariate.rb
 - lib/distribution/t.rb
 - lib/spss.rb
 - lib/statsample.rb
@@ -293,6 +281,8 @@ files:
 - lib/statsample/regression/multiple/rubyengine.rb
 - lib/statsample/regression/simple.rb
 - lib/statsample/reliability.rb
+- lib/statsample/reliability/multiscaleanalysis.rb
+- lib/statsample/reliability/scaleanalysis.rb
 - lib/statsample/resample.rb
 - lib/statsample/srs.rb
 - lib/statsample/test.rb
@@ -347,7 +337,25 @@ has_rdoc: true
 homepage: http://ruby-statsample.rubyforge.org/
 licenses: []
 
-post_install_message: 
+post_install_message: |
+  ***************************************************
+  Thanks for installing statsample.
+  
+  On *nix, you should install statsample-optimization
+  to retrieve gems gsl, statistics2 and a C extension
+  to speed some methods.
+  
+    $sudo gem install statsample-optimization
+  
+  To use it, on Ubuntu I recommend install 
+  build-essential and libgsl0-dev using apt-get and 
+  compile ruby 1.8 or 1.9 from source code first.
+  
+    $sudo apt-get install build-essential libgsl0-dev
+  
+  
+  *****************************************************
+
 rdoc_options: 
 - --main
 - README.txt