rstat 0.1.0 → 0.1.1

data/.rspec

@@ -0,0 +1,3 @@
+--tty
+--colour
+--format p

data/CHANGELOG.md

@@ -1,7 +1,15 @@
+## v0.1.1
+
+* Made rstat.rb load all files automatically instead of listing them manually.
+* Updated readme.
+* Added sample correlation coefficient.
+* Added simple linear regression.
+
 ## v0.1.0
 
 * Major reorganization! All files now in directories named after their corresponding discipline (as sort of defined by Wikipedia).
 * Added skewness.
+* Added kurtosis.
 * Mean methods now now use the sum method.
 * Minor reorganization of tests.
 * Changed Gemfile source to https.

data/Gemfile

@@ -1,4 +1,4 @@
-source "https://rubygems.org"
+source 'https://rubygems.org'
 
 # Specify your gem's dependencies in rstat.gemspec
 gemspec

data/README.md

@@ -1,12 +1,94 @@
-## Rstat
+# Rstat
 
-A very simple gem that adds some statistics methods to the core Array class.
+A very simple statistics gem.
 
-Maybe in the future it will be something more.
+## Installation
 
-Currently adds:
+    $ gem install rstat
 
-* Mean
+## Usage
+
+```ruby
+require 'rstat'
+ => true
+[1, 2, 3, 4, 5].mean
+ => 3.0
+```
+
+## Currently Adds
+
+### Arithmetic
+
+* Sum
+* Product
+* Binomial Coefficient
+
+### Descriptive Statistics
+
+#### Location
+
+* Arithmetic Mean
+* Geometric Mean
+* Harmonic Mean
+* Power Mean
 * Median
 * Mode
-* Variance / Standard Deviation
+
+#### Dispersion
+
+* Range
+* Standard Deviation
+* Coefficient of Variation
+* Percentile
+* Interquartile Range
+
+#### Shape
+
+* Variance
+* Skewness
+* Kurtosis
+* Central Moments
+* Raw Moments
+
+#### Count Data
+
+* Index of Dispersion
+
+### Regression Analysis
+
+#### Linear Regression
+
+* Simple Linear Regression
+
+## Running Tests
+
+Rstat depends upon RSpec for its tests.
+
+To run all the tests for the gem:
+
+    $ rspec spec
+
+To run one particular spec file:
+
+    $ rspec spec/rstat/array_spec.rb
+
+To run one particular line of a spec file:
+
+    $ rspec spec/rstat/array_spec.rb:177
+
+Rstat has been tested with Ruby 1.8.7 and 1.9.3. There's little reason it shouldn't work with other versions of Ruby though.
+
+## Console
+
+Rstat includes a rake task to open up an irb session with the gem source preloaded. This is particularly useful if you're working with a local development copy of the gem's source and want to experiment with your changes. To use the preloaded console:
+
+    $ rake console
+
+## Contributing
+
+If you feel like you can add something useful to Rstat then don't hesitate to send a pull request.
+
+
+## A Note of Warning
+
+This gem extends the core Array class and Fixnum class. In isolation this is usually pretty harmless. But, in combination with other gems that do the same, unpredictable behavior may result. As always, use caution, and be aware of what this gem and any others you use actually do before including it in an important project.

data/Rakefile

@@ -1,11 +1,11 @@
-require "bundler/gem_tasks"
-require "rspec/core/rake_task"
+require 'bundler/gem_tasks'
+require 'rspec/core/rake_task'
 
 RSpec::Core::RakeTask.new(:spec)
 
 task :default => :spec
 
-desc "Open an irb session preloaded with this library"
+desc 'Open an irb session preloaded with this library'
 task :console do
-  sh "irb -rubygems -I lib -r rstat.rb"
+  sh 'irb -rubygems -I lib -r rstat.rb'
 end

data/lib/rstat.rb

@@ -1,5 +1,6 @@
-require "rstat/version"
-require "rstat/core_ext"
+Dir["#{File.dirname(__FILE__)}/rstat/*.rb"].sort.each do |path|
+  require "rstat/#{File.basename(path, '.rb')}"
+end
 
 module Rstat
 end

data/lib/rstat/core_ext.rb

@@ -1,3 +1,3 @@
 Dir["#{File.dirname(__FILE__)}/core_ext/*.rb"].sort.each do |path|
-  require "rstat/core_ext/#{File.basename(path, ".rb")}"
+  require "rstat/core_ext/#{File.basename(path, '.rb')}"
 end

data/lib/rstat/core_ext/array.rb

@@ -1,3 +1,3 @@
 Dir["#{File.dirname(__FILE__)}/array/*.rb"].sort.each do |path|
-  require "rstat/core_ext/array/#{File.basename(path, ".rb")}"
+  require "rstat/core_ext/array/#{File.basename(path, '.rb')}"
 end

data/lib/rstat/core_ext/array/descriptive_statistics.rb

@@ -1,3 +1,3 @@
 Dir["#{File.dirname(__FILE__)}/descriptive_statistics/*.rb"].sort.each do |path|
-  require "rstat/core_ext/array/descriptive_statistics/#{File.basename(path, ".rb")}"
+  require "rstat/core_ext/array/descriptive_statistics/#{File.basename(path, '.rb')}"
 end

data/lib/rstat/core_ext/array/descriptive_statistics/dispersion.rb

@@ -1,3 +1,3 @@
 Dir["#{File.dirname(__FILE__)}/dispersion/*.rb"].sort.each do |path|
-  require "rstat/core_ext/array/descriptive_statistics/dispersion/#{File.basename(path, ".rb")}"
+  require "rstat/core_ext/array/descriptive_statistics/dispersion/#{File.basename(path, '.rb')}"
 end

data/lib/rstat/core_ext/array/descriptive_statistics/dispersion/percentile.rb

@@ -1,5 +1,5 @@
 class Array
-  def percentile p
+  def percentile(p)
     if p < 0 || p > 100
       nil
     else
@@ -22,4 +22,4 @@ class Array
   def middle_fifty
     self.interquartile_range
   end
-end
+end

data/lib/rstat/core_ext/array/descriptive_statistics/dispersion/range.rb

@@ -2,4 +2,4 @@ class Array
   def range
     self.max - self.min
   end
-end
+end

data/lib/rstat/core_ext/array/descriptive_statistics/index_of_dispersion.rb

@@ -22,4 +22,4 @@ class Array
   def fano_factor
     self.index_of_dispersion
   end
-end
+end

data/lib/rstat/core_ext/array/descriptive_statistics/location.rb

@@ -1,3 +1,3 @@
 Dir["#{File.dirname(__FILE__)}/location/*.rb"].sort.each do |path|
-  require "rstat/core_ext/array/descriptive_statistics/location/#{File.basename(path, ".rb")}"
+  require "rstat/core_ext/array/descriptive_statistics/location/#{File.basename(path, '.rb')}"
 end

data/lib/rstat/core_ext/array/descriptive_statistics/shape.rb

@@ -1,3 +1,3 @@
 Dir["#{File.dirname(__FILE__)}/shape/*.rb"].sort.each do |path|
-  require "rstat/core_ext/array/descriptive_statistics/shape/#{File.basename(path, ".rb")}"
+  require "rstat/core_ext/array/descriptive_statistics/shape/#{File.basename(path, '.rb')}"
 end

data/lib/rstat/core_ext/array/descriptive_statistics/shape/kurtosis.rb

@@ -2,4 +2,4 @@ class Array
   def kurtosis
     (self.central_moment(4) / self.variance ** 2) - 3
   end
-end
+end

data/lib/rstat/core_ext/array/descriptive_statistics/shape/moment.rb

@@ -1,5 +1,5 @@
 class Array
-  def central_moment k
+  def central_moment(k)
     if k < 0
       nil
     else
@@ -9,7 +9,7 @@ class Array
     end
   end
 
-  def raw_moment k
+  def raw_moment(k)
     if k < 0
       nil
     else
@@ -20,4 +20,4 @@ class Array
   # def l_moment k
   #
   # end
-end
+end

data/lib/rstat/core_ext/array/descriptive_statistics/shape/skewness.rb

@@ -2,4 +2,4 @@ class Array
   def skewness
     self.central_moment(3) / (self.variance ** 1.5)
   end
-end
+end

data/lib/rstat/core_ext/array/descriptive_statistics/shape/variance.rb

@@ -2,4 +2,4 @@ class Array
   def variance
     self.central_moment 2
   end
-end
+end

data/lib/rstat/core_ext/array/product.rb

@@ -2,4 +2,4 @@ class Array
   def product
     self.inject(:*)
   end
-end
+end

data/lib/rstat/core_ext/array/sum.rb

@@ -2,4 +2,4 @@ class Array
   def sum
     self.inject(:+)
   end
-end
+end

data/lib/rstat/core_ext/fixnum.rb

@@ -1,3 +1,3 @@
 Dir["#{File.dirname(__FILE__)}/fixnum/*.rb"].sort.each do |path|
-  require "rstat/core_ext/fixnum/#{File.basename(path, ".rb")}"
+  require "rstat/core_ext/fixnum/#{File.basename(path, '.rb')}"
 end

data/lib/rstat/core_ext/fixnum/binomial_coefficient.rb

@@ -1,5 +1,5 @@
 class Fixnum
-  def binomial_coefficient bottom
+  def binomial_coefficient(bottom)
     (self - bottom + 1 .. self).inject(1, &:*) / (2 .. bottom).inject(1, &:*)
   end
-end
+end

data/lib/rstat/regression_analysis.rb

@@ -0,0 +1,3 @@
+Dir["#{File.dirname(__FILE__)}/regression_analysis/*.rb"].sort.each do |path|
+  require "rstat/regression_analysis/#{File.basename(path, '.rb')}"
+end

data/lib/rstat/regression_analysis/linear_regression.rb

@@ -0,0 +1,3 @@
+Dir["#{File.dirname(__FILE__)}/linear_regression/*.rb"].sort.each do |path|
+  require "rstat/regression_analysis/linear_regression/#{File.basename(path, '.rb')}"
+end

data/lib/rstat/regression_analysis/linear_regression/simple_linear_regression.rb

@@ -0,0 +1,43 @@
+module Rstat
+  def self.simple_linear_regression(x, y)
+    if x.length != y.length
+      return nil
+    end
+
+    n = x.length
+
+    xy = x.zip(y).map { |a| a[0] * a[1] }
+
+    xx = x.map { |a| a * a }
+
+    sumx = x.sum
+
+    sumy = y.sum
+
+    sumxy = xy.sum
+
+    sumxx = xx.sum
+
+    slope = ((n * sumxy) - (sumx * sumy)) / ((n * sumxx) - (sumx ** 2))
+
+    intercept = (sumy - (slope * sumx)) / n
+
+    { :slope => slope, :intercept => intercept }
+  end
+
+  def self.simple_linear_regression_slope(x, y)
+    if x.length != y.length
+      return nil
+    end
+
+    Rstat.simple_linear_regression(x, y)[:slope]
+  end
+
+  def self.simple_linear_regression_intercept(x, y)
+    if x.length != y.length
+      return nil
+    end
+
+    Rstat.simple_linear_regression(x, y)[:intercept]
+  end
+end

data/lib/rstat/regression_analysis/sample_correlation_coefficient.rb

@@ -0,0 +1,20 @@
+module Rstat
+  def self.sample_correlation_coefficient(x, y)
+    if x.length != y.length
+      return nil
+    end
+
+    xm = x.mean
+    ym = y.mean
+
+    top = x.zip(y).map { |a, b| (a - xm) * (b - ym) }.inject(:+)
+
+    bottom = Math.sqrt(x.map { |a| (a - xm) ** 2 }.inject(:+)) * Math.sqrt(y.map { |b| (b - ym) ** 2 }.inject(:+))
+
+    top / bottom
+  end
+
+  def self.pearsons_correlation_coefficient(x, y)
+    Rstat.sample_correlation_coefficient x, y
+  end
+end

data/lib/rstat/version.rb

@@ -1,3 +1,3 @@
 module Rstat
-  VERSION = "0.1.0"
+  VERSION = '0.1.1'
 end

data/rstat.gemspec

@@ -1,22 +1,22 @@
 # -*- encoding: utf-8 -*-
-$:.push File.expand_path("../lib", __FILE__)
-require "rstat/version"
+$:.push File.expand_path('../lib', __FILE__)
+require 'rstat/version'
 
 Gem::Specification.new do |s|
-  s.name        = "rstat"
+  s.name        = 'rstat'
   s.version     = Rstat::VERSION
-  s.authors     = ["Sean Eshbaugh"]
-  s.email       = ["seaneshbaugh@gmail.com"]
-  s.homepage    = "http://seaneshbaugh.com/"
+  s.authors     = ['Sean Eshbaugh']
+  s.email       = ['seaneshbaugh@gmail.com']
+  s.homepage    = 'https://github.com/seaneshbaugh/rstat'
   s.summary     = %q{A Simple statistics gem.}
   s.description = %q{A Simple statistics gem.}
 
-  s.rubyforge_project = "rstat"
+  s.rubyforge_project = 'rstat'
 
   s.files         = `git ls-files`.split("\n")
   s.test_files    = `git ls-files -- {test,spec,features}/*`.split("\n")
   s.executables   = `git ls-files -- bin/*`.split("\n").map{ |f| File.basename(f) }
-  s.require_paths = ["lib"]
+  s.require_paths = ['lib']
 
-  s.add_development_dependency "rspec"
+  s.add_development_dependency 'rspec'
 end

data/spec/rstat/array_spec.rb

@@ -2,47 +2,47 @@
 require 'spec_helper'
 
 describe Rstat do
-  describe ".sum" do
-    it "calculates the sum of an array" do
+  describe '.sum' do
+    it 'calculates the sum of an array' do
       [1, 2, 3, 4, 5, 6].sum.should be(21)
     end
   end
 
-  describe ".product" do
-    it "calculates the product of an array" do
+  describe '.product' do
+    it 'calculates the product of an array' do
       [1, 2, 3, 4, 5, 6].product.should be(720)
     end
   end
 
-  context "Location" do
-    describe ".mean" do
-        it "calculates the mean of an array" do
+  context 'Location' do
+    describe '.mean' do
+        it 'calculates the mean of an array' do
           [1, 2, 3, 4, 5].mean.should be_within(0.0001).of(3.0000)
         end
     end
 
-    describe ".geometric_mean" do
-      it "calculates the geometric mean of an array" do
+    describe '.geometric_mean' do
+      it 'calculates the geometric mean of an array' do
         [1, 2, 3, 4, 5].geometric_mean.should be_within(0.0001).of(2.6052)
       end
 
-      it "calculates the geometric mean of an array with a zero element" do
+      it 'calculates the geometric mean of an array with a zero element' do
         [0, 2, 3, 4, 5].geometric_mean.should eql(0.0)
       end
 
-      it "calculates the geometric mean of an array with a negative element" do
+      it 'calculates the geometric mean of an array with a negative element' do
         [-1, 2, 3, 4, 5].geometric_mean.real?.should be_false
       end
     end
 
-    describe ".harmonic mean" do
-      it "calculates the harmonic mean of an array" do
+    describe '.harmonic mean' do
+      it 'calculates the harmonic mean of an array' do
         [1, 2, 3, 4, 5].harmonic_mean.should be_within(0.0001).of(2.1897)
       end
     end
 
-    describe ".power_mean" do
-      it "power_mean(2) should equal quadratic_mean" do
+    describe '.power_mean' do
+      it 'power_mean(2) should equal quadratic_mean' do
         p = [1, 2, 3, 4, 5].power_mean(2)
         q = [1, 2, 3, 4, 5].quadratic_mean
 
@@ -50,7 +50,7 @@ describe Rstat do
       end
     end
 
-    it "AM > GM > HM" do
+    it 'AM > GM > HM' do
       a = [1, 2, 3, 4, 5]
       am = a.arithmetric_mean
       gm = a.geometric_mean
@@ -60,142 +60,142 @@ describe Rstat do
       gm.should be > hm
     end
 
-    describe ".median" do
-      it "calculates the median of an array with an odd number of elements" do
+    describe '.median' do
+      it 'calculates the median of an array with an odd number of elements' do
         [2, 3, 5, 1, 4].median.should eql(3)
       end
 
-      it "calculates the median of an array with an even number of elements" do
+      it 'calculates the median of an array with an even number of elements' do
         [2, 3, 5, 1, 4, 6].median.should eql(3.5)
       end
     end
 
-    describe ".mode" do
-      it "calculates the mode of an empty array" do
+    describe '.mode' do
+      it 'calculates the mode of an empty array' do
         [].mode.should eql([])
       end
 
-      it "calculates the mode of an array with one element" do
+      it 'calculates the mode of an array with one element' do
         [100].mode.should eql([100])
       end
 
-      it "calculates the mode of an unimodal array" do
+      it 'calculates the mode of an unimodal array' do
         [1, 1, 1, 2, 3, 4, 5, 5, 6].mode.should eql([1])
       end
 
-      it "calculates the mode of a bimodal array" do
+      it 'calculates the mode of a bimodal array' do
         [44, 45, 46, 44, 46, 50].mode.should eql([44, 46])
       end
     end
   end
 
-  context "Dispersion" do
-    describe ".range" do
-      it "finds the range of an array" do
+  context 'Dispersion' do
+    describe '.range' do
+      it 'finds the range of an array' do
         [0, 34, 656, 400, 1000].range.should be(1000)
       end
     end
 
-    describe ".standard_deviation" do
-      it "calculates the standard deviation of an array" do
+    describe '.standard_deviation' do
+      it 'calculates the standard deviation of an array' do
         [1, 2, 3, 4, 5, 6].standard_deviation.should be_within(0.00001).of(1.70783)
       end
     end
 
-    describe ".coefficient_of_variation" do
-      it "calculates the coefficient of variation of an array" do
+    describe '.coefficient_of_variation' do
+      it 'calculates the coefficient of variation of an array' do
         [1, 2, 3, 4, 5, 6].coefficient_of_variation.should be_within(0.00001).of(0.48795)
       end
     end
 
-    describe ".percentile" do
-      it "finds the 30th percentile of an array" do
+    describe '.percentile' do
+      it 'finds the 30th percentile of an array' do
         [1, 2, 3, 4, 5, 6, 7, 8, 9, 10].percentile(30).should be(3)
       end
 
-      it "finds the 50th percentile of an array" do
+      it 'finds the 50th percentile of an array' do
         [1, 2, 3, 4, 5].percentile(50).should be(3)
       end
 
-      it "finds the 10th percentile of an array" do
+      it 'finds the 10th percentile of an array' do
         [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 15, 16, 17].percentile(10).should be(2)
       end
     end
 
-    describe ".interquartile_range" do
-      it "finds the interquartile range of an array" do
+    describe '.interquartile_range' do
+      it 'finds the interquartile range of an array' do
         [1, 2, 3, 4, 5, 6, 7, 8, 9, 10].interquartile_range.should be(5)
       end
     end
   end
 
-  context "Shape" do
-    describe ".central_moment" do
-      it "calculates the 0th central moment of an array" do
+  context 'Shape' do
+    describe '.central_moment' do
+      it 'calculates the 0th central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].central_moment(0).should be_within(0.00001).of(1.0)
       end
 
-      it "calculates the 1st central moment of an array" do
+      it 'calculates the 1st central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].central_moment(1).should be_within(0.00001).of(0.0)
       end
 
-      it "calculates the 2nd central moment of an array" do
+      it 'calculates the 2nd central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].central_moment(2).should be_within(0.00001).of(3459.44000)
       end
 
-      it "calculates the 3rd central moment of an array" do
+      it 'calculates the 3rd central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].central_moment(3).should be_within(0.00001).of(321046.00800)
       end
 
-      it "calculates the 4th central moment of an array" do
+      it 'calculates the 4th central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].central_moment(4).should be_within(0.00001).of(56245095.69920)
       end
     end
 
-    describe ".raw_moment" do
-      it "calculates the 0th central moment of an array" do
+    describe '.raw_moment' do
+      it 'calculates the 0th central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].raw_moment(0).should be_within(0.00001).of(1.0)
       end
 
-      it "calculates the 1st central moment of an array" do
+      it 'calculates the 1st central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].raw_moment(1).should be_within(0.00001).of([2, 43, 5, 65, 56, 7, 85, 3, -1, 199].mean)
       end
 
-      it "calculates the 2nd central moment of an array" do
+      it 'calculates the 2nd central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].raw_moment(2).should be_within(0.00001).of(5612.40000)
       end
 
-      it "calculates the 3rd central moment of an array" do
+      it 'calculates the 3rd central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].raw_moment(3).should be_within(0.00001).of(902497.40000)
       end
 
-      it "calculates the 4th central moment of an array" do
+      it 'calculates the 4th central moment of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].raw_moment(4).should be_within(0.00001).of(165154687.20000)
       end
     end
 
-    describe ".variance" do
-      it "calculates the variance of an array" do
+    describe '.variance' do
+      it 'calculates the variance of an array' do
         [1, 2, 3, 4, 5, 6].variance.should be_within(0.00001).of(2.91667)
       end
     end
 
-    describe ".skewness" do
-      it "finds the skewness of an array" do
+    describe '.skewness' do
+      it 'finds the skewness of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].skewness.should be_within(0.00001).of(1.57782)
       end
     end
 
-    describe ".kurtosis" do
-      it "finds the kurtosis of an array" do
+    describe '.kurtosis' do
+      it 'finds the kurtosis of an array' do
         [2, 43, 5, 65, 56, 7, 85, 3, -1, 199].kurtosis.should be_within(0.00001).of(1.69973)
       end
     end
   end
 
-  context "Count Data" do
-    describe ".index_of_dispersion" do
-      it "finds the index of dispersion of an array" do
+  context 'Count Data' do
+    describe '.index_of_dispersion' do
+      it 'finds the index of dispersion of an array' do
         [1, 2, 3, 4, 5, 6].index_of_dispersion.should be_within(0.00001).of(0.83333)
       end
     end

data/spec/rstat/fixnum_spec.rb

@@ -2,11 +2,11 @@
 require 'spec_helper'
 
 describe Rstat do
-  describe ".binomial_coefficient" do
-    it "computes the binomial coeffecient of itself and another number" do
+  describe '.binomial_coefficient' do
+    it 'computes the binomial coeffecient of itself and another number' do
       10.binomial_coefficient(4).should eq(210)
 
       50.binomial_coefficient(34).should eq(4_923_689_695_575)
     end
   end
-end
+end

data/spec/rstat/regession_analysis_spec.rb

@@ -0,0 +1,38 @@
+# encoding: utf-8
+require 'spec_helper'
+
+describe Rstat do
+  describe '.sample_correlation_coefficient' do
+    it 'computes the sample correlation coefficient of two arrays' do
+      Rstat.sample_correlation_coefficient([1, 2, 3, 4, 5], [6, 7, 8, 9, 10]).should be_within(0.00001).of(1.0)
+
+      Rstat.sample_correlation_coefficient([1, 2, 3, 4, 5], [5, 4, 3, 8, 1]).should be_within(0.00001).of(-0.24434)
+
+      Rstat.sample_correlation_coefficient([2345, 43, 346, 45, -4234], [354, 66, 100, 443, 655]).should be_within(0.00001).of(-0.59889)
+    end
+  end
+
+  describe '.simple_linear_regression' do
+    it 'should return a Hash with the slope and intercept of the linear regression equation for two arrays' do
+      result  = Rstat.simple_linear_regression([60, 61, 62, 63, 65], [3.1, 3.6, 3.8, 4, 4.1])
+
+      result.class.should eq(Hash)
+
+      result[:slope].should be_within(0.00001).of(0.18784)
+
+      result[:intercept].should be_within(0.00001).of(-7.96351)
+    end
+  end
+
+  describe '.simple_linear_regression_slope' do
+    it 'should return the slope of the linear regression equation for two arrays' do
+      Rstat.simple_linear_regression_slope([60, 61, 62, 63, 65], [3.1, 3.6, 3.8, 4, 4.1]).should be_within(0.00001).of(0.18784)
+    end
+  end
+
+  describe '.simple_linear_regression_intercept' do
+    it 'should return the slope of the linear regression equation for two arrays' do
+      Rstat.simple_linear_regression_intercept([60, 61, 62, 63, 65], [3.1, 3.6, 3.8, 4, 4.1]).should be_within(0.00001).of(-7.96351)
+    end
+  end
+end

data/spec/spec_helper.rb

@@ -1,2 +1,2 @@
 # encoding: utf-8
-require "rstat"
+require 'rstat'

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: rstat
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.1.1
   prerelease: 
 platform: ruby
 authors:
@@ -9,7 +9,7 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2012-05-16 00:00:00.000000000 Z
+date: 2013-03-04 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rspec
@@ -35,6 +35,7 @@ extensions: []
 extra_rdoc_files: []
 files:
 - .gitignore
+- .rspec
 - .rvmrc
 - CHANGELOG.md
 - Gemfile
@@ -62,12 +63,17 @@ files:
 - lib/rstat/core_ext/array/sum.rb
 - lib/rstat/core_ext/fixnum.rb
 - lib/rstat/core_ext/fixnum/binomial_coefficient.rb
+- lib/rstat/regression_analysis.rb
+- lib/rstat/regression_analysis/linear_regression.rb
+- lib/rstat/regression_analysis/linear_regression/simple_linear_regression.rb
+- lib/rstat/regression_analysis/sample_correlation_coefficient.rb
 - lib/rstat/version.rb
 - rstat.gemspec
 - spec/rstat/array_spec.rb
 - spec/rstat/fixnum_spec.rb
+- spec/rstat/regession_analysis_spec.rb
 - spec/spec_helper.rb
-homepage: http://seaneshbaugh.com/
+homepage: https://github.com/seaneshbaugh/rstat
 licenses: []
 post_install_message: 
 rdoc_options: []
@@ -94,4 +100,5 @@ summary: A Simple statistics gem.
 test_files:
 - spec/rstat/array_spec.rb
 - spec/rstat/fixnum_spec.rb
+- spec/rstat/regession_analysis_spec.rb
 - spec/spec_helper.rb