具有3个输入向量和4个输出向量的线性回归?

时间:2016-05-03 07:38:59

标签: python numpy math scikit-learn linear-regression

任务:

例如,我们有3个输入向量:

foo = [1, 2, 3, 4, 5, 6]
bar = [50, 60, 70, 80, 90, 100]
spam = [-10, -20, -30, -40, -50, -60]

此外,我们有4个输出向量,它们与输入向量具有线性相关性:

foofoo = [1, 1, 2, 2, 3, 3]
barbar = [4, 4, 5, 5, 6, 6]
spamspam = [7, 7, 8, 8, 9, 9]
hamham = [10, 10, 11, 11, 12, 12]

如何在Python中使用此数据进行线性回归?

1 个答案:

答案 0 :(得分:1)

您可以使用OLS (Ordinary Least Squares model)完成here

#imports
import numpy as np
import statsmodels.api as sm

#generate the input matrix
X=[foo,bar,spam]
#turn it into a numpy array
X = np.array(X).T
#add a constant column
X=sm.add_constant(X)

这给出了输入矩阵X

array([[   1.,    1.,   50.,  -10.],
       [   1.,    2.,   60.,  -20.],
       [   1.,    3.,   70.,  -30.],
       [   1.,    4.,   80.,  -40.],
       [   1.,    5.,   90.,  -50.],
       [   1.,    6.,  100.,  -60.]])

现在您可以调整每个所需的输出向量:

resFoo = sm.OLS(endog=foofoo, exog=X).fit()
resBar = sm.OLS(endog=barbar, exog=X).fit()
resSpam = sm.OLS(endog=spamspam, exog=X).fit()
resham = sm.OLS(endog=hamham, exog=X).fit()

result为您提供系数(对于常量,以及三列foo,bar和垃圾邮件):

>>> resFoo.params
array([-0.00063323,  0.0035345 ,  0.01001583, -0.035345  ])

您现在可以使用输入检查它:

>>> np.matrix(X)*np.matrix(resFoo.params).T
matrix([[ 0.85714286],
        [ 1.31428571],
        [ 1.77142857],
        [ 2.22857143],
        [ 2.68571429],
        [ 3.14285714]])

接近foofoo的所需输出。

有关回归的不同方法,请参阅此问题:Multiple linear regression in Python