method

fit_exponential

v1_9_3_392 - Show latest stable - 0 notes - Class: TestCase

1_8_6_287
1_8_7_72
1_8_7_330
1_9_1_378
1_9_2_180
1_9_3_125 (0)
1_9_3_392 (0)
2_1_10 (0)
2_2_9
2_4_6
2_5_5
2_6_3
What's this?

fit_exponential(xs, ys) public

To fit a functional form: y = ae^(bx).

Takes x and y values and returns [a, b, r^2].

See: http://mathworld.wolfram.com/LeastSquaresFittingExponential.html

Show source

# File lib/minitest/benchmark.rb, line 225
    def fit_exponential xs, ys
      n     = xs.size
      xys   = xs.zip(ys)
      sxlny = sigma(xys) { |x,y| x * Math.log(y) }
      slny  = sigma(xys) { |x,y| Math.log(y)     }
      sx2   = sigma(xys) { |x,y| x * x           }
      sx    = sigma xs

      c = n * sx2 - sx ** 2
      a = (slny * sx2 - sx * sxlny) / c
      b = ( n * sxlny - sx * slny ) / c

      return Math.exp(a), b, fit_error(xys) { |x| Math.exp(a + b * x) }
    end

fit_exponential

Related methods