第12章制表与绘图教学课件.ppt

/*接上页*/ axis2 label=(Degrees justify=right Centigrade) order=(-20 to 30 by 10) width=3; axis3 label=(h=3 Degrees justify=left Fahrenheit) order=(-4 to 86 by 18) width=3; proc gplot data=ResDat.minntemp; plot c2*month / caxis=red haxis=axis1 hminor=0 vaxis=axis2 vminor=1; plot2 f2*month / caxis=red vaxis=axis3 vminor=1; run; quit; 例12.21 回归结果作图。 从锐思数据库下载股票600603月持有期收益：股票——收益——持有期收益——月收益，选择600603股票并下载生成sas数据集res_sas.Monret_1。 /* 计算超额收益*/ data r600603_month; set res_sas.Monret_1; r600603=Monret-MonRfRet; r_m=Mretmc-MonRfRet; run; /* 回归与作图*/ proc reg data=r600603_month; model r600603=r_m / dw spec; output out=out1 r=r p=p l95=l u95=u; slope:test r_m=1; title2 兴业房产股票; run; /*接上页*/ goptions reset=global gunit=pct cback=white border htitle=6 htext=3 ftext=swissb colors=(black); goptions reset=symbol; proc sort data=out1 out=r_out ; by r_m; run; data regdata(keep=y_value pt_type r_m); set r_out; label pt_type=Observation Type; array regvar(4) r600603 p l u; array varlabel(4) $12. _temporary_ (Actual Predicted Lower Limits Upper Limits); do i=1 to 4; y_value=regvar(i); pt_type=varlabel(i); output; end; Run; /*接上页*/ proc gplot data=regdata; plot y_value*r_m=pt_type / haxis=axis1 vaxis=axis2 hminor=4 vminor=4; symbol1 v=* h=3.5 pct font=swissb color=black r=1; symbol2 i=join font=swissb l=1 color=blue r=1; symbol3 i=join font=swissb l=1 color=green r=1; symbol4 i=join font=swissb l=1 color=red r=1; axis1 order=(-.3 to .15 by .05); axis2 label=(angle=90 r600603) order=(-.5 to .5 by .25); title f=SIMSUN 实际值和预测值; title2 f=SIMSUN 带有上、下置信限; run; quit; 图表过程 图表过程GCHART输出高精度图表。使用GCHART过程可以制作二维或三维的柱状图和饼图。 图表过程句法 PROC GCHARTDATA=input-data-setANNOTATE=Annotate-data-set GOUT=libref.output