例4.2.sas

第一章 时间序列分析简介 第二章 时间序列的预处理 第三章 平稳时间序列分析 第四章 非平稳序列的确定性分析 第五章 非平稳序列的随机分析 第六章 多元时间序列分析

data a;
input index@@;
time=intnx('month','1jan1991'd,_n_-1);
format time year2.;
t=_n_;
t2=t**2;
cards;
130.44	133.47	120.19	113.94	114.83	137.56	143.80	178.43	180.92
218.60	259.60	292.75	313.24	364.66	381.24	445.38	1234.71	1191.19
1052.07	823.27	702.32	507.25	724.60	780.39	1198.48	1339.88	925.91
1358.78	935.48	1007.05	881.07	895.68	890.27	814.82	984.93	833.80
770.25	770.98	704.46	592.56	556.26	469.29	333.92	785.33	791.15
654.98	683.59	647.87	562.59	549.26	646.92	579.93	700.51	630.58
695.55	723.87	722.43	717.32	641.13	555.29	537.34	552.93	556.39
681.16	643.65	804.25	822.48	809.94	875.52	976.71	1032.95	917.02
964.74	1040.27	1234.62	1393.75	1285.18	1250.27	1189.76	1221.06	1097.38
1180.39	1139.63	1194.10	1222.91	1206.53	1243.01	1343.44	1411.20	1339.20
1316.91	1150.22	1242.09	1217.31	1247.42	1146.70	1134.67	1090.09	1158.05
1120.92	1279.32	1689.42	1601.45	1627.12	1570.70	1504.56	1434.97	1366.58
1535.00	1714.58	1800.22	1836.32	1894.55	1928.11	2023.54	2021.20	1910.16
1961.29	2070.61	2073.48	2065.61	1959.18	2112.78	2119.18	2213.18	2218.03
1920.32	1834.14	1764.87	1689.17
;
proc gplot;
plot index*time=1;
symbol1 c=black v=none i=join;
proc reg;
model index=t t2;
model index=t2;
output out=out p=index_cup;
proc gplot data=out;
plot index*time=1 index_cup*time=2/overlay;
symbol2 c=red v=none i=join w=2 l=3;
run;