y = xlsread('残差.xls');
Residuals = y - mean(y);
[H,pValue,stat,cValue] = archtest(Residuals,[10 15 20 ],0.05);%arch检验
%参数估算
ToEs*****l = garch(1,1);
[Es*****l,EstParamCov,logL,info] = estimate(ToEs*****l,y);
se = sqrt(diag(EstParamCov));%the estimate standard errors
%预测
rng('default')
[V,Y] = simulate(Es*****l,66);
Vf1 = forecast(Es*****l,22,'Y0',Y);
figure
plot(V,'Color',[.7,.7,.7])
hold on
plot(67:88,Vf1,'r','LineWidth',2);
title('Forecasted Conditional Variances')
legend('Observed','Forecasts with Presamples',...
'Location','NorthEast')
hold off
残差.xls
(29.5 KB)


雷达卡


京公网安备 11010802022788号







