首页 > 编程语言 >Python用GARCH、离散随机波动率模型DSV模拟和估计股

Python用GARCH、离散随机波动率模型DSV模拟和估计股

时间:2024-07-28 17:18:02浏览次数:18  
标签:Python 模型 波动 DSV 随机 GARCH 过程 模拟

原文链接: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

 

最大似然估计

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


最受欢迎的见解

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交易策略

标签:Python,模型,波动,DSV,随机,GARCH,过程,模拟
From: https://www.cnblogs.com/tecdat/p/18328487

相关文章

  • Python - Plotly Express - x 轴 xticks
    我有一个代码,我想用标签作为月份和年份的名称绘制x轴。但目前它正在x轴上绘制所有时间戳,如-代码片段是:fig=px.scatter(data_frame,x=data_frame.timestamp,y=data_frame.err_codes_integer,hover_name=data_frame.col3,......
  • 8:Python字符串的魔法属性
    test="alex"v=test.capitalize()#首字母的大写print(v)test1="alEx"v1=test1.casefold()#所有字母变小写print(v1)test2="alEx"v2=test2.center(20,"*")#设置宽度,并将内容居中,20代指总长度,*代指空白位置填充且只能填一个字符print(v2)test3="alexralexr"v3=te......
  • 如何在 Python 中创建正确显示素数的代码?
    素数是只能被自身和1整除的数。例如,数字5是素数,因为它只能被1整除和5.然而,数字6不是质数,因为它可以被整除通过2和3。编写一个名为is_prime的布尔函数,它接受一个整数作为参数如果参数是素数则返回true,否则返回false。使用程序中提示用户输入数字然后输......
  • 在 MySQL Workbench 中升级 python
    MySQLWorkbench中的python版本不断受到安全标记,因为它不是最新版本。我的电脑上有最新版本,但在MySQL文件中它似乎是旧版本。我怎样才能更新,使其与我的电脑版本相同?尝试看看是否可以更新python,但没有运气很遗憾,你无法直接升级MySQLWorkbench内置的Python版本......
  • Python科研武器库 - 文件/路径操作 - 判断路径是否存在
    使用场景:在科研中,用Python处理数据的一个核心目的是批量处理,批量处理节省了研究者大量的时间和精力,不然,还不如手动一个个地去处理。批量处理通常要求数据整体较为规整,能够进行统一的处理操作,但实际数据中总存在一些不规则的样本,甚至有些样本的命名都存在不规范,例如,整个数据集......
  • python刷题常用模板
    #=====================================素数筛Begin=====================================#MAXN=1000prime=[]isprime=[True]*(MAXN+1)defeuler():isprime[1]=Falseforiinrange(2,MAXN+1):ifisprime[i]:prime.append(i)......
  • [附开题]flask框架的全国汽车销售信息查询系统的设计与实现7m1w0(python+源码)
    本系统(程序+源码+数据库+调试部署+开发环境)带论文文档1万字以上,文末可获取,系统界面在最后面。系统程序文件列表开题报告内容研究背景随着中国汽车市场的蓬勃发展,汽车品牌的日益丰富以及消费者购车需求的多样化,汽车销售信息的准确性与时效性成为了市场关注的焦点。传统汽车......
  • [附开题]flask框架的校园停车场管理系统的设计与实现61m0e(源码+论文+python)
    本系统(程序+源码+数据库+调试部署+开发环境)带论文文档1万字以上,文末可获取,系统界面在最后面。系统程序文件列表开题报告内容研究背景随着高等教育的普及和校园规模的不断扩大,校园内车辆数量急剧增加,停车难问题日益凸显。传统的人工停车场管理模式已难以满足现代校园对高效......
  • [附开题]flask框架的校园学生管理系统s8h32(源码+论文+python)
    本系统(程序+源码+数据库+调试部署+开发环境)带论文文档1万字以上,文末可获取,系统界面在最后面。系统程序文件列表开题报告内容研究背景随着教育技术的不断进步和高校规模的不断扩大,传统的学生管理方式已难以满足现代校园管理的需求。学生数量激增、课程种类繁多、选课流程复......
  • [附开题]flask框架的校园疫情管理系统92tl0(源码+论文+python)
    本系统(程序+源码+数据库+调试部署+开发环境)带论文文档1万字以上,文末可获取,系统界面在最后面。系统程序文件列表开题报告内容研究背景随着新冠疫情的持续影响,校园作为人群密集、流动性大的场所,其疫情防控工作显得尤为重要。传统的手工记录和口头报告方式已难以满足当前复杂......