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

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

这篇文章介绍了一类离散随机波动率模型，并介绍了一些特殊情况，包括 GARCH 和 ARCH 模型。本文展示了如何模拟这些过程以及参数估计。这些实验编写的 Python 代码在文章末尾引用。

离散随机波动率模型

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

• 因此，时间实例使用非负整数进行索引 

• 获取序列的第一个 t元素 , 记

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

其中：

• Z 是 F 的噪声过程。

• φi 是实数，我假设  并且 gi ,hi 是非负值。

• fi 、gi 和 h_ihi 是过程的确定性函数。

• 过程  通常称为 _偏移_，而 σ 称为  X的_波动率。_因为σ 是一个随机过程，所以上面定义的过程 X 属于一个随机波动率模型的大家族。

• 对于噪声过程 Z，使得每个 Z_t的均值和方差都存在，我们有  和 .

案例

制定通用 DSV 模型的特化：

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

为方便起见，我设置  和 .

对于下面列表中的所有特殊情况，我假设函数 fi 、gi 和 h_i 从参数过程中选择一个元素，即 ， 和 .

GARCH 过程定义另外设置。

_GARCH(1, 1)_过程非常流行：

在 ARCH 过程中，波动性具有简化形式，对于所有 i，λi = 0，并且 。

_ARCH(1)_过程还 满足  对所有  :

模拟

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

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

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

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

---

点击标题查阅往期内容

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

左右滑动查看更多

01

02

03

04

最大似然估计

最大似然（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)

---

本文摘选《PYTHON用GARCH、离散随机波动率模型DSV模拟和估计股票收益时间序列与蒙特卡洛可视化》，点击“阅读原文”获取全文完整资料。

---

点击标题查阅往期内容

极值理论 EVT、POT超阈值、GARCH 模型分析股票指数VaR、条件CVaR：多元化投资组合预测风险测度分析

Garch波动率预测的区制转移交易策略

金融时间序列模型ARIMA 和GARCH 在股票市场预测应用

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

R语言风险价值：ARIMA，GARCH，Delta-normal法滚动估计VaR（Value at Risk）和回测分析股票数据

R语言GARCH建模常用软件包比较、拟合标准普尔SP 500指数波动率时间序列和预测可视化

Python金融时间序列模型ARIMA 和GARCH 在股票市场预测应用

MATLAB用GARCH模型对股票市场收益率时间序列波动的拟合与预测

R语言GARCH-DCC模型和DCC（MVT）建模估计

Python 用ARIMA、GARCH模型预测分析股票市场收益率时间序列

R语言中的时间序列分析模型：ARIMA-ARCH / GARCH模型分析股票价格

R语言ARIMA-GARCH波动率模型预测股票市场苹果公司日收益率时间序列

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

R语言时间序列GARCH模型分析股市波动率

R语言ARMA-EGARCH模型、集成预测算法对SPX实际波动率进行预测

matlab实现MCMC的马尔可夫转换ARMA - GARCH模型估计

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

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

R语言用多元ARMA,GARCH ,EWMA, ETS,随机波动率SV模型对金融时间序列数据建模

R语言股票市场指数：ARMA-GARCH模型和对数收益率数据探索性分析

R语言多元Copula GARCH 模型时间序列预测

R语言使用多元AR-GARCH模型衡量市场风险

R语言中的时间序列分析模型：ARIMA-ARCH / GARCH模型分析股票价格

R语言用Garch模型和回归模型对股票价格分析

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

matlab估计arma garch 条件均值和方差模型

R语言POT超阈值模型和极值理论EVT分析

R语言极值推断：广义帕累托分布GPD使用极大似然估计、轮廓似然估计、Delta法

R语言极值理论EVT：基于GPD模型的火灾损失分布分析

R语言有极值（EVT）依赖结构的马尔可夫链(MC)对洪水极值分析

R语言POT超阈值模型和极值理论EVT分析

R语言混合正态分布极大似然估计和EM算法

R语言多项式线性模型：最大似然估计二次曲线

R语言Wald检验 vs 似然比检验

R语言GARCH-DCC模型和DCC（MVT）建模估计

R语言非参数方法：使用核回归平滑估计和K-NN(K近邻算法)分类预测心脏病数据

matlab实现MCMC的马尔可夫转换ARMA - GARCH模型估计

R语言基于Bootstrap的线性回归预测置信区间估计方法

R语言随机搜索变量选择SSVS估计贝叶斯向量自回归（BVAR）模型

Matlab马尔可夫链蒙特卡罗法（MCMC）估计随机波动率（SV，Stochastic Volatility） 模型

Matlab马尔可夫区制转换动态回归模型估计GDP增长率

R语言极值推断：广义帕累托分布GPD使用极大似然估计、轮廓似然估计、Delta法

欲获取全文文件，请点击左下角“阅读原文”。

欲获取全文文件，请点击左下角“阅读原文”。