结合金融场景的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()
相关推荐
灏瀚星空20 小时前
从基础到实战的量化交易全流程学习:1.3 数学与统计学基础——概率与统计基础 | 基础概念
笔记·python·学习·金融·概率论
亿信华辰软件1 天前
大模型驱动金融数据应用的实战探索
金融
胡耀超2 天前
附1:深度解读:《金融数据安全 数据安全分级指南》——数据分类的艺术专栏系列
大数据·金融·数据治理·生命周期·数据分类·政策法规
油泼辣子多加3 天前
【风控】稳定性指标PSI
人工智能·算法·金融
oioihoii3 天前
金仓数据库 KingbaseES 产品深度优化提案:迈向卓越的全面升级
数据库·性能优化·金融·金仓数据库 2025 征文·数据库平替用金仓
财经汇报4 天前
联易融受邀参加上海审计局金融审计处专题交流座谈
金融
新讯网4 天前
马浩棋:产通链CT-Chain 破局不动产 RWA,引领数智金融新变革
金融
财经汇报4 天前
联易融出席深圳链主企业供应链金融座谈会,加速对接票交所系统
科技·金融
结冰架构5 天前
量子金融工程:蒙特卡洛算法误差压缩至0.3%
人工智能·算法·ai·金融·量子计算