第三篇 2.1-线性回归原理,推导,衡量线性回归的指标(MSE,RMSE,MAE,R^2)

本文为学习记录,写的不好请大家指正见谅,互相学习,也为自己巩固学习成果包括后续复习所用!

1.简单的线性回归

我们先来看图

线性回归就是找到一条直线,y=ax+b

加入特征值为x,即预测值为y’,实际值为y,线性回归就是找到|y-y'|的最小值,即(y-y')^2的最小值,扩展到所有样本就是找到

∑(y-y')^2的最小值 

也就是找到a和b的值,使得上述算式最小,即对a,b求导

所以得到的b=y(mean)-a*x(mean)

a=∑ (x(i)-x(mean)) * (y(i)-y(mean))  / ∑((x(i) - x(mean)) ^2

使用代码实现fit过程 即

MSE,RMSE,MAE,R^2

回到上面的线性回归,线性回归主要是找到∑(y(train)-y(train)')^2的最小值 ,这里的y指的是训练集,所以在做衡量线性回归指标的时候,

我们只要找到y测试集的∑(y(test)-y(test)')^2即可,

在这里,我们只要消除数据量m的影响即可,所以形成了

MSE(均方误差):

再根据均方误差,我们消除平方带来的特征值放大的影响,

所以有RMSE(均方根误差):

还有我们干脆直接取距离的绝对值,即:

 MAE(平均绝对误差):

封装代码:

最好的方法为R^2:

分子部分指的是我们自己的模型产生的错误,分母表示使用baseline model(基础模型)y=y(mean) 产生的错误,

两个相除再用1减去之后相当于就是我们拟合住的部分,所以对于R^2来说,越大越好

1.R^2<=1,如果等于1,表示模型完全拟合,没有任何错误。

2.当R^2=0时,表示模型和基础模型一样

3.小于0时,可能我们的模型没有线性关系了已经,可能需要更换算法了

R^2也可以写成

1-  MSE(y',y) / Var(y)  相当于分子分母同时除以一个m,这样我们在实现的时候就很简单了

 

现在我们扩展到多元的线性回归

也就是我们每一个特征X之前,都有一个theta值,相当于简单线性回归的扩展

具体的写在纸上

使用程序实现linearRegression

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,下面是实现的代码。 首先,我们需要导入所需的库和数据集: ```python import numpy as np from sklearn.datasets import load_boston from sklearn.model_selection import train_test_split from sklearn.linear_model import LinearRegression from sklearn.metrics import mean_squared_error, mean_absolute_error, r2_score # 加载波士顿房价数据集 boston = load_boston() X = boston.data y = boston.target # 将数据集划分为训练集和测试集 X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=42) ``` 然后,我们使用线性回归模型对训练数据进行拟合,并预测测试数据的房价: ```python # 使用线性回归模型对训练数据进行拟合 lr = LinearRegression() lr.fit(X_train, y_train) # 预测测试数据的房价 y_pred = lr.predict(X_test) ``` 最后,我们使用三个评估指标MSEMAER2)来评估模型的性能: ```python # 计算MSE(均方误差) mse = mean_squared_error(y_test, y_pred) print("MSE: ", mse) # 计算MAE(平均绝对误差) mae = mean_absolute_error(y_test, y_pred) print("MAE: ", mae) # 计算R2(决定系数) r2 = r2_score(y_test, y_pred) print("R2: ", r2) ``` 完整代码如下: ```python import numpy as np from sklearn.datasets import load_boston from sklearn.model_selection import train_test_split from sklearn.linear_model import LinearRegression from sklearn.metrics import mean_squared_error, mean_absolute_error, r2_score # 加载波士顿房价数据集 boston = load_boston() X = boston.data y = boston.target # 将数据集划分为训练集和测试集 X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=42) # 使用线性回归模型对训练数据进行拟合 lr = LinearRegression() lr.fit(X_train, y_train) # 预测测试数据的房价 y_pred = lr.predict(X_test) # 计算MSE(均方误差) mse = mean_squared_error(y_test, y_pred) print("MSE: ", mse) # 计算MAE(平均绝对误差) mae = mean_absolute_error(y_test, y_pred) print("MAE: ", mae) # 计算R2(决定系数) r2 = r2_score(y_test, y_pred) print("R2: ", r2) ``` 输出结果如下: ``` MSE: 24.29111947497371 MAE: 3.189091965887853 R2: 0.6687594935356313 ``` 这意味着我们的模型的平均预测误差约为3.19,R2评分为0.67,这表明模型对数据的解释力还可以,但仍有改进的空间。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值