第7章指数模型及其Python应用.doc

第7章 指数模型及其Python应用【本章精粹】本章针对均值方差模型所存在的缺点，介绍单指数模型，并讨论指数模型环境下是如何分散风险的，最后讨论指数模型的证券特征线的估计及其Python应用7.1 单指数模型在马柯维茨的均值-方差模型的讨论中，各资产间的协方差我们可以作任何假定，它们可以是由资产间存在的任意数量和种类的关系产生，而且在计算风险时所用的公式中，我们必须对所选择的资产间的协方差进行估计如果资产数目太大，我们就必须进行大量的协方差估计，使得在计算任一给定投资组合的方差时，需要花费大量时间在，公式中，如果投资者考虑的是由n种资产构成的组合，那么在求解有效资产组合时，需要掌握三个方面的基本数据：(1) 每一资产的平均收益率，共需n个；(2) 每一资产收益方差，共需n个；(3) 每一对资产之间的相关系数，共需n*(n-1)/2个总计需要2n+n*(n-1)/2个基础性数据对于每天追踪30～50种股票的投资机构来说，每天需要处理495～1325个数据；对于每天追踪150～250种股票的投资机构来说，每天需要处理11475～31625个数据；显然，这对各种投资者来说都是一件非常耗时的事情。

那么，如何使投资组合理论和方法有效实用、简便易行、真正为金融财务工作者服务，就成了金融财务经济学家极为关心的问题单指数模型能帮助我们克服这一困难，使得确定投资组合的方差计算过程变得简单在股票市场中，我们发现，当市场投资组合(如股票市场指数)的收益率显著上升或下降时，几乎所有股票的收益率都随之上升或下降虽然可能有一些股票的收益率比另一些股票的收益率上升或下降得要快，但总的来说都是呈相同趋势变化这意味着，市场投资组合收益率的变化能充分反映各种资产的共同变化趋势因此对各个资产收益率之间的协方差的计算，可以用每一资产收益率与市场投资组合收益率之间的协方差代替单指数模型就是在假定资产的收益率只受市场投资组合即单指数收益率的影响下确定投资组合的权重设资产的收益率具有简单的线性结构，即其收益率r和市场投资组合收益率rM具有关系式其中，为待估参数，为残差假定市场中有n种资产，则按上述结构，第i种资产的收益率满足，i=1,2,…,n；t=1,2,…,N在单指数模型的讨论中，假定影响各个资产收益率的因素有两类：第一类为宏观因素例如通货膨胀率、主要利率的变化、就业率等，在任何情况下，这些因素的影响都是相当大的，几乎所有企业、所有公司都不同程度地受到它们的影响，会引起资产价格总体水平的变化，再通过市场的推动，会影响到市场投资组合收益率水平，进而影响到各资产的收益率。

因此宏观因素影响整个市场的收益率第二类为微观因素例如一种新产品的推出或老产品的淘汰、局部地区或一个公司主要领导的变化，它们都只对个别企业或公司产生影响而不会影响到市场投资组合的收益率，从而使个别资产的收益率偏离市场特征线，出现残差所以微观因素仅影响个别资产的收益率其他类型的因素在单指数模型中不予考虑例如行业因素，某些事件对某一行业内的所有企业产生影响，但却不足以影响到整个经济形势或市场投资的收益率虽然这类因素也能引起残差，但我们假定残差只由微观因素所致从而我们有如下假设，对资产i,j=1,2,…,n，有同时我们还假定 (7-1) (7-2)式(7-1)说明在任一时期残差可能为正，也可能为负，但期望值为零式(7-2)说明资产残差与市场投资组合收益率不相关，即它与市场投资组合是多头或空头(销售方)无关，不因为市场投资组合为多头(购入方)而成正值，也不因为市场投资组合为空头而为负值由单指数模型结构假设和以上各项假设有 (7-3)由和式(7-3)可得即：这是指数模型的另一种假设，即任意资产的收益率由期望收益率和非期望收益率组成在下一章的套利定价理论假设中，我们要将这里的m替换成F (7-4) (7-5) (7-6)从而 (7-7)式(7-3)给出了资产i的特征方程，式(7-7)表明特征方程中的系数即模型结构中的系数恰好为资产i的风险系数。

式(7-4)给出了资产i收益率的方差，它刻画出了资产i的风险，式(7-4)右边的第一项称为资产投资的系统风险可以看作是与整个市场组合有关的风险它是由市场投资组合中各资产的风险共同作用产生的，是所有资产无法避免的风险式(7-4)右边第二项称为残差方差或非系统风险，可以看作是由微观因素所带来的风险，它仅影响到个别资产，是可以通过投资组合而消去的风险因此式(7-4)表明：资产总体风险=系统风险+非系统风险另外，系统风险本身是两项之积，第一项是资产的因子，它表示资产收益率随市场投资组合的变动而受影响的程度，第二项是市场投资组合收益率的方差，表示市场投资组合收益率的变化幅度第二项非系统风险，即残差方差，表示资产收益率由于偏离了特征线而引起的那部分方差的大小在单指数模型的假设下，资产收益率的总体方差来自两部分：一部分是特征线的变动(即系统风险)，另一部分是各点偏离特征线的程度(即非系统风险)下面考虑在单指数模型下投资组合的结构设满足单指数模型的n种资产的投资组合，则投资组合仍有单指数结构： (7-8)简写为： (7-9)由和式(7-1)，式(7-2)，有 (7-10)在单指数模型下，式(7-8)表明投资组合仍具有同类的单指数结构，式(7-9)表明投资组合的因子为各资产因子的加权平均，而式(7-10)表明投资组合的方差(风险)与单种资产类似，仍由两部分构成，第一项是由市场投资组合方差反映的系统性风险，第二项反映的是组合中各资产非系统风险的加权平均(以为权重)。

通过以上讨论，在单指数模型下，马柯维茨组合投资模型为 (7-11)根据上面的公式可知，利用单指数模型进行资产组合，所需要的估计量如下：(1) n个市场风险敏感测度；(2) n个独立的风险指标；(3) n个与市场指数无关的平均收益率；(4) 1个市场组合平均收益率；(5) 1个市场组合风险指标总计需要3n+2个基本数据这样，对于每天追踪30～50种股票的投资者，每天需要收集处理92～152个数据；对于每天追踪150～250种股票的机构投资者来说，每天仅需要收集处理452～752个数据即可这与马柯维茨组合投资模型相比，该模型所需要估计的数值大为减少，它只需要估计各资产的值、值、残差方差及市场投资组合的预期收益率和方差，这比估计各资产之间的协方差的工作量少一个数量级但该模型的精确度不如马柯维茨组合投资模型，它依赖于各资产收益率的单指数结构假设的合理性7.2 指数模型与分散化由(为超额收益率，即)可知：现在要说明的是随着投资组合数量的增加，由非市场因素引起的投资组合风险变小了，而市场因素不变以等权重投资组合为例，权重，则将与上式比较，发现：非市场成分的敏感度：投资组合对市场成分的敏感度：零均值变量： 其中取决于和，不受投资组合分散化的影响。

是公司特有成分方差的平均值，当n很大时，趋于0总之，随着分散化程度的增加，投资组合的总方差会接近系统风险7.3 指数模型的证券特征线估计的Python应用给定ABC公司超额收益率的若干历史样本，将运用于ABC公司就可以得到如下的回归方程：还可以得到此回归方程的拟合度式中：——公司ABC的系统风险；——公司ABC的非系统风险(公司的特有风险)例：有表7-1的数据样本表7-1 市场组合与证券1，2，3，4数据样本日 期市场组合证券1证券2证券3证券4rMtr1tr2tr3tr4t10.01234240.01198920.02357320.02005680.00609762-0.046799-0.037532-0.02174-0.010379-0.01069530.03502080.00560850.00911580.02535790.02127744-0.007361-0.007861-0.00546-0.011015-0.0289885-0.008848-0.025115-0.064052-0.014343-0.07553860.00171870.0092060.0038835-0.0767040.006655670.02067770.07074770.08004270.05807270.10087288-0.005059-0.017223-0.005381-0.006563-0.0227459-0.026554-0.062731-0.080492-0.111322-0.08145100.00921670.01263660.00969940.0164237-0.025275110.00400570.02145760.01722530.03030530.0152416对市场组合和证券1作回归。

先在目录G:\2glkx\data下建立tz7-1.xlsx数据文件后，取数的命令如下：import pandas as pdimport numpy as np#读取数据并创建数据表,名称为datadata=pd.DataFrame(pd.read_excel('F:\\2glkx\\data\\tz7-1.xlsx'))#查看数据表前5行的内容data.head()对表7-1中的r1t,rMt两列的数据进行回归分析，输入如下代码：import statsmodels.api as smimport pandas as pdimport numpy as npx = np.array(data[['rMt']])y = np.array(data[['r1t']])# model matrix with interceptX = sm.add_constant(x)#least squares fitmodel = sm.OLS(y, X)fit = model.fit()print fit.summary()得到如下结果： OLS Regression Results ========================================================================Dep. Variable: y R-squared: 0.557Model: OLS Adj. R-squared: 0.507Method: Least Squares F-statistic: 11.30Date: Thu, 16 Feb 2017 Prob (F-statistic): 0.00837Time: 10:。