金融计量学作业1前五次

读 bas.xls bas.xls、industry.xlsx中的数据， bas.xls（appendProcimportout=onedatafile="C:\Users\nancy\Desktop\guanghua\assignment\ DBMS=EXCELRECE; Procimportout=twoDATAFILE="C:\Users\nancy\Desktop\guanghua\assignment\ DBMS=EXCELRECE; DBMS=EXCELRECE;settwo;format keepStkcdAerA renameA0b byStkcd;byStkcd;Procsortdata=threedatabyStkcd;byIndCd1;ClassIndCd1;ESS/(k ESS(n R2(n 0.54(50F RSS/(n (TSSRSS)(k

(1R2)(k1) (41)(10.54)R&Dperformancet检验。H0R&D的系数为零，H1R&D的系(1 (0.62t统计量为t

自由度为45的t分布在10%的显著性水平下的临界值约为1.68，在5%的显著性水平下的临界值约为2.01%的显著性水平下的临界值约为2.67。而我们得到的t2.0，因此，在10%的显著性水平下我们可以原假设，即认为R&D对performance有影响。我们使用单边统计检验，t值仍2，但是自由度为45的t分布10%显著性水平下的临界值为1305%显著性水平下的临界值为1681%显著性水平下的临界值为2.40系数是相同的，并将回归的RR2记录。然后我们利用F统计量来做判断。(R2R2)/RUR2)/(50F(1RUR2)/(50Performanceൌα଴଴αଵR&ܦ଴αଶSize଴αଷAge଴αସLev଴此模型是关于α的线性模型，因此保证 无偏需要的其他条件有随机干扰项εE(R&Dsizeage,levE(0var(R&Dsizeagelev2导入个股收益率数据并拼接成一个文件，读入上证综指收益率数据DBMS=EXCELRE SHEET="TRD_Dalyr";GETNAMES=YES;MIXED=YES;DBMS=EXCELRE SHEET="TRD_Dalyr1";GETNAMES=YES; DBMS=EXCELRE SHEET="TRD_Dalyr2";GETNAMES=YES; DBMS=EXCELRE SHEET="TRD_Dalyr3";GETNAMES=YES; DBMS=EXCEL SHEET="TRD_Dalyr4";GETNAMES=YES; DBMS=EXCELRE SHEET="TRD_Dalyr";GETNAMES=YES;MIXED=YES;DBMS=EXCELRE SHEET="TRD_Dalyr1";GETNAMES=YES; DBMS=EXCELRE SHEET="TRD_Dalyr2";GETNAMES=YES; DBMS=EXCELRE SHEET="TRD_Dalyr3";GETNAMES=YES; procIMPORTOUT=trd10"F:\econometrics\hw2\hw2-1\TRD_Dalyr2-5"DBMS=EXCELRE CE;SHEET="TRD_Dalyr4";GETNAMES=YES;MIXED=YES; procIMPORTOUT=trd11"F:\econometrics\hw2\hw2-1\TRD_Dalyr3-1"DBMS=EXCELRE SHEET="TRD_Dalyr";GETNAMES=YES;MIXED=YES;USEDATE=NO;procIMPORTOUT=trd12"F:\econometrics\hw2\hw2-1\TRD_Dalyr3-2"DBMS=EXCELRE CE;SHEET="TRD_Dalyr1";GETNAMES=YES;MIXED=YES; procIMPORTOUT=trd13"F:\econometrics\hw2\hw2-1\TRD_Dalyr3-3"DBMS=EXCELRE CE;SHEET="TRD_Dalyr2";GETNAMES=YES;MIXED=YES; procIMPORTOUT=trd14"F:\econometrics\hw2\hw2-1\TRD_Dalyr3-4"DBMS=EXCELRE CE;SHEET="TRD_Dalyr3";GETNAMES=YES;MIXED=YES; procIMPORTOUT=indexDATAFILE="F:\econometrics\hw2\hw2-1\TRD_Index"DBMS=EXCELRE SHEET="TRD_Index";GETNAMES=YES;MIXED=YES;USEDATE=NO;datatrdtotal;settrd1trd2trd3trd4trd5trd6trd7trd8trd9trd10trd11trd12trd13datatrdtotal1;settrdtotal;datechar=put(trddt,$12.);year=substr(datechar,1,4);keepstkcdtrddtdretwdyear;datarettotal09;settrdtotal1;ifyear="2009";dropyear;run;datarettotal10;settrdtotal1;ifyear="2010";dropyear;dataindex1;setdatechar=put(trddt,$12.);year=substr(datechar,1,4);keeptrddtretindexyear;dataretdex09;setindex1;ifyear="2009";dropyear;run;dataretdex10;setindex1;ifyear="2010";dropyear;将个股2009年收益率数据和上证综指的收益率数据合并为一个procsortdata=rettotal09;byprocsortdata=retdex09;bytrddt;datafirstq; lengthSTKCD$12TRDDT$20DRETWD$32Retindex$18;mergeretdex09(in=x)rettotal09(in=y);bytrddt;ifx=y;run;datafirstq;setformatdretwd1retindex1best32.;keepdretwd1retindex1trddtrenamedretwd1=dretwd;rename计算系统风险和非系统procsortdata=firstq;bystkcd;run;libnameresult"F:\econometrics\hw2";run;procregoutest=result.OUTFILEEDFdata=firstqmodeldretwd=retindex;bystkcd;outputout=resout系统性风险dataresult.syst;setresult.outfile;keepstkcdretindex;非系统性风procmeansdata=resoutvarress;bystkcd;outputout=result.varresdataresult.specific;setresult.varres;keepstkcdvarres;（2）procIMPORTOUT=bsDATAFILE="F:\econometrics\hw2\hw2-1\ DBMS=EXCELRE bas";GETNAMES=YES;MIXED=YES;USEDATE=NO;databs2;setbs;if databsda;setbs2; ;keepstkcd计算个股的datarettotal10a;setrettotal10;formatdretwd1best32.;keepdretwd1trddtstkcd;renamedretwd1=dretwd;run;procmeansdata=rettotal10anoprint;VARdretwd;bystkcd;outputout=ret10a;run;dataret10b;setret10a;if_stat_="MEAN";run;dataret2010;setret10b;annualret=dretwd*_FREQ_;keepstkcdannualret;run;将资产负债率、 年收益率、系统风险和非系统风险合并为一个文procsortdata=ret2010;bystkcd;run;procsortdata=bsda;bystkcd;run;procsortdata=result.specific;bystkcd;run;procsortdata=result.syst;bystkcd;run;dataresult.secondq;mergeret2010bsdaresult.specificresult.syst;ifmissing(dar)ormissing(retindex)thendelete;run;进行回procregoutest=result.answersEDFdata=result.secondq;modelannualret=retindexvarresDAR;输出结果tPr>1111 ଴଴଴଴଴଴଴଴଴଴଴ൌ଴.଴ૠ଴଴.଴଴଴଴଴଴଴଴଴଴଴଴଴଴଴଴଴଴଴ൌ଴.଴଴଴଴.଴଴଴଴଴଴ൌ଴.଴଴଴଴଴଴଴.请解释这两个回归模型中Edu的系数的经第一个模型中，edu的系数表示在企业高管的教育水平提高一个单位后，企业业绩会提二个模型中，Edulev保持相同水平的情况高出0.002个单位。Formultipleregressionહ଴ܑ܁હ଴theconditionsunderwhichtheOLSestimatorofહ଴isunbiased,andlisttheconditionsunderwhichtheOLSestimatorofહ଴isconsistent.此模型是关于α的线性模型，因此保证 无偏需要的其他条件有1Retindex表示系统性风险，varres表示非系统性风险，dar表示资产负债率。***1%的显著性水平随机干扰项εE(R&Dsizeage,levE(0var(R&Dsizeagelev2Oneresearchusethefollowingregressiontoexaminehowfirms’R&DaffecttheirPerformanceൌα଴଴αଵR&ܦ଴αଶSize଴αଷAge଴αସLev଴wherePerformancecanbemeasuredbyROE(return‐on‐equity)orROA(return‐on‐asset),R&DcanbemeasuredbynumberofpatentsorratioofR&Demployeenumbers(tothetotalemployeenumbers),Size,Age,andLevarethesize,ageanddebt‐to‐assetratio.PleasediscussWhichmeasuretouseforperformance,ROEorWhichmeasuretouseforR&D,numberofpatentsorratioofR&Demployeenumbers?HowdoestheinclusionofAgeaffectthevarianceofestimatorofEmpiricallyhowtodecidewhethertoincludeAgeasacontrol直接由资产负债率决定的performance更合适。numberofpatents、ratioofR&DnumbersR&DR2，选择使得R2更大的一个指标作为R&D的衡量指标。 会影响，但是影响的方向不确定。var(ˆ1SST(1R2，R1R&D sizelev（在模型包含age时，也包括age）回归的R2归时，ageR12age变量后，分子变小，而分母也变小，那αଵ的方差的变化就是不确定的首先，可以先从经济意义上初步考虑，age是有必要加入的，因为不同发展阶段的公R2taget值是否显著；然后我们对不包含age的模型进行回归，再记录回归的R2。如果加入age能显著增加模Oneresearchusesthefollowingregressionmodeltostudyfirms’܀

଴ܑܖ۷଴଴ܑ଴ൌ

଴ܑܑܞ܍

଴ܑܑ܍

܁∗ܑܑܞ܍where

଴ܑܑܞܖ۷ܔ଴હ଴istheinvestmentmadebyfirmiinyeartmeasuredbytheamount,

istheinvestmentmadebyfirmiinyear

isto‐assetratiooffirmiinyear

istheageoffirmiinyeartbytheyearssinceitwasestablished,܁

isadummyvariableequaltoiffirmiisaState‐Owned‐Enterpriseinyeart.Pleasepointoutallproblemswiththisresearch对于国有公司，ROE୧ൌβ୧଴αଵInv୧଴ሺαଶ଴αସሻLev୧଴αଷAge୧଴αହlInv୧ε୧对于非国有公司，ROE୧ൌβ୧଴αଵInv୧଴αଶLev୧଴αଷAge୧଴αହlInv୧଴Inv୧4、○1234一般是介于[‐1,1]之间的数值，用多元线性回归方程不能有效的保证HW2Ret01SysRisk2IdioRisk3Lev4IdioRisk*Lev“公司报赖资”的请给出“风险不影响公司回报”的假设检验，并具体进行检验对平债，不司（1）做出的假设为：H0:4 H1:4ulibnameresult"F:\econometrics\hw2";rudataresult.hw3t5;

retindex=by= dataresult.hw3t5;setresult.hw3t5;renmedar=lv; procregnmodelannualret=sysriskidiorisklevin归回的结果为：归 计为‐0.04745，对应的t值为‐0.36，概率为72.03%，风险对公司回报的影响依赖于公司的资产负债率。 回报率的影响为：2+4lev 假设为：H0:240,H1:2,4并不一定为0.nprocregnmodelannualret=sysriskidiorisklevitestidiorisk=0,indlev=0;

归回的结果为：归F值为0.23，概率为79.82%，因此我们无法原假设没有合理的依据 procmeansdata=result.hw3t5;varlev;run;可见，平均的资产负债率为74.96193%H0274.96193%40,procreg

274.96193%4modelannualret=sysriskidiorisklevitestidiorisk+0. 归回的结果为：

上述结果显示，F值为0.35，概率为55.36%，因此我们无法 认为“对一个具有平均资产负债率的企业，风险不影响公司回报。HW2Ret01SysRisk2IdioRisk3Lev按照 将样本分为两组，检验上述模型在两组之间是否有差异test。分组SASprocrankdata=inilegroups=2out=outfile;varlev;(新数据文件outfile除了原文件infile中的变还有一个变量levrank，它将取01,1lev大于其中位数，0lev小于其中位数)procregmodelannualret=sysriskidiorisk果为：果为：我们得到回归方程为（0.01177） （0.134 模型的F值为1.12，概率为34.01%，因此我们在10%的显著性水平上无法原假设，模 procrankdata=result.hw3t5groups=2ot=result varranksprocregemodelannualret=sysriskidioriskleve回归结果为

由上述结果可以看出，levrank对应的系数的估计为‐0.06542，概率为0.22%1%的显著性水的回报率资产负债率高于中位数的公司的回报率大约比资产负债率低于中位数水平的公司的回报率低6.54% .HW2Ret01SysRisk2IdioRisk3Lev请对所有变量进行winsorize，再做回归。并解释说明情况下应该进行winsorize不进行winsorizeprocregumodelannualret=sysriskidiorisklev;ru由上述结果可知，1、 2、3的估计都是不显著的，模型整体也是不显著的（F0.3401ulibnameresult"F:\econometrics\hw2";rn;dataresult.hw4t1;setresult.hw3t5;run;u%macrowinsor1(datafile,tprocunivariatedata=&datafile.noprin;var&var.;toutputout=extremepctlpre=ppctlpts=199;procselectp1into:myp1fromselectp99into:myp99fromextreme;data&datafile.;

&var.>%sysevalf(&myp99)then&var.= and

% %

%winsor1(result.hw4t1,lev)procregumodelannualret=sysriskidiorisklev;ruwinsorize处理之后，可见3的估计和截距项的估计都是显著的，模型整体是5%的显著性水平上是显著的（F3.650.0122。winsorize处理，将大（小）于某一阀HW2Ret01SysRiskIdioRisk3Lev4Risk这里Risk是个股2009年的总风险，即个股日度回报的方差。请这个回归模案答案（1）计算2009odataresult.rettotal09;setresult.rettformatdretwd1best32.;o

keepdretwd1stkcdrenamekprocsortdata=result.rettotal09;bykprocmeansdata=result.rettotal09

eoutputout=result.rettotalvar09

twddataresult.rettotalvar09;setresult.ettotalvv vkeepstkcddretwd;rename

个报个报

nprocsortbprocsortdata=result.hw3t5b

ystkcd

;;data5mergeresult.rettotalvar09result.hw3t;bystkcd;run;5procreg modelannualret=sysriskidiorisklev

etwdvar/可见，dretwdvar350idiorisk的方差膨胀系349，dretwvar、idiorisk均与其他解释变量存在多重共线性。同时，模型整体10%的显著性水平上也是不显著的（F098。HW2Ret01SysRisk2IdioRisk3Lev，请进行异方差检验（BP方法及white方法）BPprocregdata=result.hw3t5noprint;modelannualret=sysriskidiorisklev;routpu