from sklearn.svm import LinearSVR import matplotlib.pyplot as plt data = pd.read_excel('../tmp/new_reg_data_GM11.xls') # 读取数据 data = data.set_index('Unnamed: 0') data = data.drop(index ='模型精度') feature = ['x1', 'x3', 'x4', 'x5', 'x6', 'x7', 'x8', 'x13'] #特征所在列 data_train = data.loc[range(1994, 2014)].copy() #取2014年前的数据建模 data_mean = data_train.mean() data_std = data_train.std() data_train = (data_train - data_mean) / data_std #数据标准化 x_train = data_train[feature].as_matrix() #特征数据 y_train = data_train['y'].as_matrix() #标签数据 linearsvr = LinearSVR(random_state=123) #调用LinearSVR()函数 linearsvr.fit(x_train, y_train) #预测2014 年和2015 年的财政收入,并还原结果 x = ((data[feature] - data_mean[feature]) / data std[feature]).as_matrix() data[u'y_pred'] = linearsvr.predict (x) * data_std['y'] + data_mean['y'] outputfile ='../tmp/new_reg_data_GM11_revenue.xls' data.to_excel(outputfile) print('真实值与预测值分别为: \n', data[['y', 'y_pred']]) print('预测图为: ',data[['y','y_pred']].plot(style = ['b-o','r-*'])) #画出预测结果图 plt.xlabel('年份') plt.xticks(range(1994,2015,2))

时间: 2024-02-01 14:02:07 浏览: 67
ZIP

irisdata.zip_MATLAB导入 iris_iris_iris.data.mat_iris_data.mat_iris

star5星 · 资源好评率100%
这段代码是使用线性支持向量回归(LinearSVR)进行财政收入预测的示例。首先,代码导入了所需的库,包括sklearn.svm中的LinearSVR和matplotlib.pyplot。然后,通过pd.read_excel方法读取了名为'../tmp/new_reg_data_GM11.xls'的Excel文件,并将数据设置为以'Unnamed: 0'列为索引。接下来,选择了特定的特征列,并将数据划分为训练集和测试集。然后,对训练集进行了数据标准化操作,并将特征数据和标签数据分别保存在x_train和y_train中。接着,创建了一个LinearSVR对象,并使用训练数据进行拟合。然后,对2014年和2015年的财政收入进行预测,并将结果还原至原始数据范围内。最后,将预测结果保存至Excel文件'../tmp/new_reg_data_GM11_revenue.xls',并打印出真实值和预测值。同时,代码还画出了真实值和预测值的图表,并设置了横坐标的刻度。
阅读全文

相关推荐

import numpy as np import pandas as pd from sklearn.model_selection import train_test_split, GridSearchCV from sklearn.metrics import accuracy_score, precision_score, recall_score, f1_score from sklearn.metrics import confusion_matrix import matplotlib.pyplot as plt from termcolor import colored as cl import itertools from sklearn.preprocessing import StandardScaler from sklearn.tree import DecisionTreeClassifier from sklearn.neighbors import KNeighborsClassifier from sklearn.linear_model import LogisticRegression from sklearn.svm import SVC from sklearn.ensemble import RandomForestClassifier from xgboost import XGBClassifier from sklearn.neural_network import MLPClassifier from sklearn.ensemble import VotingClassifier # 定义模型评估函数 def evaluate_model(y_true, y_pred): accuracy = accuracy_score(y_true, y_pred) precision = precision_score(y_true, y_pred, pos_label='Good') recall = recall_score(y_true, y_pred, pos_label='Good') f1 = f1_score(y_true, y_pred, pos_label='Good') print("准确率:", accuracy) print("精确率:", precision) print("召回率:", recall) print("F1 分数:", f1) # 读取数据集 data = pd.read_csv('F:\数据\大学\专业课\模式识别\大作业\数据集1\data clean Terklasifikasi baru 22 juli 2015 all.csv', skiprows=16, header=None) # 检查数据集 print(data.head()) # 划分特征向量和标签 X = data.iloc[:, :-1] y = data.iloc[:, -1] # 划分训练集和测试集 X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=42) # 6. XGBoost xgb = XGBClassifier(max_depth=4) y_test = np.array(y_test, dtype=int) xgb.fit(X_train, y_train) xgb_pred = xgb.predict(X_test) print("\nXGBoost评估结果:") evaluate_model(y_test, xgb_pred)

# 导入必要的库 import pandas as pd from sklearn.model_selection import train_test_split from sklearn.svm import SVC from sklearn.metrics import accuracy_score, confusion_matrix import seaborn as sns import matplotlib.pyplot as plt # 读入excel表格 data = pd.read_excel('E:/桌面/预测脆弱性/20230523/预测样本/预测样本.xlsx') # 定义特征变量和因变量 features = ['高程', '起伏度', '桥梁长', '道路长', '平均坡度', '平均地温', 'T小于0', '相态'] target = '交通风险' # 将特征变量和因变量分离出来,并划分训练集和验证集 x_train, x_test, y_train, y_test = train_test_split(data[features], data[target], test_size=0.2, random_state=42) # 训练支持向量机模型 svm_model = SVC(kernel='linear', C=1, gamma=1) svm_model.fit(x_train, y_train) # 计算模型精度 y_pred = svm_model.predict(x_test) accuracy = accuracy_score(y_test, y_pred) print('模型精度:', accuracy) # 生成混淆矩阵图片并保存 cm = confusion_matrix(y_test, y_pred) plt.figure(figsize=(8, 6)) sns.heatmap(cm, annot=True, cmap='Blues') plt.xlabel('Predicted Label') plt.xlabel('Predicted') plt.ylabel('Actual') plt.savefig('E:/桌面/预测脆弱性/20230523/预测样本/预测结果/交通风险预测096/支持向量机confusion_matrix.png') # 预测新的数据 new_data = pd.read_excel('E:/桌面/预测脆弱性/20230523/预测样本/预测结果/交通风险预测096.xlsx') new_x = new_data[features] new_y = svm_model.predict(new_x) new_data[target] = new_y # 输出新的excel表格 new_data.to_excel('E:/桌面/预测脆弱性/20230523/预测样本/预测结果/交通风险预测096-支持向量机结果.xlsx', index=False)运行时间很长

最新推荐

recommend-type

MATLAB-四连杆机构的仿真+项目源码+文档说明

<项目介绍> - 四连杆机构的仿真 --m3_1.m: 位置问题求解 --m2_1.m: 速度问题求解 --FourLinkSim.slx: Simlink基于加速度方程的仿真 --FourLinkSim2.slx: Simscape简化模型仿真 --FourLinkSim3.slx: Simscape CAD模型仿真 - 不懂运行,下载完可以私聊问,可远程教学 1、该资源内项目代码都经过测试运行成功,功能ok的情况下才上传的,请放心下载使用! 2、本项目适合计算机相关专业(如计科、人工智能、通信工程、自动化、电子信息等)的在校学生、老师或者企业员工下载学习,也适合小白学习进阶,当然也可作为毕设项目、课程设计、作业、项目初期立项演示等。 3、如果基础还行,也可在此代码基础上进行修改,以实现其他功能,也可用于毕设、课设、作业等。 下载后请首先打开README.md文件(如有),仅供学习参考, 切勿用于商业用途。 --------
recommend-type

ridge_regression:用于岭回归的python代码(已实现以预测下个月的CO2浓度)

ridge_regression 用于岭回归的python代码(已实现以预测下个月的CO2浓度) 资料可用性 文件 Ridge.py :标准函数和Ridge回归函数window_make.py :使用滑动窗口方法制作大小为p(窗口大小)的时间序列列表。 Final_version.ipynb :使用Co2数据对代码进行实验
recommend-type

Polygon3-3.0.8-cp35-cp35m-win_amd64.whl.rar

python whl离线安装包 pip安装失败可以尝试使用whl离线安装包安装 第一步 下载whl文件,注意需要与python版本配套 python版本号、32位64位、arm或amd64均有区别 第二步 使用pip install XXXXX.whl 命令安装,如果whl路径不在cmd窗口当前目录下,需要带上路径 WHL文件是以Wheel格式保存的Python安装包, Wheel是Python发行版的标准内置包格式。 在本质上是一个压缩包,WHL文件中包含了Python安装的py文件和元数据,以及经过编译的pyd文件, 这样就使得它可以在不具备编译环境的条件下,安装适合自己python版本的库文件。 如果要查看WHL文件的内容,可以把.whl后缀名改成.zip,使用解压软件(如WinRAR、WinZIP)解压打开即可查看。 为什么会用到whl文件来安装python库文件呢? 在python的使用过程中,我们免不了要经常通过pip来安装自己所需要的包, 大部分的包基本都能正常安装,但是总会遇到有那么一些包因为各种各样的问题导致安装不了的。 这时我们就可以通过尝试去Python安装包大全中(whl包下载)下载whl包来安装解决问题。
recommend-type

SSM Java项目:StudentInfo 数据管理与可视化分析

资源摘要信息:"StudentInfo 2.zip文件是一个压缩包,包含了多种数据可视化和数据分析相关的文件和代码。根据描述,此压缩包中包含了实现人员信息管理系统的增删改查功能,以及生成饼图、柱状图、热词云图和进行Python情感分析的代码或脚本。项目使用了SSM框架,SSM是Spring、SpringMVC和MyBatis三个框架整合的简称,主要应用于Java语言开发的Web应用程序中。 ### 人员增删改查 人员增删改查是数据库操作中的基本功能,通常对应于CRUD(Create, Retrieve, Update, Delete)操作。具体到本项目中,这意味着实现了以下功能: - 增加(Create):可以向数据库中添加新的人员信息记录。 - 查询(Retrieve):可以检索数据库中的人员信息,可能包括基本的查找和复杂的条件搜索。 - 更新(Update):可以修改已存在的人员信息。 - 删除(Delete):可以从数据库中移除特定的人员信息。 实现这些功能通常需要编写相应的后端代码,比如使用Java语言编写服务接口,然后通过SSM框架与数据库进行交互。 ### 数据可视化 数据可视化部分包括了生成饼图、柱状图和热词云图的功能。这些图形工具可以直观地展示数据信息,帮助用户更好地理解和分析数据。具体来说: - 饼图:用于展示分类数据的比例关系,可以清晰地显示每类数据占总体数据的比例大小。 - 柱状图:用于比较不同类别的数值大小,适合用来展示时间序列数据或者不同组别之间的对比。 - 热词云图:通常用于文本数据中,通过字体大小表示关键词出现的频率,用以直观地展示文本中频繁出现的词汇。 这些图表的生成可能涉及到前端技术,如JavaScript图表库(例如ECharts、Highcharts等)配合后端数据处理实现。 ### Python情感分析 情感分析是自然语言处理(NLP)的一个重要应用,主要目的是判断文本的情感倾向,如正面、负面或中立。在这个项目中,Python情感分析可能涉及到以下几个步骤: - 文本数据的获取和预处理。 - 应用机器学习模型或深度学习模型对预处理后的文本进行分类。 - 输出情感分析的结果。 Python是实现情感分析的常用语言,因为有诸如NLTK、TextBlob、scikit-learn和TensorFlow等成熟的库和框架支持相关算法的实现。 ### IJ项目与readme文档 "IJ项目"可能是指IntelliJ IDEA项目,IntelliJ IDEA是Java开发者广泛使用的集成开发环境(IDE),支持SSM框架。readme文档通常包含项目的安装指南、运行步骤、功能描述、开发团队和联系方式等信息,是项目入门和理解项目结构的首要参考。 ### 总结 "StudentInfo 2.zip"是一个综合性的项目,涉及到后端开发、前端展示、数据分析及自然语言处理等多个技术领域。通过这个项目,可以学习到如何使用SSM框架进行Web应用开发、实现数据可视化和进行基于Python的情感分析。这对于想要掌握Java Web开发和数据处理能力的学习者来说是一个很好的实践机会。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

负载均衡技术深入解析:确保高可用性的网络服务策略

![负载均衡技术深入解析:确保高可用性的网络服务策略](https://media.geeksforgeeks.org/wp-content/uploads/20240130183502/Source-IP-hash--(1).webp) # 1. 负载均衡技术概述 ## 1.1 负载均衡技术的重要性 在现代信息技术不断发展的今天,互联网应用的规模和服务的复杂性日益增长。因此,为了确保高性能、高可用性和扩展性,负载均衡技术变得至关重要。它能够有效地分配和管理网络或应用程序的流量,使得服务器和网络资源得以最优利用。 ## 1.2 负载均衡技术的基本概念 负载均衡是一种网络流量管理技术,旨
recommend-type

怎么解决头文件重复包含

解决头文件重复包含的问题主要有以下几个策略: 1. **包含 guards**:在头文件开头添加一种特殊的标识符(通常是宏),如 `#ifndef` 和 `#define` 对组合,检查某个特定宏是否已经定义过。如果没有定义,则包含内容,然后设置该宏。如果在同一文件内再次包含,由于宏已经存在,就不会再执行包含的内容,从而避免重复。 ```cpp #ifndef HEADER_NAME_H_ #define HEADER_NAME_H_ // 内容... #endif // HEADER_NAME_H_ ``` 2. **使用 extern 关键字**:对于非静态变量和函数,可以将它们
recommend-type

pyedgar:Python库简化EDGAR数据交互与文档下载

资源摘要信息:"pyedgar:用于与EDGAR交互的Python库" 知识点说明: 1. pyedgar库概述: pyedgar是一个Python编程语言下的开源库,专门用于与美国证券交易委员会(SEC)的电子数据获取、访问和检索(EDGAR)系统进行交互。通过该库,用户可以方便地下载和处理EDGAR系统中公开提供的财务报告和公司文件。 2. EDGAR系统介绍: EDGAR系统是一个自动化系统,它收集、处理、验证和发布美国证券交易委员会(SEC)要求的公司和其他机构提交的各种文件。EDGAR数据库包含了美国上市公司的详细财务报告,包括季度和年度报告、委托声明和其他相关文件。 3. pyedgar库的主要功能: 该库通过提供两个主要接口:文件(.py)和索引,实现了对EDGAR数据的基本操作。文件接口允许用户通过特定的标识符来下载和交互EDGAR表单。索引接口可能提供了对EDGAR数据库索引的访问,以便快速定位和获取数据。 4. pyedgar库的使用示例: 在描述中给出了一个简单的使用pyedgar库的例子,展示了如何通过Filing类与EDGAR表单进行交互。首先需要从pyedgar模块中导入Filing类,然后创建一个Filing实例,其中第一个参数(20)可能代表了提交年份的最后两位,第二个参数是一个特定的提交号码。创建实例后,可以打印实例来查看EDGAR接口的返回对象,通过打印实例的属性如'type',可以获取文件的具体类型(例如10-K),这代表了公司提交的年度报告。 5. Python语言的应用: pyedgar库的开发和应用表明了Python语言在数据分析、数据获取和自动化处理方面的强大能力。Python的简洁语法和丰富的第三方库使得开发者能够快速构建工具以处理复杂的数据任务。 6. 压缩包子文件信息: 文件名称列表中的“pyedgar-master”表明该库可能以压缩包的形式提供源代码和相关文件。文件列表中的“master”通常指代主分支或主版本,在软件开发中,主分支通常包含了最新的代码和功能。 7. 编程实践建议: 在使用pyedgar库之前,建议先阅读官方文档,了解其详细的安装、配置和使用指南。此外,进行编程实践时,应当注意遵守SEC的使用条款,确保只下载和使用公开提供的数据。 8. EDGAR数据的应用场景: EDGAR数据广泛应用于金融分析、市场研究、合规性检查、学术研究等领域。通过编程访问EDGAR数据可以让用户快速获取到一手的财务和公司运营信息,从而做出更加明智的决策。 9. Python库的维护和更新: 随着EDGAR数据库内容的持续更新和变化,pyedgar库也应定期进行维护和更新,以保证与EDGAR系统的接口兼容性。开发者社区对于这类开源项目的支持和贡献也非常重要。 10. 注意事项: 在使用pyedgar库下载和处理数据时,用户应当确保遵守相应的法律法规,尤其是关于数据版权和隐私方面的规定。此外,用户在处理敏感数据时,还需要考虑数据安全和隐私保护的问题。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

网络监控工具使用宝典:实时追踪网络状况的专家级技巧

![网络监控工具使用宝典:实时追踪网络状况的专家级技巧](https://docs.itrsgroup.com/docs/geneos/5.8.0/Resources/Images/netprobe/netprobe-api-images/image2_xml-rpc_api_-_user_guide_and_technical_reference.png) # 1. 网络监控工具概述 网络监控工具是确保网络稳定性和安全性不可或缺的组成部分。在这一章节中,我们将概述这些工具的重要性、工作原理以及它们如何帮助IT专业人员提高网络性能和安全性。 ## 1.1 网络监控工具的定义 网络监控工具