《股票市場(chǎng)風(fēng)險(xiǎn)指標(biāo)計(jì)算》由會(huì)員分享,可在線閱讀,更多相關(guān)《股票市場(chǎng)風(fēng)險(xiǎn)指標(biāo)計(jì)算(16頁珍藏版)》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、第7章 股票市場(chǎng)風(fēng)險(xiǎn)指標(biāo)計(jì)算,清華大學(xué)經(jīng)管學(xué)院 朱世武 Z Resdat樣本數(shù)據(jù): SAS論壇: ,計(jì)算指標(biāo)與環(huán)境,理論模型,,其中, 為股票j的超額收益; 為市場(chǎng)超額收益; 和 為參數(shù),且假定 不相關(guān)。,,,,,,將方差作為風(fēng)險(xiǎn)的度量,根據(jù)上述假設(shè)得到股票j的總風(fēng)險(xiǎn)為:,,選定時(shí)期,根據(jù)經(jīng)驗(yàn)數(shù)據(jù)可以得模型參數(shù)的最小二乘估計(jì) 以及相應(yīng)估計(jì)值 。于是,系統(tǒng)風(fēng)險(xiǎn)占該股票總風(fēng)險(xiǎn)的比例估計(jì)值為:,,,,,有,計(jì)算指標(biāo),計(jì)算1999至2005年,滬深兩地股票市場(chǎng)系統(tǒng)風(fēng)險(xiǎn)系數(shù)、個(gè)股總風(fēng)險(xiǎn)和個(gè)系統(tǒng)風(fēng)險(xiǎn)占總風(fēng)險(xiǎn)比例三個(gè)風(fēng)險(xiǎn)指標(biāo)等。即具體計(jì)算: 系統(tǒng)風(fēng)險(xiǎn)系數(shù)估計(jì)值 ; 個(gè)股總風(fēng)險(xiǎn)估計(jì)
2、值 ; 系統(tǒng)風(fēng)險(xiǎn)占總風(fēng)險(xiǎn)比估計(jì)值 ; 全市股票上述三種指標(biāo)的概況統(tǒng)計(jì)量。,,,,計(jì)算環(huán)境,計(jì)算數(shù)據(jù)集: 銀行存款利率BankIr; 基準(zhǔn)利率BchmkIr; 個(gè)股日收益Dret; 市場(chǎng)日收益DretM。 時(shí)期:2005年。 樣本:滬深兩市場(chǎng)滿足條件的2006年前上市所有股票日百分比收益。剔除1年內(nèi)交易小于50天的股票。,算法實(shí)現(xiàn),第一步:創(chuàng)建無風(fēng)險(xiǎn)利率數(shù)據(jù)集RF;2005年市場(chǎng)日收益數(shù)據(jù)。 第二步:創(chuàng)建宏文本:Stk2006.TXT。 第三步:計(jì)算2005年滬深兩市場(chǎng)A股個(gè)股風(fēng)險(xiǎn)指標(biāo)。 第四步:滬深兩市場(chǎng)個(gè)股與市場(chǎng)風(fēng)險(xiǎn)關(guān)系統(tǒng)計(jì)分析。,計(jì)算程序,創(chuàng)建無風(fēng)險(xiǎn)利率數(shù)據(jù)集: /*1998年7
3、月1日后使用七日回購利率兩周指數(shù)加權(quán)平均為基準(zhǔn)利率 */ data a; set ResDat.bchmkir; if code=B2W and date=1jun1998d; rename ir=rf; keep ir date;run; /*1998年7月1日前用一年期銀行存款利率加10%為基準(zhǔn)利率*/ data b; set ResDat.bankir; if code=d1y ; run; data c; set b; format date yymmdd10.; if enddt=. Then enddt=date(); do date=begdt to enddt; output;
4、 end; run;,data d; set c; if date<1jun1998d; keep date rf; rf=ir*1.1; /* 在一年期存款利率的基礎(chǔ)上再增加10% */ run; /*合并得到基準(zhǔn)利率表rfyr*/ data rfyr; set d a; if date<=31dec2005d; run; /*得到日基準(zhǔn)利率表rf*/ data rf; set rfyr; rf=rf/365; /* 將利率轉(zhuǎn)化為日利率 */ run;,2005年市場(chǎng)日收益數(shù)據(jù)DretM2005:總市值加權(quán)平均市場(chǎng)日收益率。 Data DretM2005; Set ResDat.DretM;
5、 If year(date)=2005 and Exchflg=0 and Mktflg=AB; /* Exchflg=0滬、深兩交易所;Mktflg=AB AB股票市場(chǎng) */ Keep date Dretmc; Run;,創(chuàng)建宏文本:Stk2006.TXT。 利用選擇的樣本股票收益數(shù)據(jù)集Dret創(chuàng)建宏。如果要計(jì)算全部股票的風(fēng)險(xiǎn)指標(biāo),可以直接利用lstkinfo創(chuàng)建宏。 Data a; Set ResDat.dret; If dif(stkcd)=0 then delete; data _null_; set a; a=%a(; b=,; c=); ; file Stk2006.txt ; p
6、ut a $ stkcd $ b $ lstknm $ c $ ; run;,%a( 000001 , S深發(fā)展A ); %a( 000002 , 萬科A ); %a( 000004 , *ST國農(nóng) ); %a( 000005 , ST星源 ); %a( 000007 , 深達(dá)聲A ); %a( 000009 , S深寶安A ); %a( 000011 , S深物業(yè)A ); %a( 000012 , 南玻A ); %a( 000016 , 深康佳A ); %a( 000017 , S*ST中華 ); %a( 200011 , 深物業(yè)B ); %a( 200012 , 南玻B ); ,計(jì)算
7、2005年滬深兩市場(chǎng)個(gè)股風(fēng)險(xiǎn)指標(biāo):模型r_s=a+beta*r_m+e /* 2005年無風(fēng)險(xiǎn)利率數(shù)據(jù)集rf2005 */ Data rf2005; Set rf; If year(date)=2005; Run; /* 2005市場(chǎng)超額收益數(shù)據(jù)集 */ Data PDretM2005; Merge DretM2005 rf2005; By date; rm= Dretmc- rf; keep date rm; run; data PDretM2005; /* 收益數(shù)據(jù)缺失時(shí),用前面交易日的數(shù)據(jù)填充 */ set PDretM2005; retain r_m; if rm=. Then r_m
8、=rm; drop rm; run; /* 2005年個(gè)股票超額收益數(shù)據(jù)集 */ Data dret; Set ResDat.Dret; Where year(date)=2005; Run;,Proc sql; Create table Pdret2005 as Select * from Dret left join rf On Dret.date=rf.date; Run; Proc sort data= Pdret2005; By stkcd date; Run; Data Pdret2005; Set Pdret2005; r_s=dret-rf; keep stkcd lstknm
9、 date r_s; run; Proc sql; Create table Capm2005 as Select * from Pdret2005 left join PdretM2005 On Pdret2005.date=PdretM2005.date; Run; Proc sort data=capm2005; By stkcd; Run;,/* 計(jì)算個(gè)股風(fēng)險(xiǎn)指標(biāo):數(shù)據(jù)集Capm2005 */ proc reg data=capm2005 outest=d cp noprint; /*輸出參數(shù)估計(jì)數(shù)據(jù)集中包括Mallow Cp估計(jì)量 */ model r_s=r_m; /* 回歸模型為
10、:r_s=a+beta*r_m+e */ by stkcd; run; data d; set d; if (_EDF_+_P_)<50 then delete; /* 刪除2005年交易日小于50天的股票,_EDF_為殘差自由度,_P_為參數(shù)個(gè)數(shù) */ data Riskl2005 (keep=sn obs _rmse_ r_m _rsq_ stkcd label=”2005年個(gè)股與市場(chǎng)風(fēng)險(xiǎn)指標(biāo)”); set d; sn=_n_; obs=(_edf_+_p_); label sn=個(gè)股序號(hào); label obs=觀測(cè)個(gè)數(shù); label _rmse_=個(gè)股風(fēng)險(xiǎn)度; label r_m=系統(tǒng)風(fēng)
11、險(xiǎn)估計(jì)值beta; label _rsq_=系統(tǒng)風(fēng)險(xiǎn)占總風(fēng)險(xiǎn)比例; run;,對(duì)數(shù)據(jù)集Riskl2005中個(gè)股的三種風(fēng)險(xiǎn)指標(biāo)作圖: %macro a(x); proc gplot data= Riskl2005; plot ,,,,,2005年股票市場(chǎng)個(gè)股風(fēng)險(xiǎn)指標(biāo)的概括統(tǒng)計(jì)量計(jì)算: proc univariate data= Riskl2005 noprint; var r_m _rmse_ _rsq_ ; output out=out n=n_r_m n_rmse_ n_rsq_ mean=mean_r_m mean_rmse_ mean_rsq_ median=median_r_m med
12、ian_rmse_ median_rsq_ min=min_r_m min_rmse_ min_rsq_ max=max_r_m max_rmse_ max_rsq_ q1=q1_r_m q1_rmse_ q1_rsq_ q3=q3_r_m q3_rmse_ q3_rsq_ range=range_r_m range_rmse_ range_rsq_; data out; set out; q3_q1_r_m =q3_r_m -q1_r_m ; q3_q1__rmse_= q3_rmse_-q1_rmse_ ; q3_q1_rsq_ = q3_rsq_-q1_rsq_ ; format _numeric_ 6.4; run;,