极值理论在股指期货 保证金设定中的应用
极值范围等问题时
详细描述
单调性分析法主要依赖于函数的导数。当函 数在某点的导数大于0时，函数在该点附近 单调递增；当导数小于0时，函数单调递减。 因此，导数由正变负或由负变正的点，即为 函数的极值点。
二次导数法
总结词
通过分析函数二次导数的符号，可以判断一 阶导数变号的个数，从而确定极值的个数及 范围。
详细描述
如果一个函数的二次导数在某区间内大于0， 那么一阶导数在这个区间内只变号一次，即
经验法
根据实际问题的背景和性质，结合经验判断 极值点。
02 极值的应用场景
经济学中的极值应用
投资组合优化
在经济学中，极值理论可以用于投资组合优化，通过确定投资组合 的风险和回报的极值范围，帮助投资者做出更明智的决策。
风险评估
在金融领域，极值理论可以用于评估极端市场情况下的风险，例如 股票价格暴涨或暴跌的可能性。
极值的类型
极大值
在某区间内，比它大的函数值只有一 点。
极小值
在某区间内，比它小的函数值只有一 点。
极值的判定方法
一阶导数判定法
如果一阶导数在某点的左右两侧由正变负或 由负变正，那么这一点就是极值点。
二阶导数判定法
如果二阶导数在某点的左右两侧由负变正或 由正变负，那么这一点就是极值点。
表格法
通过比较函数值的变化情况，确定极值点。
约束条件下的最优化问题
约束条件下的最优化问题
在某些约束条件下，寻找函数的最优解。约束条件可以是等 式或不等式形式。
约束条件下的最优解
满足约束条件的函数最优解可能是无界的，或者是局部最优 解。需要根据问题的具体情况进行分析和求解。
04 极值的计算方法
导数法
总结词
导数法是求极值最常用和最基本的方法，通过求导数并令其为零，可以找到可能的极值点。
期货投资组合交易保证金设置研究——基于Copula模型和极值理论的分析
由于 Ｃ ｐ ｌ 型 具 有 允 许 分 步 建 模 的优 点 ，因 ｏ ｕａ模 此 就 可 以 利 用 这 一 优 点 对 模 型 进 行 分 步 估 计 ，对 数 据进行 拟合 。通常 Ｃ ｐ ｌ ｏ ｕ ａ模 型 可 以运 用 两 步 极 大 似 然 估 计 法 进 行 估 计 ， 即 首 先 利 用 极 大 似 然 估 计 法 对 模 型 的 边 缘 分 布 进 行 参 数 估 计 ，然 后 再 利 用 该 方 法 对相 应 的 Ｃ ｐ ｌ ｏ ｕａ模 型 的 似 然 函 数 进 行 极 大 化 估 计 ，
正 态 Ｃ ｐ ｌ 、ｔ ｐ ｌ ｏ ｕａ － ｕａ和 阿 基 米 德 类 Ｃｏ ｕａ函 数 等 Ｃｏ ｐｌ
历史 样本 数 据容 量过 小 ，这会 导 致估计 偏 差过 大 。 上 述 研 究 在 取 得 大 量 成 果 的 同 时也 存 在 如 下 几 点
共 同 的不 足 之 处 。首 先 ，上 述 研 究 普 遍 运 用 非 参 数 的 核 估 计 方 法 对 未 来 的 对 数 收 益 率 进 行 预 测 ，但 是 该 方
标 评 价 并 选 取 模 型 。 最 后 结 合 Ｍｏ ｔ ａｌ 拟 技 术 ｎｅＣ ｒｏ模
和 极 值 理 论 方 法 来 估 计 投 资 组 合 的风 险 进 而 计 算 期 货 投资组合的保证金。
三 、 Ｃｏ ｕａ函 数 理 论 ｐｌ
（ ）Ｃｏ ｕ ａ函 数 定 义 。 一 ｐｌ
ｒ ＾
关 检 测方 法 对 Ｃｏ ｕａ模 型 的拟 合 优 度 进 行 检 验 。 ｐｌ
第 一 ，对 变 量 的 边 缘 分 布 和 连 接 它 们 的 Ｃ ｐ ｌ ｏ ｕａ
极值理论在金融风险度量中的应用
极值理论在金融风险度量中的应用作者：董兴志来源：《时代金融》2013年第14期【摘要】如何准确度量金融市场风险在金融风险管理中扮演着重要角色，而极值理论能对极端风险进行较准确的度量。
通过对沪深300指数的实证研究表明，广义帕累托分布能够很好地拟合极端日收益率数据，从而用建立在广义帕累托分布基础上的POT模型来度量投资者所面临的市场风险是合适的。
结果显示：基于正态分布假设得到的风险值小于POT模型的风险值。
【关键词】风险度量极值理论 POT模型广义帕累托分布一、引言金融风险管理的目标之一是对潜在巨大损失发生的规模和可能性大小进行准确度量。
从统计学角度来讲，巨大损失发生的规模和可能性大小分别对应着损失分布的高分位数和尾部概率。
一些参数方法和非参数方法能在有很多观测值的基础上很好地拟合经验损失分布，但对尾部的拟合效果很差，从而这些方法并不能满足风险管理者对超出观测的极端风险度量的需求。
而极值理论（Extreme Value Theory）很好地解决了这一问题，该理论是对尾部建模的一种方法，关注的是极端值而非全部数据，也就是不研究序列的整体分布情况，只关心序列的极端值分布情况，因而能对极端风险进行较准确的度量。
根据确认极值的方法不同，极值理论存在两种模型：第一种方法考虑在连续周期内的最大值，建立在此极值基础上的模型称为BMM（Block Maxima）模型，它利用广义极值分布来逼近损失分布的尾部情况；第二种方法考虑超过某一给定阈值的观察值，建立在此极值基础上的模型称为POT（Peak Over Threshold）模型，它利用广义帕累托分布（Generalized Pareto Distribution，简称GPD）来逼近损失分布的尾部情况。
由于POT模型能够更有效地使用数据、很方便地计算出在险价值（Value-at-Risk，简称VaR）和预期不足（Expected Shortfall，简称ES），因此POT模型已经成为实证研究中的首选模型。
基于极值理论方法的中国股指期货保证金设定的实证研究
第 ２ 卷第 ４期 １
２０ ｏ６年 ７月
统 计 与 信 息 论 坛
Ｓｔｔ ｔ ｓ＆ Ｉ ｆｒ ｔ ｎ Ｆ ｒ ｍ ａｉ ｉ ｓｃ ｎ ｏ ｍａ ｉ ｏ ｕ ｏ
Ｖ０ ． １ Ｎｏ ４ １２ ．
Ｊ ｌ， ０ ６ ｕｙ ２ ０
【 统计应用研究】
基于极值理论方法 的中国股指期货保证 金 设 定的实证研究
李晓渝 ， 宋 曦， 潘席龙
６０１） １０ ７ （ 西南财 经大学 中国金融研究 中心 ， ｔＩ成都 ｉ／ ￣ Ｉ
摘要： 保证金制度 的合理设定是确保股指期货 交易高效和安全进行 的保证 。文章在极值理论广义帕 累托
分布 （ Ｐ 下 ， Ｇ Ｄ） 研究 了上证 以 １０指数 为标 的的股指期货 的保证金 水平 ， 且与正态分 布假设下 的保证 金水 ８ 并 平进行 比较。建议以 Ｇ Ｄ下的损 失期 望值作为确定 中国股 指期货 保证金水平 的依据 ， Ｐ 并且 可以考虑 区分不 同头 寸设 置不同的保证金水平 ， 中国的股 指期 货交易 的保证金设定提供参 考。 为 关键词 ： 指数期货 ； 极值理论 ； 保证金水平 ； 损失期望值
（ ｎ ｒｌｅ ａｅ ｉｒ ｕｉ 下 的 在 险 价 值 Ｇｅｅａｚｄ Ｐ ｒｔ Ｄｓ ｉ ｔ ｎ） ｉ ｏ ｔｂ ｏ （ ｌｅａ Ｒｓ） 损 失 期 望 值 （ ｘｅｔｄＳ ｏｔ ｌ Ｖａｕ ｔ ｉ 和 ｋ Ｅ ｐｃｅ ｈｒａ１ ｆ ）
生产品的经验来 看 ， 股指期货是新兴市场开设金融
衍生 品交 易 的首 选 品 种 。我 国 股票 市 场 的发 展 、 机
构投资者的逐步壮大、 监管体 系和期货 市场的完善 都为股指期货开展提供 了良好 的基础。另外 ， 包括 上 证 １０指 数 、 证 ５ 数在 内 的各 种 指数 已经 或 ８ 上 ０指
湘财祈年股指期货专题二：股指期货保证金设置方法的比较
湘财祈年股指期货专题二：股指期货保证金设置方法的比较湘财祈祷研究员：刘正珺摘要：目前，中国金融期货交易所（简称：中金所）公布的沪深300股指期货保证金制度中最低保证金为15%,这个保证金水平较为保守，但在股指期货上市初期还是合理的。
因为，在股指期货推出之初，期货市场价格波动可能相对剧烈，那么相对保守的保证金能有效的防止违约风险，有助市场的健康、稳健的发展。
随着我国股指期货市场的发展，可以借鉴国外成熟市场的保证金计算方式，在不大量增加违约风险的条件下，适当降低保证金水平，可以增加市场的流动性。
本文介绍国外成熟的保证金制度中保证金设置方法，运用EWMA-VaR、GARCH-VaR 和极值理论-VaR方法，计算沪深300股指期货交易保证金水平。
●保证金制度意义：期货市场的发展需要配套的风险管理制度，保证金制度是期货是期货市场风险管理制度的核心之一，保证金水平的合理设置，能够增强市场抵御风险的能力，同时可以提高市场运作的效率。
保证金水平的设置需要在市场风险与上次流动性之间进行权衡。
●主要结论：➢在1%的显著性水平下，通过GARCH-VaR求得的保证金水平的最大值、最小值和均值均高于通过EWMA-VaR得到保证金水平的相应值。
由于EWMA-VaR方法采用固定的衰减因子，我们选用香港期货交易所对日数据选取的衰减因子0.94。
➢本文中，EWMA-VaR及GARCH-VaR对分布的假设分别基于正态和t分布。
然而，进行风险估算时，往往只需考察金融资产收益的尾部风险。
相对而言，极值理论-VaR更好的考虑到分布的尾部，得到的保证金水平比上述两种方法更加稳健。
一、国际成熟的保证金管理介绍：SPAN系统：SPAN（Standard Portfolio Analysis of Risk）是目前被广泛应用的保证金计算系统——一个基于投资组合的保证金计算与风险评估系统，它是美国金融市场对1987年的股灾反思之后，由芝加哥商业交易所(CME)根据总统顾问小组提出的加强风险控制的建议于1988年12月16日设计推出的，其核心理念与VaR具有异曲同工之妙。
基于Poisson-GP复合极值分布的股指期货保证金设置研究
律 为 Ｐ｛ ＝ｋ ｎ ｝＝Ｐ ， ＝０ １ … ， ， ， 与之 独立 的连 续
型 的建 立对 我 国股指 期货 保证 金设 置研 究有 一定 的 理 论参 考价 值 ．
目前 我 国 设 置 保 证 金 水 平 的 方 法 主 要 分 为 两 类： 第一 类 是基 于 Ｅ ＷＭＡ模 型 、 Ａ Ｃ Ｇ Ｒ Ｈ模 型 的保 证 金 水平 设置 方法 ． 类 方 法 以金 融 收益 率序 列 服 从 该
正 态分 布为 前提 ， 略 了金 融 数据 的尖 峰厚 尾性 ， 忽 存 在 低估 风 险 的可能 ． 二 类 是基 于极 值 理 论 的保 证 第
金 水平设 置 方法 ． 值理论 着 重考 虑分 布 的尾部 ， 极 能
收 稿 日期 ：０ １０ －５ ２ １ － ２ ４
作 者 简 介 ： 孜 文 ， ，９ ４ ， 士 生 ； 究 方 向 ： 货 风 险 控 制 ； — ｉ：ｏ ｏｌ 一 ｚ １３ ｃｎ 叶 男 １８ 一硕 研 期 Ｅｍａ ｎｎｎ ３ ｙｗ＠ ６ ．ｏ ｌ
卖 双方 根据 事先 约定 好 的价格 同意在未 来某 一 特定
Байду номын сангаас
祥 、 泽智 用极 值理 论 研 究 了以全 国统 一 ３０指 数 吴 ０
为标 的物 的指数 期 货 的保 证 金 水 平 ｊ结 果 表 明在 ，
违 约概 率为 １ 的情 况 下 ， 值 理 论 估 算 的保 证 金 ％ 极
水平 能 很好 的涵 盖 了 ９ ％ 的价 格 波 动风 险． ９ 目前 极
有较 大 的要求 ， 得 该类 方法 在 设 置 我 国股 指期 货 使 保 证 金水 平方 面可 能存 在较 大 的误差 ． 本 文 通过将 离散 型 随机 变量 与连续 型 随机 变量 的复合 ， 立一 个 既能反 映某 时段 内风 险发 生次 数 ， 建
极值理论在风险价值度量中的应用
极值理论在风险价值度量中的应用1、引言自20世纪70年代以来,金融市场的波动日益加剧,一些金融危机事件频繁发生,如1987年的“黑色周末”和亚洲金融危机,这使金融监管机构和广大的投资者对金融资产价值的暴跌变得尤为敏感;金融资产收益率的尖峰、厚尾现象也使传统的正态分布假定受到严重的质疑,因此如何有效地刻画金融资产收益率的尾部特征,给出其渐进分布形式,及各种风险度量模型的准确估计方法和置信区间,依此制定投资策略,确定国家监管制度,成为风险度量和管理所面临的巨大挑战;目前,对金融资产损失的估计方法主要包括历史模拟、参数方法和非参数方法;历史模拟是一种最简单的方法,它利用损失的经验分布来近似真实分布,但是该方法不能对过去观察不到的数据进行外推,更不能捕获金融资产收益序列的波动率聚类现象,而受到大量的批评;参数方法假设收益符合某种特定的分布如：正态分布、t分布等,再通过分布与样本的均值、方差的匹配对参数进行估计,或者是假设收益符合某种特定的过程如：ARMA模型、GARCH模型,该方法可以在一定程度上解释尖峰后尾现象和波动率聚类问题,具有比较好的整体拟和效果;不过参数方法只能对已经到来的灾难信息给出准确的估计,对于即将到来的灾难信息无法给出准确的预测,因此对极端事件的估计缺乏准确性;非参数方法则主要包括极值理论EVT,该理论不研究序列的整体分布情况,只关心序列的极值分布情况,利用广义帕累托分布generalized Pareto distribution或者广义极值分布generalized extreme value distribution来逼近损失的尾部分布情况;Danielsson and de Vries1997以7支美国股票构成的组合为样本比较各种模型的表现情况,发现EVT的表现比参数方法和历史模拟方法明显的好;Longin2000认为极值理论的优点在于它的没有假设特定的模型,而是让数据自己去选择,而GARCH模型作为估计风险的一种方法,它只能反映当时的波动率情况,对于没有预期到的变化缺乏准确性;不幸的是,Lee and Saltoglu2003把EVT模型应用到5个亚洲股票市场指数上,发现表现令人非常不满意,而传统的方法尽管没有一个在各个市场表现都是绝对优于其它模型的,但都比EVT模型的表现好;本人认为EVT模型之所以在亚洲市场表现不好主要是因为亚洲金融市场的数据具有很强的序列相关和条件异方差现象,不能满足EVT模型要求的独立同分布假定;另外,Jondeau and Rockinger1999,Rootzen and Kluppelberg1999,Neftci2000,Gilli and Kellezi2003和Christoffersen and Goncalves2004也分别采用极值原理和其他模型对金融数据的尾部特征进行了分析和比较;本章在传统单纯采用极值理论假设被分析数据是独立同分布的描述金融资产收益尾部特征的基础上,把ARMA－AsymmetricGARCH模型和极值理论有机的结合起来;首先利用ARMA －AsymmetricGARCH模型捕获金融数据中的序列自相关Correlation和异方差Heteroskedasticity现象,利用GMM估计参数,获得近似独立同分布的残差序列,再采用传统的极值理论对经过ARMA－AsymmetricGARCH模型筛选处理过的残差进行极值分析,在一定程度上克服了传统单纯采用极值理论时,由于金融数据序列自相关和波动率聚类现象不能满足极值理论假设所造成的估计误差;另外,本章还采用Bootstrap的方法给出了采用极值理论估计出的VaR和ES在某一置信水平 下的置信区间改进了采用似然比率法估计置信区间时,由于极值事件的小样本所造成的误差;最后,我们利用中国上证指数自1990年12月19到2004年9月30日的对数日收益率进行实证研究给出上证指数的VaR和ES值,及置信区间;2、VaR和ES的概念：VaRValue－at－Risk是一种被广泛接受的风险度量工具,2001年的巴塞耳委员会指定VaR模型作为银行标准的风险度量工具;它可以定义为在一定的置信水平p下,某一资产或投资组合在未来特定时间内的最大损失,或者说是资产组合收益损失分布函数的分位数点;假设X代表某一金融资产的收益,其密度函数为()f x,则VaR可以表示为：inf{|()(1)}p VaR x f X x p =-≤>- 1当密度函数()f x 为连续函数是也可以写作：1()p VaR F p -=-,其中1F -称为分为数函数,它被定义为损失分布()F x 的反函数;该模型计算简单,在证券组合损失X 符合正态分布,组合中的证券数量不发生变化时,可以比较有效的控制组合的风险;但是VaR 模型只关心超过VaR 值的频率,而不关心超过VaR 值的损失分布情况,且在处理损失符合非正态分布如后尾现象及投资组合发生改变时表现不稳定,会出现()()()p p p VaR X Y VaR X VaR Y +≥+ 2的现象,不满足Artzner1999提出了一致性风险度量模型的次可加性;()p ES Expected shortfull 满足Artzner1999提出的次可加性、齐次性、单调性、平移不变性条件,是一致性风险度量模型;它的定义如下：在给定的置信水平p 下,设X 是描述证券组合损失的随机变量,()[]F x P X x =≤是其概率分布函数,令1()inf{|()}F x F x αα-=≥,则()()ES X α可以表示为：11()01()()p p ES X F d p αα--=-⎰ 3 在损失X 的密度函数是连续时,()p ES 可以简单的表示为：{|()(1)}p ES E x F x p =-≤-; 本章将分别选用这两个模型来度量金融资产的风险,给出在修正过的极值模型下,其估计的方法和置信区间;3. ARMA －AsymmetricGARCH 模型ARMA －AsymmetricGARCH 模型的性质ARMA(p,q)模型：11p qt i t i j t j t i j y y μφξεε--===+++∑∑ 4其中,t ε是期望为0,方差为常数2σ的独立同分布随机变量,ARMA(p,q)模型在可逆的情况下可以表示为()AR ∞;该模型假设t y 的条件期望是可得的,条件方差为常数,通常可以用来解释时间序列的相关性,并可以对时间序列进行的短期预测;但是该模型条件方差为常数的假设,使其无法有效的解释在金融时间序列中经常被观察到的波动率聚类现象,为此,我们需要在模型中进一步引入GARCH 模型;我们令t t t z h ε=,其中t z 是期望为0,方差为常数1,的独立同分布随机变量,2t h 是t ε在t 时刻的条件方差;这里我们采用通常使用的最简单的(1,1)GARCH 模型,则条件方差可以表示为：2220111t t t h a a bh ε--=++,(1,1)GARCH 模型也可以表示成平方误2t ε的形式：22222201111()()()t t t t t t a a b b h h εεεε---=++--+- 5其中221(()|)0t t t E h F ε--=,因此(1,1)GARCH 模型本质上是平方误2t ε的ARMA(1,1);(1,1)GARCH 模型的引入不仅可以捕获到金融时间序列的波动率聚类现象,而且可以在一定程度上改善t z 尖峰后尾现象,因为44444422222222()()()()t t t t h z t t t t E Ez Eh Ez k k E Ez Eh Ez εε==≥= 6 其中4h k 和4z k 分别表示t h 和t z 的峰度,t h 的峰度明显大于等于t z 的峰度;另外,在金融序列中我们还可以明显的观察到,波动率正方向变动与收益率负方向变动的相关性大于与收益率正方向变动的相关性,一种可能的解释是收益率的负方向变动会加大波动幅度;而(1,1)GARCH 模型认为收益的正方向变动和负方向变动对波动率变动幅度有着相同的影响,为了捕获金融序列波动率变动的这一不对称性,我们引入需要Glosten et al1993提出的非对称(1,1)GARCH 模型：22220112111sgn()t t t t t h a a a bh εεε----=+++ 7其中00sgn()sgn()10t t x z x ε<⎧==⎨≥⎩,在这个模型中我们通过21sgn()t a ε-项来捕获收益率的正负变动对波动率变动的不同影响,如果收益率的波动与收益率波动率的变动像我们上面所预期的那样,则20a <;这样我们就得到了ARMA －AsymmetricGARCH 模型112222*********(sgn())()()p q t i t i j t j ti j t t t t t t t t t t y y z h a a a b b h h μφξεεεεεεεε--==----⎧=+++⎪⎪⎪=⎨⎪=+++--+-⎪⎪⎩∑∑ 8 、ARMA －AsymmetricGARCH 模型的参数估计：我们知道在条件正态分布的假设下,可以很容易的利用ARMA －AsymmetricGARCH 模型的似然函数,给出参数向量012(,,,,,,)a a a b θμφξ''=的估计值,其中1(,,)p φφφ'=,1(,,)q ξξξ'=;即使在金融收益率序列残差不满足条件正态分布的情况下,使用正态极大似然估计法,仍然可以得到参数θ的一致渐进正态非最小方差估计;但是这样我们得到的残差t z 将有很大的误差,而t z 是我们下一步进行EVT 尾部估计的输入变量,它的有效性将会直接影响我们整个的估计结果,为此我们必须寻找一个更有效的估计方法;GMMGeneralized Method of Moments 广义矩估计恰好可以满足我们的要求,它不需要假设t z 符合任何分布,只需要t z 的条件矩;在Skoglund2001“A simple efficient GMM estimator of GARCH models ”给出了该估计方法的计算过程和收敛情况;下面给去估计的步骤：首先,定义一个行向量 22[,()]t t t t r h εε=- 和广义向量 t t t g I r '=,其中t I '是工具变量,则参数θ的GMM 估计可以通过下式得到：11111min[][]T TtT t t t T g W T g θ---∈Θ=='∑∑ 9 其中11TT tt W T W -==∑是一个恰当的权重矩阵; 在Newey and McFadden1994中,我们可以知道,有效的GMM 估计可以通过另1()t t t r I θ-∂='∂∑,1()()t t t t r r W θθ-'∂∂='∂∂∑,其中11var(|)t t t r F --=∑,()t r θ∂'∂是Jacobian 行列式;把t I 和W 带入上面的目标函数9得到：1111[]()[]T T TT t t t t t t Q T g g -==='=Λ∑∑∑ 10 其中,t t t g g 'Λ=是一个含有参数θ的权重矩阵,它的元素可以表示为：其中,243[(1)]v v ∆=--,2211t t t h h c b θθ--∂∂=+∂∂,(0,0,0,0,)t t X εθ∂'=-∂,2221111112111(1,,sgn(),,2[sgn()])t t t t t t t t c h a a X εεεεε--------'=-+11(1,,,,)t t t p t t q X y y εε----=,1(|)k k t t v E z F -=通过上面对目标函数9的变化,我们得到函数T Q 是恰好可识别的,即参数θ的最优估计是使函数T Q 等于0;另外,我们要进行GMM 估计还需要一个对参数 的初始估计值和对t z 的三阶矩和四阶矩的初始估计值,而这一初始值我们可以通过对ARMA －AsymmetricGARCH 模型残差符合正态分布的情况进行最大似然估计得到;这样我们就可以得到有效的参数估计值和残差序列t z ;4、极值理论极值理论是测量极端市场条件下风险损失的一种方法,它具有超越样本数据的估计能力,并可以准确地描述分布尾部的分位数;它主要包括两类模型：BMM 模型Block Maxima Method 和POT 模型Peaks over Threshold,两类模型的主要区别有：1、极值数据的获取方法上的区别,BMM 模型通过对数据进行分组,然后在每个小组中选取最大的一个构成新的极值数据组,并以该数据组进行建模；POT 模型则通过事先设定一个阀值,把所有观测到的超过这一阀值的数据构成的数据组,以该数据组作为建模的对象,两个模型的共同点是只考虑尾部的近似表达,而不是对整个分布进行建模;2、两个模型分别采用极值理论中的两个不同的定理作为其理论依据,同时也因为获取极值数据的不同方法导致两个模型分别采用不同的分布来拟合极值数据;3、BMM 模型是一种传统的极值分析方法,主要用于处理具有明显季节性数据的极值问题上,POT 模型是一种新型的模型,对数据要求的数量比较少,是现在经常使用的一类极值模型;4、BMM 模型主要用于对未来一段较长的时间内的VaR 和ES 预测,而POT 则可以进行单步预测,给出在未来一段小的时间内VaR 和ES 的估计值;5、BMM 模型的前提条件是样本独立同分布,POT 模型的前提条件是超限发生的时间服从泊松分布,超限彼此相互独立,服从GPDgeneralized Pareto distribution 分布,且超限与超限发生的时间相互独立;样本独立同分布可以保证POT 模型的前提条件;BMM 模型的理论基础假设n M 表示我们采用BMM 方法获得的极值数据组,其中n 表示每个子样本的大小,则有下面的极限定理成立定理1：Fisher and Tippett 1928, Gnedenko 1943假设()t X 是一个独立同分布的随机变量序列,如果存在常数0n c >,n d R ∈,以及一个非退化的分布函数H ,使得 d n n nM d G c -−−→成立,则分布函数G 一定属于下面的三种标准的极值分布： Frechet ： 00()00x x x e x ααα--≤⎧⎪Φ=>⎨>⎪⎩ Weibull ： ()0()010x e x x x ααψα--⎧≤⎪=>⎨>⎪⎩Gumbel ： ()x e x e x R --Λ=∈从图1可以清楚的Frechet 分布用来描述那些极值无上界有下界的分布,Weibull 分布用来描述极值分布有上界,无下界的分布,Gumbel 分布用来描述极值无上界也无下界的分布;我们通常见到的很多分布函数都可以根据他们尾部的状况划分到上面的三种极值分布分布中去,例如：学生分布、帕累托分布Pareto distribution 、对数Gamma 分布、Cauchy distributed 根据尾部特征可以划分到Frechet 分布中去；均匀分布和Beta 分布的尾部分布可以收敛到Weibull 分布；正态分布、Gamma 分布和对数正态分布的尾部分布都收敛到Gumbel 分布;图形1：标准Frechet 、Weibull 和Gumbel 分布图但是,在实际应用中对于一个给定得极值序列,我们应该如何在这三种极值分布中做出选择呢;一种理想的方法是通过参数的形式把三种极值分布统一的表示成一个分布函数,这样我们就可以在利用最大似然估计的时候,把该参数也一块估计出来,让数据去决定它们的选择,这将极大的增加模型估计的准去性;这里我们采用 Jenkinson and Mises 的方法,把三种分布表示成如下单参数的形式：1/(1)0()0x x e e G x e ξξξξξ---+-⎧≠⎪=⎨=⎪⎩ 11其中{}|10x x x ξ∈+>,这一表达形式也被称为广义极值分布函数Generalized extreme value distribution,当1ξα-=时,表示Frechet 分布,当1ξα-=-时,表示weibull 分布,当0ξ=时表示Gumbel 分布;在定理1的基础上,对于给定一个金融资产的残差序列,我们就可以首先分组求最大值得到的极值序列记为X ;为了表达上的简洁,用u 和σ代替公式11中的n c 和n d ,则可以序列X 的近似分布函数：1,,exp((1))0()exp(exp())0u x u x u G G x u ξξσξξξσσξσ--⎧-+≠⎪-⎪==⎨-⎪--=⎪⎩ 12 其中|10x u x x ξσ-⎧⎫∈+>⎨⎬⎩⎭;然后,我们要对参数进行最大似然估计,这需要得到随机变量X 的概率密度函数,通过概率分布函数12对x 求导,我们得到随机变量X 的概率密度函数：(1)1(1),1(1)exp((1))0()1exp()exp(exp())0u x u x u g x x u x u ξξξσξξξσσσξσσσ-----⎧+-+≠⎪⎪=⎨--⎪---=⎪⎩ 13 其中|10x u x x ξσ-⎧⎫∈+>⎨⎬⎩⎭;通过似然函数就可以得到各参数的估计值： ,,,,(,,)(,,)arg max ()arg max(ln(()))u u i u u iL x g x ξσξσξσξσ=∑ 14在各参数估计值给定的基础上,我们就可以利用极值分布函数计算不同p 下的分位数值,如用p R 表示这一分位数,则在1p个周期内出现的极值收益会超过这一阀值的预期数量有且仅有一次; p R 表达形式为： ˆ1,,ˆˆ(1(ln(1)))0ˆˆ(1)ˆˆln(ln )0p u u p R G p u p ξξσσξξσξ--⎧----≠⎪=-=⎨⎪--=⎩ 15 注：关于参数的置信区间的确定我们在后面给出其计算方法; POT 模型的理论基础假设序列{}t z 的分布函数为()F x ,定义()u F y 为随机变量Z 超过门值u 的条件分布函数,它可以表示为：()(|)u F y P Z u y Z u =-≤> 0y ≥ 16根据条件概率公式我们可以得到：()()()()()1()1()()()(1())()u u F u y F u F z F u F y F u F u F z F y F u F u +--==--⇒=-+ z u ≥ 17定理2：Pickand 1975, Balkema and de Haan 1974对于一大类分布F 几乎包括所有的常用分布条件超量分布函数()u F y ,存在一个,()G y ξσ'使得：1/,/1(1)0()()10u y y F y G y u e ξξσσξξσξ--⎧-+≠⎪'≈=→∞⎨⎪-=⎩ 12当0ξ≥时,[0,)y ∈∞；当0ξ<时,[0,]y σξ∈-;分布函数,()G y ξσ'被称作广义的Pareto 分布; 图2：广义Pareto 分布在1σ=,ξ取,0,的图形从图形上我们可以看到ξ的不同取值确定了尾部的厚度,ξ越大则尾部越厚,ξ越小尾部越薄,从,()G y ξσ'函数我们还可以得到当0ξ<时,y 的最大取值为σξ-,有上界;Lee and Saltoglu2003在金融资产收益时间序列上直接使用EVT 时,由于序列的尖峰后尾,使得确定出来的ξ一定是大于零的,但是在我们的模型中,我们对残差序列进行极值分析,因此我们得到的ξ不一定大于零;根据公式12我们可以得到广义的Pareto 分布的概率密度函数,()g y ξσ'： 因此对于给定的一个样本1{,,}n z z ,对数似然函数(,|)L z ξσ可以表示为：111ln (1)ln(1)0(,|)1ln 0n i i n i i n y L y n y ξσξξσξσσξσ==⎧--++≠⎪⎪⎨⎪--=⎪⎩∑∑ 13在POT 模型中另一个重要的问题,那就是如何确定我们定理2中的阀值u ,它的确定非常重要,它是正确估计参数ξ和σ的前提;如果阀值u 选取的过高,会导致超额数据量太少,使估计出来的参数方差很大；如果阀值u 选取的过低,则不能保证超量分布的收敛性,使估计产生大的偏差;Danielsson et al1997、de Vries1997和Dupuis1998给出了对阀值u 的估计方法,一般有两种：一是根据Hill 图,令(1)(2)()n X X X >>>表示独立同分布的顺序统计量;尾部指数的Hill 统计量定义为：Hill 图定义为点1,(,)k n k H -构成的曲线,选取Hill 图形中尾部指数的稳定区域的起始点的横坐标K 所对应的数据k X 作为阀值u ;二是根据样本的超额限望图,令(1)(2)()n X X X >>>,样本的超限期望函数定义为：()()1ni i k X u e u n k =-=--∑ min{|}i k i X u => 14超限期望图为点(,())u e u 构成的曲线,选取充分大的u 作为阀值,它使得当x u ≥时()e x 为近似线性函数;另外,如果超限期望图当x u ≥时是向上倾斜的,说明数据遵循形状参数ξ为正的GPD 分布,如果超限期望图当x u ≥时是向上倾斜的,说明数据来源于尾部较短的分布,如果如果超限期望图当x u ≥时是水平的,则说明该数据来源于指数分布;这一判断方法是根据广义Pareto 分布在参数1ξ<的时候,它超限期望函数()e m 是一个线性函数;()(|)1m e m E X m X m σξξ+=->=+ 0m σξ+> 15 注：因为对于广义Pareto 分布只存在1ξ⎢⎥⎢⎥⎣⎦阶矩,如果1ξ<则存在一阶矩,否则一阶矩将不存在,就没有办法计算超限期望函数;当u 确定以后,；利用{}t z 的观测值,根据公式13进行最大似然估计得到ˆξ和ˆσ;同时,我们得到{}t z 的观测值中比阀值u 大的个数,记为u N ,根据公式17用频率代替()F u 的值,可以得到()F z 在x u >时的表达式：1/()/1/()/(1(1()))(1)()()(1())()(1)(1)1(1())010u u u z u u u u z u u N N z u N N F z F y F u F u N N e N N N z u N N e N ξσξσξσξξσξ------⎧-+-+-⎪⎪=-+=⎨⎪-+-⎪⎩⎧-+-≠⎪⎪==⎨⎪-=⎪⎩16对于给定某个置信水平p ,可以由()F z 的分布函数公式15可以得到：(()1)0ln()0(|)u p u p p p p N u p N VaR N u p N ES VaR E Z VaR Z VaR σξξσξ⎧+-≠⎪⎪=⎨⎪-=⎪⎩=+-> 17 根据GPD 的条件分布函数公式15可以得到：()111p p p p VaR u VaR u ES VaR σξσξξξξ+--=+=+--- 18 序列t z 的VaR 和ES 置信区间的估计方法：通常,对于参数置信区间的估计方法,在大样本的情况下我们可以从似然比率检验Likelihood Ratio Test 的思路中获到;似然比率检验用来检验两个同类型模型的拟和程度的好坏;两个同类型模型的似然比率符合2χ分布,它的自由度等于复杂模型中新加入的参数的个数;以POT 模型为例,要估计参数ξ和σ在给定置信水平α下的置信区间可以通过下式得到：其中,ˆξ和ˆσ为估计的最优值,(,)L x y 表示似然函数;这样我们就得到了ξ和σ的联合置信区间,如果我们希望得到p VaR 的估计值,则可以根据公式17反解出σ带入公式13得到(,)p L VaR ξ,令()max (,)p p L VaR L VaR ξξ=,p VaR 的置信区间可以通过下式得到： 但是,于超过阀值的极值数据量不会很多,使的这一估计的渐进效果可能不佳;为此,我们引入Bootstrap 方法来获得置信区间的估计;既然我们得到的序列{}t z 时独立同分布,就可以每次独立从中抽取N 个点组成新的序列,用该序列估计p VaR 和p ES ,重复这一操作,就可以得到一系列的p VaR 和p ES 估计值,求出p VaR 和p ES 的经验分布,最后根据经验分布得到p VaR 和p ES 的置信区间,并把p VaR 和p ES 的期望值作为p VaR 和p ES 的估计值;我们在这里只给出了POT 模型中p VaR 置信区间的求法,其他参数的置信区间可以类似的求得;该方法在确定置信区间的同时,也是一种检验模型稳定性的方法;5、实证分析我们采用上海证券交易所公布的日收益综合指数P 为原始数据数据来源：大智慧,样本空间选自1990年12月19日---2004年9月30日;样本容量为3391个使用Eviews 和Matlab 软件;我们定义收益为1ln ln t t t R P P -=-;我们的实证过程分为四步,1用ARMA －AsymmetricGARCH 模型对收益序列进行过虑得到近似独立同分布的残差序列t ε；2用极值理论对这一残差进行分析,给出其渐进分布,并估计出相应的VaR 和ES 值;3比较用似然比率和用Bootstrap 方法给出VaR 和ES 值的置信区间的估计;4整合第一步和第二步的结果,计算收益t R 的VaR 和ES 值;5利用BMM 模型估计VaR 值;ARMA －AsymmetricGARCH 模型形式和参数的确定首先给出收益序列t R 的描述性统计量图1,可以看到序列具有明显的尖峰后尾现象,从J －B 检验可以显着的拒绝正态性假设;对收益序列进行单位根ADF 检验见表1,因为检验的t 统计量是23.64516-,比显着性水平为1％的临界值还小,所以拒绝原假设,序列不存在单位根,是平稳序列;图1：收益序列t R 的描述性统计量ADF TestStatistic 1% Critical Value5% CriticalValue10% CriticalValueMacKinnon critical values for rejection of hypothesis ofa unit root.表1：序列的单位根检验可以进一步分析数据的自相关和偏相关见表2现象,发现滞后10期,在99％的置信水平下都不能拒绝没有自相关和偏相关的原假设,为此可以认为收益序列中不存在ARMA 现象;这样,我们就可以直接用t R 序列对常数项作最小二乘回归得到残差项t ε,然后对残差序列t ε进行ARCH 效应的LM 检验见表3,发现当q 取比较大的值8时的相伴概率仍然有0.044476p =,小于显着水平0.05,拒绝原假设,残差序列存在高阶ARCH 效应,即有GARCH 效应;表2：样本数据的自相关和偏相关表ARCH Test:F-statistic ProbabilityObsR-squared Probability表3：ARCH 效应的LM 统计量检验根据上面的分析,我们可以确定在第一步中所采用的模型公式19,并对其进行正态最大似然估计见表4;22220112111(0,1)sgn()t t t t t t t t t t t y c z h z iid h a a a bh εεεεε----=+=⎧⎪⎨=+++⎪⎩ 19 Coefficient Std. Error z-Statist ic Prob.CVariance EquationCARCH1RESID<0ARCH1GARCH1表4：公式19最大似然估计的结果从表中可以看到,正如我们所预见的那样20.12412a =-,预测不到收益的负方向变动可以导致更大的波动率出现,正方向变动会使波动率下降;1a 和b 都大于零表明过去时刻的波动对未来价格波动有着正向缓解作用,从而可以有效的解释了波动率的聚类性现象;下面我们以最大似然估计的结果为初始值按照前面所介绍的方法进行GMM 估计,其结果如下表：表5：最大似然估计和GMM 估计比较在GMM 估计值与最大似然估计值的比较中,我们可以清楚的看到,GMM 估计明显的增加了非对称项系数的绝对值,使收益的正负向变动对波动率表动的不同影响更加明显;另外,在最大似然估计中12ˆˆˆ 1.0690431a a b ++=>,这意味tε不存在有限的方差,而在GMM 估计中12ˆˆˆ0.950091a a b ++=<,保证了t ε的方差有限性;GMM 估计在没有tz 分布的情况下给出了参数的取值,并有效的降低了目标函数T Q 的取值;把GMM 估计值代入公式19,由收益序列得到残差序列t z 见图2,从图像上可以看出序列t z 变的更平稳,波动率聚类现象明显下降,更接近于独立同分布;对其进行一阶,二阶自相关和偏相关性检验和Ljung –Box 检验,结果都在很高的水平上拒绝原假设,表示残差序列t z 以没有ARMA 现象和条件异方差现象;图2：收益序列R 和残差序列zARCH Test:F-statistic ProbabilityObsR-squared Probability表6：序列t z 的ARCH 检验POT 模型的应用基于极值理论中POT 模型,我们需要利用充分大的阀值u ,对超限分布进行GPD 拟合,根据公式14,得到超限期望图见图3;发现样本的平均超限函数图在0.8u ≥时近似直线,具有明显的Pareto 分布特征;当0.8u =时数据超过阀值u 的个数369u N =；当0.9u =时287u N =；当1u =时228u N =,我们的总样本个数3390N =,在u 允许的情况下选取10％左右的数据DuMouchel1983作为极值数据组是比较合适的选择,否则可能不能抓住序列尾部分布的特征,样本内过度拟合,样本外不适用;为此,我们分别给出阀值取,的情况下,利用最大似然估计得到各参数、0.01VaR 、0.01ES 的取值和95％的置信区间见表7,以及在这些参数下的Q －Q 图和分布图见图4和图五,从图形中我们可以看到极值分布有效拟合了我们的样本分布,只有个别地方出现异常现象;且在0.8u =和0.9u =两种情况下的拟合效果没有明显的区别,为此在后面我们只给出0.9u =时的图形;估计值上界区间长度表7：参数的最大似然估计和95％置信区间图3：序列zt的超限期望图(,())u e u图4：u=时的Q－Q图0.8u=和0.9图5：0.8u=和0.9u=极值分布与经验分布的比较α=的取值范围在3到4之间,而我对于ˆξ的估计 Embrechets1999认为金融序列的1们这里计算出来的α,几乎不落在[3,4]的区域内,这主要是因为我们对金融序列用ARMA －AsymmetricGARCH 模型进行了过滤,得到的序列t z 在一定程度上消除了的尖峰后尾现象,使得ˆξ估计出来的值偏小,这与Embrechets 的结论并不矛盾;另外在Q －Q 图中,我们可以看到在的分为数之前拟合效果非常好,在后面出现了个别的异常值,这不会影响我们对0.01VaR 的估计,因为0.01VaR 只关心分为数之前的分布情况,而不受到分为数之后分布情况的影响;但是0.01ES 的估计由于受到分为数之后分布情况的影响,所以这会对0.01ES 的估计造成一定的误差,这也是为什么我们在表7中看到0.01ES 的95％估计区间明显比0.01VaR 的95％估计区间要宽的原因之一;下面我们采用Bootstrap 的方法来确定各参数的置信区间,首先在序列t z 中进行3390次重复抽取得到一个包含3390个数据的新样本,利用这些新样本估计ξ、σ、0.01VaR 和0.01ES 取值,重复上述1000次,则得到四个估计序列,其中每个序列中包含了1000关于某个参数的估计值,我们把他看作是一个样本,把这些样本与前面估计出来的参数区间相比较,如图6左其中方形区域是ξ、σ单参数确定的95％置信区间,椭圆形区域是ξ、σ的95％联合置信区间,图形中的散点表示每次估计出来的ξ、σ的值构成的点;从图形中我们可以看到大概有5％的点落在了95％的联合置信区间的外面,但是当我们考虑单参数置信区间时发现在区域以外的点大大超过了5％,这表明单参数估计的置信区间存在一定的问题,类似的现象我们还可以在ξ和0.01ES 的估计中见图6右看到,联合置信区间比较准确的捕获了数据的特性,单参数置信区间的表示方法就有较大的误差;图6：单参数和联合置信区间,以及bootstrap 的估计点图7：Bootstrap 方法得到的ξ、σ、0.01VaR 和0.01ES 的经验分布图另外,从四个参数估计序列我们可以得到四个参数的经验分布见图7,通过线性插值的。
跨品种期货套利交易最优保证金比率设计--基于Copula函数及极值理论的研究
跨品种期货套利交易最优保证金比率设计--基于Copula函数及极值理论的研究赵成珍;宋锦玲【摘要】文章主要研究期货交易中最优风险保证金的比率设定问题，通过考虑变动保证金的设定来弥补目前固定保证金不足的问题。
在期货套利交易保证金问题上，文章给出了商品期货套利交易保证金的测定理论。
在测定的过程中，引入了收益率的非对称Laplace分布和GARCH-T函数分布，在Gumbel Copula函数的基础上，应用蒙特卡洛模拟算法对两种商品的套利交易的保证金问题给出了测定，并以豆油和豆粕期货为例，选取时间序列进行了实证研究，得出结论：套利交易的保证金水平在理论上小于非套利交易下两单独品种保证金的收取之和。
此外再结合极值理论，在改进VAR的基础上，得出在极端情况下的套利交易的最优保证金比率设计。
%This paper mainly studies the setting problem of the optimal risk trading margin ratio in futures trading, by changing the variation margin setting to make up the insufficient of the fixed deposit. As for the futures arbitrage trading margin issue , this paper first presents the determination theory of commodity futures arbitrage trading margin. In the process of determination, the yield of asymmetric Laplace distribution and GARCH-T function distribution are introduced. On the base of the Gumbel Copulas function, the margin ratio determination of the two goods by the Monte Carlo Simulation Algorithm is provided. The paper takes soybean oil and soybean meal futures as an example and does empirical research on the time series of them and draws the conclusion: in theory, the arbitrage trading margin level is less than the sum of two separate varieties depositcharge in the carry trade. In addition, the paper applies the extreme value theory and concludes that in extreme cases of arbitrage, the optimal margin ratio design on the basis of the improvement of VAR.【期刊名称】《技术经济与管理研究》【年(卷),期】2014(000)012【总页数】4页(P16-19)【关键词】期货交易;套利交易;套利品种;金融收益率【作者】赵成珍;宋锦玲【作者单位】中央财经大学金融学院，北京 100081;中央财经大学金融学院，北京 100081【正文语种】中文【中图分类】F724.520世纪80年代以来，全球多次爆发金融危机，金融机构的风险控制越来越得到学术界和业界的重视。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
研究生学位论文开题报告
研究课题极值理论在股指期货保证金
设定中的应用
研究生姓名郭丽芬
学号2010304310114
专业应用数学
导师姓名孙志宾
论文工作的
起止时间2010-10——2013-05
2010年10 月16 日填写开题报告由研究生填写，一式二份（学院，研究生部）
填表说明
1．只有学籍状态为注册的研究生才允许开题，开题报告及其评审结果才能被认可。
2．硕士学位论文开题报告封面及一至七项必须用计算机输入、打印。
3．开题报告为A4大小，于左侧装订成册。
硕士研究生应逐项认真填写，各栏空格不够时请自行加页。
4．开题报告经指导教师审阅通过后，举行公开的开题报告答辩会。
答辩会专家组由三名以上高级职称的专家（包括导师）组成，成员由各学院学术委员会确定。
答辩后由专家组进行评审并给出成绩。
答辩记录及成绩交研究生部备案。
5．硕士研究生应在选题前阅读相关领域的中外文资料，并写出不少于5000字的文献综述报告，引用参考文献的篇数不得少于20篇。
文献综述报告应反映国际和国内本领域的研究历史、现状和发展趋势。
文献综述是开题报告的必要附件，开题报告通过后，一份由学院留存，一份由研究生部备案。
6．“参考文献”著录书写顺序为：
（1）文献是期刊时，其格式为：[序号]作者.文章题目.期刊名，年份，卷号，期号（若期刊无卷号，则为：年份、期号）页码。
（2）文献是图书时，其格式为：[序号]作者.书名.出版地：出版社，出版年。
一选题依据
二研究内容
三研究方案
四研究工作进度安排
五预期研究成果
六本课题创新之处
七研究基础。