statsample 0.13.1 → 0.14.0

data/po/statsample.pot

@@ -6,8 +6,8 @@
 #, fuzzy
 msgid ""
 msgstr ""
-"Project-Id-Version: statsample 0.13.0\n"
-"POT-Creation-Date: 2010-06-21 16:08-0400\n"
+"Project-Id-Version: statsample 0.13.1\n"
+"POT-Creation-Date: 2010-08-10 17:00-0400\n"
 "PO-Revision-Date: 2009-08-04 15:36-0400\n"
 "Last-Translator: FULL NAME <EMAIL@ADDRESS>\n"
 "Language-Team: LANGUAGE <LL@li.org>\n"
@@ -24,10 +24,8 @@ msgstr ""
 msgid "Mean and standard deviation"
 msgstr ""
 
-#: lib/statsample/test/t.rb:209 lib/statsample/factor/pca.rb:144
-#: lib/statsample/factor/pca.rb:149 lib/statsample/factor/pca.rb:155
-#: lib/statsample/factor/principalaxis.rb:180
-#: lib/statsample/factor/principalaxis.rb:185
+#: lib/statsample/test/t.rb:209 lib/statsample/factor/pca.rb:169
+#: lib/statsample/factor/principalaxis.rb:190
 msgid "Variable"
 msgstr ""
 
@@ -309,7 +307,6 @@ msgid "Number or factors to preserve: %d"
 msgstr ""
 
 #: lib/statsample/factor/parallelanalysis.rb:79
-#: lib/statsample/factor/pca.rb:149 lib/statsample/factor/principalaxis.rb:185
 msgid "Eigenvalues"
 msgstr ""
 
@@ -325,131 +322,207 @@ msgstr ""
 msgid "preserve?"
 msgstr ""
 
-#: lib/statsample/factor/pca.rb:43
+#: lib/statsample/factor/pca.rb:51
 msgid "Principal Component Analysis"
 msgstr ""
 
-#: lib/statsample/factor/pca.rb:143 lib/statsample/factor/principalaxis.rb:178
+#: lib/statsample/factor/pca.rb:118
+msgid "Component matrix"
+msgstr ""
+
+#: lib/statsample/factor/pca.rb:168 lib/statsample/factor/principalaxis.rb:188
 msgid "Number of factors: %d"
 msgstr ""
 
-#: lib/statsample/factor/pca.rb:144 lib/statsample/factor/principalaxis.rb:180
+#: lib/statsample/factor/pca.rb:169 lib/statsample/factor/principalaxis.rb:190
 msgid "Communalities"
 msgstr ""
 
-#: lib/statsample/factor/pca.rb:144 lib/statsample/factor/principalaxis.rb:180
+#: lib/statsample/factor/pca.rb:169 lib/statsample/factor/principalaxis.rb:190
 msgid "Initial"
 msgstr ""
 
-#: lib/statsample/factor/pca.rb:144 lib/statsample/factor/principalaxis.rb:180
+#: lib/statsample/factor/pca.rb:169 lib/statsample/factor/principalaxis.rb:190
 msgid "Extraction"
 msgstr ""
 
-#: lib/statsample/factor/pca.rb:149 lib/statsample/factor/principalaxis.rb:185
-msgid "Value"
+#: lib/statsample/factor/pca.rb:175
+msgid "Total Variance Explained"
+msgstr ""
+
+#: lib/statsample/factor/pca.rb:175
+msgid "Component"
+msgstr ""
+
+#: lib/statsample/factor/pca.rb:175
+msgid "E.Total"
+msgstr ""
+
+#: lib/statsample/factor/pca.rb:175
+msgid "%"
 msgstr ""
 
-#: lib/statsample/factor/pca.rb:155 lib/statsample/factor/principalaxis.rb:190
-msgid "Component Matrix"
+#: lib/statsample/factor/pca.rb:175
+msgid "Cum. %"
 msgstr ""
 
-#: lib/statsample/factor/principalaxis.rb:63
+#: lib/statsample/factor/pca.rb:179
+msgid "Component %d"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:64
 msgid "Variable %d"
 msgstr ""
 
-#: lib/statsample/factor/principalaxis.rb:179
+#: lib/statsample/factor/principalaxis.rb:140
+msgid "Factor Matrix"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:189
 msgid "Iterations: %d"
 msgstr ""
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:53
+#: lib/statsample/factor/principalaxis.rb:195
+msgid "Total Variance"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:195
+msgid "Factor"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:195
+msgid "I.E.Total"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:195
+msgid "I.E. %"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:195
+msgid "I.E.Cum. %"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:196
+msgid "S.L.Total"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:196
+msgid "S.L. %"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:196
+msgid "S.L.Cum. %"
+msgstr ""
+
+#: lib/statsample/factor/principalaxis.rb:203
+msgid "Factor %d"
+msgstr ""
+
+#: lib/statsample/factor/rotation.rb:35
+msgid "%s rotation"
+msgstr ""
+
+#: lib/statsample/factor/rotation.rb:125
+msgid "Rotated Component matrix"
+msgstr ""
+
+#: lib/statsample/factor/rotation.rb:142
+msgid "Component transformation matrix"
+msgstr ""
+
+#: lib/statsample/reliability/multiscaleanalysis.rb:57
 msgid "Multiple Scale analysis"
 msgstr ""
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:79
+#: lib/statsample/reliability/multiscaleanalysis.rb:85
 msgid "Scale %s"
 msgstr ""
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:110
+#: lib/statsample/reliability/multiscaleanalysis.rb:125
 msgid "Reliability analysis of scales"
 msgstr ""
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:116
+#: lib/statsample/reliability/multiscaleanalysis.rb:131
 msgid "Correlation matrix for %s"
 msgstr ""
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:121
+#: lib/statsample/reliability/multiscaleanalysis.rb:136
 msgid "PCA for %s"
 msgstr ""
 
-#: lib/statsample/reliability/multiscaleanalysis.rb:126
+#: lib/statsample/reliability/multiscaleanalysis.rb:141
 msgid "Principal Axis for %s"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:169
+#: lib/statsample/reliability/multiscaleanalysis.rb:147
+msgid "Parallel Analysis for %s"
+msgstr ""
+
+#: lib/statsample/reliability/scaleanalysis.rb:177
 msgid "Summary for %s"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:170
+#: lib/statsample/reliability/scaleanalysis.rb:178
 msgid "Items"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:171
+#: lib/statsample/reliability/scaleanalysis.rb:179
 msgid "Valid cases"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:172
+#: lib/statsample/reliability/scaleanalysis.rb:180
 msgid "Sum mean"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:173
+#: lib/statsample/reliability/scaleanalysis.rb:181
 msgid "Sum sd"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:174
+#: lib/statsample/reliability/scaleanalysis.rb:182
 msgid "Sum variance"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:175
+#: lib/statsample/reliability/scaleanalysis.rb:183
 msgid "Sum median"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:177
+#: lib/statsample/reliability/scaleanalysis.rb:185
 msgid "Item mean"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:178
+#: lib/statsample/reliability/scaleanalysis.rb:186
 msgid "Skewness"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:179
+#: lib/statsample/reliability/scaleanalysis.rb:187
 msgid "Kurtosis"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:181
+#: lib/statsample/reliability/scaleanalysis.rb:189
 msgid "Cronbach's alpha"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:182
+#: lib/statsample/reliability/scaleanalysis.rb:190
 msgid "Standarized Cronbach's alpha"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:183
+#: lib/statsample/reliability/scaleanalysis.rb:191
 msgid "Variances mean"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:184
+#: lib/statsample/reliability/scaleanalysis.rb:192
 msgid "Covariances mean"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:186
+#: lib/statsample/reliability/scaleanalysis.rb:194
 msgid "items for obtain alpha(0.8) : %d"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:187
+#: lib/statsample/reliability/scaleanalysis.rb:195
 msgid "items for obtain alpha(0.9) : %d"
 msgstr ""
 
-#: lib/statsample/reliability/scaleanalysis.rb:194
+#: lib/statsample/reliability/scaleanalysis.rb:202
 msgid "Items report for %s"
 msgstr ""
 
@@ -529,51 +602,51 @@ msgstr ""
 msgid "Y%d"
 msgstr ""
 
-#: lib/statsample/matrix.rb:110
+#: lib/statsample/matrix.rb:115
 msgid "Covariate matrix %d"
 msgstr ""
 
-#: lib/statsample/matrix.rb:152
+#: lib/statsample/matrix.rb:157
 msgid "Correlation"
 msgstr ""
 
-#: lib/statsample/matrix.rb:152
+#: lib/statsample/matrix.rb:157
 msgid "Covariance"
 msgstr ""
 
-#: lib/statsample/matrix.rb:152
+#: lib/statsample/matrix.rb:157
 msgid " Matrix"
 msgstr ""
 
-#: lib/statsample/vector.rb:614
+#: lib/statsample/vector.rb:620
 msgid "n :%d"
 msgstr ""
 
-#: lib/statsample/vector.rb:615
+#: lib/statsample/vector.rb:621
 msgid "n valid:%d"
 msgstr ""
 
-#: lib/statsample/vector.rb:616
+#: lib/statsample/vector.rb:622
 msgid "factors:%s"
 msgstr ""
 
-#: lib/statsample/vector.rb:617
+#: lib/statsample/vector.rb:623
 msgid "mode: %s"
 msgstr ""
 
-#: lib/statsample/vector.rb:618
+#: lib/statsample/vector.rb:624
 msgid "Distribution"
 msgstr ""
 
-#: lib/statsample/vector.rb:624
+#: lib/statsample/vector.rb:630
 msgid "median: %s"
 msgstr ""
 
-#: lib/statsample/vector.rb:626
+#: lib/statsample/vector.rb:632
 msgid "mean: %0.4f"
 msgstr ""
 
-#: lib/statsample/vector.rb:627
+#: lib/statsample/vector.rb:633
 msgid "sd: %0.4f"
 msgstr ""
 

data/test/test_factor.rb

@@ -1,6 +1,31 @@
 require(File.dirname(__FILE__)+'/helpers_tests.rb')
 
 class StatsampleFactorTestCase < MiniTest::Unit::TestCase
+  def test_parallelanalysis
+    pa=Statsample::Factor::ParallelAnalysis.with_random_data(305,8,150)
+    assert_in_delta(1.2454, pa.ds_eigenvalues['ev_00001'].mean, 0.01)
+    assert_in_delta(1.1542, pa.ds_eigenvalues['ev_00002'].mean, 0.01)
+    assert_in_delta(1.0836, pa.ds_eigenvalues['ev_00003'].mean, 0.01)
+  end
+  def test_map
+    fields=%w{height arm.span forearm lower.leg weight bitro.diameter chest.girth chest.width}
+    m=Matrix[ 
+          [ 1, 0.846, 0.805, 0.859, 0.473, 0.398, 0.301, 0.382],
+          [ 0.846, 1, 0.881, 0.826, 0.376, 0.326, 0.277, 0.415],
+          [ 0.805, 0.881, 1, 0.801, 0.38, 0.319, 0.237, 0.345],
+          [ 0.859, 0.826, 0.801, 1, 0.436, 0.329, 0.327, 0.365],
+          [ 0.473, 0.376, 0.38, 0.436, 1, 0.762, 0.73, 0.629],
+          [ 0.398, 0.326, 0.319, 0.329, 0.762, 1, 0.583, 0.577],
+          [ 0.301, 0.277, 0.237, 0.327, 0.73, 0.583, 1, 0.539],
+          [ 0.382, 0.415, 0.345, 0.365, 0.629, 0.577, 0.539, 1]
+    ]
+    map=Statsample::Factor::MAP.new(m)
+    assert_in_delta(map.minfm, 0.066445,0.00001)
+    assert_equal(map.number_of_factors, 2)
+    assert_in_delta(map.fm[0], 0.312475,0.00001)
+    assert_in_delta(map.fm[1], 0.245121,0.00001)
+
+  end
   # Tested with SPSS and R
   def test_pca
       a=[2.5, 0.5, 2.2, 1.9, 3.1, 2.3, 2.0, 1.0, 1.5, 1.1].to_scale
@@ -44,7 +69,7 @@ class StatsampleFactorTestCase < MiniTest::Unit::TestCase
   def test_principalaxis
       matrix=::Matrix[
       [1.0, 0.709501601093587, 0.877596585880047, 0.272219316266807],  [0.709501601093587, 1.0, 0.291633797330304, 0.871141831433844], [0.877596585880047, 0.291633797330304, 1.0, -0.213373722977167], [0.272219316266807, 0.871141831433844, -0.213373722977167, 1.0]]
-      fa=Statsample::Factor::PrincipalAxis.new(matrix,:m=>1)
+      fa=Statsample::Factor::PrincipalAxis.new(matrix,:m=>1, :max_iterations=>50)
       
       cm=::Matrix[[0.923],[0.912],[0.507],[0.483]]
       
@@ -56,12 +81,12 @@ class StatsampleFactorTestCase < MiniTest::Unit::TestCase
       }
       eigen1=2.175
       assert_in_delta(eigen1, fa.eigenvalues[0],0.001)
-
+      assert(fa.summary.size>0)
       fa=Statsample::Factor::PrincipalAxis.new(matrix,:smc=>false)
+            
       assert_raise RuntimeError do
         fa.iterate
       end
-      assert(fa.summary.size>0)
 
   end
 
@@ -82,6 +107,7 @@ class StatsampleFactorTestCase < MiniTest::Unit::TestCase
     assert(!varimax.h2.nil?, "H2 shouldn't be empty")
     
     _test_matrix(expected,varimax.rotated)
+    assert(varimax.summary.size>0)
   end
   def _test_matrix(a,b)
     a.row_size.times {|i|

data/test/test_matrix.rb

@@ -2,6 +2,8 @@ require(File.dirname(__FILE__)+'/helpers_tests.rb')
 
 
 class StatsampleMatrixTestCase < MiniTest::Unit::TestCase
+  
+  
   def test_covariate
     a=Matrix[[1.0, 0.3, 0.2], [0.3, 1.0, 0.5], [0.2, 0.5, 1.0]]
     a.extend Statsample::CovariateMatrix

data/test/test_reliability.rb

@@ -122,52 +122,52 @@ class StatsampleReliabilityTestCase < MiniTest::Unit::TestCase
           }
         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
+      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