结合金融场景的Scipy模块编程

结合金融场景的Scipy模块编程

数据链接:https://pan.baidu.com/s/1VMh8-4IeCUYXB9p3rL45qw

提取码:c6ys

python 复制代码
import numpy as np
import pandas as pd
import statsmodels
import matplotlib.pyplot as plt
from pylab import mpl
mpl.rcParams['font.sans-serif']=['FangSong']
mpl.rcParams['axes.unicode_minus']=False
python 复制代码
P_BoComm=pd.read_excel('data/交通银行A股和H股每日收盘价数据(2018年至2020年).xlsx',sheet_name='Sheet1',header=0,index_col=0)

P_BoComm.plot(figsize=(9,6),grid=True, fontsize=13) # 股票可视化
python 复制代码
R_BoComm=np.log(P_BoComm/P_BoComm.shift(1)) # 计算交通银行股票的对数收益率
R_BoComm.describe()
python 复制代码
R_BoComm.isnull().any()
python 复制代码
Index_dropna = R_BoComm.dropna()  # 直接删除
Index_dropna
python 复制代码
import statsmodels.api as sma

Y = Index_dropna.iloc[:,-1]
X = Index_dropna.iloc[:,0]

X_addcons = sma.add_constant(X)

model = sma.OLS(endog=Y, exog=X_addcons)

result=model.fit()
result.summary()
python 复制代码
plt.figure(figsize=(9,6))
plt.scatter(X,Y,c='b',marker='o')
plt.plot(X,result.params[0]+result.params[1]*X,'r-',lw=2.5)
plt.xticks(fontsize=13)
plt.xlabel('交通银行A股日收益率',fontsize=13)
plt.yticks(fontsize=13)
plt.ylabel('交通银行H股日收益率',fontsize=13)
plt.title('交通银行A股日收益率的散点图和线性拟合',fontsize=13)
plt.grid()
plt.show()
相关推荐
Ashlee_code12 小时前
美联储降息趋缓叠加能源需求下调,泰国证券交易所新一代交易系统架构方案——高合规、强韧性、本地化的跨境金融基础设施解决方案
java·算法·金融·架构·系统架构·区块链·需求分析
万粉变现经纪人2 天前
如何解决pip安装报错ModuleNotFoundError: No module named ‘tkinter’问题
python·beautifulsoup·pandas·pip·策略模式·httpx·scipy
quant_19863 天前
如何通过 WebSocket 接口订阅实时外汇行情数据(PHP 示例)
开发语言·网络·后端·websocket·网络协议·金融·php
boyedu5 天前
区块链应用场景深度解读:从金融革命到社会治理的全方位革新
金融·区块链
前端开发与ui设计的老司机5 天前
从UI设计到数字孪生实战应用:构建智慧金融的风险评估与预警平台
ui·金融
超龄超能程序猿6 天前
(5)机器学习小白入门 YOLOv:数据需求与图像不足应对策略
人工智能·python·机器学习·numpy·pandas·scipy
马特说6 天前
基于随机森林的金融时间序列预测系统:从数据处理到实时预测的完整流水线
算法·随机森林·金融
超龄超能程序猿7 天前
使用 Python 对本地图片进行图像分类
开发语言·人工智能·python·机器学习·分类·数据挖掘·scipy
Gauss松鼠会8 天前
GaussDB应用场景全景解析:从金融核心到物联网的分布式数据库实践
数据库·分布式·物联网·金融·database·gaussdb
马特说8 天前
金融时间序列机器学习训练前的数据格式验证系统设计与实现
python·机器学习·金融