计算方差-协方差矩阵.ppt

lecture 10 FINANCIALMODELING 金融建模 第10章计算方差 协方差矩阵 要计算有效投资组合 我们就必须计算股票收益数据的方差 协方差矩阵 本章中 我们将讨论在Excel中怎样实现这个计算 其中最显而易见的计算为样本方差 协方差矩阵 这是直接由历史收益计算而得的矩阵 我们介绍几种计算方差 协方差的方法 包括在电子表中用超额收益矩阵直接计算 VBA实现该方法计算 即使样本方差 协方差矩阵看起来像一个很明显的选择 但我们将用大量的文字说明它也许不是方差与协方差最好的估计 样本方差 协方差矩阵有两个不尽人意的缺陷 一是它常使用不现实的参数 二是它难以用于预测 这些将主要在10 5和10 6节中讨论 作为样本矩阵的替换 10 9和10 10节将讨论用于优化方差 协方差矩阵估计的 压缩 方法 在开始本章之前 你应先阅读第34章数组函数的内容 里面有一些Excel函数 其参数是向量和矩阵 它们的实施与标准Excel函数略有不同 本章重点讨论这些数组函数Transpose 和MMult 还有 自制 的数组函数的使用 10 1引言 我们用我们的数字例子来说明计算方差 协方差矩阵的矩阵方法 我们通过减去资产各自的平均收益 得到超额收益矩阵 接下来的电子表中的42 52行 在55 61行中我们计算样本方差 协方差矩阵 10 2 1一个稍微更有效率的替代方法正如你所期望那样 的确存在其他计算方差 协方差矩阵可选方法 这里所介绍的方法跳过了超额收益的计算 并且直接使用单元格B71 G76中的公式进行计算 它通过使用数组函数 MMULT TRANSPOSE B23 G33 B35 G35 B23 G33 B35 G35 10 通过写入B23 G33 B35我们直接将每项收益减去均值得到超额收益向量 10 3我们应该除以M还是M 1 Excel与统计量 在前面的计算中我们除以M 1而非M 以此得到无偏的方差和协方差的估计 不过这个选择看起来几乎没有多大影响 我们引用主流的教科书 对于为什么要用M 1取代M这儿有一段很长的历史 如果你从来没有听说过 你可以参考任何一本好的统计教材 这里我们主要想提醒你 如果你在计算一个分布的方差时 这个分布存在已知的先验的均值 而不需要从历史数据估计的时候 那么M 1应该变回M 我们同样想说关于在分母上用M 1替代M上 我们认为对你是已知的 但这却是对你不负责任的 例如 试图用图例说明去证明一个充满疑问的假设 Excel本身某程度上在除以M还是M 1这个问题上也有些混乱 在下面的电子表中我们给出几种计算均值 方差 标准差和协方差的方法 Excel区分总体方差 Varp 除以M 样本方差 Var 除以M 1 以及总体和样本标准差 分别为Stdevp和Stdev 但是Excel并没有在协方差函数Covar中作此区分 你可以看到B30中Covar除以M 和Varp一样 如果你想得到一个相应的除以M 1的协方差函数 那么你得像单元格B33那样用Covar乘以 或者你需要使用像单元格B32那样的数组函数 MMULT TRANSPOSE B3 B13 B16 C3 C13 C16 10 如果Excel是完全合理的 它应该有两个函数 Covarp 它除以M 对应Varp或Stdevp 以及Covar 它除以M 1 对应Var或Stdev 困惑了吗 没关系 正如该部分开始的教科书引用指出的那样 它不是一个至关重要的问题 10 4计算方差 协方差矩阵的其他方法在这一节中 我们介绍两种替代计算方差 协方差矩阵的方法 第一种是使用一个VBA数组函数 它可以直接计算出样本方差 协方差矩阵 第二种是使用Excel的Offset函数 10 4 1一个计算方差 协方差矩阵的VBA函数我们的第一种替代方法是用一个VBA函数 我们的第一种替代方法是用一个VBA函数 FunctionVarCovar rng Range AsVariantDimiAsIntegerDimjAsIntegerDimnumColsAsIntegernumCols rng Columns CountDimmatrix AsDoubleReDimmatrix numCols 1 numCols 1 Fori jTonumColsForj 1TonumColsmatrix i 1 j 1 Application WorksheetFunction Covar rng Columns i rng Columns j NextjVarCovar matrixEndFunction 这个函数是一个数组函数 意味着它必须使用 Ctrl Shift Enter 下面有一个例子 应注意以上6只股票的GMVP包含两个空头 BA和MSFT 且在GE及IBM上有一个非常大正向头寸 这是一个潜在的不主张使用样本方差 协方差矩阵计算GMVP的例子 一个寻找最小方差的投资者是不会将77 的BA和10 的MSFT的卖空得资金 投放到投资组合的61 的GE和100 的IBM上去的 这些用样本方差 协方差矩阵计算产生的不合理的投资组合使得我们必须寻找各种其他方法来计算该矩阵 它们将在10 7 10 9节中讨论 但在进入这些主题之前 我们还是先讨论有效投资组合的计算 该电子表中的行16给出了每只股票在数据的最后一天的市场价值 从2004年1月开始 行17计算出以每项资产的市场权重为投资比例的投资组合 市场权重和有效投资组合权重之间有显著差别 该电子表显示了另外两个利用样本方差 协方差矩阵优化配置投资组合所产生的问题 第一个问题是 这个最优的投资组合包含了一些非常大的卖空和一些的不切实际的非常大的多头 在运行模拟运算表计算不同c值下的包络线投资组合时我们发现 整条有效前沿均含有非常大的空头和多头投资组合 被标注的部分显示出由c变化而导致的变动 只有JNJ在所有包络线组合中保持着正的权数 如果我们认为通过这样的方式获得投资组合的比例是不合理的 那么在该优化过程中 不是方差 协方差矩阵出了问题 就是均值向量出了问题 甚至两者都是 我们将在本章的10 7 10 10节和第13章的Black Litterman模型中再回到该主题上 这个例子揭示的第二个问题与该样本方差 协方差矩阵不像真实的相关系数有关 在接下来的电子表中我们给出相关系数矩阵 在调查期内 最大的相关系数 通用汽车公司和波音公司之间 为0 89 检查该矩阵显示存在一些非常大且不像真实的相关系数 有6个相关系数大于0 5 这是难以置信的 最小的相关系数 波音公司和家乐氏之间 为 0 1 这也许也是一个问题 当波音公司的收益上升时 谷类食品销量会下降 使用Excel中的数组函数 见第35章 我们可以计算SIM方差 协方差函数 这个数组函数需要被单独提取出来 它包含了一个If函数 我们知道它包括三部分 IF condition answerifconditionisture answerifconditionisfalse 条件B24 G24 A25 A30表示在行B24 G24的 跨越输入 它包含了一行股票名字 是和列向量A25 A30的输入相等 为了让大家更清楚地看到这个条件 我们在下面的电子表中插入了公式IF B24 G24 A25 A30 1 0 可以发现它生成了一个对角线上为1 而其余为0的电子表 首先 在我们的SIM方差 协方差矩阵中的第一个IF条件为IF B24 G24 A25 A30 B18 G18 这说明了假如我们在对角线上 我们应该放在资产的收益方差中 其次 IF条件说明当不在对角线上时 我们应该输入MMULT TRANSPOSE B20 G20 B20 G20 H18 公式MMULT TRANSPOSE B20 G20 B20 G20 生成了一个 矩阵 该公式用S P500的方差乘以该矩阵 在H列给出 在电子表中我们通过该6只股票样本相关系数均值来确定该常数相关系数 单元格B19 当然这是允许变化的 假如我们觉得股票平均相关系数在未来是0 3 那么可以从收益数据中直接估计出该方差 协方差矩阵 目前很少理论涉及怎么选择合适的收缩估计 我们建议选择一个收缩因子使得GMVP所有权数均为正 详见下一节 GMVP的平均收益为12 73 单元格F17 标准差为7 73 单元格F19 有一个更加简洁的计算方法 这是一个更加简洁的方法 它借用了数组函数IF A14 A19 A14 A19 1 0 来计算1的向量 10 10 2用替代方差