用于执行线性或非线性最小二乘近似的Ruby库?

时间:2011-05-25 00:57:39

标签: ruby linear-regression least-squares

是否有一个Ruby库允许我对一组数据进行线性或非线性最小二乘近似。

我想做的是以下内容:

  • 给出一系列[x,y]数据点
  • 针对该数据生成线性或非线性最小二乘近似
  • 库不需要弄清楚它是否需要进行线性或非线性近似。图书馆的来电者应该知道他们需要什么类型的回归

我不想尝试移植一些C / C ++ / Java库来获得这个功能,所以我希望有一些我可以使用的现有Ruby库。

3 个答案:

答案 0 :(得分:7)

尝试使用' statsample'宝石。您可以使用下面提供的示例执行对数,指数,幂或任何其他转换。我希望这有帮助。

require 'statsample'

# Independent Variable
x_data = [Math.exp(1), Math.exp(2), Math.exp(3), Math.exp(4), Math.exp(5)]

# Dependent Variable
y_data = [3, 5, 7, 9, 11]

# Logarithmic Transformation of X data 
# Math.log in Ruby has the base of Euler's number 'e' ~= '2.71828', 
# instead of the base '10'. Just a note.
log_x_data = x_data.map { |x| Math.log(x) }

# Linear Regression using the Logarithmic Transformation
x_vector=log_x_data.to_vector(:scale)
y_vector=y_data.to_vector(:scale)
ds={'x'=>x_vector,'y'=>y_vector}.to_dataset
mlr=Statsample::Regression.multiple(ds,'y')
mlr.summary

# Provides the value of the y-intercept 
#p mlr.constant

# Lists the coefficients of each casual variable. In this case, we have only one--'x'.
#p mlr.coeffs

# The regression output produces the line y = 1 + 2*x, but 
# considering that we transformed x earlier, it really produces
# y = 1 + 2*ln(x).

答案 1 :(得分:5)

我使用this snippet来计算一些回归。第一个参数是包含x坐标的数组,第二个参数是包含y坐标的数组,最后一个是您要查找的多项式的次数。不确定这是否是您正在寻找的,但希望它有所帮助。

答案 2 :(得分:0)

我正在维护一个C库,用于非线性最小二乘最小化http://apps.jcns.fz-juelich.de/lmfit,它附带Ruby的swig文件。