Hi, I am trying to use proc model to run a time series regression yt = A*xt by 10 group ids and then use the regression result to forecast y(t+1) at x(t+1). The coefficients came out with 10 groups. However the forecast result came out just as one group. I wonder how can I generate forecast by groups as well. My code is as follows. data shock1; *generate x1 to xt then forecast y at x(t+1) and x(t+2) r = 0; * 5 lagged values assumed be zero; do i=1 to 5; output; end; i=6; * impulses; r = 1; output shock1; r = .; * future values; do i=7 to 8; output; end; run; proc model data=regdata; by groupID; r = A1*lag5(1,r)+A2*lag5(2,r)+A3*lag5(3,r)+A4*lag5(4,r)+A5*lag5(5,r); fit r/ sur outs=vcov (drop=_nused_) outest=est out=res; solve / forecast data=shock1 out=prediction; The out=prediction dataset only generate one group of forecast without groupID. I need the forecast to be also by groupID. I would appreciate it very much if anyone could help me on this.
... View more