VAR模型及其在投资组合中的应用.doc

二〇一五 年 七 月
VAR模型及其在投资组合中的应用

内容提要
20世纪90年代以来，随着金融衍生产品市场的迅猛发展，加剧了金融市场的波动，2008年的金融危机使得大量的金融机构和投资者破产，风险管理再一次成为金融活动的核心内容。基于VaR的风险管理理论也在巴塞尔协议II的推广下开始广泛地被金融机构所运用，成为目前市场上主流的风险管理工具。本文将VaR及其延伸概念边际VaR和成分VaR的风险管理理论运用到证券市场的投资组合风险调整过程中，选取能够覆盖多数行业的40只个股构成一个投资组合，运用蒙特卡洛法分别计算投资组合在95%的置信水平和持有期为1天的条件下组合的VaR，以此来分析投资组合的风险分布及单只个股的风险贡献度；同时将VaR运用均值-VaR的组合优化理论确定投资组合的最小VaR投资组合，对比调整前后的损益走势图来说明VaR在投资组合风险调整优化过程中的有效性。

【关键词】 投资组合 风险管理 VaR 均值-VaR 组合优化理论

一、序言
（一）研究背景及意义
20 世纪 90 年代以来，随着世界金融市场在业务范围和产品规模上的急剧扩张，使得世界各国经济体之间的一体化和联动性不断增强，近些年的金融危机在国家之间的传导也更为迅速，往往带来整个行业的衰退和大量金融机构的破产。08 年的全球金融危机最初只是美国房地产市场上的次债危机，但由于涉及大量金融衍生产品如 CDO、MBO 和全球范围内的大量机构投资者，使得次债危机最终演变为全球范围内的金融危机，雷曼兄弟等众多金融机构破产倒闭，全球经济也迅速进入衰退周期。
因此可以总结出：世界经济一体化和联动性的增强在横向上扩大了金融风险影响的范围。对此，以巴塞尔委员会为首的全球金融监管机构开始重新制定金融风险管理标准，风险管理再次成为金融活动的核心内容。尤其对于证券公司、基金公司来说，他们持有的不再是单一的一种资产，而是众多资产组成的一揽子投资组合，如何运用一种有效的风险管理标准全面地衡量组合的风险，成为他们首要考虑的问题，VaR 正是在这种 背景下产生并快速发展起来的。
早期的VaR只是作为一种衡量风险的方式，便于向管理层和决策者汇报，是一种消极被动的运用；随后管理者发现可以运用VaR进行主动的风险调控和绩效评估，为优化资源配置提供依据，此时VaR已经演变成为一种主动的积极的管理策略。 目前，VaR作为风险管理领域的主流工具，广泛地被银行、保险公司、机构投资者、非金融机构及监管层机构所运用，应用的范围不仅限于单个的资产或者项目，还包括投资组合、衍生金融工具如理财产品定价、信用风险的度量等方面。
而我国的资本市场起步晚，但是在规模和数量上却发展迅速。在全球经济联动性增强、我国资本市场开放程度不断加大的趋势下，投资者面临的风险将会更加复杂、国际化、多样化，这对投资者的管理能力和风险控制能力提出了更高的要求。尤其是对于管理资金庞大的基金管理人来说，任何细微的失误都会造成重大的损失。因此，VaR风险衡量法的推广在我国资本市场上具有很大的意义。
首先，对于证券市场上的投资者或是基金管理人来说，随着投资组合中的股票数量逐渐增多，投资者希望了解组合整体的风险水平，VaR作为风险控制依据，基金公司可以为每个交易员设定VaR数额限制，能够有效地约束交易员的过度投机行为，避免一些重大的损失。 同时，VaR 可以作为基金业绩评估标准，在投资活动中风险和收益呈正向关系，高收益往往伴随着高风险，因此目前基金业绩评估指标中不再简单地以收入高低来评价业绩，而是开始将风险因素考虑到绩效评估中，防止基金管理人过度追求高收益而忽略对风险的防范。
（二）文献综述
1． VaR研究现状
关于VaR的研究，最早由 JP.Morgan 推出的 VaR（Value-at-Risk）模型，之后发展成为“信用风险估价”（Credit Value at Risk）模型，主要是在正态分布的假定下用RiskMetrics 计算 VaR。随后其他学者将VaR的风险管理理论运用到投资组合、衍生金融工具如理财产品定价、信用风险的度量等各个方面，并在此基础上延伸出 CVaR、MVaR 等概念。 Rachel Campbell 和 Ronald Huisman、Kees Kodeijk 在 2001 年通过实证研究用历史模拟法和VaR风险管理模型对资产组合进行选择，然后同基于收益率正态分布假设的均值—方差模型资产组合的结果进行对比，得出传统的均值—方差模型会低估风险资产组合所面临的市场风险的结论。
Giuseppe AlesiiT 在 2005 年认为现金流的管理在实值期权的风险管理中具有重要地位，因此将VaR引入到现金流风险管理中，用马尔科夫链的蒙特卡洛模拟法对现金流的净现值建模，定义未来每一个时期的现金流 CFs，从而对代表波动性的 VaR进行估算，在此基础上考虑实值期权的最优决策问题，结论是基于净现值的VaR不仅能够保值，还可以降低操作风险。
Chonghui Jiang, Yongkai Ma, Yunbi An 在 2013 年将VaR的风险管理思想运用到保险策略中，提出了基于VaR的保险组合策略（VBPI）。结合中国保险市场的分析，假定风险资产符合几何布朗运动，通过把 VBPI 策略和传统的买入持有策略（B&H）、固定投资比例 CPPI 策略的对比，用组合收益率表现来说明 VBPI 法的优越性。结果显示在考虑交易成本的条件下 CPPI 策略只能维持最低价值，而 VBPI 能在很大程度上解释组合保护条约的内涵；同时两种策略都能够对冲风险下行带来的损失，保险价值和置信水平越高，则限制风险下行的效果越明显。
在国内，VaR 作为一种新的衡量风险的方法，主要运用在资本市场中。彭寿康在 2003 年利用上证 A 股指数、上证 30 指数收益率，用VaR的历史模拟法对股价指数进行了考察，结果表明我国股价指数收益率存在明显的尖峰厚尾特征，用历史模拟法和 Iosistic 分布模型比较适合度量股价指数的市场风险。
目前，基于VaR度量金融风险已成为国外大多数金融机构广泛采用的衡量金融风险大小的方法。VaR 模型提供了衡量市场风险和信用风险的大小，不仅有利于金融机构进行风险管理，而且有助于监管部门有效监管。
2. 投资组合优化问题的研究现状
投资组合优化理论最早源于马克维茨的组合选择理论，目的在于通过多样化的投资来分散风险。目前学术界以均值-方差组合优化模型为基础，衍生出一系列组合优化模型，如考虑 VaR、CVaR 等因素，在国际上的研究进展有： Gordon J. Alexandera, Alexandre M. Baptistab 在 2002 年就将VaR运用到投资组合选择中，通过对均值-VaR 模型四个方面的研究来证明其优越性。第一，对比了均值-方差模型和均值-VaR 模型的有效前沿的变化；第二，怎么将均值-VaR 与期望效用函数最大求解相结合；对比机构运用方差和VaR分别代替风险时的最优化结果，进行实际经济含义的验证。
Robert J. Elliott、Tak Kuen Siu 和 Alex Badescu 在 2010 年提出了一种基于马尔科夫链主导控制下的 BS 经济考虑下的均值方差组合优化模型的解决方式。他们认为主流的马克维茨的均值方差模型是基于均值和方差这两个静态变量的数学模型，只考虑了单一时期经济内的组合最优化，并且这种假设只有在收益率分布符合正态假定，同时经济体的效用函数是二次函数时才有效。在基于马尔科夫链的模型中连续时间和马尔科夫链的假定暗含着经济体的不同状态，通过分离定理和随即最大化原则，可以放宽马尔科夫链的限制，为均值方差模型提供一种更直接详细的解。
投资组合优化理论在国内主要的运用是结合沪深股市的股票组合，对组合的风险进行衡量，以达到降低风险的目的；或是结合交易费用、卖空等因素下的最有投资策略解，国内主要的研究现状如下：
王波、高岳林在 2008 年将基于VaR的条件风险价值 CVaR 运用到中国沪深两市的组合风险管理中，因为 CVaR 可以度量置信水平下的平均损失，可以很好地处理厚尾问题。在实证中选取沪深股市的 16 只股票构成投资组合，考虑市场不允许卖空和整手买入的约束机制下建立 CVaR 投资组合，运用差分进化法进行求解，通过计算不同收益阈值下的买入量、损失值、收益值、资金投入量等有效地验证了 CVaR的有效性。
蒋翠侠、许启发、张世英在 2013 年提出由于金融资产收益多数具有的非正态性和厚尾分布，同时消费者的效用函数可能是二阶以上函数时，需要考虑更高阶矩的时变特性，为此建立基于多目标优化技术和效用理论的高阶距动态投资模型。实证中通过对全球几个主要股票市场的研究发现：金融市场收益率存在高阶矩、并且具有时变性，对组合投资决策有显著的影响。

二、VaR理论概述
风险管理的首要任务是选定合适的风险度量方法。市场上的风险度量方法很多，主要有资产收益率的标准差σ法、β系数、判定系数R、及在险价值VaR法。而VaR凭借其独特优势成为国际上风险管理的主流方法，下面我们将详细介绍VaR的理论及其优越性。
（一）VaR的定义
VaR（Value at Risk）即在险价值，衡量投资者对某项金融资产在Δt的持有期内，给定置信水平c的条件下，投资组合P的最大损失值是多少，用公式表示为：

其中，ΔP=P−P为在时间内的损益函数；P0为期初价值；Pt为期末价值。
VaR的定义中首要涉及持有期和置信水平的设定。常用的置信水平是99%、95%、90%，风险管理部门会根据自己的风险偏好来选择置信水平。比如社保基金、养老基金等机构对风险比较敏感，就会要求相对较高的置信水平，而如股票型基金等追求高风险高收益的机构，则会选择相对低一点的置信水平。
持有期一般与投资组合波动率的大小呈正比，持有期越长，波动率越大。期货市场及衍生品市场对风险波动比较敏感，适合以每日为周期计算VaR，其他一些期限较长的头寸如养老基金等可以每周作为计算周期。
（二）VaR的计算原理
假定投资组合的期初价值为P0，在Δt的投资期限内收益率为R，则期末价值为P=P0(1+R)。在置信度为c的条件下，投资组合的最低价值P*=P0(1+R*)，R*假定为持有期内的最低收益率。
此时，可以算出投资组合最低价值状态下相对于其均值的风险差值，即为相对VaR，公示表示为：

也可以计算出最低价值相对于期初状态时的风险差值，即为绝对VaR，用公式表示为：

由定义可以看出，求解VaR的实质就是求解一段时期内在一定的置信水平下，投资组合的最低收益值P* 或者最低收益率 R *。
假定投资组合的年收益率 是随机变量，服从均值为μ 和波动率为σ 的分布。同时假定投资组合的年收益率与投资期限无关，则在 的持有期限内，组合的收益率和方差分别 μΔt 和Δt σ2。此时相对VaR为：

所以，只要求出在置信水平c下的R* 或者P* ，VaR 就迎刃而解。在这里，我们要考虑到收益率R* 所服从的分布，分为一般分布和正态分布两种。
（1）一般分布条件下
假定未来投资收益P 的概率密度函数是 f ( p )，在给定置信水平c，投资组合未来的最低价值P* 可以表示为：

在一般分布条件下，对于分布的离散或是连续、瘦尾或是厚尾都没有限制，任何分布状况都适用。
（2）正态分布条件下：
由于正态分布的很多特性可以将VaR的计算大大简化，是一种更优的计算方法。假定标准正态分布密度函数为 ϕ (ε)，投资组合的最低价值，此时一般为负数，则

其中， a ( a > 0)是标准正态分布条件下置信水平 c 的分位数。
同时，推出最低收益率R* 表示为：

则可以计算出持有期限 内的绝对和相对 VaR：

（三）VaR的计算方法
由于VaR的计算涉及到收益率的分布假定问题，目前计算方法主要有两大不同的类别：参数法和非参数法。参数法是假定收益率服从随机独立的正态分布为前提的，包括方差-协方差法和蒙特卡洛模拟法；非参数法则不对收益率的分布做任何假定，主要为历史模拟法。下面将主要介绍下历史模拟法和蒙特卡洛模拟法的主要思想。
1. 历史模拟法
历史模拟法对于市场风险因子的统计分布假定很少，是一种独立的非参数化的方法，主要理想是通过模拟市场风险因子的历史变化来构造投资组合的未来收益的概率分布，进而通过置信水平确定的分位数确定 VaR。主要的步骤如下：
a. 确定风险因子，并用风险因子表示投资组合的价值估价公式；
b. 根据风险因子在过去 T+1 个时期内的历史表现数据，可以计算出 T 个实际价值变化率。然后结合投资者设定的基准日的风险因子的实际值就可以得到下一个持有期 N 个市场风险因子的可能值；
c. 利用估值的公式，用得到的 T 个可能的市场风险因子的数值，计算 T 个下一持有期假定的组合价值，减去当前组合的价值，即为 T 个潜在的损益；
d. 将 T 个损益值按从大到小排列，选择置信水平如 95%，倒数（1-95%）*N即为VaR值。
由于历史模拟法是根据模型设定的市场因子的实际历史价格数据作为样本进行模拟，并以此代替未来的价格分布，因此对数据数量的要求比较高，只适用于大样本数据。
2. 蒙特卡洛模拟法
蒙特卡洛（Monte Carlo）又被称为随机抽样法，其主要的思想是概率原理，当需要判断某个结果发生的概率，或是某个随机变量的期望值时，就通过计算机模拟实验过程，模拟大量实验得到事件发生的频率，或是随机变量的均值。计算时假定 资产价格的变动是一种随机过程，通过计算机模拟一段时间内资产价格的随机变动路径，使得最后产生的价格路径近似等于真实过程，以此计算出未来收益率，进一步计算出 VaR。一般情况下，模拟价格随即路径时假定符合几何布朗运动（GBM），因此是一种参数与半参数相结合的方法。 蒙特卡洛模拟法计算VaR是最有效的，不仅能够处理线性分布，还可以解释非线性价格风险、波动性风险、肥尾及极端事件等。但是由于其生成的随机数列可能是伪随机的，使得结果可能产生误差。
（四）VaR成为风险管理的主流工具
相对于其他的风险管理方法，VaR 最大的优势在于，它能够将金融机构面临的风险敞口在一定的概率下，以直观的数字展现出来，以便于更好的理解分析。同时，基于VaR的衍生指标 CVaR 和 MVaR 能够更加深入地衡量投资组合的风险分布和单只个股对组合风险的影响。目前，随着VaR风险管理历年的转变，它已经成为金融机构资本配置、管理风险的重要依据，也成为监管者提供了更好的风险管理标准。
VaR应用理念的转变：有消极管理到积极管理 随着 2004 年的巴塞尔协议对风险管理提出新的标准以来，各大金融机构开始运用VaR进行风险管理，它的应用理念也在慢慢转变。最早的时候VaR只是作为一种衡量风险的方式，便于向管理层和决策者汇报，是一种消极被动地运用；随后，管理者发现可以运用VaR进行主动的风险调控，进行绩效评估或是为优化资源配置提供依据，此时VaR已经演变成为一种主动的积极的管理策略。
1. 作为风险控制的工具，调整风险头寸结构
VaR 把风险量化和金融机构传统的风险控制方式如总额限制、头寸限制、止损限制等方法结合在一起，通过对预期收益和风险的调整进行风险控制。这一方式不仅适用于金融机构、也适用于一般的公司。 例如在基金公司运作体系中，基金管理人运用VaR风险控制，将会很好地改善基金的风险状况。对于基金管理人来说，绩效决定着基金产品的市场销售状况，越来越多的基金管理人开始选择量化地风险分析模型来衡量基金投资组合中的风险分布。VaR 的延伸指标边际 VaR、成分VaR进一步把组合整体的风险细化到单个成分 股的风险贡献中，成为未来投资组合风险管理的一种趋势。
目前在国际上，一些大投资银行机构如花旗、德意志银行在自身运用VaR风险管理的同时，也开始研发一些风险管理产品，客户可以直接按照自己的条件和需求运用软件计算 VaR。
2. 作为风险预算与资源配置的依据
VaR 在深入分析头寸结构的风险分布后，可以为管理者的投资项目的资产配置提供优化的依据。比如像养老基金、社保基金等对风险要求必要严格的机构，在进行投资之前都要对项目进行风险评估、根据结果合理地分配资金，使得风险和收益达到最优的结构。
目前市场上在涉及到多个资产的投资组合的风险管理中，管理者往往会通过马克维茨的均值-方差组合优化模型计算出最小方差投资组合，以达到风险一定下的收益最大化，或是收益一定的条件下方差最小化。随着VaR广泛运用，考虑将VaR代替方差作为衡量风险下行的标准，形成均值-VaR 优化模型，通过非线性的优化模型，求得最小VaR投资组合，这也将是本文研究的主体部分。
（五）均值-VaR 组合优化模型
在马克维茨的均值-方差模型中，用方差来刻画投资组合的风险，但是只能说明投资组合风险的偏离程度，并不能准确地说明偏离的方向，而投资者最关心的是资产价格向下波动时带来的风险损失。而VaR正是基于概率原理刻画了资产价格下行带来的最大损失值，因而可以用VaR代替方差来描述投资组合的风险，由此形成了均值-VaR 组合优化模型：

其中： 是给定的目标期望收益率; 是标准正态分布对应的置信水平c下的分位数。
将VaR的表达式化简为：

带入到均值-方差模型确定的有效前沿中：

推出均值-VaR 模型的有效前沿方程:

三、实证数据
在介绍完论文涉及的VaR理论和计算方法之后，我们开始进入实证分析部分。具体的思路是将VaR的理论运用到证券市场的股票投资组合风险管理中，通过分析VaR在投资组合风险分析中的应用来说明VaR在现实市场中的适用性和有效性。具体过程如下：首先通过参数法和蒙特卡洛模拟法来计算股票投资组合的整体VaR、成分VaR及边际VaR，以此来作为依据深入分析投资组合的风险分布；其次运用均值-VaR投资组合优化模型来计算最小VaR投资组合，并通过对比原始权重和最小VaR权重下的损益走势图来说明VaR在投资组合风险结构调整过程中的有效性，对实际市场操作尤其是基金管理人对基金组合的风险管理具有重要的指导意义。首先在这一章我们将对实证数据进行简单地统计分析。
为了使我们的实证过程能和实际投资活动结合得更紧密，我们将选取股票市场上的开放式基金实际持有的部分成分个股作为我们的证券投资组合（以下简称投资组合），构成比例依据基金公告公布的持仓比例为准。
本文选取的样本是开放式基金—000011华夏大盘精选，基金的管理人是华夏基金管理有限公司，投资目标是追求基金资产的长期增值，属于稳健性的投资风格。华夏大盘精选是我国证券市场上第一支开放式基金，具有广泛地代表性。同时由于华夏大盘精选的成分股票有100多只，为了便于我们的分析，我们选择2015年3月31日基金半年报中公布的持仓明细中前40只股票构成一个投资组合，其权重将按照基金中的权重为准，经过调整后折合为投资组合的新权重。
样本数量的时间统计区间选为2013年7月1日-2015年6月30日的交易时间，节假日及周末不计算在内共有487个交易数据。置信水平设c为95%，持有期设ΔT为1天，我们立足于在2015年6月30日，目标是为了计算投资组合7月1日的VaR值。
表1：2015 年6月30日华夏大盘投资组合头寸结构
代码
名称
持股数（万股）
市值（万元）
权重占比
板块名称
600256
广汇股份
3107.96
73907.27
13.64%
房地产业
600068
葛洲坝
4000.68
47968.17
8.85%
建筑业
600086
东方金钰
1620.68
31959.75
5.90%
其他制造业
600050
中国联通
6000.01
31500.03
5.81%
信息技术业
000822
山东海化
2790
29574
5.46%
石油、化学、塑胶、
600019
宝钢股份
4500.36
27137.14
5.01%
金属、非金属
600316
洪都航空
800.72
25807.09
4.76%
机械、设备、仪表
600831
广电网络
2000.97
21370.35
3.94%
传播与文化产业
000659
珠海中富
1700.01
17952.08
3.31%
石油、化学、塑胶、
000888
峨眉山Ａ
895
16289
3.01%
社会服务业
600107
美尔雅
1050.06
15666.93
2.89%
纺织、服装、皮毛
000090
深天健
1000
11550
2.13%
建筑业
601328
交通银行
2000
11080
2.05%
金融、保险业
600657
信达地产
1700.01
10829.03
2.00%
房地产业
600962
国投中鲁
621.19
10877.05
2.01%
食品、饮料
600683
京投银泰
1200
10152
1.87%
房地产业
000069
华侨城Ａ
1200.01
9492.06
1.75%
社会服务业
002033
丽江旅游
459.99
9264.22
1.71%
社会服务业
000046
泛海建设
1000
8780.03
1.62%
房地产业
000002
万科Ａ
1000
8450
1.56%
房地产业
600499
科达机电
472.58
7566.01
1.40%
机械、设备、仪表
000718
苏宁环球
900.91
7207.29
1.33%
房地产业
600511
国药股份
450.01
6831.09
1.26%
批发和零售贸易
000978
桂林旅游
600.92
6646.22
1.23%
社会服务业
600000
浦发银行
650
6396
1.18%
金融、保险业
600597
光明乳业
800.07
6344.59
1.17%
食品、饮料
600392
太工天成
250.08
6104.49
1.13%
信息技术业
600239
云南城投
455
5851.33
1.08%
房地产业
600997
开滦股份
350
5586
1.03%
采掘业
600811
东方集团
800.09
5504.6
1.02%
综合类
002238
天威视讯
300.03
5373.45
0.99%
传播与文化产业
000430
ST张家界
650
5342.96
0.99%
社会服务业
600325
华发股份
450
5049
0.93%
房地产业
600812
华北制药
412.66
4753.81
0.88%
医药、生物制品
600246
万通地产
946.67
4733.33
0.87%
房地产业
600075
新疆天业
400
4707.95
0.87%
石油、化学、塑胶、
000024
招商地产
250.01
4575.12
0.84%
房地产业
000709
河北钢铁
1000
4520
0.83%
金属、非金属
600135
乐凯胶片
288.25
4511.09
0.83%
石油、化学、塑胶、
601166
兴业银行
336.49
4535.95
0.84%
金融、保险业
由统计数据可以看出投资组合基本涵盖了股票市场上的所有行业，主要侧重于房地产行业、建筑行业、传统制造业、服务行业、信息技术产业及金属行业，能够很好地反映市场整体风险状况。这样涉及大多数行业的资产配置也可以反映出投资组合的管理理念是稳健的风格，通过不同行业的资产配置在追求大盘指数的平均收益的基础上追求资本增值。

四、基于投资组合的实证研究
在完成对实证数据的收集整理工作后，下面将开始分析VaR在组合风险结构调整过程中的应用。首先分别运用参数法和蒙特卡洛模拟法分别计算投资组合的VaR，运用失败率检验法验证计算结果的有效性；然后运用均值-VaR优化模型求出最小VaR投资组合，对比最优组合下的权重和原始权重的损益走势来说明VaR对投资组合风险调整过程的指导意义。
（一）Monte carlo法计算投资组合VaR
参数法下我们使用的是算数收益率，而在蒙特卡洛法先我们使用对数收益率，
在收益率很小的条件下，两者是等价的关系。在对数收益率的假定下：
表2：蒙特卡洛法计算结果统计
序号
代码
名称
单个VaR
序号
代码
名称
单个VaR
1
600256
广汇股份
-22990958.923
21
600499
科达机电
-2884481.633
2
600068
葛洲坝
-19890000.239
22
000718
苏宁环球
-1611511.639
3
600086
东方金钰
-12228408.990
23
600511
国药股份
-2255484.396
4
600050
中国联通
-8581608.962
24
000978
桂林旅游
-2072020.039
5
000822
山东海化
-16437688.195
25
600000
浦发银行
-1566248.529
6
600019
宝钢股份
-5320651.498
26
600597
光明乳业
-1759107.925
7
600316
洪都航空
-9528867.410
27
600392
太工天成
-2051708.431
8
600831
广电网络
-8654294.471
28
600239
云南城投
-2872603.695
9
000659
珠海中富
-4588114.030
29
600997
开滦股份
-1636284.060
10
000888
峨眉山Ａ
-4399413.895
30
600811
东方集团
-1438279.483
11
600107
美尔雅
-7628899.166
31
002238
天威视讯
-1873771.928
12
000090
深天健
-4209944.636
32
000430
ST张家界
-1971028.100
13
601328
交通银行
-1547954.990
33
600325
华发股份
-1568263.141
14
600657
信达地产
-4181343.028
34
600812
华北制药
-1378475.909
15
600962
国投中鲁
-4752731.141
35
600246
万通地产
-1108840.603
16
600683
京投银泰
-4629536.271
36
600075
新疆天业
-1416442.126
17
000069
华侨城Ａ
-2732503.399
37
000024
招商地产
-1616876.737
18
002033
丽江旅游
-3345336.756
38
000709
河北钢铁
-1383434.609
19
000046
泛海建设
-1384373.604
39
600135
乐凯胶片
-1526080.021
20
000002
万科Ａ
-1932464.483
40
601166
兴业银行
-1142947.557
注：在持有期 ΔT 为一天,c=95%

以上我们运用两种方法来计算投资组合的VaR，参数法主要是运用算数收益率来计算VaR，而蒙特卡洛则运用的对数收益率来计算，两者之间的不同在于是否考虑时间的连续性。结果显示截止到2015年6月30日为止，投资组合的总市值为54.17亿元，在持有期为1天，置信水平为95%的条件下，蒙特卡洛模拟法计算出的组合在7月1日有5%的可能性发生1.15亿的损失，占总市值的比重为2.12%；参数法计算的的投资组合在7月1日的VaR为1.24亿，占总市值的2.29%，两种方法相差0.17%，因此我们认为计算的VaR是有效的。
（二）参数法下组合VaR与历史实际损益对比-失败率检验法
以上我们是通过两种不同的方法计算投资组合的VaR，由于股票的收益率并不非足够小，因而两种方法计算的结果相差1.17%，在可接受的误差范围内。为了更精确地衡量VaR的有效用，下面我们将采用失败率检验法来验证。失败率检验法主要是运用历史数据来验证VaR对实际损失的覆盖程度，通过衡量一定置信水平下，实际损失超过VaR的频率与模型设定的上限值是否接近。如果失败率在置信度1-C附近则可认为模型有效；反之则模型误差比较大。我们基于统计的487个历史日期，在每个时间点根据历史数据计算出未来持有一天的VaR值，将会有487个VaR的预测值；同时计算历史数据期内的每天的实际损益值，将两者之间的关系用图像表示。结果显示：实际损益超过VaR值得下限次数为27次，占总观察周期的6.1%；而实际损益超过VaR上界的次数为8次，占总观察周期的1.8%。我们所设的置信水平是5%，因此可以认为我们模型估计的VaR结果是有效的。

图1：投资组合历史真实损益与VaR带对比图
（三）基于均值-VaR 约束下的最小VaR投资组合
通过以上的分析可知，可以通过CVaR衡量成分股对投资组合的风险贡献度，还可以根据成分股MVaR来调整个股权重，进而降低投资组合的整体风险水平。但是随着投资组合中成分股数量的增多，成分股之间的相关性对组合的风险和收益具有很大影响。如何在考虑相关性和MVaR的条件下求得最小的VaR投资组合，下面我们将运用均值-VaR组合优化模型来计算。4.5.1 计算投资组合的最小VaR投资组合。
均值-VaR 优化模型的目标是在期望收益率给定的条件下，使得组合的VaR最小的投资组合，在不同的期望收益率下会形成投资组合的有效前沿。在下面饿实证中，我们用样本数据历史表现的平均收益率来表示期望收益率，模型用公式表示为：

其中，E是给定的目标期望收益率，根据给定的期望收益率就可以求得最优投资组合权重向量x及相应的VaR；Ζc是标准正态分布对应的置信水平C下的分位数。
由2.6节的分析结果可以得出有效前沿为：

结合历史数据，E为我们给定的历史收益率的均值，即为我们设定的期望收益率，Ζc为95%的置信水平下正态分布的分位数1.648，持有期ΔT为一天，计算结果如下：
a=8.27e+004，b=-47.08，c=0.4336，d=3.3662e+004，=3.23%，VaR=2023400
最终求得的权重x是个40*1的矩阵，x矩阵表示为：
X=（ -7.77%,2.09%,-4.79%,5.29%,7.76%,-6.47%,-7.92%,15.39%,-14.56%,29.61%,-2.39%,-27.66%,113.36%,-5.72%,11.05%,6.66%,-12.69%,-2.63%,7.78%,18.4%,16.09%,13.45%,-3.65%,-10.00%,-28.49%,30.87%,9.91%,-6.68%,-6.43%,3.94%,-3.33%,-7.53%,7.41%,2.74%,-17.93%,-9.77%,2.22%,8.14%,-34.92%）

图2：均值-VaR模型投资组合有效前沿

将调整前后的权重比进行对比结果如下：
表3：投资组合优化前后权重变化
序号
股票
调整前x
调整后x
序号
股票
调整前x
调整后x
1
广汇股份
13.64%
-7.77%
21
科达机电
1.40%
16.09%
2
葛洲坝
8.85%
2.09%
22
苏宁环球
1.33%
13.45%
3
东方金钰
5.90%
-4.79%
23
国药股份
1.26%
-3.65%
4
中国联通
5.81%
5.29%
24
桂林旅游
1.23%
-10.00%
5
山东海化
5.46%
7.76%
25
浦发银行
1.18%
-28.49%
6
宝钢股份
5.01%
-6.47%
26
光明乳业
1.17%
30.87%
7
洪都航空
4.76%
-7.92%
27
太工天成
1.13%
9.91%
8
广电网络
3.94%
15.39%
28
云南城投
1.08%
-6.68%
9
珠海中富
3.31%
-14.56%
29
开滦股份
1.03%
-6.43%
10
峨眉山Ａ
3.01%
29.61%
30
东方集团
1.02%
3.94%
11
美尔雅
2.89%
-2.39%
31
天威视讯
0.99%
-3.33%
12
深天健
2.13%
-27.66%
32
ST张家界
0.99%
-7.53%
13
交通银行
2.05%
113.36%
33
华发股份
0.93%
7.41%
14
信达地产
2.00%
-5.72%
34
华北制药
0.88%
2.74%
15
国投中鲁
2.01%
11.05%
35
万通地产
0.87%
-17.93%
16
京投银泰
1.87%
6.66%
36
新疆天业
0.87%
-9.77%
17
华侨城Ａ
1.75%
-12.69%
37
招商地产
0.84%
9.17%
18
丽江旅游
1.71%
-2.63%
38
河北钢铁
0.83%
2.22%
19
泛海建设
1.62%
7.78%
39
乐凯胶片
0.83%
8.14%
20
万科Ａ
1.56%
18.40%
40
兴业银行
0.84%
-34.92%
来源：Matlab结果统计

在置信水平为95%的条件下，均值-VaR优化结果显示：调整后最优组合的VaR为2023400，相比调整前的124234968，风险水平大大降低，说明新的权重分配有效地降低了组合的整体风险水平。
（四）权重调整前后投资组合损益对比分析
对比最优组合点和原始权重时，由于允许卖空操作，因而权重分配具有很大的变化，为了更好地说明最优投资点对组合风险管理的优势。我们基于2015年6月30日到2015年10月20日共75天地交易数据，分别在原始权重和最优权重的条件下计算每日的损益。将计算结果绘制成表显示如下：
图3：原始权重和最优权重的条件下每日损益比较表

由上图可以发现，粗线为最小VaR投资组合权重下的损益曲线，细线为原始权重下的损益曲线，粗线的波动性明显小于细线，即最优投资组合能有效地控制组合的风险水平，表现为损益曲线的相对平缓，而原始权重下的损益具有很大的波动性。主要原因是由于模型中的VaR主要依据方差为基础计算而来，VaR的最优化相当于方差的最优化选择，因此最优组合下组合的波动性比较平缓。
同时我们注意到在2015年7月22日、8月9日及9月6日的交易日附近，随着大盘横盘后开始向下调整，个股进入普跌阶段，投资组合收益率陡然下降，比较调整前后的损益水平，可以明显看出最小VaR投资组合有效地限制了损失程度，而此时原始权重会给投资组合带来大幅的损失