|
data example5_3; input x@@; t=_n_; cards; 10.77 13.30 16.64 19.54 18.97 20.52 24.36 23.51 27.16 30.80 31.84 31.63 32.68 34.90 33.85 33.09 35.46 35.32 39.94 37.47 35.24 33.03 32.67 35.20 32.36 32.34 38.45 38.17 32.14 39.70 49.42 47.86 48.34 62.50 63.56 67.61 64.59 66.17 67.50 76.12 79.31 78.85 81.34 87.06 86.41 93.20 82.95 72.96 61.10 61.27 71.58 88.34 98.70 97.31 97.17 91.17 80.20 85.12 81.40 70.87 57.75 52.35 67.50 87.95 85.46 84.55 98.16 102.42 113.02 119.95 122.37 126.96 122.79 127.96 139.20 141.05 140.87 137.08 145.53 145.59 134.36 122.54 106.92 97.23 110.39 132.40 152.30 154.91 152.69 162.67 160.31 142.57 146.54 153.83 141.81 157.83 161.79 142.07 139.43 140.92 154.61 172.33 191.78 199.27 197.57 189.29 181.49 166.84 154.28 150.12 165.17 170.32 ; proc gplot data=example5_3; plot x*t=1; symbol1 c=black i=join v=star; proc autoreg data=example5_3; model x=t/nlag=5 dwprob archtest; model x=t/nlag=2 noint garch=(p=1,q=1); output out=out p=xp; proc gplot data=out; plot x*t=2 xp*t=3/overlay; symbol2 v=star i=none c=black; symbol3 v=none i=join c=red w=2; run;
|