浅谈基于Copula

【摘要】由于许多金融资产收益率的波动可能存在非对称性,本文通过对收益率采用garch或者egarch建模,对存在杠杠效应的金融资产采用egarch建模。并将copula函数和monte carlo技术应用于融通深证100基金,计算得到其前十大重仓的股票及其投资组合的风险值var以及es,结果表明,基于copula-es的风险度量方法可以为我国基金管理公司评估和管理市场风险,从而控制和减少资产损失提供参考。

1.引言

近年来金融工具及其衍生物越来越多元化,其带来的不确定因素也越来越大,因而金融市场的风险也就越来越高。金融市场间的关系更是变得日趋复杂,更多的呈现出非线性、非对称和厚尾的特性,金融波动和危机的频繁出现使聚合风险管理和金融市场间相依关系分析成为国内外关注的焦点。

现阶段最常使用的风险度量指标是20世纪90年代j.p.morgan和g30[1]集团提出的var(value at risk)方法,var旨在一定的置信水平下,估计金融资产或组合受市场因子波动影响,而在未来特定的一段时间内的最大可能损失。严格的说,var描述了在一定的目标期间内收益和损失的预期分布的分位数。如果c代表置信水平,var对应的是较低的尾部水平1-c。可表示为:

其中,表示某事件的概率,表示资产或资产组合在持有期内的损失,c为置信度水平。在最近这些年var作为金融风险度量工具得到了广泛的应用,然而,研究发现var不具有次可加性和一致性的风险度量,后来针对这一问题,acerbi等[2-3]提出了期望损失es(expected shortfall)的定义。

假设r为持有期内资产或者资产组合的损益,并假设r的累积分布函数f(r)(cdf)是连续的,那么对于置信水平,var也可以用如下定义:

式中,表示r的分布在给定显著水平的下侧分位数。假设表示r的概率密度函数,那么置信水平为1-c下的es可以定义为:

式中,为示性函数。es实质上是将资产价值r乘以权重的从-到0的积分,这样它就把超过var水平的损失部分考虑进去了。从经济意义上讲,es就是指当损失超过var时的平均损失。由于它同时具有了次可加性和一致性,是一个较好的风险度量工具。rappoport(1993)[4]第一次在金融行业中用它来做风险管理的一个工具,后来acerbi(1997)[3]等人证明了该方法是一个一致性的风险度量工具。

同时,通过引入copula函数度量资产组合集成风险的方法已经越来越成熟。schweizer和sklar最早提出copula函数的概念及其它的一些性质。后来sklar指出了copula函数可以把具有不同类型边缘分布函数连接起来,并且能抓住它在风险管理应用中的本质特征(例如:尖峰厚尾性);emberchts等第一次把该方法引入到金融类相关研究之中。许多研究学者在他们的基础上做了很多有意义的研究。例如:breymann等人研究表明了学生t-copula的经验拟合比高斯copula优越很多;ceske,hemandez(1999)提出可以将copula函数与montecarlo技术结合计算相关损失;matteis对archimedean copula做了很好的总结。

在我国,copula函数方法在金融上的应用才刚刚起步,且其中绝大多数文献做的是介绍性、引入性的研究。最早见的是张尧庭(2002)提出copula函数在金融风险领域大有可为;史道济利用copula函数研究外汇组合的相关性;司继文(2004,2005)分别将copula函数应用于国内外的股票市场和期货市场;韦艳华、张世英(2004)将garch模型应用于copula函数,来度量金融时间序列的自相关结构。前人的研究主要集中在利用copula函数对股市或资产组合的相关性研究。而韦艳华(2004)利用garch模型拟合正态copula函数的边缘分布,然后运用monte carlo仿真技术计算投资组合的var。

本文创新一是采用garch或者egarch模型来拟合t-copula函数的边缘分布,克服了传统garch模型不能处理特定非对称金融时间序列的局限性。对此,本文也比较分析了单独使用garch下和本文采取的方法下的风险值,研究表明本文提出改进的思路对风险估计更为准确。改进二在于对于风险衡量的指标不是仅仅采用var,而是利用var与es双监管的风险度量方法,克服了传统风险度量var的很多缺陷(不具有次可加性、正齐次性等)。最后通过度量我国股票型开放式基金的市场风险为出发点,运用monte carlo仿真技术计算投资组合的var以及es,具有很强的实用价值。

2.garch或者egarch收益率分布模型

对于某一金融资产,投资者最想知道的是将来某个时刻该资产收益率的信息。由于金融资产收益率的尖峰厚尾性、条件异方差性、波动聚集性和杠杆效应等,普通的garch模型对对前三者能较好的刻画,但是对于杠杠效应garch模型不能刻画出,因此,本文对不存在杠杆效应的收益率序列采用garch模型拟合资产收益率的特征,对存在杠杠效应的收益率序列采用egarc模型刻画。该模型是glosten jagannathan和runkle在engle提出自回归条件异方差(arch)模型和bofloerselev提出的广义自回归条件异方差(garch)模型的基础上改进提出的,它考虑了坏消息和好消息对波动性不同的影响。

假设投资组合中有d种金融资产,对于资产i,直接根据最近的n期历史收益率数据(t=1,2,…,n)运用garch其中eagch模型中条件方差采用自然对数形式,意味着非负,且杠杠效应为指数型的。模型中引入了一个重要参数,若,说明信息作用非对称,存在杠杠效应。为第i个资产收益率序列;为的条件均值项;v是t分布的自由度。为待估参数;随机扰动项服从均值为0,方差为1的标准t分布。即:

它的形式使得garch或者egarch模型能够较好地描述收益率序列的各种特性(如条件异方差性、波动聚集性)。

假定利用观察资产收益率历史数据样本可以得到,在估计出参数后,可以得到下一时刻收益率的条件分布:

其中,是自由度为v的t分布函数,是到时刻t为止的信息集。利用matlab可以很方便地由样本观测数据估计出garch模型的各个参数,从而由式(2)得到给定t时刻前信息集的条件概率分布。

3.copula函数的参数估计

在资产分配、衍生品定价和风险管理等许多金融领域中,有关风险因素之间的相关性度量及其重要。许多文献中常采用多元正态函数,然而各种金融资产的收益率之间一般并不符合多元正态分布的假设,为此,本文使用copula函数来解决这个问题。

由sklar定理可知,对于一个具有边际分布函数为()的金融资产的联合分布函数f,一定能找到一个copula函数c,使得:

如果所有的边际分布函数都连续则从上式定义的copula函数是唯一的。从上式可以计算得出copula:

其中,

文献表明,t-copula能更好地刻画各个金融资产的尾部相关性,本文研究的是t-copula连接函数分布。

其中,表示相关系数矩阵为r,自由度为v的维标准t分布函数,表示自由度为v的单变量标准t分布函数的反函数。copula模型为:

参数v为t分布的自由度。为d维t-copula分布,表示均值为0,方差为1,自由度参数为的正规化t分布函数,即:

式中是伽马函数。由于t-copula的密度函数对任意维数都不是一个简单的形式,本文根据t-copula函数形式使用matlab工具估计其参数,过程如下:

(1)把资产收益率数据通过概率积分变换转化为一致分布;

(2)用对数似然函数法估计学生t-copula的参数:

此处的copula函数c为公式(4)给出的;

(3)令,此处是单变量累积标准正态分布函数;

4.利用模特卡罗模拟资产组合的var和es

根据embrechts关于利用t-copula函数模拟随机变量的方法,多次模拟资产组合资产收益率的随机扰动项。具体模拟步骤为:

(1)由上述估计出的随机扰动项的相关系数矩阵r,对其进行cholesky分解。

(2)根据标准正态分布,模拟d个相互***的随机变量,。

(3)产生与y相互***的变量e,服从分布。

(4)令。

(5)令,则x为服从自由度为v的t分布。

(6)计算得到。

(7)根据得到联合分布为,连接函数为的d维随机扰动项。

(8)根据egarch或garch模型,得到金融资产收益率的条件均值和条件方差,然后根据随机波动方程,得到资产组合的资产收益率向量。

(9)给定资产在投资组合中的权重,计算投资组合收益r的值。

(10)重复上述过程5000次,模拟得到其经验分布,容易求出var和es的值。

5.实证研究

(1)数据的选取和边缘分布的估计

本文选取融通深证100基金2010年3月31日公布的前10大重仓股票作为观测样本,如表1所示。

本文采用从2004年7月1日至20010年3月31日共937个(对空缺数据已做处理)交易日的收盘价作为原始数据,计算出每只股票的对数收益率,并根据公式(1)、(2)估计出每只股票的边际分布,利用aic准则、sc准则和杠杆系数检验可以得出万科a、苏宁电器、中兴通讯、格力电器、吉林敖东存在杠杠效应,适宜采用egarch模型建模(其余采用garch建模拟合效果更好),下面以第三、四只股票五粮液(000858)和苏宁电器(002024)为例,分析其边缘分布函数的估计和拟合效果的评价。参数估计结果见表2和3:

其中,表2、3中括号的数据表示相应t统计量的估计值。从表2、3中给出的k-s相伴概率可知,利用cml方法对原序列做概率积分变换后,序列服从[0,1]的均匀分布。由此可以说明本文提出的模型可以较好地描述相关资产的边缘分布。同样的方法检验了其他8只股票的边缘分布拟合效果,都说明了garch或者egarch模型能较好拟合各自的边缘分布,因此用本文的模型描述收益率序列的边缘分布是充分的。

(2)copula函数参数估计和monte carlo模拟var和es

此处为了对比分析采用egarch或garch拟合边缘分布与仅仅采用garch拟合的效果,根据上面估计得各个股票收益序列的边缘分布,利用文中第四部分的估计copula函数参数的方法,估计得出其t分布的自由度dof=7.5848和各自的相关矩阵(表4、表5)。

由上述结果可知,仅仅采用garch拟合边缘分布使得各个资产之间的相关系数整体性的变小。从而可以推断出可能会导致风险的低估,从而对准确度量基金风险存在一定的影响。进一步的风险值比较分析可见表6、表7。

按照表1的投资比例,假设投资者处于t时刻,这里的t时刻指的是样本时间段的最后一天,即2010年3月31日,t时刻的投资组合价值为:

假设资产持有期从t时刻到t+1时刻,置信水平选择95%和99%进行估计。根据本文的copula-garch(egarch)模型,运用t-copula函数的模特卡罗模拟仿真模拟5000次,可以得出t+1时刻各个股票的收益率序列,进一步可以计算出t+1时刻各股票的损失序列,给定置信水平,容易得出t时刻到t+1时刻相应的投资组合var和es值,表6、7分别给出了本文提出的采用egarch或garch拟合边缘分布和传统方法仅仅采用garch模型拟合的情况下各个股票和证券投资组合的var和es值。

表6、表7清晰的显示本文所提出的方法对准确估计风险更进了一步,比较而言传统方法只使用garch模型拟合边缘分布导致了单个资产和资产组合的风险值都偏小。对于文中提出使用egarch模型拟合万科a、苏宁电器、中兴通讯、格力电器、吉林敖东,从表6与7中可以看出,var与es的风险值都比其他个股风险值偏离程度更大,说明采用egarch模型针对特定(存在杠杠效应)金融序列拟合效果更好,风险估计更为准确。

单独分析表6可以看出,在投资额一定的情况下,基金的风险值要比单个资产进行投资的风险值小,可见该基金选择的各个股票之间的相关关系有较大差别,说明投资组合可以大大降低投资风险。从var和es的风险值看,es都大于var,说明es比var度量风险更为保守,也说明了var在度量风险上存在计算风险值偏低的现象,这样对基金控制风险和减少资产损失极为不利,特别是当极端事件发生时,资产组合的风险值与实际值就会发生偏差。另外从var与es的差值可出看出,置信水平越高,投资组合降低风险的程度也就越大,但是由于var不具有次可加性,从es的差值能很明显看出。

6.结论

本文为了描述特定资产具有非对称性的特征,通过对资产收益率进行egarch建模,对不存在杠杠效应的资产仍使用传统的garch模型,这与copula可以连接具有不同边际分布的函数的相关关系相符,同时考虑到var度量风险的不足,引入了es一致性风险度量,通过t-copula函数和monte carlo模拟计算出了证券投资组合的var以及es的值。最后文章对融通深证100基金风险度量的实证研究可以看出,有的金融资产收益率序列可能存在杠杠效应;而且var的确存在低估风险的不足;同时也得出了风险值var或者es在置信度越高,它们的差值越为明显,说明了本文copula-egarch(garch)模型能较好地刻画投资组合二中不同资产间非正态非线性非对称的相关性。基于该模型的风险度量方法可以为我们基金管理公司评估和管理资产组合的市场风险,从而有利于公司控制和减少资产损失提供一定的参考作用。

参考文献

[1]张尧庭.连接函数(copula)技术与金融风险分析明[j].统计研究,2002(4):48-51.

[2]史道济,邱男.关于外汇组合风险相关性的分析[j].系统工程,2005,23(6):90-94.

[3]司继文,蒙坚玲,龚朴.国内外期货市场相关性研究田[j].华中科技大学学报,2004,21(1):16-19.

[4]韦艳华,张世英.金融市场的相关性分析copula-garch模型及其应用[j].系统工程,2004(4):7-12.

资助项目:国家自然科学基金项目(71071071,70871058);***人文社会科学研究规划项目(09yja790100);清华大学经济管理学院

转载请注明出处学文网 » 浅谈基于Copula

学习

小小棋童成长阶梯

阅读(15)

学棋,能促进孩子的智力发展,更能培养孩子的意志品质。PART1为什么要学棋不是每个孩子都喜欢下棋,尤其是女孩,他们一开始常常是被动地接触棋子,后来才慢慢喜欢上的。关于“为什么学棋”,我们常常可以听到三种不同的声音。孩子说:“就是喜欢嘛!可

学习

中西古典首饰对比

阅读(19)

本文为您介绍中西古典首饰对比,内容包括新古典主义首饰特征,中国古典珠宝首饰设计。首饰就像天上的星星一样,给人们带来闪耀和美丽的同时,所映衬的是当时社会背景和文化。谈起中国首饰是离不开玉文化的,自古有:国无玉不盛,家无玉不富,人

学习

这个世界有没有圣诞老人

阅读(23)

早晨,阳光很好,我开始收拾家务,用一块干净的毛巾擦拭电脑,中断了妞在电脑上看的动画片。她要求在我的笔记本上接着看,我答应了。她开机,浏览,忽然沉默了,抬起头望着我:“妈妈,你说,这世界上根本没有圣诞老人是不是?每年的圣诞老人都是你装的。”我一

学习

花盖梨优质丰产栽培关键措施

阅读(23)

本文为您介绍花盖梨优质丰产栽培关键措施,内容包括花盖梨功效和用处以及禁忌,东北花盖梨包邮多少钱一斤。通过介绍建园、高光效树形、科学施肥、灌水前促后控、合理控制负载量、植物生长调节剂应用、病虫害防治、适时采收等技术,总结了花

学习

文言虚词初探

阅读(26)

文言虚词是相对文言实词而言的,虽不表示实际意义,但在表达上却各自具有不同的作用。从某种意义上说,不能掌握古代汉语的虚词,就很难读懂古文。所以,我们要重视文言虚词的学习,并根据不同的语境,分析不同虚词的意义和用法,这样不仅能在考试中游刃

学习

阻击因病致贫

阅读(23)

“有什么别有病”“病来如山倒”,这些老百姓耳熟能详的俗语,从一个侧面证明了疾病尤其是大病重病的可怕。一场大病不啻为一场灾难,沉重的经济负担压垮了不少家庭,使得因病致贫、因病返贫屡屡上演。重大疾病已经成为横亘在贫困人口脱贫路上最

学习

流星飞镖——咻作文600字

阅读(22)

在我家的装饰架上,摆着一个小和尚,它们是爸爸从四川给我买回来的,我特别喜欢它。这个小和尚拿着一个流星镖,只有拳头大小,正在练功呢!真是栩栩如生啊。只见这个小和尚脑袋剃的光溜溜的,摸上去非常光滑。它上身的袈裟是橙黄色的,下身穿着红色的练

学习

外地人容易误解的扬州方言特色词语研究

阅读(20)

摘要:扬州方言是扬州地方文化的一个重要标志和载体,在扬州地域文化传承中起着重要的作用。但是,扬州方言和外来人员所讲语言的差异,造成了扬州本地人和来扬人员交流的不便与误解。在此背景下,如何在保护和保持扬州方言的前提下,拉近扬州

学习

迟子建《世界上所有的夜晚》中的死亡主题解读

阅读(28)

本文为您介绍迟子建《世界上所有的夜晚》中的死亡主题解读,内容包括迟子建写世界上所有的夜晚,迟子建作品世界上所有的夜晚。《世界上所有的夜晚》是当代女作家迟子建的代表作品,这部小说与作家带有符号化的北极村温馨的童话故事明显不同

学习

郑州古玩市场激战在所难免

阅读(21)

中原古玩城、天下收藏、郑州珍宝城、碧波园珍宝大世界……一时间,郑州新建古玩市场遍地开花,而且大多是2万平方米以上的大店。加上郑州古玩城、北古玩城,郑州的大型古玩城不下6座。郑州古玩市场的容量究竟有多大?古玩行当新兵老将大打商战,最

学习

汇率风险管理

阅读(86)

本文为您介绍汇率风险管理,内容包括汇率风险管理工作成效,汇率风险管理工作情况的报告。一、保险公司汇率风险的类型及特点汇率风险,又称外汇风险,是指经济主体在涉外业务中因汇率波动而蒙受损失的可能性。保险公司所面临的汇率风险通常是

学习

高中化学教学中关于原电池与电解池的综合分析

阅读(27)

原电池和电解池是氧化还原反应原理的具体应用,探究电极反应发生的原因,是理解及应用原电池、电解池的前提和关键。因此,在分析、研究原电池或电解池的问题时,要紧紧抓住氧化还原反应这一主线索。关键词:原电池电解池电极电极反应综合应用一、

学习

湾流G280体验飞行

阅读(23)

本文为您介绍湾流G280体验飞行,内容包括湾流g280一年养护费用,私人飞机湾流g280价格。G280公务机是美国湾流公司与以色列航宇工业公司共同研制的最新产品,与同级别的其它公务机相比,拥有舒适的客舱、较大的航程、较高的速度、更强的动力等

学习

浅谈对蓝海战略的认识

阅读(20)

【摘要】蓝海战略一经提出,在世界范围内获得了很大的反响。本文通过对蓝海战略内容的介绍,提出对蓝海战略的一个理性的认识,在此基础上得出企业正确实施蓝海战略的方法。【关键词】蓝海战略;认识;实施一、蓝海战略的主要内容(1)蓝海战略的内涵

学习

优降糖致低血糖症临床分析

阅读(22)

2型糖尿病的发病率在全球呈上升趋势,我国的糖尿病患者已达4000万人。由于经济状况等原因,优降糖在我国仍是使用最广泛的降糖药,从而导致低血糖症的发生率居高不下。优降糖所致的低血糖症具有表现多样化,持续时间长,反复性强的特点,重者

学习

ITH钻孔技术应用分析

阅读(23)

本文通过普通凿岩机与ITH凿岩设备钻孔方式的应用对比,分析了ITH钻孔技术在凿岩钻孔应用中的优点和缺点,从而为操作使用者选择使用普通凿岩设备或ITH凿岩设备提供应用指导依据。关键词:ITH;钻孔技术;应用分析InTheHole(简称ITH)钻孔技术将钻杆、

学习

歌剧咏叹调《林间小鸟》演唱分析

阅读(34)

奥芬巴赫一生的轻歌剧创作中仅仅留下一部正歌剧《霍夫曼故事》,这也是奥芬巴赫获得极大成功的一部歌剧。“霍夫曼”讲述了三段悲伤的爱情故事,其中由女主人公演唱的这三首咏叹调是全剧的经典,《林间小鸟》则是其中最具趣味性的咏叹调。关键

学习

对建筑节能设计中的遮阳技术分析

阅读(25)

本文通过对建筑遮阳的技术原理、遮阳系统的界定及分类状况的分析介绍,并结合目前国内外遮阳技术的研究状况,详细阐述了遮阳设计中应该考虑的因素与方法。关健词:建筑节能设计;遮阳技术;原理;系统界定;分类状况1建筑遮阳原理夏季南向窗接受太阳

学习

板式换热器不锈钢板片腐蚀失效原因分析

阅读(28)

【摘要】板式换热器凭借着诸多优势,成为了食品、冶金、石油化工等领域的主导换热设备。但由于换热器的工作环境复杂,板式换热器容易出现因腐蚀而导致设备失效。本文结合应用实例,就板式换热器不锈钢板片腐蚀失效原因进行分析,可为此类故障防

学习

跨文化视角下中国人文纪录片的符号分析

阅读(32)

【内容摘要】本文试从文化折扣、高低级语境以及文化中心论的跨文化视角,选取中央电视台出品的《故宫》以及美国国家地理频道的《故宫》国际精编版——《解密紫禁城》进行比较,从语言符号、画面符号以及文化符号等三个具体层面,分析中国人文

学习

《红楼梦》中贾迎春角色分析

阅读(24)

本文为您介绍《红楼梦》中贾迎春角色分析,内容包括红楼梦中贾迎春的特点,红楼梦中贾迎春的结局。对文学材料的内部人物进行分析是文学评论的一个重要部分,也是业内学者经常进行的一项工作,本文作者将针对《红楼梦》中的贾迎春这一角色进行

学习

重庆东水门大桥主桥设计计算分析

阅读(29)

本文为您介绍重庆东水门大桥主桥设计计算分析,内容包括重庆东水门大桥设计说明,如何区别东水门大桥和千厮门大桥。摘要:以重庆东水门大桥为例,通过电算与手算相结合来复核斜拉桥的受力是否合理。本文通过电算方法,来验算斜拉桥特殊截面的