结合金融场景的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()
相关推荐
uncle_ll2 小时前
金融科技白皮书:2022-2023年度回顾与前瞻
科技·金融
Hotchip华芯邦科技1 天前
MEMS硅麦克风应用电子烟雾化产业稳步爬升,耐高温、 防油、防酸、防腐蚀等性能优势和可实现自动化贴片及极高的一致性等特性使其必将成为主流
科技·单片机·金融·生活·社交电子·健康医疗·制造
Leuanghing2 天前
使用Python生成卡方分布表并导出为Excel文件
python·excel·pandas·scipy·卡方分布表
gooxi_hui2 天前
Gooxi受邀参加海通证券AI+应用生态大会,助力数智金融高质量发展
人工智能·金融
君臣Andy2 天前
金融量化框架选择
金融
AIGC方案2 天前
免费下载 | 2024年金融数字化转型白皮书
金融
中国云报2 天前
金融数据中心容灾“大咖说” | Veritas的“高举高打”之道
金融
lyw_YTU_Sussex2 天前
人工智能在金融领域的创新与应用
人工智能·金融
NiNg_1_2342 天前
Python中Tushare(金融数据库)入门详解
数据库·python·金融
中国云报2 天前
金融数据中心容灾“大咖说” | 英方软件的“分而治之”之道
金融