关于【对数收益率是什么】:对数收益率是什么，今天向乾小编给您分享一下，如果对您有所帮助别忘了关注本站哦。
- 内容导航：
- 1、对数收益率是什么
- 2、资产收益率的非平稳性——为何机器学习预测效果不佳？
1、对数收益率是什么
对数收益率是2个阶段资产净值取对数后差值，即财产好几个阶段的对数收益率相当于其各阶段对数收益率之和。
考量股票价格股权投资收益的水准指标值具体有股利分配年化收益率与持有期收益率和拆股后持有期收益率等。
1、股利收益率
股利收益率，又被称为获利率，就是指股份有限公司以现钱方式发放的股利或分红与股市价钱的比例其计算方法为：
该年化收益率可以用测算已得的股利分配年化收益率，也可以用以预知未来很有可能的股利分配年化收益率。
2、持有期收益率
持有期收益率指投资人持有股票期内的股利收益和买卖交易价差之和与股票价格买入价的比例。其计算方法为：
股票价格都还没到期还款日的，投资人持有股票时间较短则几日、长则为多年，持有期收益率便是反应投资人在相应持有期里的所有股利收入及其资本利得占项目投资本钱的比例。
2、资产收益率的非平稳性——为何机器学习预测效果不佳？
02 数据获取
使用tushare获取股票指数收盘价数据，计算日对数收益率。
#先引入后面可能用到的包（package）import pandas as pd import numpy as npimport matplotlib.pyplot as pltimport seaborn as sns sns.set() %matplotlib inline #正常显示画图时出现的中文和负号from pylab import mplmpl.rcParams['font.sans-serif']=['SimHei']mpl.rcParams['axes.unicode_minus']=False
import tushare as tsdef get_price(code,start,end): df=ts.get_k_data(code,start,end) df.index=pd.to_datetime(df.date) return df.closecodes=['sh','sz','cyb','zxb','hs300']names=['上证综指','深证综指','创业板指','中小板指','沪深300']end_day = pd.to_datetime('2020-10-30')start_day = end_day - 10 * 252 * pd.tseries.offsets.BDay()start=start_day.strftime('%Y-%m-%d')end=end_day.strftime('%Y-%m-%d')#指数收盘价数据df = pd.DataFrame({name:get_price(code, start, end) for name,code in zip(names,codes)})
(df/df.iloc[0]-1).plot(figsize=(12,7))plt.title('指数累计收益率\n2011-2020',size=15)plt.show()
计算股票指数的日对数收益率。实际上在经济金融上，采用对数收益率已经是约定俗成了，当然这样处理主要是基于对数处理的统计特性比较适合建模，如为了使数据更加平滑，克服数据本身的异方差等。
rs=(np.log(df/df.shift(1))).dropna()rs.plot(figsize=(12,5))plt.title('指数日对数收益率',size=15)plt.show()
03 收益率平稳性可视化分析
所谓的时间序列平稳性，简单的理解是，序列的均值和方差不应该随着时间的推移而改变。时间序列中收益之间的协方差也不应是时间的函数。
下面将演示和分析在预测时间序列时当数据不满足模型假设条件的挑战。在开始之前，我们要知道，证券收益率序列往往不满足平稳性的要求。
import scipy.stats as statsdef add_mean_std_text(x, **kwargs): mean, std = x.mean(), x.std() mean_tx = f"均值: {mean:.4%}\n标准差: {std:.4%}" txkw = dict(size=14, fontweight='demi', color='red', rotation=0) ymin, ymax = plt.gca().get_ylim() plt.text(mean+0.025, 0.8*ymax, mean_tx, **txkw) returndef plot_dist(rs, ex): plt.rcParams['font.size'] = 14 g = (rs .pipe(sns.FacetGrid, height=5,aspect=1.5) .map(sns.distplot, ex, kde=False, fit=stats.norm, fit_kws={ 'lw':2.5, 'color':'red','label':'正态分布'}) .map(sns.distplot, ex, kde=False, fit=stats.laplace, fit_kws={'linestyle':'--','color':'blue', 'lw':2.5, 'label':'拉普拉斯分布'}) .map(sns.distplot, ex, kde=False, fit=stats.johnsonsu, fit_kws={'linestyle':'-','color':'green','lw':2.5, 'label':'约翰逊分布'}) .map(add_mean_std_text, ex)) g.add_legend() sns.despine(offset=1) plt.title(f'{ex}收益率',size=15) return
plot_dist(rs, '上证综指')
从上图我们可以看出正态分布对数据集的拟合程度是很差，而拉普拉斯分布和约翰逊分布在拟合收益率方面相对较好。统计分布图表明，在样本时间段内，股票指数收益率不是正态分布。
下面进一步用分位数图来比较参数的分位数分布，考察与正态分布的偏离情况。
plt.rcParams['figure.figsize'] = 10,7def quantile_plot(x, **kwargs): res = stats.probplot(x, fit=True, plot=plt) _slope, _int, _r = res[-1] ax = plt.gca() ax.get_lines()[0].set_marker('s') ax.get_lines()[0].set_markerfacecolor('r') ax.get_lines()[0].set_markersize(13.0) ax.get_children()[-2].set_fontsize(22.) txkw = dict(size=14, fontweight='demi', color='r') r2_tx = "r^2 = {:.2%}\nslope = {:.4f}".format(_r, _slope) ymin, ymax = ax.get_ylim() xmin, xmax = ax.get_xlim() ax.text(0.5*xmax, .8*ymin, r2_tx, **txkw) return
quantile_plot(rs['上证综指'])
分位数图清晰地显示出，上证综指收益率表现出“肥尾”的现象（概率分布图也可以看出是“尖峰后尾”）。这意味着，极端回报(无论是正回报还是负回报)出现的频率，远远高于正态分布预测的频率。
但是，我们还没有真正进入讨论时间序列平稳性的问题。
检验收益率是否平稳的一种方法是考察每年收益率的均值和标准差的变动情况。如果这两个变量不是时间的函数，那么它们每年的均值应该是非常相似的。
def plot_facet_hist(rs, ex): plt.rcParams['font.size'] = 12 df = rs.assign(year=lambda df: df.index.year) g = (sns.FacetGrid(df, col='year',col_wrap=2, height=4, aspect=1.2) .map(sns.distplot, ex, kde=False, fit=stats.norm, fit_kws={ 'lw':2.5,'color':'red', 'label':'正态分布'}) .map(sns.distplot, ex, kde=False, fit=stats.laplace, fit_kws={'linestyle':'--','color':'blue', 'lw':2.5, 'label':'拉普拉斯分布'}) .map(sns.distplot, ex, kde=False, fit=stats.johnsonsu, fit_kws={'linestyle':'-', 'color':'green','lw':2.5, 'label':'约翰逊分布'}) .map(add_mean_std_text, ex)) g.add_legend() g.fig.subplots_adjust(hspace=.20) sns.despine(offset=1) return
plot_facet_hist(rs, '上证综指')
(rs.groupby(rs.index.year)['上证综指'] .agg(['mean', 'std']) .plot(marker='o', subplots=True))plt.show()
显然，直方图中拟合的正态密度图在每一年中差异很大，某些年份近似服从正态分布，而某些年份表现出较强的“尖峰后尾”。这是否意味着我们之前关于收益率不是正态分布的结论只是一个时间尺度的问题?下面通过分位数图进一步考察这个问题。
def plot_facet_qq(rs, ex): df = rs.assign(year=lambda df: df.index.year) g = (df .pipe(sns.FacetGrid, col='year',col_wrap=2, height=7,aspect=1.3) .map(quantile_plot, ex) .fig.subplots_adjust(hspace=0.2)) sns.despine(offset=1, trim=True) return
plot_facet_qq(rs, '上证综指')
历年的收益率QQ图也显示出，在某些年份收益率分布近似满足正态分布；而在其他年份，则不满足正态分布。所以可以认为，通过视觉观察，收益率序列的均值和方差是时间的函数，即收益率序列的均值和方差会随着时间的变动而变化。
模拟正态分布样本和累积分布。
norm = stats.normRANDOM_STATE=888def generate_norm_rvs(ser, N=None): if not N: N = ser.shape[0] return norm.rvs(ser.mean(), ser.std(), size=N, random_state=RANDOM_STATE)def generate_norm_pdf(ser, N=None): if not N: N = ser.shape[0] _min, _max = ser.min(), ser.max() x = np.linspace(_min, _max, N) y = norm.pdf(x, ser.mean(), ser.std()) return x, ydef generate_norm_cdf(ser, N=None): if not N: N = ser.shape[0] _min, _max = ser.min(), ser.max() x = np.linspace(_min, _max, N) y = norm.cdf(x, ser.mean(), ser.std()) return x, y
最后，基于Kolmogorov-Smirnov对收益率分布进行统计检验。
使用收益率的累积分布（CDF），并使用kolmogorov-smirnov的二样本检验(scipy.stats.kstest)与拟合CDF或经验CDF进行比较。该检验可以验证两个样本的CDF是否来自相同的分布，例如正态分布。kstest输出ks统计量和一个p值。p值小于0.05意味着我们可以拒绝两个系列来自同一分布的原假设。p值越小，我们就越确信它们来自不同的分布。
首先，我们将分析整个收益率序列，然后比较各个年份的检验结果。
def plot_cdf(ser, **kwds): g = sns.kdeplot(ser, cumulative=True, lw=3, color='blue') x, y = generate_norm_cdf(ser) # 生成正态分布CDF g.plot(x, y, color='red', lw=3, label='fitted normal CDF') ks, p = stats.kstest(ser, 'norm', args=(ser.mean(), ser.std())) xmin,xmax=plt.gca().get_xlim() ymin,ymax=plt.gca().get_ylim() txkw = dict(size=14, fontweight='demi', color='red', rotation=0) tx_N = ser.shape[0] tx_args = (ks, p, tx_N, ser.shape[0]) tx = 'ks 2 sample test\nks: {:.4f}, p: {:.4f}\nrvs sample N: {:.0f}\nreturn sample N: {:.0f}'.format(*tx_args) plt.text(xmax*0.2, 0.2*ymax, tx, **txkw) sns.despine(offset=1) (plt.legend(frameon=True, prop={'weight':'demi', 'size':12}) .get_frame()) return
plot_cdf(rs['上证综指'])
def plot_facet_cdf(rs, ex, **kwds): df=rs.assign(year=lambda df:df.index.year) g = (df .pipe(sns.FacetGrid, col='year', col_wrap=2, height=5, aspect=1.3) .map(plot_cdf, ex, **kwds)) g.add_legend() g.fig.subplots_adjust(hspace=.20) sns.despine(offset=1) return
plot_facet_cdf(rs, '上证综指')
总体收益率的自相关系数和偏相关系数图
关于时间序列的自\偏相关系数图代码可参见《【Python量化基础】时间序列的自相关性与平稳性》。
历年日收益率的自相关系数和偏相关系数图
上述图形清晰表明，上证综指收益率在均值和方差上是非平稳的，并且显示了变量自相关的时间段。此外，我们看到上证综指收益率的总体分布非正态分布，但在各种不同的时间期间，如每年，有一段时期是近似正态分布的。那么，当我们试图使用过去的收益率近似于正态分布的统计指标(均值、方差)来预测未来的收益分布时，会出现什么问题呢?
为了探究这个问题，下面使用机器学习包sklearn的 TimeSeriesSplit方法考察样本的分布情况，该方法提供了交叉验证的前向walf-forward形式，使用先前的数据点来序列地预测下一个周期，从而保留了与时间相关的信息。
from sklearn.model_selection import TimeSeriesSplitxx = rs['上证综指'].copy() _base = 252 # 1年测试样本_max_train_sizes = [_base*1, _base*2, _base*3, _base*5] _n_split=5 gs = gridspec.GridSpec(_n_split, len(_max_train_sizes), wspace=0.0)fig = plt.figure(figsize=(20,25))rows = []for j, max_size in enumerate(_max_train_sizes): tscv = TimeSeriesSplit(n_splits=_n_split, max_train_size=max_size) for i, (train, test) in enumerate(tscv.split(xx)): tmp_train = xx.iloc[train] tmp_test = xx.iloc[test] min_train_dt, max_train_dt = tmp_train.index.min(), tmp_train.index.max() min_test_dt, max_test_dt = tmp_test.index.min(), tmp_test.index.max() ks, p = stats.ks_2samp(tmp_train, tmp_test) # 获得ks检验统计指标 df_row = (max_size, ks, p, min_train_dt.date(), max_train_dt.date(), min_test_dt.date(), max_test_dt.date()) rows.append(df_row) tmp_ax = plt.subplot(gs[i, j]) if i in [0,1,2,3,4] and j != 0: tmp_ax.set_yticks([]) sns.kdeplot(tmp_train, cumulative=True, lw=3, color='blue', ax=tmp_ax, label='train') sns.kdeplot(tmp_test, cumulative=True, lw=3, color='red', ax=tmp_ax, label='test') plt.title('max train size: {}, ks: {:.4f}, p: {:.4f}\ntrain dates: {}_{}\ntest dates: {}_{}' .format(max_size, ks, p, min_train_dt.date(), max_train_dt.date(), min_test_dt.date(), max_test_dt.date()), fontsize=11.) plt.subplots_adjust(top=1.03) plt.tight_layout()
从图中不难看出，在不同的回溯期间，训练样本和测试样本CDF的变化程度，包括它们准确同步的时间。为了量化我们的观察结果，可以创建一个dataframe，其中包含ks测试的输出以及回溯期和数据拆分（即训练集和测试集）的信息。
# create data table with ks test results and tss informationcols = ['max_size', 'ks', 'p', 'min_train_dt', 'max_train_dt', 'min_test_dt', 'max_test_dt']tss_ks_df = pd.DataFrame(rows, columns=cols)obj_cols = tss_ks_df.select_dtypes(include=['object']).columns接下来，可以检查在0.05和0.01的置信水平下，拒绝alpha值的原假设中训练/测试对中所占的百分比。
get_pct_reject = lambda df, col, alpha: df.query(f'{col} < {alpha}').shape[0] / df.shape[0]a1 = 0.05a2 = 0.01get_pct_reject(tss_ks_df, 'p', a1), get_pct_reject(tss_ks_df, 'p', a2)
(0.5, 0.4)
结果似乎不太理想，50%的训练/测试对拒绝原假设（原假设是训练和测试数据来自相同的分布）。这表明，如果我们用训练数据的均值/方差参数来预测测试数据的收益率分布，那么在很多时候都会出错。这从另一个角度表明，数学建模（如机器学习或深度学习）在预测金融资产收益率上往往效果不佳。
04 结语
对金融时间序列进行预测是量化建模分析的重要组成部分，不少量化交易策略正是基于历史数据的预测而构建的。本文以指数收益率为例（上证综指），基于图形观察的视角，展示了收益率的统计特征和“尖峰厚尾”现象，揭示出其不满足通常条件下的平稳性假设条件，同时利用机器学习的方法，考察数学建模对收益率进行预测可能面临的问题——训练集和测试集来自不同的分布。本文也从侧面反映了盲目照搬照套数学模型（含机器学习和深度学习）进行时间序列预测可能效果不佳。
参考资料：
Mixture Model Trading (Part 1 - Motivation)
http://www.blackarbs.com/blog/mixture-model-trading-part-1/1/16/2018
本文关键词：对数收益率是什么意思，对数收益率和简单收益率，收益率和对数收益率，对数收益率和简单收益率的关系，使用对数收益率的优点。这就是关于《对数收益率是什么，对数收益率和简单收益率（资产收益率的非平稳性——为何机器学习预测效果不佳）》的所有内容，希望对您能有所帮助！更多的知识请继续关注《赛仁金融》百科知识网站：http://yzsryq.com/！
还没有评论，来说两句吧...