关于R语言下tGarch-copula-CoVaR 系统风险程序与教学

tulipsliu

34281

收藏 2019-12-04

教学编号：2019-002课中国古语：授人以鱼不如授人以渔；

第一节

      1.1 我的Rsudio信息
   先展示我的Rsudio 版本信息，输入命令 sessionInfo()
> sessionInfo()
R version 3.6.1 (2019-07-05)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 17763)

Matrix products: default

locale:
[1] LC_COLLATE=Chinese (Simplified)_China.936
[2] LC_CTYPE=Chinese (Simplified)_China.936
[3] LC_MONETARY=Chinese (Simplified)_China.936
[4] LC_NUMERIC=C
[5] LC_TIME=Chinese (Simplified)_China.936

attached base packages:
[1] stats    graphics  grDevices utils    datasets  methods base

loaded via a namespace (and not attached):
[1] compiler_3.6.1 tools_3.6.1
>

附件列表

原图尺寸 38.56 KB

原图尺寸 42.7 KB

原图尺寸 50.78 KB

原图尺寸 50.49 KB

原图尺寸 50.91 KB

原图尺寸 62.4 KB

原图尺寸 49.88 KB

原图尺寸 90.1 KB

RCoVaRCopula-Tulipsliu.zip

大小:13.42 MB

只需: 1 个论坛币马上下载

R版本计算systemic risk 的CoVaR

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

tulipsliu

2019-12-4 11:51:45

不好意思，本来把主程序 demoCoVaR 的代码发论坛的，字数超过论坛限制的字数，最后不得不删除一半的内容。导致大家看不到程序代码。

可以下载后打开 demoCoVaR.R程序文件

对比论坛我的课程了解哪些主要行的意思。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Bright66

2019-12-4 13:46:18

tulipsliu 发表于 2019-12-4 11:15
教学编号：2019-002课中国古语：授人以鱼不如授人以渔；

前言：这个工具包比较复杂需要大量的 ...

感谢楼主分享

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tmdxyz

2019-12-5 08:54:57

感谢楼主分享

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

向阳的阳

2019-12-8 15:53:10

感谢楼主分享！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

leeyaya

2019-12-10 16:34:50

感谢指导！感谢分享！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

npyyh

2019-12-12 20:09:13

谢谢楼主分享！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

npyyh

2019-12-12 22:18:12

我想问一下模型估计出来的结果：
Family
------
No: 2
Name:  t

Parameter(s)
------------
par:  0.08  (SE = 0.02)
par2: 10  (SE = 1.59)
Dependence measures
-------------------
Kendall's tau: 0.05 (empirical = 0.05, p value < 0.01)
Upper TD:       0.01
Lower TD:       0.01

Fit statistics
--------------
logLik:  16.83
AIC: -29.67
BIC: -17.35

这里fitCopula输出的par和par2分别是t-Copula的哪个参数？还有就是这里的t-copula在CDVine包的PDF文件解释中不能进行method=‘itau’，否则会报错，为什么正常算出了结果。（除csv数据改动外，其他未进行改动。）

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tulipsliu

2019-12-14 22:53:43

npyyh 发表于 2019-12-12 22:18
我想问一下模型估计出来的结果：
Family
------

好像你加我QQ了吧？这个问题前几天加我QQ好友的人问过。

首先这里用的不是 CDVine packages ;
用的是  VineCopula package ;

方法上来说 method="mle"
很少用你说的 method="itau"

这些参考手册就可以，最近确实没时间看金融时间计量模型，在做 DSGE 模型的学习。和使用 LATEX 编辑数学公式文件。

只有 copula 的参数；不是所有 family 都有两个参数
有的只有  par1
有的有两个 par1 par2

你得下载张世英、韦艳华《Copula在金融建模上的应用》忘记书名；你自己去看；

实际上示例里我还给了一个  BB7 copula  ；在国内的教材张世英韦艳华都好像还没写。

应该是一个加拿大的学者有出书籍多元copula ；
里面就有  BB1  BB6  BB7  等 family ；

非常适合 VineCopula :
最近几天都没来论坛，确实忙其他呢。所以今天才回复。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

15736076576

2019-12-18 16:05:51

tulipsliu 发表于 2019-12-4 11:15
教学编号：2019-002课中国古语：授人以鱼不如授人以渔；

前言：这个工具包比较复杂需要大量的 ...

感谢楼主分享

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

15736076576

2019-12-18 16:07:29

tulipsliu 发表于 2019-12-4 11:15
教学编号：2019-002课中国古语：授人以鱼不如授人以渔；

前言：这个工具包比较复杂需要大量的 ...

能加下楼主QQ吗

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tulipsliu

2019-12-19 18:30:10

15736076576 发表于 2019-12-18 16:07
能加下楼主QQ吗

我是群主，你下载程序入群就加了我。

好吧，QQ 280201722 潇湘旭

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Firedragon1

2019-12-27 12:14:15

感谢楼主！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

life1sbrilliant

2020-2-5 16:40:04

感谢楼主分享

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

-->

感谢楼主！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

宁柠凝

2020-5-18 19:59:50

感谢楼主分享

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tulipsliu

2020-5-19 06:15:55

宁柠凝发表于 2020-5-18 19:59
感谢楼主分享

谢谢，但是这个代码还要升级。后面几个月可能是一定的升级。

以后也许可能会开发为R 工具包。
package.seleton('fincovar")

这个程序即使我发论坛，要能用也有点难上手的，你可以按照程序里的群号码加群。

不过群名字更改了，中国青年金融学者基地。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

翔翔星星

2020-5-31 22:32:26

这个代码里面有个SCJ copula ，但在vine copula 里找不到对应的family，怎么办呢？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tulipsliu

2020-6-1 18:21:31

翔翔星星发表于 2020-5-31 22:32
这个代码里面有个SCJ copula ，但在vine copula 里找不到对应的family，怎么办呢？

你问了同样的问题，我记得应该有的，你百度 R vinecopula ，找到 cran 网页。

下载 VineCopula.pdf 文件手册。里面的 family 就是copula 的种类。

在里面应该有。我最近忙其他的呢很久没有认真的再看这个。

***喜欢对你有帮助***

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

juefu65505

2020-10-12 23:07:26

楼主你好，请问我用原始数据运行代码第80行通过rugarch工具包的 pdist函数做均匀变换的时候报错
Error in pdist("sstd", z, mu = 0, sigma = 1, skew = coef(fittemp)["skew"], :
参数没有用(z, mu = 0, sigma = 1, skew = coef(fittemp)["skew"], shape = coef(fittemp)["shape"])
请问是哪里参数设置错误

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

juefu65505

2020-10-12 23:09:01

请问楼主运行代码第80行pdist函数做均匀变换时报错是什么原因？
Error in pdist("sstd", z, mu = 0, sigma = 1, skew = coef(fittemp)["skew"], :
参数没有用(z, mu = 0, sigma = 1, skew = coef(fittemp)["skew"], shape = coef(fittemp)["shape"])

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

whitebait301

2020-10-25 10:03:56

谢谢楼主好人

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tulipsliu

2020-10-25 10:28:51

whitebait301 发表于 2020-10-25 10:03
谢谢楼主好人

我只不过增加了作者程序里没提供的 rugarch 工具包估计的 t-garch 模型代码。
作者的代码在 www.github.com 。原作者更是好人，是一个意大利学者。哈哈

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tulipsliu

2020-10-25 10:35:53

juefu65505 发表于 2020-10-12 23:09
请问楼主运行代码第80行pdist函数做均匀变换时报错是什么原因？
Error in pdist("sstd", z, mu = 0, sigma ...

非常抱歉，最近几个月非常忙，我都没看到你的提问；
这里应该是你的 garch 模型边缘分布拟合不好。

这里要看数据，作者说自己的模型符合 TGARCH.
实际应用中，比如你可以看 rugarch 工具包的PDF 文件。
我给你一个链接地址：
https://cran.r-project.org/web/packages/rugarch/
这里可以下载作者的帮助文件手册PDF 文件。

其实这个工具包可以拟合很多种类的GARCH，比如 AGARCH ，GJR,EGARCH。
这里具体到你的数据，应该是不符合 TGARCH 模型。因为拟合的数据在去除条件异方差，也就是要做 PIT 概率积分变换，获得 copula 建模必须有的无条件方差数据。这个程序无法拟合出这样的边缘分布数据，所以在进行 PIT 概率积分变换转换为均匀分布的 U[0,1] 分布数据的时候。会出现你提到的这个错误。

这个地方就是非常重要的地方，也是金融学硕士、博士研究生在拟合这个GARCH-COPULA 建模的时候经常遇到的问题，这里的功底在于理解和加深对 GARCH 模型建模的技巧，模型选择的评估手段和方法。

实际上，一个解决方案是：
第一， TGARCH 模型边缘分布拟合的时候判断是否是 TGARCH(1,1) 好，还是 TGARCH(2,1) ,OR TGARCH(1,2), TGARCH(2,3) 等不同阶数的TGARCH 模型是否更适合自己的数据；
第二，不用TGACH ，测试一下 GJR 或者 EGARCH 是否更好？
我将把你的提问置顶为优秀回帖，以提示后面的建模者，在边缘分布GARCH 模型建模上对细节的把握和谨慎的态度。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

小hr

2020-11-13 22:50:42

楼主你好，我研究的是Vine-Copula，数据没有自相关性和ARCH效应，在做边缘分布那还可以用ARMA(0，1)-GARCH(1，1)拟合吗？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

黄欣文

2021-1-14 10:59:16

感谢楼主分享！楼主好人！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

YxyXM

2021-2-22 00:53:22

楼主，想请问一下，程序中没有时变量，为什么CoVaR计算出来后是动态的呢，不应该是5%分位数下的一个值吗

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

余博士伦

2021-2-26 17:10:03

太贵了，压缩文件下不起了，能发邮箱吗？564279942@qq.com，谢谢！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tulipsliu

2021-2-26 20:13:11

小hr 发表于 2020-11-13 22:50
楼主你好，我研究的是Vine-Copula，数据没有自相关性和ARCH效应，在做边缘分布那还可以用ARMA(0，1)-GARCH( ...

不可以；今年国际上 SCI 撤了很多国内的论文；
数据修改过；

有 ARCH 效应你就建立 GARCH 模型，没有就做其他模型；

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tulipsliu

2021-2-26 20:17:04

YxyXM 发表于 2021-2-22 00:53
楼主，想请问一下，程序中没有时变量，为什么CoVaR计算出来后是动态的呢，不应该是5%分位数下的一个值吗

这个代码是国际上的，是国际上一个学者的全英文论文的代码；

是动态的，我知道你的疑惑；

国内的前几年做出来过，算出来的是一个值；

MATLAB 官方的某个 VaR ， value at risk 也是一个值；

但是国际上另外非官方的 MATLAB 程序计算出来的也是动态 VaR ，是和收益率的时间一样长的数值。

计算的方法和逻辑都不一样，这里是你们没理解而已。

我知道国内计算的 CoVaR 是一个数值居多；甚至很多是公开发表的博士论文；博士毕业论文。

这里不但计算出来了动态的 CoVaR ，动态的 VaR 也计算出来了；

如果已经下载了这个资料，程序文件夹里： demoCoVaR.R 程序文件里，开头就给了英文论文的名字；

最近很少看这个帖子，不知道已经有留言。今天才回复。