options ls=72 ps=50; data arima;input t x y; x1=dif(x);y1=dif(y); cards; 1 108.377 106.415 2 110.823 106.972 3 113.964 108.935 4 116.134 110.576 5 119.565 111.196 6 122.050 111.065 7 123.751 113.136 8 125.643 113.416 9 128.457 114.757 10 130.071 114.819 11 131.685 117.026 12 133.828 118.144 13 135.770 119.107 14 138.502 119.792 15 140.825 122.468 16 144.173 122.074 17 147.591 124.914 18 149.641 122.471 19 152.275 126.401 20 155.434 125.789 21 157.352 127.960 22 160.852 127.094 23 164.692 129.393 24 169.367 131.400 25 173.796 131.925 26 177.584 131.996 27 179.534 133.563 28 180.687 133.188 29 180.365 133.410 30 180.209 136.611 31 177.917 136.652 32 177.130 137.308 33 176.063 139.098 34 175.594 138.667 35 175.831 141.188 36 176.569 143.312 37 177.178 143.851 38 176.153 143.349 39 176.368 145.479 40 176.861 146.053 41 178.001 145.198 42 179.204 147.253 43 181.523 147.642 44 183.492 149.644 45 186.408 150.171 46 189.361 150.979 47 191.226 154.276 48 191.927 154.378 49 191.410 153.612 50 190.127 155.565 51 189.356 155.494 52 191.509 155.153 53 194.482 159.418 54 198.870 159.103 55 200.991 159.605 56 201.980 161.858 57 202.588 163.555 58 203.738 163.133 59 204.037 163.813 60 203.971 165.099 61 204.233 165.756 62 205.146 166.266 63 206.590 168.721 64 208.199 168.784 65 210.815 170.185 66 213.223 171.790 67 215.111 170.919 68 217.959 171.843 69 221.914 172.919 70 226.284 174.903 71 229.555 175.518 72 232.527 176.331 73 233.870 175.443 74 234.579 176.312 75 235.410 177.676 76 237.470 180.410 77 239.555 182.109 78 242.617 182.461 79 245.283 181.727 80 247.481 185.294 81 247.774 185.397 82 247.944 184.576 83 247.392 188.606 84 247.011 186.442 85 248.064 186.860 86 248.154 189.173 87 248.929 189.841 88 248.741 190.755 89 250.057 193.146 90 250.023 192.801 91 251.518 194.451 92 251.733 194.274 93 252.369 197.181 94 252.894 196.462 95 252.227 200.087 96 253.721 199.964 97 253.556 199.187 98 254.269 199.752 99 255.058 199.834 100 256.168 203.643 ; proc plot;plot x*t;run; proc plot;plot x1*t;run; proc plot;plot y*t;run; proc plot;plot y1*t;run; proc arima data=arima; identify var=x; identify var=x1; estimate p=1 method=ml; forecast out=a1 back=0 lead=10 id=t;run; proc plot data=a1(firstobs=90); plot x1*t='*' forecast*t='f' l95*t='l' u95*t='u'/overlay;run; proc arima data=arima; identify var=y; identify var=y1; estimate q=1 method=ml; forecast out=a2 back=0 lead=10 id=t;run; proc plot data=a2(firstobs=90); plot y1*t='*' forecast*t='f' l95*t='l' u95*t='u'/overlay;run;