结合金融场景的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()
相关推荐
青云交6 小时前
Java 大视界 -- Java 大数据机器学习模型在金融衍生品定价中的创新方法与实践(166)
java·大数据·金融·数据采集·机器学习模型·java 大数据·金融衍生品定价
NEET_LH1 天前
金融数据分析(Python)个人学习笔记(6):安装相关软件
python·金融·数据分析
程序猿chen1 天前
《JVM考古现场(十六):太初奇点——从普朗克常量到宇宙弦的编译风暴》
jvm·git·后端·程序人生·金融·java-ee·量子计算
Yvsanf2 天前
交易引擎中的设计模式
设计模式·金融
CryptoPP4 天前
springboot 对接马来西亚数据源API等多个国家的数据源
spring boot·后端·python·金融·区块链
进击的雷神4 天前
质量金融:如何用SLA数据发行债券?
金融·质量工程师
人类群星闪耀时4 天前
区块链技术如何重塑金融衍生品市场?
金融·区块链
CryptoPP4 天前
深入实践:基于WebSocket的全球化金融数据实时对接方案。 马来西亚、印度、美国金融数据API
websocket·网络协议·金融
雷焰财经5 天前
中和农信:让金融“活水”精准浇灌乡村沃土
金融
银河金融数据库5 天前
历史分钟高频数据
数据库·金融