用jupyter写一个用BP预测多维股票价格的代码,其中评价指标要有mape,mae,rmse,r2,并画出拟合效果图。

时间: 2024-04-10 16:28:14 浏览: 19
以下是使用BP神经网络进行多维股票价格预测的代码,并计算了MAPE、MAE、RMSE和R2这四个评价指标,并绘制了拟合效果图。请注意,这只是一个示例代码,实际使用中可能需要根据具体情况进行调整。 ```python import numpy as np import pandas as pd from sklearn.preprocessing import MinMaxScaler from sklearn.model_selection import train_test_split from sklearn.metrics import mean_absolute_percentage_error, mean_absolute_error, mean_squared_error, r2_score import matplotlib.pyplot as plt import tensorflow as tf from tensorflow.keras.models import Sequential from tensorflow.keras.layers import Dense # 加载数据 data = pd.read_csv('stock_prices.csv') # 假设数据保存在stock_prices.csv文件中 # 数据预处理 scaler = MinMaxScaler() scaled_data = scaler.fit_transform(data) # 划分训练集和测试集 train_data, test_data = train_test_split(scaled_data, test_size=0.2, shuffle=False) # 准备训练集和测试集的输入和输出数据 train_X = train_data[:, :-1] train_y = train_data[:, -1] test_X = test_data[:, :-1] test_y = test_data[:, -1] # 构建BP神经网络模型 model = Sequential() model.add(Dense(64, activation='relu', input_dim=train_X.shape[1])) model.add(Dense(32, activation='relu')) model.add(Dense(1, activation='linear')) # 编译模型 model.compile(optimizer='adam', loss='mean_squared_error') # 训练模型 model.fit(train_X, train_y, epochs=100, batch_size=32, verbose=0) # 使用模型进行预测 train_pred = model.predict(train_X) test_pred = model.predict(test_X) # 反归一化 train_pred = scaler.inverse_transform(train_pred) train_y = scaler.inverse_transform([train_y]) test_pred = scaler.inverse_transform(test_pred) test_y = scaler.inverse_transform([test_y]) # 计算评价指标 mape_train = mean_absolute_percentage_error(train_y[0], train_pred[:, 0]) mape_test = mean_absolute_percentage_error(test_y[0], test_pred[:, 0]) mae_train = mean_absolute_error(train_y[0], train_pred[:, 0]) mae_test = mean_absolute_error(test_y[0], test_pred[:, 0]) rmse_train = np.sqrt(mean_squared_error(train_y[0], train_pred[:, 0])) rmse_test = np.sqrt(mean_squared_error(test_y[0], test_pred[:, 0])) r2_train = r2_score(train_y[0], train_pred[:, 0]) r2_test = r2_score(test_y[0], test_pred[:, 0]) print("Train MAPE: %.4f" % mape_train) print("Test MAPE: %.4f" % mape_test) print("Train MAE: %.4f" % mae_train) print("Test MAE: %.4f" % mae_test) print("Train RMSE: %.4f" % rmse_train) print("Test RMSE: %.4f" % rmse_test) print("Train R2: %.4f" % r2_train) print("Test R2: %.4f" % r2_test) # 绘制拟合效果图 plt.figure(figsize=(10, 6)) plt.plot(train_y[0], label="Actual train prices") plt.plot(train_pred[:, 0], label="Predicted train prices") plt.plot(len(train_y[0]) + np.arange(len(test_y[0])), test_y[0], label="Actual test prices") plt.plot(len(train_y[0]) + np.arange(len(test_y[0])), test_pred[:, 0], label="Predicted test prices") plt.xlabel("Time") plt.ylabel("Stock Price") plt.title("Stock Price Prediction") plt.legend() plt.show() ``` 请确保将实际的股票价格数据保存在名为"stock_prices.csv"的文件中,并根据需要对代码进行适当的调整。

相关推荐

最新推荐

recommend-type

预测问题评价指标:MAE、MSE、R-Square、MAPE和RMSE

以上是对于预测问题的评价指标。 1.平均绝对误差(Mean Absolute Error, MAE) 误差越大,该值越大。 2.均方误差(Mean Squared Error, MSE) 误差越大,该值越大。 SSE(和方差)与MSE之间差一个系数n,即SSE = n *...
recommend-type

python之MSE、MAE、RMSE的使用

我就废话不多说啦,直接上代码吧! target = [1.5, 2.1, 3.3, -4.7, -2.3, 0.75] prediction = [0.5, 1.5, 2.1, -2.2, 0.1, -0.5] error = [] for i in range(len(target)): error.append(target[i] - prediction...
recommend-type

防雷及短路计算软件.zip

防雷及短路计算软件
recommend-type

电线穿管选用小软件.zip

电线穿管选用小软件
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

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

实现实时数据湖架构:Kafka与Hive集成

![实现实时数据湖架构:Kafka与Hive集成](https://img-blog.csdnimg.cn/img_convert/10eb2e6972b3b6086286fc64c0b3ee41.jpeg) # 1. 实时数据湖架构概述** 实时数据湖是一种现代数据管理架构,它允许企业以低延迟的方式收集、存储和处理大量数据。与传统数据仓库不同,实时数据湖不依赖于预先定义的模式,而是采用灵活的架构,可以处理各种数据类型和格式。这种架构为企业提供了以下优势: - **实时洞察:**实时数据湖允许企业访问最新的数据,从而做出更明智的决策。 - **数据民主化:**实时数据湖使各种利益相关者都可
recommend-type

用matlab绘制高斯色噪声情况下的频率估计CRLB,其中w(n)是零均值高斯色噪声,w(n)=0.8*w(n-1)+e(n),e(n)服从零均值方差为se的高斯分布

以下是用matlab绘制高斯色噪声情况下频率估计CRLB的代码: ```matlab % 参数设置 N = 100; % 信号长度 se = 0.5; % 噪声方差 w = zeros(N,1); % 高斯色噪声 w(1) = randn(1)*sqrt(se); for n = 2:N w(n) = 0.8*w(n-1) + randn(1)*sqrt(se); end % 计算频率估计CRLB fs = 1; % 采样频率 df = 0.01; % 频率分辨率 f = 0:df:fs/2; % 频率范围 M = length(f); CRLB = zeros(M,1); for
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。
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。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依