最优投资组合选择

第11章最优投资组合选择,清华大学经管学院朱世武ZhushwResdat样本数据：SAS论坛：,用线性规划选择投资组合,用线性规划求解最优投资组合步骤：用means过程计算股票收益；用data步生成procLp的输人数据集；用procLp求解最优投资组合权；用procLp进行灵敏度分析；用data步，根据最优权及投资组合的大小计算投资于每只股票的金额。,创建数据集,创建收益数据集return。数据集Return包括1995年2005年A股市场月持有期收益，及其它8只股票月持有期收益数据。return变量解释：Stkcd为股票代码；Date为日期；Mretmc为A股市场的月持有期收益；Monret为个股的月持有期收益。,计算期望收益,在使用PROCLP解决线性规划问题之前，要先估计期望收益这里用过去的平均收益来估计期望收益。用PROCMEANS语句计算股票收益。procmeansdata=returnnoprint;bystkcd;varmonret;outputout=m_out;datam_out1a;setm_out;where_stat_=mean;keepstkcdmonret;run;,风险度量,常用风险度量的指标有两种：股票收益标准差；资本资产定价模型（CAPM）的。,股票收益标准差,标准差是统计学对波动性的度量，而股票收益的波动性正是持有股票的风险所在。计算所选8只股票的标准差，保存在数据集m_out1b中。datam_out1b;setm_out;where_stat_=std;keepstkcdmonret;renamemonret=std;labelmonret=月收益率标准差;run;,CAPM的,第二种风险度量指标是CAPM的，它代表系统风险。,根据CAPM模型，资产i的期望收益和市场期望收益之间有如下关系：,投资组合p的CAPM：,.,投资组合权重为时（i=1,2,N），参数和可以表示成单个股票的线性组合：,总风险定义为收益的标准差，根据CAPM，资产i收益的方差为：,对于投资组合p，收益方差为：,当投资组合分散时，第二项会减小，分散好的投资组合，该项接近0，其收益的方差近似为：,因此，对于分散好的投资组合，就是要最小化,计算最优投资组合权重,计算出股票期望收益和风险水平后，就可以用PROCLP来找出在最大可接受风险的前提下收益最大的投资组合权。该线性规划问题为以下形式：,其中：A是一个的系数阵（rhs）。b是一个常数向量。c是一个的价格系数向量.x是一个的结构变量向量是的一个下界。PROCLP中缺省下界是0。是的一个上界。,下面的例中用PROCLP来最大化以下目标函数以求解投资组合权重：,.,该目标函数有以下约束条件：投资组合风险是1.2，即选择投资组合权重时要使权重与各股票的乘积之和是1.2。,权重和为l：,.,投资组合权重的上界是1。,下面的DATA步创建一个名为WEIGHT1的数据集，按PROCLP的格式说明线性规划问题。dataweight1;input_id_:$10.r000002r000007r000011r000016r600601r600604r600651r600653_type_$_rhs_;cards;exp_return0.02470.01270.01240.0080.02160.00680.02630.0144maxbeta1.11850018911.34701230111.31857492331.02606541291.35094999650.90389584431.25591093351.2961143173eq1.2sum_wts1.01.01.01.01.01.01.01.0eq1.0available11111111upperbd.available00000000lowerbd.;run;,灵敏度分析,可以通过灵敏度分析来考察投资组合权重（线性规划问题的主要解）的稳定性质。proclpdata=weight2primalin=lp_out2rangepricerangerhs;title2灵敏度分析;run;quit;,限定组合最高风险水平,风险约束也可以是不等式，即当投资组合的小于或等于一个最大可接受值时计算出的最优投资组合权。下面考察允许投资组合小于等于1.2时上面问题主要解的变化情况。dataweight3;setweight2;if_id_=betathen_type_=le;run;proclpdata=weight3primalout=lp_out3;run;quit;,用整数规划产生投资组合,例11.1用PROCLP求解整数规划。/*创建价格数据集price*/dataprice;mergeresdat.qttndiststk8(in=samp);bystkcd;ifsampanddate=30dec2005d;keepstkcdclpr;run;dataweight4;input_id_:$10.r000002r000007r000011r000016r600601r600604r600651r600653_type_$_rhs_;cards;exp_eger12345678integer.;/*求解整数规划*/proclpdata=weight4imaxit=200primalout=lp_out4dualout=d_out;run;,/*整理求解输出数据集*/datalp_out4a(keep=_var_value_rename=_var_=name);setlp_out4;if_n_8thendelete;run;/*打印主要解，列出每只股票应购买的手数*/procprintdata=lp_out4a;title整数规划;title2购买手数;run;注意，在PROCLP输出中，LOTS约束条件（限制投资金额）是一个最大值为$100000的不等式约束。在这$100000中，只有71975用于创建投资组合。,/*打印对偶解数据集*/procprintdata=d_out;title整数规划;title2收益、风险和投资金额;run;对偶解数据集包含以下信息：投资组合收益；最佳投资金额；最佳投资组合风险水平。注意，这里投资组合收益和风险度量与前面线性规划问题中的度量不一样，前面的投资组合收益和风险是每只股票收益和风险按最优投资组合权重的线性组合，因此不能直接比较。,用非线性规划产生投资组合,在Markowitz模型中，股票收益指的是平均收益（也叫算术平均收益）。股票风险是股票收益的标准差。证券投资组合收益是每只股票的加权收益。在不允许卖空，即组合权重在0和1之间，且和为1的条件下，证券投资组合的收益为：,两只股票组合的风险定义为：,用DATA步和PROCCORR产生投资组合,考虑从以下三只股票中选择两只构造投资组合，万科A（000002），深达声（000007），S深物业A（000011）。试图找出在既定风险水平下收益最大的投资组合。,计算股票收益的均值和方差,/*从RETURN数据集中抽出3只股票*/datareturn1;setreturn;wherestkcdin(000002000007000011);/*对数据集进行处理，将3只股票的数据合为一个数据集中的三个变量，以便用PROCCORR求3只股票月收益的协方差*/proctransposedata=return1out=return1(drop=_name_label_);bystkcd;varmonret;proctransposedata=return1out=return1(drop=_name_);proccorrdata=return1covoutp=cov_out1nosimple;varcol1col2col3;titleMarkowitz模型;quit;,计算股票投资组合的平均收益,一般说来，两种股票投资组合的收益用以下方法计算：,/*数据集COV_OUTl1中，增加变量X，表示权重。例中，X从0到l，步长为0.05。*/datacov_out2(drop=_name_);setcov_out1;if_type_neMEANthendelete;dox=0to1by.05;output;end;renamecol1=r000002col2=r000007col3=r000011;labelx=投资组合的权重;run;,/*计算投资组合的收益*/datamean1;setcov_out2;pfol_m1=x*r000002+(1-x)*r000007;pfol_m2=x*r000002+(1-x)*r000011;pfol_m3=x*r000007+(1-x)*r000011;run;procprintdata=mean1;titleMarkowitz模型;title2两种股票投资组合的收益;run;,图中，组合1用符号*，组合2用符号，组合3用+。由图可以看出，投资组合1在所有组合权重水平中都有最高收益。,计算投资组合的风险,投资组合的风险可通过下式计算：,程序参考教材,含有两只股票投资组合的风险水平Obsxpfol_r1pfol_r2pfol_r310.000.157650.157610.1576120.050.153890.154570.1550630.100.150350.151710.1527540.150.147040.149030.1506850.200.143990.146550.1488660.250.141210.144270.1473170.300.138710.142210.1460280.350.136520.140360.1450190.400.134640.138760.14429,从列表可以看出风险水平随投资组合权重变化。例中，投资组合1在观测14处风险最小（0.13042），即65的万科A（000002）和35的深达声（000007）。此外，对于每种权重，投资组合1都有最低风险。结合上面分析，组合1是最佳选择。,识别有效证券组合,有效证券组合在任何既定风险水平下都有最大收益。可以通过GPLOT过程画出收益对风险的图形，从而形象地识别出有效证券组合。,图中，竖轴表示收益，横轴表示风险。曲线的端点（靠右边）是不分散投资组合（一种股票占100，另一种为0），分散投资降低了风险水平。投资组合的有效边界是在既定风险水平下具有最大收益的组合，也即最左上的组合。组合3远不如组合1和2好。,先去掉部分没有优势的投资组合，然后画出投资组合的一个子集。下面把权重小于0.4的投资组合删去，因为从前面图中可看出它们的收益较小而风险较高，同时去掉组合3。,由图可以看出，投资组合的有效边界是权重从0.65到1的组合1。其次是权重从0.7到1.0的组合2。,用PROCNLP产生投资组合,NLP过程可以解决大量线性和非线性最优化问题。用PROCNLP过程时，还可以加入边界条件和一般线性条件。,下面仍以本章选择的8只股票为例，说明使用PROCNLP求解投资组合的过程。权重计为XlX8，要求在0与1之间，和为1。最后要求期望收益最少为1.9。,datareturn2;setreturn;proctransposedata=return2out=return2(drop=_name_label_);bystkcd;varmonret;proctransposedata=return2out=return2(drop=_name_);datareturn2;setreturn2;renamecol1=x1col2=x2col3=x3col4=x4col5=x5col6=x6col7=x7col8=x8;proccorrdata=return2covoutp=cov_out3nosimple;varx1-x8;titleMarkowitz模型;quit;,procnlpoutest=nlp_out4noprint;parmsx1-x8;var=x1*x1*0.0194244422+x2*x2*0.0248545995+x3*x3*0.024841303+x4*x4*0.0150937098+x5*x5*0.029767682+x6*x6*0.0156948233+x7*x7*0.0295971799+x8*x8*0.0227559726+2*x1*x2*0.012652434+2*x1*x3*0.0150048423+2*x1*x4*0.0085234643+2*x1*x5*0.0087603573+2*x1*x6*0.0043113555+2*x1*x7*0.0078285695+2*x1*x8*0.0068767269+2*x2*x3*0.0164591308+2*x2*x4*0.0102511593+2*x2*x5*0.0125957989+2*x2*x6*0.0091172983+2*x2*x7*0.0140708542+2*x2*x8*0.0126957816+2*x3*x4*0.0084318348+2*x3*x5*0.0097300933+2*x3*x6*0.0075174348+2*x3*x7*0.0084425639+2*x3*x8*0.0090521457+2*x4*x5*0.009197054+2*x4*x6*0.0043473766+2*x4*x7*0.0070145728+2*x4*x8*0.0091100859+2*x5*x6*0.0101411356+2*x5*x7*0.0175577636+2*x5*x8*0.0167131136+2*x6*x7*0.0104334122+2*x6*x8*0.0083393987+2*x7*x8*0.0136299731;risk=sqrt(var);minrisk;bounds0=x1-x8=1;lincon1=x1+x2+x3+x4+x5+x6+x7+x8,.019=0.0246951391*x1+0.0127497123*x2+0.0124444582*x3+0.0079714507*x4+0.0216395959*x5+0.0067692931*x6+0.0263324389*x7+0.0143749164*x8;run;,计算投资金额,dataamt_out3;setnlp_out4;if_TYPE_nePARMSthendelete;arrayxa(8)x1-x8;arraya(8)a1-a8;doi=1to8;a(i)=round(xa(i)*100000,.01);end;total=sum(ofa1-a8);run;dataamt_out4;setamt_out3;renamea1=r000002a2=r000007a3=r000011a4=r000016a5=r600601a6=r600604a7=r600651a8=r600653;run;procprintdata=amt_out4;varr000002r000007r000011r000016r600601r600604r600651r600653total;title2每只股票投资的数额;run;,使用INEST=输入数据集,在PROCNLP中使用INEST=输入收据集，可以使其余PROCNLP语句简化，同时也简化诊断检查。INEST=数据集为TYPE=EST中的数据。所有股票的初始值均取为0.125。如果这些初始值不可行，PROCNLP将计算可行的初始值。下界是0，上界是1。收益不等式约束的右边值设为0.019（即1.9个百分点）。,datain_nlp4(type=est);input_type_$8.x1-x8_rhs_;cards;parms.25.125.lowerbd00000000.upperbd11111111.eq111111111ge0.02469513910.01274971230.01244445820.00797145070.02163959590.00676929310.02633243890.01437491640.019,procnlpinest=in_nlp4outest=nlp_out5;parmsx1-x8;var=x1*x1*0.0194244422+x2*x2*0.0248545995+x3*x3*0.024841303+x4*x4*0.0150937098+x5*x5*0.029767682+x6*x6*0.0156948233+x7*x7*0.0295971799+x8*x8*0.0227559726+2*x1*x2*0.012652434+2*x1*x3*0.0150048423+2*x1*x4*0.0085234643+2*x1*x5*0.0087603573+2*x1*x6*0.0043113555+2*x1*x7*0.0078285695+2*x1*x8*0.0068767269+2*x2*x3*0.0164591308+2*x2*x4*0.0102511593+2*x2*x5*0.0125957989+2*x2*x6*0.0091172983+2*x2*x7*0.0140708542+2*x2*x8*0.0126957816+2*x3*x4*0.0084318348+2*x3*x5*0.0097300933+2*x3*x6*0.0075174348+2*x3*x7*0.0084425639+2*x3*x8*0.0090521457+2*x4*x5*0.009197054+2*x4*x6*0.0043473766+2*x4*x7*0.0070145728+2*x4*x8*0.0091100859+2*x5*x6*0.0101411356+2*x5*x7*0.0175577636+2*x5*x8