原文链接：http://tecdat.cn/?p=20678

预测股价曾经受到了投资者，政府，企业和学者宽泛的关注。然而，数据的非线性和非平稳性使得开发预测模型成为一项简单而具备挑战性的工作。在本文中，我将解释如何将 GARCH，EGARCH和 GJR-GARCH 模型与Monte-Carlo 模仿联合应用，以建设无效的预测模型。金融工夫序列的峰度，稳定率和杠杆效应特色证实了GARCH的合理性。工夫序列的非线性特色用于查看布朗运动并钻研工夫演变模式。非线性预测和信号分析方法因其在特征提取和分类中的鲁棒性而在股票市场上越来越受欢迎。

动力学零碎能够用一组时变（间断或离散）变量来形容，这些变量形成信号剖析非线性办法的根底。如果工夫的以后值和状态变量能够精确地形容下一时刻的零碎状态，则能够说这样的零碎是确定性的。另一方面，如果工夫和状态变量的以后值仅形容状态变量的值随工夫变动的概率，则将动力学零碎视为随机零碎。

因而，在应用GARCH 建模办法之前，我将采纳 分形维数（FD），重定范畴和 递归量化剖析（RQA）的 数据建模 技术来总结数据的非线性动力学行为并实现钻研指标。

办法

Hurst系数 （H） 是长期依赖的特征参数，与 FD （_FD + H = 2_）无关。 R / S剖析是数据建模的外围工具。教训钻研表明，与同类中的其余办法相比， _R / S_带来了更好的后果，例如自相干，光谱合成的剖析。它是度量工夫序列差别的度量，该工夫序列的定义是给定持续时间_（T）_的均值范畴，除以该持续时间的标准偏差 [ R / S = k * T（H） ]； _ķ_ 是一个取决于工夫序列的常数。H度量工夫序列的长期记忆，将其表征为均值回复，趋势或随机游走。

H <0.5示意均值回复
H> 0.5示意趋势序列，并且
H = 0.5示意随机游走。

我将展现如何应用 GARCH 模型进行危险评估。

GARCH 模型的一个要害限度是对其参数施加非负束缚，以确保条件方差的正性。这样的束缚会给预计GARCH 模型带来艰难。

因而，提出了 非对称GARCH 模型，即俗称的 GJR-GARCH 模型，以解决对称GARCH 模型的局限性。更重要的是，指数 GARCH 或 EGARCH 模型绝对于传统的GARCH 模型具备潜在的改良。

数据挖掘

查看数据。

在过来的几十年中，原油价格呈现出较大的稳定，尤其是在2008年左右。能够看出，随着屡次回升和降落，价格放弃在绝对较低的程度。从自相干图能够看出原始数据中显著的自相干。QQ和PP图的形态表明该过程靠近正态，然而_重尾散布_。

简略收益率的罕用模式为：r（t）= {p（t）— p（t-1）} / p（t-1），对数收益率= ln（pt / p（t-1），pt每日原油价格，r（t）是每日收益。

对数收益率在这里被视为本文的每日收益率。原始价格和对数收益率的直观显示分明地证实了以简直恒定的均值，应用对数收益率是正当的。

收益率序列图显示了高和低变动周期。在图中能够看到一个随机且集中在零左近的过程。大幅度稳定的正收益和负收益都减少了风险投资和治理的难度。每日收益率的平均值基本上在零程度程度左近，并且具备显著的波动性聚类，表明存在异方差性。ACF很小，然而高度相干。QQ和PP图的形态没有显著变动。

sns.distplot(df.returns, color=’blue’) #密度图# 汇总统计数据print(df.returns.describe())

收益率的偏度（-0.119）和向右偏离表明，收益率反比负收益率高，峰度（7.042）反映了油价稳定大。

规范正态分布的偏度和峰度别离为0和3。

Jarque-Bera测验的值表明，传统的正态分布假如不适用于原油收益的实在散布。

ADF = ADF(df.returns)print(ADF.summary().as_text())kpss = KPSS(df.returns)print(kpss.summary().as_text())

进行了VR测验，以测试对数收益率序列是否是纯正的随机游走，以及是否具备肯定的可预测性。我在这里比拟了1个月和12个月的对数收益率，并且回绝了该系列为纯随机游走的空值。用负测验统计量VA（-11.07）回绝零示意在工夫序列中存在序列相关性。用ADF，KPSS，DFGLS，PP和ZA统计量对单位根和平稳性进行的测验均显示出显着性，表明应用 GARCH型模型来拟合收益序列是适合的。

非线性动力学

应用_Hurst_对平稳性的钻研。

# 计算最近价格的Hurst系数tau = \[sqrt(std(subtract(closes\_recent\[lag:\], closes\_recent\[:-lag\]))) for lag in lags\]m = polyfit(log(lags), log(tau), 1)hurst = m\[0\]*2

_H_（0.531）示意具备长期依赖性的随机静止工夫序列。证实了本钻研中GARCH模型的合理性。

computation = RQAComputation.create(settings,verbose=True)result = computation.run()result.min\_diagonal\_line_length = 2

此处，低 R 示意较低的周期性和随机行为。此外，较低的 DET 值示意不确定性。这证实了应用GARCH 办法的合理性。

GARCH模型

在估算GARCH类型的模型之前，将收益率乘以100。因为稳定率截距与模型中其余参数十分靠近，因而这有助于优化程序进行转换。

X = 100* df.returns

让咱们拟合一个 ARCH 模型并绘制平方残差以查看自相关性。

def getbest(TS):best_aic = np.inffor i in pq_rng:for d in d_rng:for j in pq_rng:try:tmp_mdl = smt.ARIMA(TS, order=(i,d,j)).fit(#aic: 22462.01 | order: (2, 0, 2)

gam = arch_model(Model.resid, p=2, o=0, q=2, dist=’StudentsT’)gres = gam.fit(update_freq=5, disp=’off’)print(gres.summary())tsplot(gres.resid**2, lags=30)

咱们能够看到平方残差具备自相干的根据。让咱们拟合一个GARCH模型并查看其性能。我将依照以下步骤进行操作：

通过ARIMA（p，d，q） 模型的组合进行迭代，以拟合最优工夫序列。
依据具备最低AIC的ARIMA模型抉择 GARCH模型。
将 GARCH（p，q） 模型拟合到工夫序列。
查看模型残差和平方残差进行自相干

因而，咱们在这里发现，最好的模型是 ARIMA（2,0,2）。当初，咱们对残差进行绘图，以确定它们是否具备条件异方差。

 arch_model(X, p=2, q=2, o=1,power=2.0, vol=’Garch’, dist=’StudentsT’)

am = arch_model(X, p=2, q=2, o=1,power=2.0, vol=’Garch’, dist=’StudentsT’)

所有3个GARCH 模型的输入都以表格格局显示。 _（）_ 是白噪声，alpha和beta是模型的参数。此外， _[1] +[1] <1_ 示意稳固的模型。 EGARCH 仿佛是最好的三个这模型。

最好在训练/测试中拆分数据并取得MSE / MAE / RMSE后果以比拟最佳模型拟合。

通过将残差除以条件稳定率来计算标准化残差。

std\_resid =  resid /  conditional\_volatilityunit\_var\_resid =  resid /  resid.std()

标准化残差和条件稳定图显示了一些误差，但幅度不大。

sns.kdeplot(squared_resid, shade=True)sns.kdeplot(std_resid, shade=True)sns.kdeplot(unit\_var\_resid, shade=True)

还标绘了标准化残差以及非标准化的残差。残差的平方在核心更加尖峰，表明散布的尾部比规范残差的尾部更重。让咱们检查一下ACF图。

plot\_acf(std\_resid)

看起来有些尖峰超出了暗影的相信区。让咱们查看残差平方。

残差平方显示数据点在蓝色暗影的置信度区域（95％）内，示意模型拟合较好。

res = am.fit()fig = res.hedgehog_plot(type=’mean’)

图显示了整个2019年的预测。橙色线示意在不同工夫区间的预测。

基于模仿的预测

这里应用基于仿真的办法从EGARCH 模仿中取得预测稳定率的置信区间。要从EGARCH 模型取得稳定预测，该模型是从拟合模型的最初一次察看中模仿得出的。反复此过程很屡次，以取得稳定率预测。预测点是通过对模仿求平均值来计算的，别离应用模仿散布的2.5％和97.5％的分位数来计算95％的置信区间。思考均匀收益率（mu）为0.0292，年稳定率（vol）为（26.48）* sqrt 252 = 37.37％。

#定义变量T = 252 #交易天数mu = 0.0622 #收益vol = 0.3737 #稳定率daily_returns=np.random.normal((1+mu)**(1/T),vol/sqrt(T),T)#生成图-价格序列和每日收益的直方图plt.plot(price_list)plt.hist(daily_returns-1, 100)

在最上方的图上，依据遵循正态分布的随机每日收益，显示了一个交易年度（252天）内潜在价格序列演变的模仿。第二个图是一年中这些随机每日收益的直方图。然而，能够通过运行成千上万的模仿来取得洞察，每次模仿都基于雷同的特色（价格交易量）产生一系列不同的潜在价格演变。

#设置一个空列表来保留咱们每个模仿价格序列的最终值result = \[\]S = df.Price\[-1\] #起始股票价格（即最初可用的理论股票价格）T = 252 #交易天数mu = 0.0622 #收益率vol = 0.3737 #稳定率#抉择要模仿的运行次数-我抉择了10,000for i in range(10000):#应用随机正态分布创立每日收益表daily_returns= np.random.normal((1+mu)**(1/T),vol/sqrt(T),T)#设定起始价格，并创立由上述随机每日收益生成的价格序列#将每次模仿运行的完结值增加到咱们在开始时创立的空列表中result.append(price_list\[-1\])

因为这些是对每日收益的随机模仿，因而此处的后果会略有不同。因为每个模仿中蕴含的门路，平均值趋向于“ mu”应用的均匀收益。上面的直方图显示了价格散布的两个分位数，以理解高收益率或低收益率的可能性。

很显然，原油价格有5％的机会最终跌破29.72元，有5％的机会高于101.75美元。

概括

在原油价格高稳定的背景下，我钻研并提出了混合时变长记忆 GARCH 和基于模仿的预测模型，该模型思考了诸如非对称性和异方差，时变危险，长记忆和重尾散布等稳定事实。教训证据表明，具备布朗运动的原油数据往往在其工夫动静方面显示出肯定水平的可预测性。这项钻研思考了2000年至2019年的数据，过后股市经验了几次金融危机和危机后阶段。应用此期间的数据训练的模型无望具备杰出的预测能力。

当解决长时间稳定的原油价格的工夫序列数据时，GARCH （2,2）模型预计了方差的持久性。进行了蒙特卡洛剖析，以查看后果的稳健性。蒙特卡洛 模仿的输入表明，即便在管制了无关因素之后，后果依然是牢靠的。因而，这些发现提供了杰出的混合 EGARCH 和 蒙特卡洛 模仿的的预测模型，其中思考了波动性特色，如波动性聚类和不对称性，时变危险和重尾散布，来掂量原油价格。

最受欢迎的见解

1.HAR-RV-J与递归神经网络（RNN）混合模型预测和交易大型股票指数的高频稳定率

2.R语言中基于混合数据抽样(MIDAS)回归的HAR-RV模型预测GDP增长

3.稳定率的实现：ARCH模型与HAR-RV模型

4.R语言ARMA-EGARCH模型、集成预测算法对SPX理论稳定率进行预测

5.GARCH（1,1），MA以及历史模拟法的VaR比拟

6.R语言多元COPULA GARCH 模型工夫序列预测

7.R语言基于ARMA-GARCH过程的VAR拟合和预测

8.matlab预测ARMA-GARCH 条件均值和方差模型

9.R语言对S＆P500股票指数进行ARIMA + GARCH交易策略