1、统计分析

统计分析

1、计算回报率

​ 通常在金融行业会用日回报率来计算贝塔值，以衡量一家公司的市场风险。我们可以通过价格数据来计算回报率（也称作收益率），有时会需要将日收益每周或者每月收益率，所以接下来主要介绍的就是收益率的计算。计算回报率最重要的就是要明白这些价格的排序。

import pandas as pd
import tushare as ts

data = ts.get_k_data('000001')  # 获取平安的k线数据
close = np.array(data.close)  # 转换为ndarray数组，以便做批量计算
ret  = (close[1:]-close[:-1])/close[:-1]
ret[0:2]
运行结果：
array([0.00775724, 0.00441476])

# 验证回报率：
> data.date[0:3]
>:
0    2016-08-05
1    2016-08-08
2    2016-08-09
Name: date, dtype: object
        
> data.close[0:3]
>:
0    8.766
1    8.834
2    8.873
Name: close, dtype: float64
        
> (8.834-8.766)/8.766
> 0.007757243896874243

1.1、从日回报率计算月回报率

$P_0$为上个月最后一个交易日的收盘价，$P_1$是这个月第一个交易日的收盘价。

月回报率：
$$R_{monthly}=\frac{P_{20}-P_0}{P_0}$$
月对数回报率：
$$logretur{n}_{monthly}=log(\frac{P_{20}}{P_0})$$
日对数回报率：
$$logretur{n}_i^{daily}=log(\frac{P_i}{P_{i-1}})$$
百分比回报率和对数回报率之间的关系：
$$R_{monthly}=exp(logreturn)-1$$
对数回报率的总和：
$$logretur{n}_{monthly}=log(\frac{P_{20}}{P_0})=\sum _{i=1}^{20}logretur{n}_i^{daily}$$

# 对数回报率
logret = np.log(close[1:]/close[:-1])

# 处理时间格式
dd = []
for dt in data.date:
    # 字符串格式的时间对象只有先转换为datetime格式才能进行转换到date格式
    d0 = dt.strptime(i,'%Y-%m-%d')  # 转换为datetime.datetime格式
    d1 = d0.date()  # 转换为datetime.date格式
    dd.append(d1)
# 将新列表转换为ndarray数组
da = np.array(dd)

# 月回报率
ym = []
for i in range(0,np.size(logret)):
    ym.append(''.join([da[i].strftime("%Y"),da[i].strftime("%m")]))
# 通过pandas进行批量计算
y = pd.DataFrame(logret,ym,columns=['ret_monthly'])
ret_monthly = y.groupby(y.index).sum()
ret_monthly
运行结果:
    	ret_monthly
201608	0.044402
201609	-0.035544
201610	0.002146
201611	0.049096
201612	-0.046842
...
201811	-0.022410
201812	-0.141794
201901	0.197798
201902	0.130401
201903	-0.053110

以上程序依旧通过获取的平安的K线数据来进行计算

1.2、从日回报率计算年回报率

方法类似

# 年回报率
date = []
for i in range(0,np.size(logret)):
    date.append(''.join([da[i].strftime("%Y")]))
y = pd.DataFrame(logret,date,columns=['ret_annual'])
ret_annual = np.exp(y.groupby(y.index).sum())-1
运行结果:
    	ret_annual
2016	0.013347
2017	0.517730
2018	-0.318350
2019	0.316649

2、T-检验

金融领域常用t-统计量来检验某个零假设。当零假设成立时，t-统计量服从t分布。

以下的代码生成1000个服从标准正态分布的随机数，然后进行两个检验。第一个检验的零假设是均值为0.5，第2个检验的零假设是均值为0.

from scipy import stats
import scipy as sp
import tushare as ts
import numpy as np

sp.random.seed(1235)
x = stats.norm.rvs(size=1000)
plt.figure(figsize=(7,4))

# 分别为两条数据添加图例
plt.plot(x,lw=1.0,label='1st')  
plt.grid(True)  # 网格设置
plt.legend(loc=0)  # 图例标签位置设置
plt.axis('tight')
plt.xlabel('index')
plt.ylabel('value')

图例:

print(stats.ttest_1samp(x,0.5))
print(stats.ttest_1samp(x,0))

运行结果:
Ttest_1sampResult(statistic=-16.755214405175575, pvalue=1.0140316700193316e-55)
Ttest_1sampResult(statistic=-1.1953686865174769, pvalue=0.23222669336297136)

测试中国平安的平均日回报率是否为0

data = ts.get_k_data('000001')
pa_data_close = np.array(data.close)
ret = pa_data_close[1:]/pa_data_close[:-1] - 1

print(round(sp.mean(ret),5),stats.ttest_1samp(ret,0))

运行结果:
0.00065 Ttest_1sampResult(statistic=0.9167715496440845, pvalue=0.35960876859507807)

结果显示中国平安的平均日回报率为0.065%，t值是0.92而p值是0.36。所以，统计检验的结果不支持平均日回报率等于0.

3、基于52周最高价和最低价的交易策略

​ 当某只股票今天的价格接近过去52周的最高价格，卖出该股票；当其今天得价格接近过去52周达到的最低价格，买入该股票。

import tushare as ts
from datetime import datetime
import numpy as np

wk_data = ts.get_k_data("000002",start='2018-01-01',end='2018-12-30')  # 获取万科地产2018一年的K线数据
wk_data = np.array(wk_data)  # 转换为ndarray数组
wk_data_1 = wk_data[-1]

y = np.array(wk_data.tolist())[:,-3]  # 最低价
aa = y.astype(np.float)  
high = max(aa)  #  算出一年中最高的
low = min(aa)  #  算出一年中最低的

print("""
today         price    high      low       from low
%s    %s    %s    %s    %s
"""%(wk_data_1[0],wk_data_1[1],high,low,round((wk_data_1[-3]-low)/(high-low)*100,2)))

运行结果：

today         price    high      low       from low
2019-04-19    30.61    38.939    9.964    70.74

根据过去52周最高价和最低价的交易策略

ata_1[1],high,low,round((wk_data_1[-3]-low)/(high-low)*100,2)))

运行结果：

today price high low from low
2019-04-19 30.61 38.939 9.964 70.74

根据过去52周最高价和最低价的交易策略