Python用GARCH、离散随机波动率模型DSV模拟和估计股票收益时间序列与蒙特卡洛可视化

这篇文章介绍了一类离散随机波动率模型。

由Kaizong Ye,Coin Ge撰写

介绍了一些特殊情况,包括 GARCH 和 ARCH 模型。

本文展示了如何模拟这些过程以及参数估计。本文为这些实验编写的 Python 代码在文章末尾引用。

离散随机波动率模型

是一个随机基,有一个完整的  的可测量子集  , 一个概率测量和一个过滤


视频

时间序列分析模型 ARIMA-ARCH GARCH模型分析股票价格数据

探索见解

去bilibili观看

探索更多视频


视频

马尔可夫链蒙特卡罗方法MCMC原理与R语言实现

探索见解

去bilibili观看

探索更多视频


视频

随机波动率SV模型原理和Python对标普SP500股票指数时间序列波动性预测

探索见解

去bilibili观看

探索更多视频

  • 因此,时间实例使用非负整数进行索引 
  • 获取序列的第一个 t元素 , 记号.

_离散随机波动率_(   DSV) 模型  是一个实值 stochastic process (一系列随机变量)满足以下方程:

其中:

  • Z 是 F 的噪声过程。
  • φi 是实数,我假设  并且 gi ,hi 是非负值。
  • fi 、gi 和 h_ihi 是过程的确定性函数。
  • 过程 通常称为 _偏移_,而 σ 称为  X的_波动率。_因为σ 是一个随机过程,所以上面定义的过程 X 属于一个随机波动率模型的大家族。
  • 对于噪声过程 Z,使得每个 Z_t的均值和方差都存在,我们有  和 .

案例

制定通用 DSV 模型的特化:

后移算子 ,对于 ,产生其参数过程的滞后版本,即 , 和 , 如果 . 例如

为方便起见,我设置  和 .

对于下面列表中的所有特殊情况,我假设函数 f\_i、gi 和 hi 从参数过程的历史中选择一个元素,即 , 和 .

GARCH 过程定义另外设置


Python使用GARCH,EGARCH,GJR-GARCH模型和蒙特卡洛模拟进行股价预测

阅读文章


模拟

离散随机波动率模型通常用于对观察到的时间序列的对数收益进行建模。因此,为了模拟原始时间序列的路径,我们需要模拟其对数收益并计算 .

由带参数的高斯噪声驱动的 GARCH(1,1) 过程的样本路径 :

path( \[0.001, 0.2, 0.25\])
cumprod* repeat.reshape
plt.subplots


随时关注您喜欢的主题


注意 σ 过程为  不能低于 ≈0.0353

最大似然估计

最大似然(ML)参数估计是所有讨论模型的选择方法,因为转换密度,即给定过去信息的 X_t 的密度是明确已知的。因此,过程样本路径 x 的对数似然函数由下式给出

其中,而 是 Z的密度。将上述对数似然函数最小化可得到 的最大似然估计

.

蒙特卡罗研究

为了测试 ML 参数估计过程,我进行了以下蒙特卡罗实验。

  • 使用参数 (0.001, 0.2, 0.25) 模拟长度为 5000 的 2500 个独立 GARCH(1,1) 过程路径。我使用了高斯噪声,即 .
  • 将这些路径中的每一个都输入到 ML 估计并获得估计的参数向量 .
  • 此优化过程中参数的搜索范围限制为 [1e-8, 1]。
  • 将原始与估计的 进行比较。
  • 使用参数向量模拟 GARCH(1,1),计算均值和标准差,并将它们与“真实”均值和标准差(分别为 5.098 和 1.084)进行比较。

正如期望的那样,估计量 非常不准确,并且在大多数情况下,甚至不接近真实向量。特别是,估计的通常设置为零(参见下面的直方图)。

ps = \[0.001, 0.2, 0.25\]

cumprod * repeat
print, np.std

另一方面,来自估计的的过程均值和标准偏差要准确得多。这是一件好事,因为我们通常更关心恢复未知数据生成过程的特征,而不是模型的真实参数值。

    mes, stvs, esms 
  
    ax\[1\].hist
    fig.tight_layout

噪声过程  不必归一化为均值 0 和方差 1。

柯西噪音

实际上,我们只需要确保随机变量 Zt 的分布具有密度即可。如果是这种情况,过程模拟和 ML 估计都可以按照描述的方式工作。

那么如何用从柯西分布中采样的噪声替换高斯噪声呢?在许多概率论书籍中,柯西分布被用作反例,因为它具有许多“病态”特性。例如,它没有均值,因此也没有方差。

我不知道柯西分布中的不稳定样本是什么样子的。看一下带有参数向量的 GARCH(1,1) 过程的示例路径 :

如果使用路径生成函数的时间足够长,甚至可能会生成溢出异常。因此,我用来生成上面显示的直方图的 Python 函数失败了。为了了解原因,让我们使用来自柯西分布的样本生成一些直方图:

柯西分布具有分位数函数

评估 给出

这意味着,例如,在 0.0001 的概率下,采样值大于 3183.10。为了比较,让我们计算标准正态分布的相应分位数:

norm.ppf(0.99)
norm.ppf(0.999)
norm.ppf(0.9999)


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。在此对他对本文所作的贡献表示诚挚感谢,他在上海财经大学完成了统计学专业的硕士学位,专注人工智能领域。擅长Python.Matlab仿真、视觉处理、神经网络、数据分析。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。

​非常感谢您阅读本文,如需帮助请联系我们!

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498

关注有关新文章的微信公众号


永远不要错过任何见解。当新文章发表时,我们会通过微信公众号向您推送。

技术干货

最新洞察

This will close in 0 seconds