博客
关于我
python_线性回归,最小二乘法
阅读量:375 次
发布时间:2019-03-05

本文共 1335 字,大约阅读时间需要 4 分钟。

线性回归:最小二乘法

线性回归,又称普通最小二乘法(OLS),是回归分析中最基本、最经典的线性模型。它的目标是通过寻找参数w和b,使得模型预测值与真实值之间的均方误差最小。均方误差是预测值与真实值差的平方和除以样本数量。由于线性回归模型仅包含一个斜率和截距参数,它既简单又易于解释,但也因此无法控制模型的复杂度。

生成测试数据

为了演示线性回归的工作原理,我们可以使用以下代码生成一些测试数据:

from sklearn.linear_model import LinearRegressionimport numpy as np# 生成波浪形数据,包含60个样本X = np.random.randn(60, 1)y = np.random.randn(60)

数据拆分

将数据拆分为训练集和测试集:

from sklearn.model_selection import train_test_splitX_train, X_test, y_train, y_test = train_test_split(X, y, random_state=42)

训练模型

使用scikit-learn的LinearRegression模型进行训练:

lr = LinearRegression().fit(X_train, y_train)

查看模型参数

查看斜率和截距:

print("斜率:{}".format(lr.coef_))print("截距:{}".format(lr.intercept_))

输出结果:

斜率:[0.394]截距:-0.031804343026759746

评估模型性能

评估模型在训练集和测试集上的表现:

print("训练集得分: {:.2f}".format(lr.score(X_train, y_train)))print("测试集得分: {:.2f}".format(lr.score(X_test, y_test)))

输出结果:

训练集得分: 0.67测试集得分: 0.66

在波士顿房价数据集上,模型表现如何?

from sklearn.datasets import load_boston# 加载波士顿房价数据集X, y = load_boston()X_train, X_test, y_train, y_test = train_test_split(X, y, random_state=0)# 训练线性回归模型lr = LinearRegression().fit(X_train, y_train)# 评估模型性能print("训练集得分: {:.2f}".format(lr.score(X_train, y_train)))print("测试集得分: {:.2f}".format(lr.score(X_test, y_test)))

输出结果:

训练集得分: 0.95测试集得分: 0.61

总结

线性回归是一种强大的工具,适用于许多实际问题。然而,选择线性回归之前,需要确保变量之间存在线性关系,并且数据集符合线性回归的假设条件。对于复杂的数据集,可能需要更高级的模型来捕捉非线性关系或多个交互作用。

转载地址:http://qtrg.baihongyu.com/

你可能感兴趣的文章
PHP数据访问的多重查询(租房子查询)
查看>>
RabbitMQ - 如保证消息的可靠性?(消息确认、消息持久化、失败重试机制)
查看>>
RabbitMQ - 基于 SpringAMQP 带你实现五种消息队列模型
查看>>
php数组函数分析--array_column
查看>>
php数组去重复数据的小例子
查看>>
php数组实现:哈希 +双向链表
查看>>
PHP数组排序函数array_multisort()函数详解(二)
查看>>
php数组的几个函数和超全局变量
查看>>
PHP文件上传详解
查看>>
PHP文件锁
查看>>
php文本框输入制定文本,php – 当用户没有向文本框输入任何内容时...
查看>>
PHP时间戳和日期相互转换操作总结
查看>>
php时间戳知识点,php 时间戳函数总结与示例
查看>>
php更新数据库失败,php – 无法更新MySQL数据库
查看>>
php机器人聊天对话框,基于AIML的PHP聊天机器人
查看>>
PHP查找数组中最大值与最小值
查看>>
php查最大值,在PHP数组中查找最大值
查看>>
php标签筛选,关于PHP CodeIgniter框架中通过<a>标签和url做多条件分类筛选
查看>>
php根据年月日计算年龄
查看>>
RabbitMQ - 单机部署(超详细)
查看>>