sklearn调包侠之线性回归

  1. 云栖社区>
  2. 博客>
  3. 正文

sklearn调包侠之线性回归

是罗罗攀啊 2018-06-28 11:46:00 浏览405
展开阅读全文
img_745f3c3f47200dd2ea7987ae88e25fec.png

线性回归原理

如图所示,这是一组二维的数据,我们先想想如何通过一条直线较好的拟合这些散点了?直白的说:尽量让拟合的直线穿过这些散点(这些点离拟合直线很近)。

img_f9d60858152aef8453ffcfdc4dc3992e.png
目标函数(成本函数)

要使这些点离拟合直线很近,我们需要用数学公式来表示:

img_25d94181bc646e8869b5dc8fd2e7c20b.png
梯度下降法

之前在讲解回归时,是通过求导获取最小值,但必须满足数据可逆,这里通常情况下使用梯度下降法,也就是按着斜率方向偏移。详细可看这篇文章(https://www.jianshu.com/p/96566542b07a)。
tips:这篇文章讲解了梯度上升法,梯度下降法类似。

img_934720a23364477c7220565869d167ac.png

实战——房价预测

数据导入

该数据使用sklearn自带的数据集,通过sklearn.datasets导入我们的boston房价数据集。

from sklearn.datasets import load_boston
boston = load_boston()

通过DESCR属性可以查看数据集的详细情况,这里数据有14列,前13列为特征数据,最后一列为标签数据。

print(boston.DESCR)
img_70621c004595a63115e30471f7dfa485.png

boston的data和target分别存储了特征和标签:

img_937d430fd46f23dc29022b1cfbd77bcd.png
切分数据集
from sklearn.model_selection import train_test_split
X_train, X_test, y_train, y_test = train_test_split(boston.data, boston.target, test_size = 0.2, random_state=2)
数据预处理

普通的线性回归模型太简单,容易导致欠拟合,我们可以增加特征多项式来让线性回归模型更好地拟合数据。在sklearn中,通过preprocessing模块中的PolynomialFeatures来增加特征多项式。
其重要参数有:

  • degree:多项式特征的个数,默认为2
  • include_bias:默认为True,包含一个偏置列,也就是 用作线性模型中的截距项,这里选择False,因为在线性回归中,可以设置是否需要截距项。
from sklearn.preprocessing import PolynomialFeatures
poly = PolynomialFeatures(degree=2,include_bias=False)
X_train_poly = poly.fit_transform(X_train)
X_test_poly = poly.fit_transform(X_test)
模型训练与评估

线性算法使用sklearn.linear_model 模块中的LinearRegression方法。常用的参数如下:

  • fit_intercept:默认为True,是否计算截距项。
  • normalize:默认为False,是否对数据归一化。

简单线性回归

from sklearn.linear_model import LinearRegression

model2 = LinearRegression(normalize=True)
model2.fit(X_train, y_train)
model2.score(X_test, y_test)

# result
# 0.77872098747725804

多项式线性回归

model3 = LinearRegression(normalize=True)
model3.fit(X_train_poly, y_train)
model3.score(X_test_poly, y_test)

# result
# 0.895848854203947
总结

多项式的个数的不断增加,可以在训练集上有很好的效果,但缺很容易造成过拟合,没法在测试集上有很好的效果,也就是常说的:模型泛化能力差。

网友评论

登录后评论
0/500
评论
是罗罗攀啊
+ 关注