statsample 0.6.0 → 0.6.1

data/History.txt

@@ -1,3 +1,6 @@
+=== 0.6.1 / 2010-02-08
+* Bug fix on DominanceAnalysis summary for Ruby1.9
+* Some extra documentation
 === 0.6.0 / 2010-02-05
 * New Statsample::Factor module. Include classes for extracting factors (Statsample::Factor::PCA and  Statsample::Factor::PrincipalAxis) and rotate component matrix  ( Statsample::Factor::Rotation subclasses). For now, only orthogonal rotations
 * New Statsample::Dataset.crosstab_with_asignation, Statsample::Dataset.one_to_many

data/lib/statsample.rb

@@ -108,7 +108,7 @@ end
 # * Dataset: An union of vectors.
 #
 module Statsample
-  VERSION = '0.6.0'
+  VERSION = '0.6.1'
   SPLIT_TOKEN = ","
   autoload(:Database, 'statsample/converters')
   autoload(:Anova, 'statsample/anova')

data/lib/statsample/bivariate/tetrachoric.rb

@@ -27,11 +27,11 @@ module Statsample
     # See http://www.john-uebersax.com/stat/tetra.htm for extensive 
     # documentation about tetrachoric correlation.
     #
-    # This class uses algorithm AS116 from Applied Statistics(1977) 
-    # vol.26, no.3.
-    #
-    # You can see FORTRAN code on http://lib.stat.cmu.edu/apstat/116
+    # This class uses Brown(1977) algorithm. You can see FORTRAN code on http://lib.stat.cmu.edu/apstat/116
     # 
+    # == References:
+    # * Brown, MB. (1977) Algorithm AS 116: the tetrachoric correlation and its standard error. _Applied Statistics, 26_, 343-351.
+    #
     # <b>Usage</b>.
     # With two variables x and y on a crosstab like this:
     # 

data/lib/statsample/dominanceanalysis.rb

@@ -1,269 +1,321 @@
 require 'statsample/dominanceanalysis/bootstrap'
 module Statsample
-    class DominanceAnalysis
-        include GetText
-        bindtextdomain("statsample")
-        def initialize(ds,y_var, r_class = Regression::Multiple::RubyEngine)
-            @y_var=y_var
-            @dy=ds[@y_var]
-            @ds=ds
-            @r_class=r_class
-            @ds_indep=ds.dup(ds.fields-[y_var])
-            @fields=@ds_indep.fields
-            create_models
-            fill_models
-        end
-        def fill_models
-            @models.each{|m|
-                @fields.each{|f|
-                    next if m.include? f
-                    base_model=md(m)
-                    comp_model=md(m+[f])
-                    base_model.add_contribution(f,comp_model.r2)
-                }
-            }
-        end
-        def dominance_for_nil_model(i,j)
-            if md(i).r2>md(j).r2
-                1
-            elsif md(i).r2<md(j).r2
-                0
-            else
-                0.5
-            end           
-        end
-        # Returns 1 if i D k, 0 if j dominates i and 0.5 if undetermined
-        def total_dominance_pairwise(i,j)
-            dm=dominance_for_nil_model(i,j)
-            return 0.5 if dm==0.5
-            dominances=[dm]
-            @models_data.each{|k,m|
-              if !m.contributions[i].nil? and !m.contributions[j].nil?
-                if m.contributions[i]>m.contributions[j]
-                    dominances.push(1)
-                elsif m.contributions[i]<m.contributions[j]
-                    dominances.push(0)
-                else
-                  return 0.5
-                    #dominances.push(0.5)
-                end
-              end
-            }
-            final=dominances.uniq
-            final.size>1 ? 0.5 : final[0]
-        end
-        
-        # Returns 1 if i cD k, 0 if j cD i and 0.5 if undetermined
-        def conditional_dominance_pairwise(i,j)
-            dm=dominance_for_nil_model(i,j)
-            return 0.5 if dm==0.5
-            dominances=[dm]
-            for k in 1...@fields.size
-                a=average_k(k)
-                if a[i]>a[j]
-                    dominances.push(1)
-                elsif a[i]<a[j]
-                    dominances.push(0)
-                else
-                  return 0.5
-                    dominances.push(0.5)
-                end                 
-            end
-            final=dominances.uniq
-            final.size>1 ? 0.5 : final[0]            
+  # Dominance Analysis is a procedure based on an examination of the R2 values
+  # for all possible subset models, to identify the relevance of one or more 
+  # predictors in the prediction of criterium.
+  #
+  # See Budescu(1993) and Azen & Budescu (2003) for more information.
+  #  Use:
+  #  a=1000.times.collect {rand}.to_scale
+  #  b=1000.times.collect {rand}.to_scale
+  #  c=1000.times.collect {rand}.to_scale
+  #  ds={'a'=>a,'b'=>b,'c'=>c}.to_dataset
+  #  ds['y']=ds.collect{|row| row['a']*5+row['b']*3+row['c']*2+rand()}
+  #  da=Statsample::DominanceAnalysis.new(ds,'y')
+  #  puts da.summary
+  #  ==>
+  #  Resultado del Analisis de Dominancia de a, b, c en y
+  # 
+  #  ----------------------------------------------------------------
+  #  |                      | r2    | sign  | a     | b     | c     |
+  #  ----------------------------------------------------------------
+  #  |  Modelo 0            |       |       | 0.637 | 0.260 | 0.115 |
+  #  ----------------------------------------------------------------
+  #  |  a                   | 0.637 | 0.000 | --    | 0.239 | 0.109 |
+  #  |  b                   | 0.260 | 0.000 | 0.617 | --    | 0.103 |
+  #  |  c                   | 0.115 | 0.000 | 0.632 | 0.249 | --    |
+  #  ----------------------------------------------------------------
+  #  |  k=1 Promedio        |       |       | 0.624 | 0.244 | 0.106 |
+  #  ----------------------------------------------------------------
+  #  |  a*b                 | 0.877 | 0.000 | --    | --    | 0.098 |
+  #  |  a*c                 | 0.746 | 0.000 | --    | 0.229 | --    |
+  #  |  b*c                 | 0.363 | 0.000 | 0.612 | --    | --    |
+  #  ----------------------------------------------------------------
+  #  |  k=2 Promedio        |       |       | 0.612 | 0.229 | 0.098 |
+  #  ----------------------------------------------------------------
+  #  |  a*b*c               | 0.975 | 0.000 | --    | --    | --    |
+  #  ----------------------------------------------------------------
+  #  |  Promedios generales |       |       | 0.624 | 0.244 | 0.106 |
+  #  ----------------------------------------------------------------
+  #  
+  #  De a pares
+  #  
+  #  ----------------------------
+  #  |  Pares | T   | C   | G   |
+  #  ----------------------------
+  #  |  a - b | 1.0 | 1.0 | 1.0 |
+  #  |  a - c | 1.0 | 1.0 | 1.0 |
+  #  |  b - c | 1.0 | 1.0 | 1.0 |
+  #  ----------------------------
+  #
+  # == References:
+  # * Budescu, D. V. (1993). Dominance analysis: a new approach to the problem of relative importance of predictors in multiple regression. _Psychological Bulletin, 114_, 542-551.
+  # * Azen, R. & Budescu, D.V. (2003). The dominance analysis approach for comparing predictors in multiple regression. _Psychological Methods, 8_(2), 129-148.
+  class DominanceAnalysis
+    include GetText
+    bindtextdomain("statsample")
+    # Creates a new DominanceAnalysis object
+    # Params:
+    # * ds: A Dataset object
+    # * y_var: Name of dependent variable
+    # * r_class: Class to generate the regressions. Could be any subclass of 
+    #            Statsample::Regression::Multiple::BaseEngine
+    # 
+    def initialize(ds,y_var, r_class = Regression::Multiple::RubyEngine)
+    @y_var=y_var
+    @dy=ds[@y_var]
+    @ds=ds
+    @r_class=r_class
+    @ds_indep=ds.dup(ds.fields-[y_var])
+    @fields=@ds_indep.fields
+    create_models
+    fill_models
+    end
+    def fill_models
+      @models.each do |m|
+        @fields.each do |f|
+          next if m.include? f
+          base_model=md(m)
+          comp_model=md(m+[f])
+          base_model.add_contribution(f,comp_model.r2)
         end
-        # Returns 1 if i gD k, 0 if j gD i and 0.5 if undetermined        
-        def general_dominance_pairwise(i,j)
-          ga=general_averages
-          if ga[i]>ga[j]
-              1
-          elsif ga[i]<ga[j]
-              0
+      end
+    end
+    def dominance_for_nil_model(i,j)
+      if md([i]).r2>md([j]).r2
+        1
+      elsif md([i]).r2<md([j]).r2
+        0
+      else
+        0.5
+      end           
+    end
+    # Returns 1 if i D k, 0 if j dominates i and 0.5 if undetermined
+    def total_dominance_pairwise(i,j)
+      dm=dominance_for_nil_model(i,j)
+      return 0.5 if dm==0.5
+      dominances=[dm]
+      @models_data.each do |k,m|
+        if !m.contributions[i].nil? and !m.contributions[j].nil?
+          if m.contributions[i]>m.contributions[j]
+              dominances.push(1)
+          elsif m.contributions[i]<m.contributions[j]
+              dominances.push(0)
           else
-              0.5
-          end                 
-        end
-        def pairs
-          @models.find_all{|m| m.size==2}
-        end
-        def total_dominance
-            pairs.inject({}){|a,pair|
-              a[pair]=total_dominance_pairwise(pair[0], pair[1])
-              a
-            }
-        end
-        def conditional_dominance
-            pairs.inject({}){|a,pair|
-              a[pair]=conditional_dominance_pairwise(pair[0], pair[1])
-              a
-            }
-        end
-        def general_dominance
-            pairs.inject({}){|a,pair|
-              a[pair]=general_dominance_pairwise(pair[0], pair[1])
-              a
-            }
-        end
-        
-        def md(m)
-            @models_data[m.sort]
-        end
-        # Get all model of size k
-        def md_k(k)
-          out=[]
-          models=@models.each{|m| out.push(md(m)) if m.size==k }
-          out
-        end
-        
-        # For a hash with arrays of numbers as values
-        # Returns a hash with same keys and 
-        # value as the mean of values of original hash
-        
-        def get_averages(averages)
-          out={}
-          averages.each{|key,val| out[key]=val.to_vector(:scale).mean }
-          out
-        end
-        # Hash with average for each k size 
-        # model
-        def average_k(k)
-          return nil if k==@fields.size
-          models=md_k(k)
-          averages=@fields.inject({}) {|a,v| a[v]=[];a}
-          models.each do |m|
-            @fields.each do |f|
-              averages[f].push(m.contributions[f]) unless m.contributions[f].nil?
-            end
+            return 0.5
+              #dominances.push(0.5)
           end
-          get_averages(averages)
         end
-        def general_averages
-          if @general_averages.nil?
-            averages=@fields.inject({}) {|a,v| a[v]=[md(v).r2];a}
-            for k in 1...@fields.size
-                ak=average_k(k)
-                @fields.each{|f|
-                    averages[f].push(ak[f])
-                }
-            end
-            @general_averages=get_averages(averages)
-          end
-          @general_averages
+      end
+      final=dominances.uniq
+      final.size>1 ? 0.5 : final[0]
+    end
+    
+    # Returns 1 if i cD k, 0 if j cD i and 0.5 if undetermined
+    def conditional_dominance_pairwise(i,j)
+      dm=dominance_for_nil_model(i,j)
+      return 0.5 if dm==0.5
+      dominances=[dm]
+      for k in 1...@fields.size
+        a=average_k(k)
+        if a[i]>a[j]
+            dominances.push(1)
+        elsif a[i]<a[j]
+            dominances.push(0)
+        else
+          return 0.5
+            dominances.push(0.5)
+        end                 
+      end
+      final=dominances.uniq
+      final.size>1 ? 0.5 : final[0]            
+    end
+    # Returns 1 if i gD k, 0 if j gD i and 0.5 if undetermined        
+    def general_dominance_pairwise(i,j)
+      ga=general_averages
+      if ga[i]>ga[j]
+        1
+      elsif ga[i]<ga[j]
+        0
+      else
+        0.5
+      end                 
+    end
+    def pairs
+      @models.find_all{|m| m.size==2}
+    end
+    def total_dominance
+      pairs.inject({}){|a,pair| a[pair]=total_dominance_pairwise(pair[0], pair[1])
+      a
+      }
+    end
+    def conditional_dominance
+      pairs.inject({}){|a,pair|
+      a[pair]=conditional_dominance_pairwise(pair[0], pair[1])
+      a
+      }
+    end
+    def general_dominance
+      pairs.inject({}){|a,pair|
+      a[pair]=general_dominance_pairwise(pair[0], pair[1])
+      a
+      }
+    end
+    
+    def md(m)
+      @models_data[m.sort]
+    end
+    # Get all model of size k
+    def md_k(k)
+      out=[]
+      models=@models.each{|m| out.push(md(m)) if m.size==k }
+      out
+    end
+    
+    # For a hash with arrays of numbers as values
+    # Returns a hash with same keys and 
+    # value as the mean of values of original hash
+    
+    def get_averages(averages)
+      out={}
+      averages.each{|key,val| out[key]=val.to_vector(:scale).mean }
+      out
+    end
+    # Hash with average for each k size model.
+    def average_k(k)
+      return nil if k==@fields.size
+      models=md_k(k)
+      averages=@fields.inject({}) {|a,v| a[v]=[];a}
+      models.each do |m|
+        @fields.each do |f|
+          averages[f].push(m.contributions[f]) unless m.contributions[f].nil?
         end
-        def create_models
-          @models=[]
-          @models_data={}
-          for i in 1..@fields.size
-              c=Statsample::Combination.new(i,@fields.size)
-              c.each{|data|
-                  convert=data.collect {|i1|
-                      @fields[i1]
-                  }
-                  @models.push(convert)
-                  ds_prev=@ds.dup(convert+[@y_var])
-                  modeldata=ModelData.new(convert,ds_prev, @y_var, @fields, @r_class)
-                  @models_data[convert.sort]=modeldata
-              }
+      end
+      get_averages(averages)
+    end
+    def general_averages
+      if @general_averages.nil?
+        averages=@fields.inject({}) {|a,v| a[v]=[md([v]).r2];a}
+        for k in 1...@fields.size
+          ak=average_k(k)
+          @fields.each do |f|
+            averages[f].push(ak[f])
           end
         end
-        def summary(report_type=ConsoleSummary)
-            out=""
-            out.extend report_type
-            out << _("Summary for Dominance Analysis of %s on %s\n") % [@fields.join(", "),@y_var]
-                t=Statsample::ReportTable.new
-                t.header=["","r2","sign"]+@fields
-                row=[_("Model 0"),"",""]+@fields.collect{|f|
-                    sprintf("%0.3f",md(f).r2)
-                }
-                t.add_row(row)
-                t.add_horizontal_line
-            for i in 1..@fields.size
-                mk=md_k(i)
-                mk.each{|m|
-                    t.add_row(m.add_table_row)
-                }
-                # Report averages
-                a=average_k(i)
-                if !a.nil?
-                    t.add_horizontal_line
-                    row=[_("k=%d Average") % i,"",""] + @fields.collect{|f|
-                        sprintf("%0.3f",a[f])
-                    }
-                    t.add_row(row)
-                    t.add_horizontal_line
-                    
-                end
-                
-            end
-            
-            g=general_averages
-	     t.add_horizontal_line
-            
-            row=[_("Overall averages"),"",""]+@fields.collect{|f|
-                        sprintf("%0.3f",g[f])
-            }
-            t.add_row(row)
-            out.parse_table(t)
-                    
-            out.nl
-            out << _("Pairwise")+"\n"
-            td=total_dominance
-            cd=conditional_dominance
-            gd=general_dominance
-            t=Statsample::ReportTable.new([_("Pairs"),"T","C","G"])
-            pairs.each{|p|
-                name=p.join(" - ")
-                row=[name, sprintf("%0.1f",td[p]), sprintf("%0.1f",cd[p]), sprintf("%0.1f",gd[p])]
-                t.add_row(row)
-            }
-            out.parse_table(t)
-            return out
-        end
-            class ModelData
-            attr_reader :contributions
-            def initialize(name,ds,y_var,fields,r_class)
-                @name=name
-                @fields=fields
-                @contributions=@fields.inject({}){|a,v| a[v]=nil;a}
-                r_class=Regression::Multiple::RubyEngine if r_class.nil?
-                @lr=r_class.new(ds,y_var)
-            end
-            def add_contribution(f,v)
-                @contributions[f]=v-r2
-            end
-            def r2
-                @lr.r2
-            end
-            def add_table_row
-                begin
-                sign=sprintf("%0.3f", @lr.significance)
-                rescue RuntimeError
-                    sign="???"
-                end
-                [@name.join("*"), sprintf("%0.3f",r2), sign] + @fields.collect{|k|
-                    v=@contributions[k]
-                    if v.nil?
-                        "--"
-                    else
-                    sprintf("%0.3f",v)
-                    end
-                }
-            end
-            def summary
-                out=sprintf("%s: r2=%0.3f(p=%0.2f)\n",@name.join("*"),r2,@lr.significance,@lr.sst)
-                out << @fields.collect{|k|
-                    v=@contributions[k]
-                    if v.nil?
-                        "--"
-                    else
-                    sprintf("%s=%0.3f",k,v)
-                    end
-                }.join(" | ") 
-                out << "\n"
-                
-                return out
-            end
-            end
+        @general_averages=get_averages(averages)
+      end
+      @general_averages
+    end
+    def create_models
+      @models=[]
+      @models_data={}
+      for i in 1..@fields.size
+      c=Statsample::Combination.new(i,@fields.size)
+      c.each  do |data|
+        convert=data.collect {|i1| @fields[i1] }
+        @models.push(convert)
+        ds_prev=@ds.dup(convert+[@y_var])
+        modeldata=ModelData.new(convert,ds_prev, @y_var, @fields, @r_class)
+        @models_data[convert.sort]=modeldata
+      end
+      end
+    end
+    def summary(report_type=ConsoleSummary)
+    out=""
+    out.extend report_type
+    out << _("Summary for Dominance Analysis of %s on %s\n") % [@fields.join(", "),@y_var]
+      t=Statsample::ReportTable.new
+      t.header=["","r2","sign"]+@fields
+      row=[_("Model 0"),"",""]+@fields.collect{|f|
+        sprintf("%0.3f", md([f]).r2)
+      }
+      t.add_row(row)
+      t.add_horizontal_line
+    for i in 1..@fields.size
+      mk=md_k(i)
+      mk.each{|m|
+        t.add_row(m.add_table_row)
+      }
+      # Report averages
+      a=average_k(i)
+      if !a.nil?
+          t.add_horizontal_line
+          row=[_("k=%d Average") % i,"",""] + @fields.collect{|f|
+              sprintf("%0.3f",a[f])
+          }
+          t.add_row(row)
+          t.add_horizontal_line
+          
+      end
+      
     end
     
+    g=general_averages
+    t.add_horizontal_line
+    
+    row=[_("Overall averages"),"",""]+@fields.collect{|f|
+              sprintf("%0.3f",g[f])
+    }
+    t.add_row(row)
+    out.parse_table(t)
+          
+    out.nl
+    out << _("Pairwise")+"\n"
+    td=total_dominance
+    cd=conditional_dominance
+    gd=general_dominance
+    t=Statsample::ReportTable.new([_("Pairs"),"T","C","G"])
+    pairs.each{|p|
+      name=p.join(" - ")
+      row=[name, sprintf("%0.1f",td[p]), sprintf("%0.1f",cd[p]), sprintf("%0.1f",gd[p])]
+      t.add_row(row)
+    }
+    out.parse_table(t)
+    return out
+    end
+    class ModelData
+      attr_reader :contributions
+      def initialize(name,ds,y_var,fields,r_class)
+        @name=name
+        @fields=fields
+        @contributions=@fields.inject({}){|a,v| a[v]=nil;a}
+        r_class=Regression::Multiple::RubyEngine if r_class.nil?
+        @lr=r_class.new(ds,y_var)
+      end
+      def add_contribution(f,v)
+        @contributions[f]=v-r2
+      end
+      def r2
+        @lr.r2
+      end
+      def add_table_row
+        begin
+        sign=sprintf("%0.3f", @lr.significance)
+        rescue RuntimeError
+            sign="???"
+        end
+        [@name.join("*"), sprintf("%0.3f",r2), sign] + @fields.collect{|k|
+          v=@contributions[k]
+          if v.nil?
+              "--"
+          else
+          sprintf("%0.3f",v)
+          end
+        }
+      end
+      def summary
+        out=sprintf("%s: r2=%0.3f(p=%0.2f)\n",@name.join("*"),r2,@lr.significance,@lr.sst)
+        out << @fields.collect{|k|
+          v=@contributions[k]
+          if v.nil?
+              "--"
+          else
+            sprintf("%s=%0.3f",k,v)
+          end
+        }.join(" | ") 
+        out << "\n"
+        return out
+      end
+    end # end ModelData
+  end # end Dominance Analysis
 end

data/lib/statsample/dominanceanalysis/bootstrap.rb

@@ -1,5 +1,8 @@
 module Statsample
 class DominanceAnalysis
+  # Generates Bootstrap sample to identity the replicability of a Dominance Analysis. See Azen & Bodescu (2003) for more information.
+  # References: 
+  # * Azen, R. & Budescu, D.V. (2003). The dominance analysis approach for comparing predictors in multiple regression. _Psychological Methods, 8_(2), 129-148.
     class Bootstrap
         include GetText
         include Writable

data/lib/statsample/regression/multiple.rb

@@ -3,18 +3,44 @@ module Statsample
   module Regression
     # Module for Linear Multiple Regression Analysis.
     # 
-    # You can call Regression::Multiple.listwise or Regression::Multiple.pairwise or instance directly the engines.
+    # You can call Statsample::Regression::Multiple.listwise,  Statsample::Regression::Multiple.pairwise or instance directly the engines.
     # 
-    # Example.
+    #  Use:.
     #
     #  require 'statsample'
-    #  @a=[1,3,2,4,3,5,4,6,5,7].to_vector(:scale)
-    #  @b=[3,3,4,4,5,5,6,6,4,4].to_vector(:scale)
-    #  @c=[11,22,30,40,50,65,78,79,99,100].to_vector(:scale)
-    #  @y=[3,4,5,6,7,8,9,10,20,30].to_vector(:scale)
-    #  ds={'a'=>@a,'b'=>@b,'c'=>@c,'y'=>@y}.to_dataset
-    #  lr=Statsample::Regression::Multiple.listwise(ds,'y')        
-    #  #<Statsample::Regression::Multiple::AlglibEngine:0x7f21912e4758 @ds_valid=#<Statsample::Dataset:69891073182680 @fields=[a,b,c,y] labels={"a"=>nil, "b"=>nil, "y"=>nil, "c"=>nil} cases=10, @lr=#<Alglib::LinearRegression:0x7f21912df118 @model=#<Alglib_ext::LinearModel:0x7f21912df708>, @ivars=3, @cases=10, @report=#<Alglib_ext::LrReport:0x7f21912df168>>, @y_var="y", @ds=#<Statsample::Dataset:69891073182680 @fields=[a,b,c,y] labels={"a"=>nil, "b"=>nil, "y"=>nil, "c"=>nil} cases=10, @fields=["a", "b", "c"], @lr_s=nil, @dep_columns=[[1, 3, 2, 4, 3, 5, 4, 6, 5, 7], [3, 3, 4, 4, 5, 5, 6, 6, 4, 4], [11, 22, 30, 40, 50, 65, 78, 79, 99, 100]], @ds_indep=#<Statsample::Dataset:69891073180060 @fields=[a,b,c] labels={"a"=>nil, "b"=>nil, "c"=>nil} cases=10, @dy=Vector(type:scale, n:10)[3,4,5,6,7,8,9,10,20,30]>
+    #  a=1000.times.collect {rand}.to_scale
+    #  b=1000.times.collect {rand}.to_scale
+    #  c=1000.times.collect {rand}.to_scale
+    #  ds={'a'=>a,'b'=>b,'c'=>c}.to_dataset
+    #  ds['y']=ds.collect{|row| row['a']*5+row['b']*3+row['c']*2+rand()}
+    #  lr=Statsample::Regression::Multiple.listwise(ds,'y')
+    #  puts lr.summary
+    #  Summary for regression of a,b,c over y
+    #  *************************************************************
+    #  Engine: Statsample::Regression::Multiple::AlglibEngine
+    #  Cases(listwise)=1000(1000)
+    #  r=0.986
+    #  r2=0.973
+    #  Equation=0.504+5.011a + 2.995b + 1.988c
+    #  ----------------------------
+    #  ANOVA TABLE
+    #  --------------------------------------------------------------
+    #  |  source     | ss       | df  | ms      | f         | s     |
+    #  --------------------------------------------------------------
+    #  |  Regression | 2979.321 | 3   | 993.107 | 12040.067 | 0.000 |
+    #  |  Error      | 82.154   | 996 | 0.082   |           |       |
+    #  |  Total      | 3061.475 | 999 |         |           |       |
+    #  --------------------------------------------------------------
+    #  Beta coefficientes
+    #  -----------------------------------------------
+    #  |  coeff    | b     | beta  | se    | t       |
+    #  -----------------------------------------------
+    #  |  Constant | 0.504 | -     | 0.030 | 16.968  |
+    #  |  a        | 5.011 | 0.832 | 0.031 | 159.486 |
+    #  |  b        | 2.995 | 0.492 | 0.032 | 94.367  |
+    #  |  c        | 1.988 | 0.323 | 0.032 | 62.132  |
+    #  -----------------------------------------------
+    # 
     module Multiple
         # Creates an object for listwise regression. 
         # Alglib is faster, so is prefered over GSL

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification 
 name: statsample
 version: !ruby/object:Gem::Version 
-  version: 0.6.0
+  version: 0.6.1
 platform: ruby
 authors: 
 - Claudio Bustos
@@ -9,7 +9,7 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2010-02-06 00:00:00 -03:00
+date: 2010-02-08 00:00:00 -03:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency