data resdat.garch_port8(keep=date);
set resdat.idx000001;
where 2004<=year(date)<=2005;
%macro a(x);
data a(keep=date r_1);
set resdat.stk&x;
where 2004<=year(date)<=2005;
adjclpr=clpr*mcfacpr;
r_1=(adjclpr-lag(adjclpr))/lag(adjclpr); /*日百分比收益*/
data r_port8;
set a;
if r_1=. then r_1=0;
else r_1=r_1;
/*单个股票波动率计算:*/
proc autoreg data=r_port8;
model r_1=/ nlag=1 garch=(q=1,p=1,tr);
output out=out cev=cev;
data out(keep=date cev_1);
set out;
cev_1=sqrt(cev);
data resat.garch_port8(rename=(cev_1=cev&x));
merge resat.garch_port8 out;
by date;
%mengd a;
%include"random1995.txt";
run;



雷达卡



京公网安备 11010802022788号







