statsample 0.15.0 → 0.15.1

data/History.txt

@@ -1,3 +1,10 @@
+=== 0.15.1 / 2010-10-20
+* Statsample::Excel and Statsample::PlainText add name to vectors equal to field name 
+* Statsample::Dataset.delete_vector accept multiple fields.
+* Statsample::Dataset.dup_only_valid allows duplication of specific fields 
+* ScaleAnalysis doesn't crash on one-item scales 
+* Updated references
+
 === 0.15.0 / 2010-09-07
 * Added class Statsample::Reliability::ICC for calculation of Intra-class correlation (Shrout & Fleiss, 1979; McGraw & Wong, 1996). Tested with SPSS and R values.
 * References: Updated and standarized references on many classes. Added grab_references.rb script, to create a list of references for library

data/Rakefile

@@ -22,6 +22,7 @@ end
 task :release do
 system %{git push origin master}
 end
+
 desc "Update pot/po files."
 task "gettext:updatepo" do
   require 'gettext/tools'
@@ -48,14 +49,14 @@ h=Hoe.spec('statsample') do
 ***************************************************
 Thanks for installing statsample.
 
-On *nix, you should install statsample-optimization
+On *nix, you could install statsample-optimization
 to retrieve gems gsl, statistics2 and a C extension
 to speed some methods.
 
   $ sudo gem install statsample-optimization
 
 On Ubuntu, install  build-essential and libgsl0-dev 
-using apt-get and  compile ruby 1.8 or 1.9 from 
+using apt-get. Compile ruby 1.8 or 1.9 from 
 source code first.
 
   $ sudo apt-get install build-essential libgsl0-dev
@@ -106,4 +107,5 @@ task :publicar_docs => [:clean, :docs] do
   }
   sh %{rsync #{h.rsync_args} #{local_dir}/ #{host}:#{remote_dir}}
 end
+
 # vim: syntax=Ruby

data/lib/statsample.rb

@@ -115,7 +115,7 @@ module Statsample
     @@has_gsl
   end
   
-  VERSION = '0.15.0'
+  VERSION = '0.15.1'
   SPLIT_TOKEN = ","
   autoload(:Database, 'statsample/converters')
   autoload(:Anova, 'statsample/anova')

data/lib/statsample/converters.rb

@@ -76,18 +76,15 @@ module Statsample
   class SpreadsheetBase
     class << self
       def extract_fields(row)
-=begin
-        fields=[]
-        row.to_a.collect {|c|
+        i=0;
+        fields=row.to_a.collect{|c|
           if c.nil?
-            break
+            i+=1
+            "var%05d" % i 
           else
-            fields.push(c)
-          end
+            c.to_s.downcase
+          end        
         }
-=end
-raise "Should'nt be empty headers: [#{row.to_a.join(",")}]" if row.to_a.find_all {|c| c.nil?}.size>0
-        fields=row.to_a.collect{|c| c.downcase}
         fields.recode_repeated
       end
                                          
@@ -132,6 +129,9 @@ raise "Should'nt be empty headers: [#{row.to_a.join(",")}]" if row.to_a.find_all
           end
           convert_to_scale_and_date(ds,fields)
           ds.update_valid_data
+          fields.each {|f|
+            ds[f].name=f
+          }
           ds
         end
       end
@@ -220,6 +220,9 @@ raise "Should'nt be empty headers: [#{row.to_a.join(",")}]" if row.to_a.find_all
         end
         convert_to_scale_and_date(ds, fields)
         ds.update_valid_data
+        fields.each {|f|
+          ds[f].name=f
+        }
         ds
       end
     end

data/lib/statsample/dataset.rb

@@ -64,7 +64,7 @@ module Statsample
     # Ordered ids of vectors
     attr_reader :fields
     # Name of dataset
-    attr_accessor:name
+    attr_accessor :name
     # Number of cases
     attr_reader :cases
     # Location of pointer on enumerations methods (like #each)
@@ -157,27 +157,30 @@ module Statsample
 
     # Creates a copy of the given dataset, deleting all the cases with
     # missing data on one of the vectors
-    def dup_only_valid
-      if @vectors.any?{|field,vector| vector.has_missing_data?}
-        ds=dup_empty
-        each_array { |c|
-          ds.add_case_array(c) unless @fields.find{|f| @vectors[f].data_with_nils[@i].nil? }
+    # @param array of fields to include. No value include all fields
+    def dup_only_valid(*fields_to_include)
+      if fields_to_include.size==1 and fields_to_include[0].is_a? Array
+        fields_to_include=fields_to_include[0]
+      end
+      fields_to_include=@fields if fields_to_include.size==0
+      if fields_to_include.any? {|f| @vectors[f].has_missing_data?}
+        ds=Dataset.new(fields_to_include)
+        fields_to_include.each {|f| ds[f].type=@vectors[f].type}
+        each {|row|
+          unless fields_to_include.any? {|f| @vectors[f].has_missing_data? and !@vectors[f].is_valid? row[f]}
+            row_2=fields_to_include.inject({}) {|ac,v| ac[v]=row[v]; ac}
+            ds.add_case(row_2)
+          end
         }
-        ds.update_valid_data
       else
-        ds=dup()
+        ds=dup fields_to_include
       end
       ds
     end
-    # Returns an array with the fields from first argumen to last argument
-    def from_to(from,to)
-      raise ArgumentError, "Field #{from} should be on dataset" if !@fields.include? from
-      raise ArgumentError, "Field #{to} should be on dataset" if !@fields.include? to
-      @fields.slice(@fields.index(from)..@fields.index(to))
-    end
+    
     # Returns a duplicate of the Database
     # If fields given, only include those vectors.
-    # Every vector will be dup.
+    # @param array of fields to include. No value include all fields    
     def dup(*fields_to_include)
       if fields_to_include.size==1 and fields_to_include[0].is_a? Array
         fields_to_include=fields_to_include[0]
@@ -192,10 +195,20 @@ module Statsample
       }
       Dataset.new(vectors,fields)
     end
+    
+    
+    # Returns an array with the fields from first argumen to last argument
+    def from_to(from,to)
+      raise ArgumentError, "Field #{from} should be on dataset" if !@fields.include? from
+      raise ArgumentError, "Field #{to} should be on dataset" if !@fields.include? to
+      @fields.slice(@fields.index(from)..@fields.index(to))
+    end
+    
     # Returns (when possible) a cheap copy of dataset.
     # If no vector have missing values, returns original vectors.
     # If missing values presents, uses Dataset.dup_only_valid
-    #
+    # @param array of fields to include. No value include all fields
+
     def clone_only_valid(*fields_to_include)
       if fields_to_include.size==1 and fields_to_include[0].is_a? Array
         fields_to_include=fields_to_include[0]
@@ -209,6 +222,7 @@ module Statsample
     end
     # Returns a shallow copy of Dataset.
     # Object id will be distinct, but @vectors will be the same.
+    # @param array of fields to include. No value include all fields
     def clone(*fields_to_include)
       if fields_to_include.size==1 and fields_to_include[0].is_a? Array
         fields_to_include=fields_to_include[0]
@@ -267,7 +281,7 @@ module Statsample
       }
       Matrix.rows(rows)
     end
-    # We have the same datasets if vectors and fields are the same
+    # We have the same datasets if +vectors+ and +fields+ are the same
     def ==(d2)
       @vectors==d2.vectors and @fields==d2.fields
     end
@@ -336,10 +350,17 @@ module Statsample
       @fields.each{|f| @vectors[f].set_valid_data}
       check_length
     end
-    # Delete vector named <tt>name</tt>.
-    def delete_vector(name)
-      @fields.delete(name)
-      @vectors.delete(name)
+    # Delete vector named +name+. Multiple fields accepted.
+    def delete_vector(*args)
+      if args.size==1 and args[0].is_a? Array
+        names=args[0]
+      else
+        names=args
+      end
+      names.each do |name|
+        @fields.delete(name)
+        @vectors.delete(name)
+      end
     end
     
     def add_vectors_by_split_recode(name_,join='-',sep=Statsample::SPLIT_TOKEN)
@@ -493,6 +514,7 @@ module Statsample
         raise DatasetException.new(self, e)
       end
     end
+    
     # Returns each case as hash and index
     def each_with_index # :yield: |case, i|
       begin
@@ -530,11 +552,13 @@ module Statsample
       }
       @i=nil
     end
-    # Set fields order. If you omit one or more vectors,  
+    # Set fields order. If you omit one or more vectors, they are
+    # ordered by alphabetic order.
     def fields=(f)
       @fields=f
       check_order
     end
+    
     def check_order
       if(@vectors.keys.sort!=@fields.sort)
         @fields=@fields&@vectors.keys
@@ -640,16 +664,18 @@ module Statsample
     def to_multiset_by_split_one_field(field)
       raise ArgumentError,"Should use a correct field name" if !@fields.include? field
       factors=@vectors[field].factors
-      ms=Multiset.new_empty_vectors(@fields,factors)
+      ms=Multiset.new_empty_vectors(@fields, factors)
       each {|c|
         ms[c[field]].add_case(c,false)
       }
       #puts "Ingreso a los dataset"
       ms.datasets.each {|k,ds|
         ds.update_valid_data
+        ds.name=@vectors[field].labeling(k)
         ds.vectors.each{|k1,v1|
           #        puts "Vector #{k1}:"+v1.to_s
           v1.type=@vectors[k1].type
+          v1.name=@vectors[k1].name
         }
       }
       ms
@@ -675,7 +701,10 @@ module Statsample
       
       ms.datasets.each do |k,ds|
         ds.update_valid_data
-        ds.vectors.each{|k1,v1| v1.type=@vectors[k1].type }
+        ds.vectors.each{|k1,v1| 
+          v1.type=@vectors[k1].type
+          v1.name=@vectors[k1].name
+        }
       end
       ms
       

data/lib/statsample/reliability.rb

@@ -6,9 +6,11 @@ module Statsample
       def cronbach_alpha(ods)
         ds=ods.dup_only_valid
         n_items=ds.fields.size
+        return nil if n_items<=1
         s2_items=ds.vectors.inject(0) {|ac,v|
         ac+v[1].variance }
         total=ds.vector_sum
+        
         (n_items.quo(n_items-1)) * (1-(s2_items.quo(total.variance)))
       end
       # Calculate Chonbach's alpha for a given dataset
@@ -33,6 +35,7 @@ module Statsample
       # from +r+ current reliability, achieved with
       # +n+ items
       def n_for_desired_reliability(r,r_d,n=1)
+        return nil if r.nil?
         (r_d*(1-r)).quo(r*(1-r_d))*n
       end
       

data/lib/statsample/reliability/icc.rb

@@ -4,6 +4,17 @@ module Statsample
     # According to Shrout & Fleiss (1979, p.422): "ICC is the correlation 
     # between one measurement (either a single rating or a mean of 
     # several ratings) on a target and another measurement obtained on that target"
+    # == Usage
+    #   require 'statsample'
+    #   size=1000
+    #   a = size.times.map {rand(10)}.to_scale
+    #   b = a.recode{|i|i+rand(4)-2}
+    #   c  =a.recode{|i|i+rand(4)-2}
+    #   d = a.recode{|i|i+rand(4)-2}
+    #   ds={'a'=>a,'b'=>b,'c'=>c,'d'=>d}.to_dataset
+    #   # Use :type attribute to set type to summarize
+    #   icc=Statsample::Reliability::ICC.new(ds, :type=>:icc_1_k)
+    #   puts icc.summary
     # 
     # == Reference
     # * Shrout,P. & Fleiss, J. (1979). Intraclass Correlation: Uses in assessing rater reliability. Psychological Bulletin, 86(2), 420-428

data/lib/statsample/reliability/scaleanalysis.rb

@@ -15,14 +15,26 @@ module Statsample
       attr_reader :ds,:mean, :sd,:valid_n, :alpha , :alpha_standarized, :variances_mean, :covariances_mean, :cov_m
       attr_accessor :name
       def initialize(ds, opts=Hash.new)
-        @ds=ds.dup_only_valid
-        @k=@ds.fields.size
+        @dumped=ds.fields.find_all {|f|
+          ds[f].variance==0
+        }
+        
+        @ods=ds
+        @ds=ds.dup_only_valid(ds.fields - @dumped)
+        
+        
+        @k=@ds.fields.size        
         @total=@ds.vector_sum
+
+        @o_total=@dumped.size > 0 ? @ods.vector_sum : nil
+        
         @vector_mean=@ds.vector_mean
         @item_mean=@vector_mean.mean
         @item_sd=@vector_mean.sd
+        
         @mean=@total.mean
         @median=@total.median
+        
         @skew=@total.skew
         @kurtosis=@total.kurtosis
         @sd = @total.sd
@@ -37,8 +49,8 @@ module Statsample
         @variances_mean=@variances.mean
         @covariances_mean=(@variance-@variances.sum).quo(@k**2-@k)
         begin
-          @alpha = Statsample::Reliability.cronbach_alpha(ds)
-          @alpha_standarized = Statsample::Reliability.cronbach_alpha_standarized(ds)
+          @alpha = Statsample::Reliability.cronbach_alpha(@ds)
+          @alpha_standarized = Statsample::Reliability.cronbach_alpha_standarized(@ds)
         rescue => e
           raise DatasetException.new(@ds,e), "Error calculating alpha"
         end
@@ -62,7 +74,7 @@ module Statsample
         end
         total.each do |f,var|
           var.each do |tot,v|
-            out[f][tot]=out[f][tot].to_f / total[f][tot]
+            out[f][tot]=out[f][tot].quo(total[f][tot])
           end
         end
         out
@@ -164,7 +176,7 @@ module Statsample
         @sif||=stats_if_deleted_intern
       end
       def stats_if_deleted_intern # :nodoc:
-        
+        return Hash.new if @ds.fields.size==1
         @ds.fields.inject({}) do |a,v|
           cov_2=@cov_m.submatrix(@ds.fields-[v])
           #ds2=@ds.clone
@@ -182,12 +194,33 @@ module Statsample
       end
       def report_building(builder) #:nodoc:
         builder.section(:name=>@name) do |s|
+          
+          
+          if @dumped.size>0
+            s.section(:name=>"Items with variance=0") do |s1|
+              s.table(:name=>_("Summary for %s with all items") % @name) do |t|
+                t.row [_("Items"), @ods.fields.size]
+                t.row [_("Sum mean"),     "%0.4f" % @o_total.mean]
+                t.row [_("S.d. mean"),     "%0.4f" % @o_total.sd]
+
+              end
+            
+              s.table(:name=>_("Deleted items"), :header=>['item','mean']) do |t|
+                @dumped.each do |f|
+                  t.row(["#{@ods[f].name}(#{f})", "%0.5f" % @ods[f].mean])
+                end
+              end
+            end
+          end
+          
+          
           s.table(:name=>_("Summary for %s") % @name) do |t|
-          t.row [_("Items"), @ds.fields.size]
+            t.row [_("Valid Items"), @ds.fields.size]
+          
           t.row [_("Valid cases"), @valid_n]
           t.row [_("Sum mean"),     "%0.4f" % @mean]
           t.row [_("Sum sd"),       "%0.4f" % @sd  ]
-          t.row [_("Sum variance"), "%0.4f" % @variance]
+#          t.row [_("Sum variance"), "%0.4f" % @variance]
           t.row [_("Sum median"),   @median]
           t.hr
           t.row [_("Item mean"),    "%0.4f" % @item_mean]
@@ -196,26 +229,51 @@ module Statsample
           t.row [_("Skewness"),     "%0.4f" % @skew]
           t.row [_("Kurtosis"),     "%0.4f" % @kurtosis]
           t.hr
-          t.row [_("Cronbach's alpha"), "%0.4f" % @alpha]
-          t.row [_("Standarized Cronbach's alpha"), "%0.4f" % @alpha_standarized]
+          t.row [_("Cronbach's alpha"), @alpha ? ("%0.4f" % @alpha) : "--"]
+          t.row [_("Standarized Cronbach's alpha"), @alpha_standarized ? ("%0.4f" % @alpha_standarized) : "--" ]
           t.row [_("Mean rpb"), "%0.4f" % mean_rpb]
           
           t.row [_("Variances mean"),  "%g" % @variances_mean]
           t.row [_("Covariances mean") , "%g" % @covariances_mean]
           end
-          s.text _("Items for obtain alpha(0.8) : %d" % Statsample::Reliability::n_for_desired_reliability(@alpha, 0.8, @ds.fields.size))
-          s.text _("Items for obtain alpha(0.9) : %d" % Statsample::Reliability::n_for_desired_reliability(@alpha, 0.9, @ds.fields.size))          
+          
+          if (@alpha)
+            s.text _("Items for obtain alpha(0.8) : %d" % Statsample::Reliability::n_for_desired_reliability(@alpha, 0.8, @ds.fields.size))
+            s.text _("Items for obtain alpha(0.9) : %d" % Statsample::Reliability::n_for_desired_reliability(@alpha, 0.9, @ds.fields.size))          
+          end
+          
           
           sid=stats_if_deleted
           is=item_statistics
           itc=item_total_correlation
           
-          
           s.table(:name=>_("Items report for %s") % @name, :header=>["item","mean","sd", "mean if deleted", "var if deleted", "sd if deleted"," item-total correl.", "alpha if deleted"]) do |t|
             @ds.fields.each do |f|
-              t.row(["#{@ds[f].name}(#{f})", sprintf("%0.5f",is[f][:mean]), sprintf("%0.5f",is[f][:sds]), sprintf("%0.5f",sid[f][:mean]), sprintf("%0.5f",sid[f][:variance_sample]), sprintf("%0.5f",sid[f][:sds]),  sprintf("%0.5f",itc[f]), (sid[f][:alpha].nil?) ?  "--" : sprintf("%0.5f",sid[f][:alpha])])
+              row=["#{@ds[f].name}(#{f})"]
+              if is[f]
+                row+=[sprintf("%0.5f",is[f][:mean]), sprintf("%0.5f", is[f][:sds])]
+              else
+                row+=["-","-"]
+              end
+              if sid[f]
+                row+= [sprintf("%0.5f",sid[f][:mean]), sprintf("%0.5f",sid[f][:variance_sample]), sprintf("%0.5f",sid[f][:sds])]
+              else
+                row+=%w{- - -}
+              end
+              if itc[f]
+                row+= [sprintf("%0.5f",itc[f])]
+              else 
+                row+=['-']
+              end
+              if sid[f] and !sid[f][:alpha].nil?
+                row+=[sprintf("%0.5f",sid[f][:alpha])]
+              else
+                row+=["-"]
+              end
+              t.row row
             end # end each
           end # table
+          
         end # section
       end # def
     end # class

data/lib/statsample/vector.rb

@@ -105,12 +105,14 @@ module Statsample
     private :check_type
 
     # Return a vector usign the standarized values for data
-    # with sd with denominator n-1
+    # with sd with denominator n-1. With variance=0, returns nil
+    # 
     
     def vector_standarized(use_population=false)
       check_type :scale
       m=mean
       sd=use_population ? sdp : sds
+      return nil if sd==0.0
       @data_with_nils.collect{|x|
         if !x.nil?
           (x.to_f - m).quo(sd)
@@ -119,6 +121,7 @@ module Statsample
         end
       }.to_vector(:scale)
     end
+    
     alias_method :standarized, :vector_standarized
     # Return a vector with values replaced with the percentiles
     # of each values
@@ -623,7 +626,7 @@ module Statsample
         s.table(:name=>_("Distribution")) do |t|
           frequencies.sort.each do |k,v|
             key=labels.has_key?(k) ? labels[k]:k
-            t.row [key,v, ("%0.2f%%" % (v.quo(n_valid)*100))]
+            t.row [key, v , ("%0.2f%%" % (v.quo(n_valid)*100))]
           end
         end
         s.text _("median: %s") % median.to_s if(@type==:ordinal)

data/references.txt

@@ -2,7 +2,7 @@ References
 * Azen, R. & Budescu, D.V. (2003). The dominance analysis approach for comparing predictors in multiple regression. <em>Psychological Methods, 8</em>(2), 129-148.
 * Azen, R. & Budescu, D.V. (2006). Comparing predictors in Multivariate Regression Models: An extension of Dominance Analysis. <em>Journal of Educational and Behavioral Statistics, 31</em>(2), 157-180.
 * Budescu, D. V. (1993). Dominance analysis: a new approach to the problem of relative importance of predictors in multiple regression. <em>Psychological Bulletin, 114</em>, 542-551.
-* Cochran(1972)
+* Cochran, W.(1972). Sampling Techniques [spanish edition].
 * Cohen et al. (2003). Applied Multiple Reggression / Correlation Analysis for the Behavioral Sciences
 * Dinneen, L., & Blakesley, B. (1973). Algorithm AS 62: A Generator for the Sampling Distribution of the Mann- Whitney U Statistic. <em>Journal of the Royal Statistical Society, 22</em>(2), 269-273
 * Dziuban, C., & Shirkey E. (1974). When is a correlation matrix appropriate for factor analysis? Some decision rules. Psychological Bulletin, 81(6), 358-361.

data/test/test_dataset.rb

@@ -340,6 +340,10 @@ class StatsampleDatasetTestCase < MiniTest::Unit::TestCase
     expected=Statsample::Dataset.new({'v1'=>[1,4].to_vector(:scale), 'v2'=> [5,8].to_vector(:scale), 'v3'=>[9, 12].to_vector(:scale)})
     assert_equal(expected,ds2)
     assert_equal(expected.vectors.values,Statsample::only_valid(v1,v2,v3))
+    expected_partial=Statsample::Dataset.new({'v1'=>[1,3,4].to_vector(:scale), 'v3'=>[9, 11,12].to_vector(:scale)})
+    assert_equal(expected_partial, ds1.dup_only_valid(%w{v1 v3}))
+    
+    
   end
   def test_filter
     @ds['age'].type=:scale

data/test/test_reliability_icc.rb

@@ -105,12 +105,20 @@ class StatsampleReliabilityIccTestCase < MiniTest::Unit::TestCase
         assert_in_delta(0.676, @icc.icc_3_k_ci[0], 0.001)
         assert_in_delta(0.986, @icc.icc_3_k_ci[1], 0.001)
       end
+      should "incorrect type raises an error" do
+        assert_raise(::RuntimeError) do 
+          @icc.type=:nonexistant_type
+        end
+      end
     end
     
     begin
       require 'rserve'
       require 'statsample/rserve_extension'
       context "McGraw and Wong" do
+        teardown do
+          @r=$reliability_icc[:r].close
+        end
         setup do
           if($reliability_icc.nil?)
             size=100
@@ -135,15 +143,17 @@ class StatsampleReliabilityIccTestCase < MiniTest::Unit::TestCase
               icc_a_k=icc(ds,'t','a','a'))
               ")
             @iccs=@r.eval('iccs').to_ruby
-            $reliability_icc={ :icc=>@icc, :iccs=>@iccs
+            $reliability_icc={ :icc=>@icc, :iccs=>@iccs, :r=>@r
             }
+            
           end
           @icc=$reliability_icc[:icc]
           @iccs=$reliability_icc[:iccs]
-          
+          @r=$reliability_icc[:r]
+
         end
         [:icc_1, :icc_k, :icc_c_1, :icc_c_k, :icc_a_1, :icc_a_k].each do |t|
-          context "ICC Type #{t}" do
+          context "ICC Type #{t} " do
             should "value be correct" do
               @icc.type=t
               @r_icc=@iccs[t.to_s]
@@ -174,8 +184,10 @@ class StatsampleReliabilityIccTestCase < MiniTest::Unit::TestCase
               @icc.type=t
               @r_icc=@iccs[t.to_s]
               assert_in_delta(@r_icc['lbound'],@icc.lbound)
-              assert_in_delta(@r_icc['ubound'],@icc.ubound)
-              
+              assert_in_delta(@r_icc['ubound'],@icc.ubound)  
+            end
+            should "summary generated" do
+              assert(@icc.summary.size>0)
             end
           end
         end

data/test/test_rserve_extension.rb

@@ -9,6 +9,9 @@ class StatsampleRserveExtensionTestCase < MiniTest::Unit::TestCase
     setup do
       @r=Rserve::Connection.new
     end
+    teardown do
+      @r.close
+    end
     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

data/test/test_vector.rb

@@ -30,10 +30,10 @@ class StatsampleTestVector < MiniTest::Unit::TestCase
       @correct_b=@correct_b.to_scale
       
       @common=lambda  do |av,bv|
-        assert_equal(@correct_a,av)
-        assert_equal(@correct_b,bv)
-        assert(!av.has_missing_data?)
-        assert(!bv.has_missing_data?)        
+        assert_equal(@correct_a, av, "A no es esperado")
+        assert_equal(@correct_b, bv, "B no es esperado")
+        assert(!av.has_missing_data?, "A tiene datos faltantes")
+        assert(!bv.has_missing_data?, "b tiene datos faltantes")        
       end
     end
     should "return correct only_valid" do
@@ -285,6 +285,10 @@ class StatsampleTestVector < MiniTest::Unit::TestCase
     assert_equal(0,vs.mean)
     assert_equal(1,vs.sds)
   end
+  def test_vector_standarized_with_zero_variance
+    v1=100.times.map {|i| 1}.to_scale
+    assert(v1.standarized.nil?)
+  end
 
   def test_add
     a=Statsample::Vector.new([1,2,3,4,5], :scale)

metadata

@@ -5,8 +5,8 @@ version: !ruby/object:Gem::Version
   segments: 
   - 0
   - 15
-  - 0
-  version: 0.15.0
+  - 1
+  version: 0.15.1
 platform: ruby
 authors: 
 - Claudio Bustos
@@ -35,7 +35,7 @@ cert_chain:
   rpP0jjs0
   -----END CERTIFICATE-----
 
-date: 2010-09-07 00:00:00 -04:00
+date: 2010-10-20 00:00:00 -03:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -181,8 +181,8 @@ dependencies:
         segments: 
         - 2
         - 6
-        - 1
-        version: 2.6.1
+        - 2
+        version: 2.6.2
   type: :development
   version_requirements: *id011
 description: |-
@@ -371,14 +371,14 @@ post_install_message: |
   ***************************************************
   Thanks for installing statsample.
   
-  On *nix, you should install statsample-optimization
+  On *nix, you could install statsample-optimization
   to retrieve gems gsl, statistics2 and a C extension
   to speed some methods.
   
     $ sudo gem install statsample-optimization
   
   On Ubuntu, install  build-essential and libgsl0-dev 
-  using apt-get and  compile ruby 1.8 or 1.9 from 
+  using apt-get. Compile ruby 1.8 or 1.9 from 
   source code first.
   
     $ sudo apt-get install build-essential libgsl0-dev