example5_3.sas

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

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;