当前位置: 首页 > news >正文

机器学习线性归回实战(单因子和多音字分别建立预测房价模型)

以面积为输入变量,建立单因子模型,评估模型表现,预测合理房价

import pandas as pd
import numpy as npdata = pd.read_csv('usa_housing_price.csv')
data.head()

在这里插入图片描述

from matplotlib import pyplot as pltfig = plt.figure(figsize=(10,10)) #图的尺寸
fig1 = plt.subplot(231) #2行3列第一个图
plt.scatter(data.loc[:,'Avg.Area Income'],data.loc[:,'Price'])
plt.title('Price VS Income')fig2 = plt.subplot(232) 
plt.scatter(data.loc[:,'Avg.Area House Age'],data.loc[:,'Price'])
plt.title('Price VS House Age')fig3 = plt.subplot(233) 
plt.scatter(data.loc[:,'Avg.Area Number of Rooms'],data.loc[:,'Price'])
plt.title('Price VS Number of Rooms')fig4 = plt.subplot(234) 
plt.scatter(data.loc[:,'Area Population'],data.loc[:,'Price'])
plt.title('Price VS Area Population')fig5 = plt.subplot(235) 
plt.scatter(data.loc[:,'size'],data.loc[:,'Price'])
plt.title('Price VS size')plt.show()

在这里插入图片描述

# 定义X,Y
X = data.loc(:,'size') # 面积
Y = data.loc(:,'Price') # 价格X = np.array(X).reshape(-1,1) #维度转换 print(X.shape)import sklearn.linear_model import LineearRegression
LR1 = LinearRegression()
LR1.fit(X,y) #模型训练y_predict_1 = LR1.predict(X)
print(y_predict_1) #计算预测的价格和面积from sklean.metrics import mean_squared_error,r2_score # 评估模型
mean_squared_error_1 = mean_squared_error(y,y_predict_1)
r2_score_1 = r2_score(y,y_predict_1) 
print(mean_squared_error_1,r2_score_1) #r2_score_1越接近1越好

在这里插入图片描述

fig6 = plt.figure(figsize=(8,5))
plt.scatter(X,y)
plt.plot(X,y_predict_1,'r') #红色方式可视化

在这里插入图片描述

多因子模型

#重新定义X
X_multi = data.drop(['Price'],axis=1)
x_Multi 

在这里插入图片描述

# 建立线性回归的模型
LR_multi = LinearRegression()
LR_multi.fit(X_multi,y)# 模型预测
y_predict_multi = LR_multi.predict(X_multi)
print(y_predict_multi)mean_squared_error_multi = mean_squared_error(y,y_predict_multi)
r2_score_multi = r2_score(y,y_predict_multi)
print(mean_squared_error_multi,r2_score_multi)

在这里插入图片描述
在这里插入图片描述

# 可视化
fig7 = plt.figure(figsize=(8,5))
plt.scatter(y,y_predict_multi)
plt.show()

在这里插入图片描述

预测Income=65000,House Age=5,Number of Rooms=5,Population=30000,size=200合理房价

X_test = [65000,5,5,30000,200]

X_test =  [65000,5,5,30000,200]
X_test = np.array(X_test).reshape(1,-1) #转换成数组,维度转化为一行若干列y_test_predict = LR_multi.predict(X_test) #调用模型进行预测
print(y_test_predict)
http://www.lryc.cn/news/615919.html

相关文章:

  • 一个基于 Next.js 和 Puppeteer 的 Markdown 转图片服务,支持 Docker 部署和 API 集成
  • Node.js面试题及详细答案120题(01-15) -- 基础概念篇
  • python | numpy小记(十):理解 NumPy 中的 `np.random.multinomial`(进阶)
  • Stlink识别不到-安装驱动
  • 医防融合中心-智慧化慢病全程管理医疗AI系统开发(下)
  • 整数规划-分支定界
  • Docker Compose 部署高可用 MongoDB 副本集集群(含 Keepalived + HAProxy 负载均衡)
  • AI编程插件对比分析:CodeRider、GitHub Copilot及其他
  • 给AI装上“翻译聚光灯”:注意力机制的机器翻译革命
  • 【精彩回顾·成都】成都 User Group×柴火创客空间:开源硬件驱动 AI 与云的创新实践!
  • 打卡day34
  • openpnp - 顶部相机如果超过6.5米影响通讯质量,可以加USB3.0信号放大器延长线
  • Spark执行计划与UI分析
  • AutoCAD 2026 的主要功能
  • 变量详解:创建初始化与内存管理
  • lesson34:深入理解Python线程:从基础到实战优化
  • XGBoost算法在机器学习中的实现
  • Android Camera 打开和拍照APK源码
  • Android 开发问题:Invalid id; ID definitions must be of the form @+id/ name
  • Android 16 KB页面大小适配的权威技术方案总结
  • Ubuntu 安装 Kibana
  • 神经机器翻译(NMT)框架:编码器-解码器(Encoder-Decoder)结构详解
  • 支持selenium的chrome driver更新到139.0.7258.66
  • 去除Edge微软浏览器与Chrome谷歌浏览器顶部出现“此版本的Windows不再支持升级Windows 10”的烦人提示
  • Elasticsearch QueryDSL 教程
  • Linux操作系统从入门到实战(十八)在Linux里面怎么查看进程
  • 三、k8s 1.29 之 安装1网络 / ikuai路由器虚拟机安装
  • Linux810 shell 条件判断 文件工具 ifelse
  • JavaWeb(苍穹外卖)--学习笔记18(Apache POI)
  • 【QT】常⽤控件详解(七)容器类控件 GroupBox TabWidget 布局管理器 Spacer