第第3章章 金融市场风险的度量金融市场风险的度量 1学习目标学习目标 通过本章学习，您可以了解或掌握：1. 金融市场风险度量方法的发展与演变；2. 灵敏度方法的基本原理及应用；3. 波动性方法的基本原理及应用；4. VaR方法的基本原理及应用；5. 基于历史模拟法的VaR计算方法；6. 基于Monte Carlo模拟法的VaR计算方法；7. 基于Delta、Gamma灵敏度指标的VaR计算方法；8. 压力试验和极值理论。2主要内容主要内容第一节 金融市场风险度量方法的演变第二节 灵敏度方法第三节 波动性方法第四节 VaR方法第五节 基于历史模拟法的VaR计算第六节 基于Monte Carlo模拟法的VaR计算第七节 基于Delta、Gamma灵敏度指标的VaR计算第八节 压力试验第九节 极值理论3第一节第一节金融市场风险度量方法的演变金融市场风险度量方法的演变4一、名义值度量法一、名义值度量法 1. 名义值度量法名义值度量法(Notional Amounts)的基本思想： 将资产组合的价值作为该组合的市场风险值。2.方法评价优点：方便简单缺点：只是粗略估计，一般会高估市场风险的大小5二、灵敏度二、灵敏度方法方法1. 灵敏度方法灵敏度方法(Sensitivity Measures)的基本思想可以通过基于Taylor展示式的资产组合价值随市场因子变化的二阶形式来展现： 6三、波动性方法三、波动性方法 1. 波动性方法波动性方法(Volatility Measure)的基本思想： 利用因市场风险因子变化而引起的资产组合收益的波动程度来度量资产组合的市场风险。2. 波动性方法实则统计学中方差或标准差的概念在风险度量中的应用。 7四、四、VaR方法方法 1. VaR(Value at Risk)的定义：指市场处于正常波动的状态下，对应于给定的置信度水平，投资组合或资产组合在未来特定的一段时间内所遭受的最大可能损失。2.VaR的应用领域金融风险度量确定内部经济资本需求设定风险限额绩效评估金融监管8五、压力试验和极值理论五、压力试验和极值理论1. 压力试验压力试验(Stress Testing)的核心思想：通过构造、 模拟一些极端情景，度量资产组合在极端情景发生时的可能损失大小。2. 极值理论极值理论(Extreme Value Theory)的核心思想： 应用极值统计方法来刻画资产组合价值变化的尾部统计特征，进而估计资产组合所面临的最大可能损失。9六、集成风险或综合风险度量六、集成风险或综合风险度量 1. 集成风险或综合风险集成风险或综合风险的定义：在各种风险“共同作用” 下金融机构所面临的整体风险。2.集成风险或综合风险的度量基于Copula函数的度量方法，其基本思想和步骤简要介绍如下：(1) 将引致集成风险的所有不同类型的风险驱动因子组成一个联合随机向量；(2) 得到单个风险因子的边缘分布函数；(3) 引入Copula函数，利用边缘分布函数获得随机向量的联合分布函数；(4) 基于联合分布函数，运用VaR等方法度量集成风险。 10第二节第二节灵敏度方法灵敏度方法11一、简单缺口模型一、简单缺口模型 1. 简单缺口模型简单缺口模型(Simple Gap Model)主要考察经营者 所持有的各种金融产品的缺口或净暴露情况以及市 场因子变动的幅度。 几个相关概念 正暴露：正暴露：有可能获得额外收益的金融产品的暴露； 负暴露：负暴露：有可能遭受损失的金融产品的暴露； 净暴露：净暴露：正暴露与负暴露之差的绝对值。12一、简单缺口模型一、简单缺口模型( (续续) ) 2.简单缺口模型的评价： 没有考虑期限对风险的影响，或者说没有考 虑正暴露和负暴露的期限结构对风险的影响。13二、到期日缺口模型二、到期日缺口模型1.利用到期日缺口模型度量金融风险的基本公式： GRSGR其中， GRSG：敏感性总缺口敏感性总缺口 R：某市场因子的变动幅度14二、到期日缺口模型二、到期日缺口模型(续续) 2.评价 (1)优点计算简单，便于实施。 (2)缺点没有考虑资产和负债所面临的市场风险；以经营者的资产负债表为基础，不能体现表外项目的市场风险；考察期的划分不可避免地存在着误差。15三、久期三、久期(一)久期的概念1.债券定价的基本公式 (3.2.1)2. 一阶泰勒展式 (3.2.2) 16三、久期三、久期 (一一)久期的概念久期的概念(续续)3. Macaulay久期由(3.2.1)式和(3.2.2)式，得Macaulay久期久期 (3.2.4) 17三、久期三、久期 （一）久期的概念（一）久期的概念(续续)4. 离散形式的久期公式 （3.2.5） 18三、久期三、久期 （一）久期的概念（一）久期的概念(续续)5. 调整久期调整久期或修正久期修正久期 （3.2.6） 19三、久期三、久期 （一）久期的概念（一）久期的概念( (续续) )6. 有效久期(Effective Duration)针对结构更为复杂的产品，提出有效久期有效久期的概念，定义如下： 20三、久期三、久期( (续续) )(二) 久期的性质 性质性质1 零息债券的久期是其到期期限，息票债券 久期的上限是相应的永久债券的久期。 性质性质2 息票债券的久期与息票率之间呈反向关系。 性质性质3 久期与贴现率之间呈反向关系。 性质性质4 债券到期日与久期之间呈正向关系。 性质性质5 债券组合的久期是该组合中各债券久期的 加权平均。21三、久期三、久期( (续续) )(三) 久期的缺陷1.对不同期限的现金流采用了相同贴现率，这与实际常常不符；2.仅仅考虑了收益率曲线平移对债券价格的影响，没有考虑不同期限的贴现率变动的不同步性； 3.仅仅考虑了债券价格变化和贴现率变化之间的线性关系，只适用于贴现率变化很小的情况。 22四、久期缺口四、久期缺口模型模型(一) 基本公式 (3.2.9) 其中， 称为久期缺口久期缺口(Duration Gap)。23四、久期缺口四、久期缺口模型模型( (续续) )(二) 评价1. 优点：考虑了每笔现金流量的时间价值，避免了到期日缺口模型中因时间区间划分不当而有可能带来的的误差，从而比到期日缺口模型更加精确。2.缺点： 计算较为复杂，对小规模的金融机构可能不够经济； 作为模型基础的久期概念存在一些不足。24五、凸性五、凸性(一)凸性的定义1. 二阶泰勒展式 （3.2.10）2.结合二阶泰勒展式和久期公式，得 其中 称为凸性凸性 25五、凸性五、凸性 (一一)凸性的定义凸性的定义(续续)3. 有效凸性对于内含期权以及其他现金流不确定的利率衍生产品，可以定义有效凸性有效凸性如下：26五、凸性五、凸性(续续)(二) 凸性的性质 性质性质1 贴现率增加会使得债券价格减少的幅度 比久期的线性估计值要小，而贴现率减少会使得债券价格增加的幅度比久期值估计值要大；而且凸性越大，上述效应越明显。 性质性质2 收益率和久期给定时，息票率越大，债券的凸性越大。 27五、凸性五、凸性 (二二)凸性的性质凸性的性质(续续) 性质性质3 通常债券的到期期限越长，债券的凸性越大，并且债券凸性增加的速度随到期期限的增加越来越快。 性质性质4 债券组合的凸性是组合内各种债券凸性的加权平均。 28六、六、系数和风险因子敏感系数系数和风险因子敏感系数 (一) 系数与资本资产定价模型1. 系数的公式表示 根据CAPM，在证券市场处于均衡状态时， （3.2.13）其中， 即为系数系数。 29六、六、系数和风险因子敏感系数系数和风险因子敏感系数 (一一) 系数与资本资产定价模型系数与资本资产定价模型(续续)2. 系数的理解 i系数实际上反映了证券i的超额期望收益率对市场组合超额期望收益率的敏感性；当系数取正值时，说明所考察的证券与市场组合的走势刚好一致，反之则反是；系数满足可加性。30六、六、系数和风险因子系数和风险因子敏感系数敏感系数(续续) (二) 风险因子敏感系数和套利定价模型 1. 风险因子敏感系数来源于Ross于1976年提出的套利定价理论（APT）。2. 套利定价理论的一般形式 (3.2.15) 其中， 称为第 k 个风险溢价因子 的风险风险 因子敏感系数因子敏感系数。 31七、金融衍生品的灵敏度测量七、金融衍生品的灵敏度测量 1. 金融衍生品的价格F 可以表示成下面的形式 F = F（S, t, r, ） （3.2.16） 其中：S表示标的物资产的当前价格，t表示当前时间，r表示无风险利率， 表示标的物资产价格的波动率。 32七、金融衍生品的灵敏度测量七、金融衍生品的灵敏度测量( (续续) )2. 金融衍生品定价公式的泰勒展式 （3.2.17） 33灵敏度指标灵敏度指标 公式公式 含义含义 (Delta)反映金融衍生品价格对其标的物资产价格的线性敏感性 (Gamma)反映灵敏度系数对标的物资产价格S的灵敏性 (Theta)反映金融衍生品价格对时间变化的敏感性 (Vega)反映衍生证券价格对其标的物资产价格波动率的线性敏感性（Rho）反映金融衍生品价格关于利率的线性敏感性七、金融衍生品的灵敏度测量七、金融衍生品的灵敏度测量( (续续) )3. 金融衍生品灵敏度指标的含义解析 34灵敏度指标灵敏度指标无收益资产组合的远期合约无收益资产组合的远期合约不付红利的欧式看涨期权不付红利的欧式看涨期权 Delta1Gamma0ThetaVega0Rho相互关系七、金融衍生品的灵敏度测量七、金融衍生品的灵敏度测量( (续续) )4.远期合约和期权的灵敏度指标 35八、灵敏度度量法评述八、灵敏度度量法评述 1.主要特点：简明直观；应用方便；最适合于由单个市场风险因子驱动的金融工具且市场因子变化很小的情形。36八、灵敏度度量法评述八、灵敏度度量法评述( (续续) ) 2. 不足：可靠性难以保证；难以定义受多个市场风险因子影响的资产组合的灵敏度指标；无法对不同市场因子驱动的风险大小进行横向比较；不能给出资产组合价值损失的具体数值；一阶灵敏度方法一般不考虑风险因子之间的相关性。37第三节第三节波动性方法波动性方法38一、单种资产风险的度量一、单种资产风险的度量 1.假设某种金融资产收益率r为随机变量，该资产的风险可用收益率标准差即波动系数来度量。 越大说明该资产面临的市场风险越大，反之则反是。39一、单种资产风一、单种资产风险的度量险的度量( (续续) ) 2. 当无法准确知道资产收益率的概率分布时，可利用随机变量r的若干个历史样本观测值来估计 r的数学期望和标准差：期望：标准差：40二、资产组合风险的度量二、资产组合风险的度量 (一)基本思路 用收益率的方差或标准差来度量资产组合的风险。(二)相关的计算公式1.数学期望 (3.3.3)2.方差 （3.3.4）3.相关系数 (3.3.5) 41三、特征风险、系统性风险与风险分三、特征风险、系统性风险与风险分散化散化(一)资产组合收益率方差 令 ,且所有单个资产的风险相同,则可得 资产组合收益率的方差为(二)讨论 1. 若 ,则 ，从而 。 2. 若 ,则 42四、波动性方法的优缺点评述四、波动性方法的优缺点评述 1. 优点：含义清楚，应用也比较简单。2. 缺点：对资产组合未来收益概率分布的准确估计比较困难;仅描述资产组合未来收益的波动程度，并不能说明资产组合价值变化的方向;无法给出资产组合价值变化的具体数值 。43第四节第四节 VaR方法方法44一、一、VaR方法的基本概念方法的基本概念(一) VaR的定义 指市场处于正常波动的状态下，对应于给定的置信度水平，投资组合或资产组合在未来特定的一段时间内所遭受的最大可能损失。用数学语言可表示为 (3.4.1) 45一一、VaR方法的基本概念方法的基本概念( (续续) )(二) VaR的基本特点：1.仅在市场处于正常波动的状态下才有效，无法准确度量极端情形时的风险；2.VaR值是一个概括性的风险度量值；3.VaR值具有可比性(Comparable) ；4.时间跨度越短，假定收益率服从正态分布计算的VaR值越准确、有效；5.置信度和持有期是影响VaR值的两个基本参数。 46一、一、VaR方法的基本概念方法的基本概念( (续续) )(三) 置信度和持有期的选择和设定1.持有期的选择和设定需考虑以下因素: (1) 考虑组合收益率分布的确定方式; (2) 考虑组合所处市场的流动性和所持组合头 寸交易的频繁性。47一一、VaR方法的基本概念方法的基本概念 (三三) 置信度和持有期选择和设定置信度和持有期选择和设定(续续) 2. 置信度的选择和设定需考虑以下因素： (1) 考虑历史数据的可得性、充分性； (2) 考虑VaR的用途； (3) 考虑比较的方便。48二、二、VaR的计算的计算(一) VaR的计算方法概括 1.计算VaR值的核心问题是估计资产组合未来损益P的概率分布。2.计算VaR的一般步骤 (1) 建立映射关系； (2) 建模； (3) 给出估值模型和VaR值。 49二二、VaR的计算的计算 (一一) VaR的计算方法概括的计算方法概括( (续续) ) 3. VaR计算方法的分类(根据P 分布确定方法划分)(1) 收益率映射估值法：收益率映射估值法：直接应用组合中资产的投资收益率来确定P分布。(2)风险因子映射估值法：风险因子映射估值法：将组合价值表示成风险因子的函数，然后通过风险因子的变化来估计组合的未来损益分布。进一步分为： 风险因子映射估值模拟法 风险因子映射估值分析法50二、二、VaR的计算的计算( (续续) )(二) 基于收益率映射估值法的VaR计算1. 绝对绝对VaR和相对相对VaR的概念(1) 以组合的初始值为基点考察持有期内组合的价值变化即 PA=P-P0= P0R (3.4.3) 由此求得的VaR称为绝对绝对VaR，记为VaRA 。(2) 以持有期内组合的预期收益为基点考察持有期内组合的价值变化，即PR=P-E(P)= P0(R-) (3.4.4) 由此求得的VaR称为相对相对VaR ，记为VaRR 。51二、二、VaR的的计算计算 (二二) 基于收益率映射估值法的基于收益率映射估值法的VaR计算计算(续续) 2. 组合的投资收益率服从正态分布的日VaR计算假设初始价值为P0，日投资收益率R服从正态分布，期望收益率与波动率分别为和，于是在置信度c下分别得到日绝对VaRA 和日相对VaRR：(3.4.6)(3.4.7) 52二二、VaR的计算的计算 (二二) 基于收益率映射估值法的基于收益率映射估值法的VaR计算计算(续续) 3. 组合中资产的投资收益率服从正态分布的日VaR计算 假设组合由n种资产构成，组合中n种资产的日投资收益率向量服从n维正态分布 ，则该组合的日绝对VaRA为其中，53二、二、VaR的计算的计算 (二二) 基于收益率映射估值法的基于收益率映射估值法的VaR计算计算(续续) 4. 关于资产组合的VaR计算 资产组合的初始价值 ，在置信度c下资产组合的日绝对VaR和日相对VaR分别为：日绝对VaR：日相对VaR： 54二、二、VaR的的计算计算 (二二) 基于收益率映射估值法的基于收益率映射估值法的VaR计算计算(续续) 5. 关于VaR的时间加总问题(1)基本思路：当求出1单位的VaR，可直接利用时间加总公式求出持有期为t的VaR。(2)计算公式根据独立同分布随机变量和的分布特征可知，组合在t日的投资收益率服从正态分布 ，于是t日的绝对VaR和相对VaR分别为 绝对VaR: 相对VaR:55二、二、VaR的计算的计算 (二二) 基于收益率映射估值法的基于收益率映射估值法的VaR计算计算(续续) 6.收益率映射估值法优缺点评述(1)优点：原理简单，尤其在正态分布假设下应用更加方便。 (2) 缺点：组合中金融工具之间相关系数的确定常常比较困难，计算量大;正态分布的假设常常与实际中的尖峰厚尾现象不符合。 56三、边际三、边际VaR、增量、增量VaR和成分和成分VaR(一) 边际VaR(Marginal VaR，简记为M-VaR) 设资产组合 ，所谓的边际边际VaR是指资产组合中资产的头寸变化而导致的组合VaR的变化，即 (3.3.14) 57三、边际三、边际VaR、增量、增量VaR和成分和成分VaR(续续)(二)增量VaR(Incremental VaR，简记为I-VaR) 假设在原来资产组合 的基础上，新增加另一个资产组合 ，并将调整后的资产组合的VaR记为VaR(w+dw)。 于是，dw的VaR，即增量增量VaR被定义为 I-VaR(dw)VaR(w+dw)VaR(w) (3.4.17)58三、边际三、边际VaR、增量、增量VaR和成分和成分VaR(续续)(三) 成分VaR(Component VaR，简记为C-VaR)1.定义 若资产组合 中资产 i 的VaR(记为 )满足 ， 则称 为该资产i 的成分成分VaR。 59三、边际三、边际VaR、增量、增量VaR和成分和成分VaR (三三) 成分成分VaR ( (续续) )2.C-VaR的特性：(1)组合中所有资产的成分VaR之和恰好等于组合的VaR；(2)资产i的成分VaR恰好为资产i对组合VaR的贡献份额；60三、边际三、边际VaR、增量、增量VaR和成分和成分VaR (三三) C-VaR (续续) (3)若某资产的成分VaR为负，则该资产可对冲组合其余部分的风险，且对冲量为成分VaR；(4)当资产组合的 n 维收益率向量 R 服从 n 维正态分布 时，资产 i 的成分VaR为： 61四、四、VaR方法的优缺方法的优缺点评述点评述 1. VaR方法的优点： VaR方法可以测量不同风险因子、不同金融工具构成的复杂资产组合以及不同业务部门所面临的总体风险；VaR方法提供了一个概括性的且具有可比性的风险度量值；VaR方法能更加体现出投资组合分散化对降低风险的作用。 62四、四、VaR方法的优方法的优缺点评述缺点评述( (续续) )2.VaR方法的局限性： (1) 决定组合价值变化的风险因子在未来的发展变化同过去的行为一致的隐含假定与实际不符； (2) 正态性假设不能准确刻画资产收益率分布经常出现的尖峰、厚尾、非对称等分布特征； (3) 基于同样的历史数据，运用不同方法所计算的VaR值往往差异较大； (4) 不能准确度量金融市场处于极端情形时的风险；63四、四、VaR方法的方法的优缺点评述优缺点评述( (续续) ) (5) 可能不满足次可加性(Sub-additive)； (6) 对组合损益的尾部特征的描述并不充分，从而对风险的刻画也不完全； (7) VaR方法得到的是统计意义上的结论，对个体所得结论并不确定； (8) 计算VaR时对历史数据的搜集、处理一般比较繁杂，而且有时还无法获得相应的历史数据；同时，计算复杂，计算量也比较大。 64第五节第五节基基于历史模拟法的于历史模拟法的VaR计算计算65引言引言1. 历史模拟法历史模拟法(Historical Simulation，有时简记为HS)包括：标准历史模拟方法标准历史模拟方法(Standard Historical Simulation，简记为SHS)加权历史模拟法加权历史模拟法(Weighted Historical Simulation)滤波历史模拟法滤波历史模拟法(Filtered Historical Simulation) 66一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的的基本原理和实施步骤基本原理和实施步骤(一) 基本原理 1. 将各个风险因子在过去某一时期上的变化分布或变化情景准确刻画出来，作为该风险因子未来的变化分布或变化情景；2. 在上述基础上通过建立风险因子与资产组合价值之间的映射表达式模拟出资产组合未来可能的损益分布；3. 计算给定置信度下的VaR。67一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的的基本原理和实施步骤基本原理和实施步骤( (续续) )(二) 一般计算步骤1. 识别风险因子变量，建立证券组合价值与风险因子变量之间的映射关系；2. 用历史数据模拟风险因子未来可能取值；3. 计算证券组合未来价值水平或损益分布；4. 基于损益分布计算置信度c下的VaR。68一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的的基本原理和实施步骤基本原理和实施步骤( (续续) )(三) 举例1. 假设某美国公司于1998年12月31日持有一份三个月后到期、以16.5百万美元交换10百万英镑的远期合约。用标准历史模拟法计算这家美国公司于1998年12月31日持有该合约在c=95%置信度下的日VaR值。 692. 定义以下符号：S ：以美元表示的英镑的即期价格；K ：货币远期合约中的约定价格，K=1.65；f ：远期合约的市场价值；r ：用年化的百分率表示的3个月的美元利率；r*：用年化的百分率表示的3个月的英镑利率； ：合约的到期期限，=92/365年； ：3个月的美元折现因子； ：3个月的英镑折现因子。 一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )703.分四步计算(1)第一步确定风险因子，分别为即期汇率S、美元利率r以及英镑利率r*；再建立远期合约的市场价值与上述市场风险因子之间的函数表达式，即 一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )71(2) 第二步，选取市场风险因子从1998年8月10日至拟考察的日期12月31日之间101个交易日的连续历史数据，并对应地计算出即期汇率S的值，见后面的表格。 一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )72t日期（1998年）r（%/年）r*（%/年）S（美元/英镑）012月31日4.93755.96881.6637112月30日4.9686.9671.6108212月29日5.0076.9891.6087312月28日4.9947.0161.60855412月25日4.9746.9821.61005一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )73(3) 第三步计算S、r和r*在1999年1月4日的100个可能取值，此时T=100；并对应计算出远期合约价值和损益值在1999年1月4日的100个可能取值，具体的计算结果见后面的表格。 一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )74t风险因子可能值r（%/年）r*（%/年）S（美元/英镑）14.9074.97061.716624.89855.94681.665834.95055.94181.6638544.95756.00281.662254.93755.94881.66315一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )75t风险因子可能值r（%/年）r*（%/年）S（美元/英镑）14.9074.97061.716624.89855.94681.665834.95055.94181.6638544.95756.00281.662254.93755.94881.66315r(1)r(0) r(0)r(-1)= 4.9375+ 4.9375- 4.968一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )76t风险因子可能值r（%/年）r*（%/年）S（美元/英镑）14.9074.97061.716624.89855.94681.665834.95055.94181.6638544.95756.00281.662254.93755.94881.66315S(4)S(0) S(-3)S(-4)=1.6637+1.60855-1.61005一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )77(4)第四步将远期合约在1999年1月4日的100个损益值的可能取值从大到小排列，可得到远期合约在1999年1月4日的损益分布；计算出95%置信度下的分位数为Tc= 10095% =95，则第Tc+1=96个数值26408.2977美元，即为美国公司持有该合约在95%置信度下的日VaR值。具体结果见后面的表格。 一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )78t风险因子可能值远期合约价值f的可能取值（美元）远期合约损益值f的可能取值（美元）r（%/年）r*（%/年）S（美元/英镑）14.9074.97061.7166655265.6353561684.635324.89855.94681.6658113592.00320011.00334.95055.94181.6638596672.741763091.74175544.95756.00281.662278235.15927-15345.840754.93755.94881.6631588969.37713-4611.62287f2f2f0=113592.003-93581一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )79Kt风险因子可能值远期合约价值f的可能取值（美元）远期合约损益值f的可能取值（美元）r（%/年）r*（%/年）S（美元/英镑）934.87655.90781.661672894.64051-20686.3595944.93455.99581.661671680.40446-21900.5955954.91856.05581.661971572.65537-22008.3446964.90855.92281.6609567172.70234-26408.2977974.86155.98281.6613566802.61306-26778.3869f从大到小排列一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )80Kt风险因子可能值远期合约价值f的可能取值（美元）远期合约损益值f的可能取值（美元）r（%/年）r*（%/年）S（美元/英镑）934.87655.90781.661672894.64051-20686.3595944.93455.99581.661671680.40446-21900.5955954.91856.05581.661971572.65537-22008.3446964.90855.92281.6609567172.70234-26408.2977974.86155.98281.6613566802.61306-26778.3869 95%置信度下的分位数：Tc=10095% =95VaR=f (kTc+1) =f (96)一、基于一、基于标准历史模拟法计算标准历史模拟法计算VaR的基本原的基本原 理和实施步骤理和实施步骤 (三三)举例举例( (续续) )81二、计算二、计算VaR的标准历史模拟法评述的标准历史模拟法评述 (一) 优点 1.直观、简单、便于理解，计算过程容易掌握；2.非参估计，减少参数估计风险和模型风险；3.不用假定市场风险因子未来变化的分布形式，可以处理非对称和尖峰厚尾等问题；4.能够处理非线性问题；5.原理简单而实用，所以容易与计算VaR的其他方法相融合，从而也容易被改进和推广。 82二、计算二、计算VaR的标准历史模拟法评述的标准历史模拟法评述( (续续) )(二) 不足1.风险因子的未来变化等同于历史数据变化的基本假设与现实不符；2.风险因子历史数据在未来时刻等概率出现的假设，与现实也经常不符；3.获取大量连续历史数据并非易事；4.得到的VaR值的波动性较大，稳健性较差；5.第三节所言有关VaR方法的一些缺陷仍然存在。 83三、计算三、计算VaR的标准历史模拟法的标准历史模拟法的修正及扩展的修正及扩展 (一) 时间加权历史模拟法 1.假设风险因子在过去第t期的变化值fi (-t)、可能价值V(t)和可能损益值V(t)在未来出现的可能性(权数)都是2. 根据置信度c计算分位数时，即求满足 的最大值m，于是损益分布中所对应的第m个值V(km)即为置信度c下的VaR。84三、计算三、计算VaR的标准历史模拟法的标准历史模拟法的修正及扩展的修正及扩展( (续续) )(二) 波动率加权历史模拟法1.根据历史数据建立风险因子时间序列模型；2.模拟风险因子在历史数据选用区间中的波动率以及未来时期的波动率 和 ；3.用下式对历史数据权重加以调整，再选择标准历史模拟法或者时间加权历史模拟法计算VaR。 85第六节第六节基于基于Monte Carlo模拟法的模拟法的VaR计算计算86一、一、Monte Carlo模拟模拟法法 (一) 基本原理与实施步骤1. 解决问题时如果没有实际数据，则无法借助随机抽样统计分析方法对总体进行推断，怎么办?在美国研制原子弹的“曼哈顿计划”中，需要计算中子进入反应堆屏障的随机性运动，但无法获得实际数据。解决办法：运用计算机产生随机数。冯诺伊曼（Von Neumann）借用赌城Monte Carlo来为这种方法命名。87一、一、Monte Carlo模拟法模拟法 (一一) 基本原理与实施步骤基本原理与实施步骤( (续续) )2. Monte Carlo模拟法的应用领域：(1)求解确定性问题积分的数值计算；各类方程的求解等。(2)求解随机性问题运筹学中的库存问题；随机服务系统中的排队问题；金融资产价格的变化问题等。883. 本文通过Monte Carlo模拟法计算资产组合VaR所涉及到的有关金融问题几乎都是随机性的，求解随机性问题的Monte Carlo模拟法的成功实施主要取决于三个基本要素：用以模拟随机变量未来变化路径的随机模型的准确性；每次模拟的独立性；足够多的模拟次数。 一、一、Monte Carlo模拟法模拟法 (一一) 基本原理与实施步骤基本原理与实施步骤( (续续) )89一、一、Monte Carlo模拟法模拟法( (续续) )(二) 单变量资产价格的随机模拟与随机数的产生1.单变量资产价格的随机模拟(1)几何布朗运动：(2)确定股票初始价格St，并估计出参数和；(3)利用计算机生成n个相互独立的标准正态随机数，代入上式递推得到股价的时间序列 ，得到股票价格的一条样本轨道；(4)重复得到N条样本轨道，及股价的变化分布。902. 单变量随机数的产生(1)第一步基于0,1上均匀分布的随机数的产生借助于在计算机上设立的所谓“随机数发生器”来实现 ；通过迭代算法生成大量的“伪随机数”。(2) 第二步，通过累积密度函数(或分布函数)的逆函数，把第一步产生的0,1上均匀分布的随机数转化为特定概率分布的随机数。 一、一、Monte Carlo模拟法模拟法 (二二) 单变量单变量 资产价格的随机模拟与随机数的产生资产价格的随机模拟与随机数的产生( (续续) )91一、一、Monte Carlo模拟法模拟法( (续续) )(三) 多变量资产价格的随机模拟与随机数的产生1. 若风险因子不相关按单变量的方法分别模拟每个风险因子变量。2. 若风险因子相关基于Cholesky因子分解法模拟资产价格以及产生随机数。92二、基于二、基于Monte Carlo模拟法的模拟法的计算计算VaR的基本步骤的基本步骤 1.识别风险因子变量，建立资产组合价值与风险因子变量之间的映射关系；2.对风险因子未来变化进行随机模拟，得到各个风险因子变量未来变化的一条样本轨道；3.利用第1步给出的映射关系计算组合价值及组合价值的变化值；4.不断重复第二与第三步，直至达到模拟要求的次数；5.基于损益分布计算置信度c下的VaR。 93三、基于三、基于Monte Carlo模拟法计算模拟法计算VaR的应用举例的应用举例o2016年11月10日，某投资者买入10000股招商银行股票，价格为17.80元/股。o请问：n假定招商银行的股价变动服从几何Brown运动，请用连续型方法模拟出招商银行股价在未来两天内变化的一条样本路径。n下一个交易日，在99%置信度水平下，该投资者的最大可能损失是多大？94三、基于三、基于Monte Carlo模拟法计算模拟法计算VaR的应用举例的应用举例o2016年11月10日，某投资者构造了一个投资组合，将100万元按50%的权重分别投资到了招商银行和新华保险这两个股票。请基于Monte Carlo模拟法分析该投资者在下一个交易日的损益情况。o请问：n该投资头寸的期望值是多少？n在95%置信度水平下，该投资者的最大可能损失是多大？95三、基于三、基于Monte Carlo模拟法计算模拟法计算VaR的应用举例的应用举例第一步 以历史数据估计风险因子的分布特征风险因子协方差矩阵S（美元/英镑）r（%/年）r*（%/年）S3.02E-05-1.89E-05-0.000515r-1.89E-050.0013050.000145r*-0.0005150.0001450.011152风险因子样本均值Srr*0.00055-0.000395-0.01013296第二步 风险因子协方差矩阵的Cholesky分解a0.005495453下三角矩阵Tb-3.44E-03d-9.37E-02c0.035960699e-0.00493042f0.048429394三、基于三、基于Monte Carlo模拟法计算模拟法计算VaR的应用举例的应用举例(续续)97第三步 利用Monte Carlo模拟方法生成三个风险因子的样本第四、第五步 估值并计算VaR风险因子的随机数风险因子未来变化的可能取值远期合约价值和损益的可能取值e1e2e3S（美元/英镑）r（%/年） r*（%/年） （美元）12.35 -0.510.351.6771643154.91E+005.757937795233680.3677140099.367721.711.42-0.371.6736472254.98E+005.773520937201256.9875107675.987531.760.760.551.6739219974.96E+005.816645037201250.0586107669.0586961.992.36-0.311.6533140495.03E+006.118482105-11153.012-104734.01299-2.13 0.340.811.6525446854.95E+006.19580045-24732.5699-118313.57100 -2.39 -0.770.341.6511158674.92E+006.20287341-40656.5788-134237.579三、基于三、基于Monte Carlo模拟法计算模拟法计算VaR的应用举例的应用举例(续续)98四、基于四、基于Monte Carlo模拟法模拟法VaR计算的评述计算的评述 1. 优势：(1) 与历史模拟法相比，该法结果更精确。可以产生大量关于风险因子未来取值的模拟样本；可以最大程度地将风险因子未来变化的各种可能情景模拟出来；不必受到历史数据在数量与质量等方面所存在的种种制约 。 99四、基于四、基于Monte Carlo模拟法模拟法VaR计算的评述计算的评述( (续续) ) (2) 是一种完全估值法，可以处理非线性、非正态问题。 (3) 通过建立随机模型并对模型中相关参数进行估计和修正，从而使得对风险因子变化的模拟更加贴近于现实 。 (4) 可以借助计算机来完成，从而大大提高有效性和精确性。 100四、基于四、基于Monte Carlo模拟法模拟法VaR计算的评述计算的评述( (续续) )2. 不足：(1) 依赖于随机模型以及估计模型参数的历史数据，存在模型风险和参数估计误差；(2) 伪随机数可能导致模拟错误和失效；(3) 收敛速度慢、计算效率低，计算量大；(4) 随机模拟次数不够多时，方差比较大；(5) 对小区间内变化用静态方法处理会产生相应偏差，分割数过少会大大地加剧该偏差。101五、五、Monte Carlo模拟法的模拟法的改进与扩展介绍改进与扩展介绍 1. 下文有时会根据需要将前文介绍的Monte Carlo模拟法，称为传统Monte Carlo模拟法。 2.对传统Monte Carlo模拟法的改进和修正主要体现在三方面：(1)降低伪随机数的集聚性，合理减少风险因子数量，以提高计算效率和收敛速度；(2)降低样本的方差，以提高计算的准确性；(3)引入Markov过程，以降低用静态方法处理时所产生的偏差。 102五、五、Monte Carlo模拟法的模拟法的改进与扩展介绍改进与扩展介绍( (续续) )(一) 对收敛速度和计算效率的改进1. 拟拟Monte Carlo方法方法(QuasiMonte Carlo Simulation)(1) 用预先设定的确定性方法在空间中产生一些低偏差(Low Discrepancy)的拟随机数。(2) 优点：能够更加均匀地分布在间隔域中；收敛速度更快，从而计算精度也更高。 103五、五、Monte Carlo模拟法的改进与扩展介绍模拟法的改进与扩展介绍 (一一) 对收敛速度和计算效率的改进对收敛速度和计算效率的改进( (续续) )2. 情景情景Monte Carlo模拟法模拟法(QuasiMonte Carlo Scene Simulation)(1) 先采用主成分分析方法从众多风险因子中提取少数几个主成分，再进行下一步的模拟。 (2)优点：由于各个主成分的所有可能不同组合的个数有限 ，因此情景Monte Carlo模拟方法会提高计算的速度和效率。 104五、五、Monte Carlo模拟法的模拟法的改进与扩展介绍改进与扩展介绍( (续续) )(二) 对降低样本方差的改进1. 镜像变量法镜像变量法(Antithetic Variables)抽取标准正态随机变量的样本时，取其相反数-为另一个样本。2. 控制变量法控制变量法(Control Variates) 假设有两种金融衍生工具A和B， 如果B 存在着封闭的定价公式，则可以利用工具B的定价误差来对工具A的定价结果进行调整。1053. 重要抽样法重要抽样法(Importance Sampling)通过变换随机样本的概率测度，以适当加大对我们所研究的问题具有重要影响的样本出现的可能性。4. 分层抽样法分层抽样法(Stratified Sampling) 为避免不能抽取到随机变量在某些取值范围内的样本，可以设法将拟要抽取的样本比较均匀地分布在随机变量的取值空间中。五、五、Monte Carlo模拟法的改进与扩展模拟法的改进与扩展介绍介绍 (二二) 对降低样本方差的改进对降低样本方差的改进( (续续) )1065. 矩匹配法矩匹配法(Moment Matching) 在模拟生成某个分布的样本时，可以对已经生成的样本进行一个变换，使得变换之后样本的某些矩与被模拟分布理论上的矩保持一致；然后，再将变换后的样本运用到定价或者是VaR估计中去。 五、五、Monte Carlo模拟法的改进与扩展模拟法的改进与扩展介绍介绍 (二二) 对降低样本方差的改进对降低样本方差的改进( (续续) )107五、五、Monte Carlo模拟法的模拟法的改进与扩展介绍改进与扩展介绍( (续续) )(三) 马尔可夫链Monte Carlo模拟法马尔可夫链马尔可夫链(Markov链链)Monte Carlo模拟法模拟法， (MCMC Simulation)，将Markov过程引入到传统Monte Carlo模拟法之中，从而实现动态模拟的目的，即抽样分布能够随着模拟的进行而不断改变。108第七节第七节基于基于Delta、Gamma灵敏度灵敏度指标的指标的VaR计算计算109引言引言 1. Delta类方法用Taylor一阶展式近似资产组合的价值2. Delta-Gamma类方法用Taylor二阶展式近似资产组合的价值110111一、基于一、基于Delta类方法的类方法的VaR计算计算 1. Delta类方法主要包括：Delta-正态方法Delta-加权正态方法Delta-混合正态方法Delta-GARCH方法 一、基于一、基于Delta类方法的类方法的VaR计算计算(续续)(一) 基于Delta正态方法的VaR计算 1. 基本原理：资产组合的损益近似为 若风险因子x =(x1，xn)的收益率向量服从正态分布 ，则资产组合的损益 也近似地服从正态分布；估计出近似损益分布的均值和方差参数；基于近似损益分布计算VaR。112一、基于一、基于Delta类方法的类方法的VaR计算计算 (一一) 基于基于Delta-正态方法的正态方法的VaR计算计算(续续)2.基于Delta - 正态方法计算VaR的具体步骤：(1) 识别风险因子变量，建立组合价值与风险 因子之间的映射关系；(2) 估计风险因子收益率的协方差矩阵；(3) 计算灵敏度系数Delta； (4) 估计资产组合未来损益的近似分布； (5) 计算VaR。113一、基于一、基于Delta类方法的类方法的VaR计算计算 (一一) 基于基于Delta-正态方法的正态方法的VaR计算计算(续续)3. 假设某美国公司于1998年12月31日持有一份三个月后到期、以16.5百万美元交换10百万英镑的远期合约。基于Delta正态方法的VaR计算持有该远期合约在95置信度下的VaR。 (1) 写出远期价值和风险因子之间的函数关系 ：3个月的美元折现因子； ：3个月的英镑折现因子。 114115一、基于一、基于Delta类方法的类方法的VaR计算计算 (一一) 基于基于Delta-正态方法的正态方法的VaR计算计算(续续) (2) 计算远期合约在未来一天内价值变化的一阶近似式 (3) 估计三个风险因子收益率 S/S、P*/P*以及P/P 的协方差矩阵； (4) 计算风险因子S、P*以及P的协方差矩阵； (5) 根据资产组合价值对每个风险因子的Delta值计算出组合价值变化 f 的标准差； (6) 计算得到在95%的置信度下持有该远期合约的VaR为$98,150.135 。一、基于一、基于Delta类方法的类方法的VaR计算计算 (一一) 基于基于Delta-正态方法的正态方法的VaR计算计算(续续)风险因子收益率协方差矩阵S/SP*/P*P/PS/S1.17E-057.88E-072.90E-08P*/P*7.88E-076.76E-088.83E-10P/P2.90E-088.83E-107.96E-09风险因子暴露向量SP*SP*（-KP）16392392.7216392392.72-16298811.5组合损益的方差3.56E+09组合损益的标准差59665.73556VaR（95%置信度）$98150.135一、基于一、基于Delta类方法的类方法的VaR计算计算 (一一) 基于基于Delta-正态方法的正态方法的VaR计算计算(续续)4. Delta-正态方法评述 (1) 优点：计算简单、操作方便；当风险因子变化很小时，计算误差较小。 (2) 缺陷：不能反映非线性风险；多元正态分布假设不尽合理；使用简单移动平均方法估计收益率协方差矩阵易导致预测结果失真;历史数据长度对协方差估计的可靠性也会产生影响。117一、基于一、基于Delta类方法的类方法的VaR计算计算(续续)(二) 基于Delta加权正态方法的VaR计算1.Delta加权正态模型，又称“RiskMetrics”方法方法：由JP Morgan公司提出；解决Delta-正态方法采用简单移动平均方法估计协方差矩阵时容易出现失真的问题。 2. 基本原理：除采用对历史数据赋权重的方法估计风险因子收益率向量r的协方差矩阵 的方法外，其它原理均与Delta-正态方法相同。118一、基于一、基于Delta类方法的类方法的VaR计算计算(续续)(三) 基于Delta-GARCH方法的VaR计算1.与Delta-正态方法、Delta-加权正态方法比较，Delta-GARCH方法更擅长于刻画：金融时间序列的厚尾分布；金融时间序列的波动性聚集特征。2. 基本原理：除利用GARCH模型来估计风险因子收益率向量 r 的协方差矩阵 外，该法的基本原理、思路和计算步骤完全类似于Delta-正态方法。 119120一、基于一、基于Delta类方法的类方法的VaR计算计算(续续)(四) 基于Delta-EGARCH-GED方法的VaR计算1. 与Delta-正态方法、Delta-加权正态方法和Delta-GARCH方法比较， Delta-EGARCH-GED方法更加擅长刻画：金融时间序列的厚尾分布；金融时间序列的波动性聚集特征；杠杆效应。2. Delta-EGARCH-GED方法用广义误差分布(Generalized Error Distribution，GED)描述具有厚尾特征的资产组合价值的未来变化。一、基于一、基于Delta类方法的类方法的VaR计算计算( (续续) )(五) Delta-混合正态模型1. Delta-混合正态模型用混合正态分布来描述风险因子收益率或者是资产组合价值变化的厚尾分布特征。 121二、基于二、基于Delta-Gamma类方法类方法的的VaR计算计算 1. Delta-Gamma类方法主要包括： Delta-Gamma-GARCH方法 Delta-Gamma-Wilson方法 Gamma-CF方法 Gamma-Johnson方法 122二、基于二、基于Delta-Gamma类方法类方法的的VaR计算计算(一) 基于Delta-Gamma正态方法计算VaR1. Delta-Gamma正态方法正态方法的基本原理：采用资产组合价值的变化或损益 关于风险因子变化向量的二阶Talyor展式作为 的近似； 仍假设风险因子的收益率服从正态分布；仿照Delta-正态方法的实施步骤来进行。123二、基于二、基于Delta-Gamma类方法的类方法的VaR计算计算(二) 其他重要的Delta-Gamma类方法1. Delta-Gamma-GARCH方法方法的基本原理为：采用资产组合价值的损益 关于风险因子变化向量的二阶Talyor展式作为 的近似； 风险因子的收益率变化服从多元GARCH模型。1242. Delta-Gamma-Wilson方法方法的基本思路：将VaR计算归结为一个最优化问题来求解 s.t. 其中， 为标准正态分布下对应于 的分位数。 二、基于二、基于Delta-Gamma类方法的类方法的VaR计算计算 (二二) 其他重要的其他重要的Delta-Gamma类方法类方法1253. Gamma-CF方法方法和Gamma-Johnson方法方法产生的主要背景：若采用资产组合价值的损益 关于风险因子变化向量的二阶Talyor展式作为 的近似，即使风险因子的收益率服从正态分布， 也不服从正态分布，此时计算VaR就比较麻烦。二、基于二、基于Delta-Gamma类方法的类方法的VaR计算计算 (二二) 其他重要的其他重要的Delta-Gamma类方法类方法1264.Gamma-CF方法的基本思想：通过非标准正态分布 的偏度和峰度对标准正态分布的分位数进行调整，从而得到 自身的分位数近似计算公式。5.Gamma-Johnson方法的基本思想：对非正态分布 做一个单调变换，使得变换之后的随机变量服从标准正态分布。 二、基于二、基于Delta-Gamma类方法的类方法的VaR计算计算 (二二) 其他重要的其他重要的Delta-Gamma类方法类方法127三、基于三、基于Hull-White正态变换方法正态变换方法1. Hull-White正态变换方法，最早由Hull和White (1998)提出，核心思想是：(1) 利用变换 将风险因子收益率 rit (非正态随机变量，分 布函数为Git) 变换为标准正态随机变量 fit；(2) 利用Cholesky分解方法，模拟生成关于多元 正态随机向量(f1t , f2t , , fnt)的样本；128三、基于三、基于Hull-White正态变换方法正态变换方法(续续)(3) 利用逆变换 得到风险因子收益率随机向量 (r1t, r2t,rnt) 的样本；(4) 在上述样本基础上计算VaR。129第八节第八节压力试验压力试验130引言引言1. 现实市场中的非正常波动或者极端波动的事件和情景时有发生，金融风险因子或金融资产价值的变化分布往往呈现出明显的“厚尾”特征，此时继续运用经典的VaR方法度量厚尾分布事件的风险将有可能产生较大的估计偏差。131引言引言( (续续) )2.对于厚尾分布厚尾分布，通常有两种理解：(1)一种是与正态分布比较，把峰度比正态分布高的分布称为厚尾分布。包括t-分布、对数正态分布、广义误差分布、混合正态分布等。本文厚尾分布皆是这种意义的。(2)另一种是Ramazan Gencay的定义，即满足 的分布F（x）称为厚尾分布。 按该定义，上述的分布都不是厚尾的。 1323. 压力试验和极值理论是目前度量厚尾分布事件风险的两种基本方法。本节介绍压力试验，下一节介绍极值理论。4. 压力试验压力试验，是在模拟或构造未来可能出现的极端情景的基础上，对极端情景及其影响下的资产组合的价值变化做出评估和判断。5. 压力试验的两种主要方法：(1) 情景分析法(2) 系统化压力试验 引言引言( (续续) )133一、情景分析法一、情景分析法(一) 情景分析法情景分析法是最常用的压力试验方法，主要用于评估一个或几个市场风险因子突然从当前市场情景变化到某些极端情景或事件的过程中对资产组合价值变化的影响程度。各种情景分析法的实施主要包含两个关键步骤：情景构造；情景评估。134(二) 情景构造情景构造1. 市场情景和极端情景：(1) 对影响资产组合价值 P 的n个市场风险因子r1,r2,rn 依次赋值，得到的一个取值组合，称为一个市场情景市场情景，记为 。 (2) 对未来某个极端情景s下的n个风险因子进行评估 赋值，此时的市场情景就是我们所构造的极端情景极端情景，记为 。 一、情景分析法一、情景分析法1352. 情景构造包括：典型情景构造典型情景构造；历史情景构造历史情景构造；VaR情景构造情景构造；特殊事件假定法特殊事件假定法；Monte Carlo情景构造情景构造。一、情景分析法一、情景分析法 (二二) 情景构造情景构造(续续)1363. 典型情景构造典型情景构造是通过对市场风险因子的未来变化进行构造模拟，进而生成极端情景，又称为标准压力情景法标准压力情景法。衍生产品政策集团(Derivatives Policy Group，简称DPG)针对银行业所提出的标准压力情景法已为许多银行所广泛应用。典型情景构造法的优势是可比性。一、情景分析法一、情景分析法 (二二) 情景构造情景构造(续续)1374. 历史情景构造(1)主要根据历史上曾发生过的极端事件来模拟、构造未来的极端情景；或者说选择历史上发生极端事件时风险因子的时间序列数据，来模拟、构造未来的极端情景。(2)分类不考虑风险因子之间相关性的历史情景构造考虑风险因子之间相关性的历史情景构造一、情景分析法一、情景分析法 (二二) 情景构造情景构造(续续)1385.VaR情景构造构造与持有期相对应的各个风险因子的波动性以及风险因子之间相关性的极端情景。6.特殊事件假定法通过设想未来可能发生的一次突发性事件（例如可能发生的自然灾害、突发性的政治事件等）对市场风险因子可能产生的影响，进而生成极端情景。一、情景分析法一、情景分析法 (二二) 情景构造情景构造(续续)1397. Monte Carlo模拟情景构造(1) 基本思想：用Monte Carlo法生成压力情景。(2) 具体做法：假设市场风险因子变化服从特定分布，通过历史数据估计出分布函数中的参数；模拟生成市场风险因子变化的大量样本；根据事先确定的压力情景的判断标准，在已经生成的大量情景中，筛选出压力情景。一、情景分析法一、情景分析法 (二二) 情景构造情景构造(续续)140一、情景分析法一、情景分析法( (续续) )(三) 情景评估1. 情景评估情景评估：对极端情景发生时资产组合价值未来的可能变化进行评估。2. 情景评估的具体做法：根据资产定价理论和映射关系，得到市场极端情景rs下资产组合的价值以及资产组合的当前市场价值；前者减去后者得到极端情景发生时资产组合价值在未来的一个变化量。 141一、情景分析法一、情景分析法( (续续) )(四) 情景分析法的评价1. 优点：可考察市场风险因子的极端变动的影响；VaR情景分析法还能够用来评估市场风险因子波动率和相关系数发生极端变化时对资产组合价值的可能影响；可以使得金融机构的高层和风险管理部门能较为准确地评估和把握极端事件的影响，提高风险管理策略的有效性和可靠性。 142一、情景分析法一、情景分析法 (四四) 情景分析法的评价情景分析法的评价(续续)2. 不足：压力情景的构造往往具有较大的主观性；不可能考虑到所有可能的压力情景；仅能指出损失，而不能给出可能性；根据历史构造的压力情景，未必是未来的最坏情景；基于历史数据构造的压力情景很少考虑到资产组合自身的风险特征。 143二、系统化压力试验二、系统化压力试验(一) 基本原理 1. 系统化压力试验系统化压力试验是在一定条件下，对影响资产组合价值P的风险因子r1，r2，rn，采用数学或者统计的方法生成大量的市场情景，然后评估这些情景对资产组合价值变化的影响，从中搜寻资产组合的最坏情景，即导致资产组合价值损失最大的压力情景。 144二、系统化压力试验二、系统化压力试验 (一一) 基本原理基本原理 (续续)2. 系统化压力试验与情景分析法的主要区别：前者针对一系列不同的压力情景进行压力试验，因此更加彻底、系统化；前者既考虑了风险因子在历史上的极端变动，又考虑到未来潜在的所有可能压力情景，因而本质上是一种前瞻性的情景分析法；前者考虑了资产组合的风险特征在确定其最坏情景过程中的作用。 145二、系统化压力试验二、系统化压力试验( (续续) )(二) 搜寻最坏情景的基本思路和步骤 1. 所谓最坏情景最坏情景（Worst-case Scenarios）是指导致资产组合在未来发生最大损失的市场情景。2. 搜寻最坏情景的基本思路：在一些合理性条件下，搜寻资产组合最小价值对应的市场情景。 (1) 合理性条件合理性条件：在搜寻最坏情景时，不可能保证所有市场情景都会被搜寻，所以必须考虑不同市场情景的合理性，即发生的可能性。合理性条件记为C 。146二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)(2) 满足合理性条件C的所有市场情景的集合，称为允许域允许域（Admissibility Domain），记为D。 (3) 系统化搜寻最坏情景就转化为下面的一个最优化问题：r* 其中P(r1,r2,rn)是资产组合的价值函数；表示允许域D中的一个市场情景；r*指的是该资产组合在给定的合理性条件C下的最坏情景。 1473. 系统化压力试验的一般步骤:(1) 确定合理性条件和允许域；不考虑风险因子之间相关性的情形考虑风险因子之间相关性的情形(2) 搜寻最坏情景。因素推动方法网格搜寻方法Monte Carlo法和拟Monte Carlo模拟法二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)148合理性条件合理性条件C 价值函数价值函数P（r）允许域允许域D最坏情景最坏情景r* 二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)1494.合理性条件以及允许域的确定 (1)不考虑风险因子之间相关性的情形:二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)r1r2150 (2) 考虑风险因子之间相关性的情形: D = r : （r0 - r） -1（r0 - r）T k2 其中为风险因子的变化r - r0 的协方差矩阵。二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)1515. 搜寻最坏情景 (1) 因素推动方法因素推动方法确定市场因子的最不利变化方向；让各风险因子朝最不利的变化方向移动一个给定的数值，使资产组合的价值减少最多；并据此确定资产组合的最坏情景以及最坏情景下的资产组合价值。二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)152(2) 网格搜寻方法网格搜寻方法先把前面得到的“长方体”允许域进一步分割为若干小的长方体，每个小长方体的顶点就对应着一个压力情景；再根据资产组合的价值函数表达式，可计算出上述每个压力情景处的资产组合价值；最后，通过比较不同情景处资产组合价值的大小，可近似得到资产组合的最坏情景。 二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)153r0r1r2二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)154(3)Monte Carlo法和拟Monte Carlo方法确定一个从n维单位长方体到n维空间中的长方体允许域的变换T；采用Monte Carlo或者拟Monte Carlo方法生成N个随机向量，可得到N个压力情景 T(xi)；分别计算所对应的资产组合价值，其中那个最小值所对应的压力情景就是近似的最坏情景。二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)155r变化的幅度r*因子变动-0.1445-0.10837-0.0723-0.0361200.0361250.072250.108370.1445-0.422412112558.8127110.3141655156201.4170741.1185280.2199817214350228881-0.3168169840.684392.0898936.9113483.2128022.9142562157098171632186163-0.2112127145.4441696.9256241.870788.0685327.7399866.81114403128937143468-0.1056-15531.55-980.06613564.828111.0842650.7457189.8371726.386259.81007910-58182.29-43630.8-29086-14539.71.93E-0514539.0829075.643609.158140.40.105603-100812.1-86260.6-71716-57169.4-42629.8-28090.7-13554979.315510.60.211206-143419.7-128868-114323-99777-85237.4-70698.3-56162-41628.3-270970.316809-186005.1-171454-156909-142362-127823-113284-98747-84213.8-69682.50.4224121-228568.5-214017-199472-184926-170386-155847-141311-126777-112246二、系统化压力试验二、系统化压力试验(二二) 搜寻最坏情景的基本思路和步骤搜寻最坏情景的基本思路和步骤 (续续)156三、压力试验方法评述三、压力试验方法评述 (一) 良好特性：压力试验能为风险管理部门提供比VaR模型更多的有价值信息：压力试验对于度量资产波动率的变化所产生的影响也具有重要作用；运用压力试验还可以识别风险因子之间相关性的变化对资产组合风险的影响。157(二) 不足之处压力试验受限于实施压力试验者的主观判断及其经验；全值估值法计算量往往很大，从而需要花费很大的人力、物力。 压力试验没有给出每一种结果发生的可能性。 三、压力试验方法评述三、压力试验方法评述(续续) 158第九节第九节极值理论极值理论159160引言引言1.极值理论主要运用统计（特别是极值统计）理论和方法测度厚尾分布事件或极端事件所导致的风险损失。2.极值理论主要包括：分块样本极值理论分块样本极值理论 对大量独立同分布样本分块后的极大值进行建模POT (Peaks Over Threshold) 对超过给定阈值的样本数据进行建模160161一、分块样本极值理论一、分块样本极值理论(一) 分块样本极值理论的基本思想1.关键概念：X1，X2，Xn，：与总体 X 的独立同分布的非退化随机样本序列，这里表示金融风险因子的随机样本序列；样本极大值样本极大值： 极大值随机样本序列极大值随机样本序列：样本极小值样本极小值：极小值随机样本序列极小值随机样本序列：161一、分块样本极值理论一、分块样本极值理论 (一一) 分块样本极值理论的基本思想分块样本极值理论的基本思想(续续)2.分块样本极值理论的基本思想在于：由于总体 X 的分布函数F(x)与极值 和 的分布函数满足 (3.8.4) (3.8.5) 所以，要获得总体 X 的分布函数F(x)，只要得到极大值 或极小值 的分布函数即可； 分块样本极值理论考察总体X的尾部特征，即x充分大时的 。 162162一、分块样本极值理论一、分块样本极值理论(续续)(二) 极大值的分布估计基于Fisher-Tippett定理1. 由于难获得极大值 的精确分布，所以一般情况下用极大值 的渐近分布极值分布进行近似。 极值分布极值分布的定义：假设存在实数序列cn和dn，n=1,2, cn 0 和某个非退化分布函数 ，使得当 时有： 则称 为一个极值渐近分布，简称为极值分布。同时称总体X的分布函数F(x)处于 的最大吸引场中，记为FMDA(H)。 163163一、分块样本极值理论一、分块样本极值理论 (二二) 极大值的分布估计极大值的分布估计(续续)2. 估计极值分布的预备工作广义极值分布(1) 设： 则称 为广义极值分布(Generalized Extreme Value Distribution)，其中 ；为使 , 均有 ， 应满足 ；当 时， 定义为 。 164164一、分块样本极值理论一、分块样本极值理论 (二二) 极大值的分布估计极大值的分布估计(续续)(2) 广义极值分布的分类按 的取值范围分类： 0，Frchet分布 时，称为标准广义极值分布，记为 。165165一、分块样本极值理论一、分块样本极值理论(续续)3. 基于Fisher-Tippett 定理估计极值分布的基本思路： Fisher-Tippett 定理证明了，样本极大值 的渐 近分布，即极值分布 ，必为标准广义极值分布 ，或者必是标准的Weibull分布、Gumbel 分布和 Frchet 分布中的一种； 估计出 cn、dn和 之后，可以用 进一 步估计X 的分布 F(x) 。166166二、二、POT模型模型(一) POT模型的基本思想1. 分块样本极值理论实际操作中存在的主要问题： 在对极值分布 H(x) 进行估计时经常没有充足的样本数据； 已有的不能用以刻画尾部特征的数据对极值分布拟合的意义不大。167167二、二、POT模型模型 (一一) POT模型的基本思想模型的基本思想(续续)2. POT模型的基本思想：以超过某个充分大阈值u的所有观测值序列为样本，估计总体 X (金融风险因子)的超阈值 u 的尾部分布，即当x充分大时的 ，其中的 即下文定义的超额分布函数。超额分布函数超额分布函数：总体 X 的分布函数为F(x)，阈值 u 为常数，定义F(x)的右端点为 。对阈值 u xF，定义 为总体 X 对阈值 u 的超额分布函数超额分布函数(Excess Distribution Function)。168168二、二、POT模型模型(二) 超额分布函数的估计基于P-B-H定理1. 预备定义 广义帕累托分布： 其中 ；当 0时， ；当 0时有根据Pickands-Balkama-de Haan定理，超额分布函数可以用广义帕累托分布来逼近。 三、三、广义极值分布与广义帕雷托分布对广义极值分布与广义帕雷托分布对厚尾分布的拟合及参数、分位数的估计厚尾分布的拟合及参数、分位数的估计(一) 分块样本极值理论中广义极值分布的估计1.主要估计方法包括：Hill估计Hall试算法三、广义极值分布与广义帕雷托分布对厚尾分三、广义极值分布与广义帕雷托分布对厚尾分布的拟合及参数、分位数的估计布的拟合及参数、分位数的估计(续续) (一一)分块样本极值理论中广义极值分布的估计分块样本极值理论中广义极值分布的估计(续续)2.估计思路： (1) 估计广义极值分布的参数 、标准化或正规化参数cn和dn ； (2) 将上述估计结果代入如下两个式子的右边 从而得到F (x)的尾部估计，及X 的p-分位数 。 三、三、广义极值分布与广义帕雷托分布对厚广义极值分布与广义帕雷托分布对厚尾分布的拟合及参数、分位数的估计尾分布的拟合及参数、分位数的估计(续续)(二) POT模型中的广义帕累托分布的估计1. 估计思路： (1) 估计 根据 Hill 图或超额均值函数 的图形估计阈值u；用超过阈值的样本估计 。三、三、广义极值分布与广义帕雷托分布对厚尾分广义极值分布与广义帕雷托分布对厚尾分布的拟合及参数、分位数的估计布的拟合及参数、分位数的估计 (二二) POT模型中的广义帕累托分布的估计模型中的广义帕累托分布的估计(续续)(2) 估计估计广义帕累托分布的参数 与 ；将上述估计结果代入如下式子的右边(3) 将 和 的估计结果代入如下式的右边 , 从而得到F (x)的尾部估计，以及X 的p-分位数 。 三、三、广义极值分布与广义帕雷托分布对厚尾分广义极值分布与广义帕雷托分布对厚尾分布的拟合及参数、分位数的估计布的拟合及参数、分位数的估计 (二二) POT模型中的广义帕累托分布的估计模型中的广义帕累托分布的估计(续续)176四、基于极值理论四、基于极值理论VaR计算计算1. 利用POT模型计算VaR的步骤： 根据经验超额均值函数确定资产组合的收益率分布是否具有厚尾特征； 根据前文所述思路和步骤，估计用来逼近超额分布函数的广义帕累托分布的参数； 给定充分大的置信度c ，计算资产组合的VaR：176五、基于极值理论计算五、基于极值理论计算VaR方法评述方法评述1. 运用极值理论计算VaR的优点： 可以较为准确地描述厚尾分布的尾部特征； 具有计算分位数或VaR的具体公式，便于操作。2. 弱点： 只适于描述资产组合价值变化分布的尾部特征； 有关历史数据较少，这会加大该法的应用难度，和影响应用效果； 某些条件在现实市场中并不一定能得到满足。177177思考题思考题1. 简述金融市场风险度量方法的演变过程。2. 请简要介绍现代金融市场风险度量的主要方法，并对各种方法的应用步骤、适用范围、优缺点进行讨论。3. 试对久期、凸性以及灵敏度方法中所涉及的灵敏度系数的含义、特性、应用范围、局限性等进行剖析和比较。4. 简述在VaR方法中选择和设定置信度和持有期时应考虑的基本要素。178思考题思考题(续续)5. 简述VaR方法的一般计算步骤并讨论该法存在的局限性。6. 请运用历史模拟法和股价的历史数据，计算中信证券公司的VaR。7. 请运用Monte Carlo模拟法并选择相应的股价历史数据对招商银行进行VaR计算。8. 以本章第五节中远期合约为例，请运用Delta-Gamma方法计算资产组合的VaR。9. 试析压力试验与极值理论方法较传统的VaR方法有哪些异同和改进。179思考题思考题(续续)10. 受金融风暴的影响，从2007年第四季度开始，美国花旗集团连续五个季度亏损，2008年全年亏损达到187亿美元，股价从2007年5月的近50美元狂泻至最低不到1美元（0.97美元），仅2008年11月19日-21日的3天之内就下跌了55%，11月20日单日跌幅达到了26%。为使花旗集团免于破产，美国财政部通过问题资产处置计划（简称为TARP）于2008年10-11月间两次向花旗集团注资合计450亿美元，同时提供了高额的资产担保。试运用压力试验与极值理论，揭示美国财政部向花旗集团注资的有效性和合理性，并给出具体的分析步骤。180181谢谢观看/欢迎下载BYFAITHIMEANAVISIONOFGOODONECHERISHESANDTHEENTHUSIASMTHATPUSHESONETOSEEKITSFULFILLMENTREGARDLESSOFOBSTACLES.BYFAITHIBYFAITH