|
|
地板
楼主 |
发表于 2008-8-1 11:26:03
|
只看该作者
Re: 用SAS预测时间序列的问题
sas编程并不麻烦的,你可以参考一些时间序列书的,
一般时间序列编程,通常会调用一个过程函数,比如说arima模型,reg模型等,
你可以通过编程把数据导入,也可以通过其它方式导入, 然后编程,具体编程语言要看你对数据进行时间序列预处理后,它显现出的规律来作了.
你好象只是对销售额这一个变量做时间序列分析,它用到的程序会简单一些的,如果是几个变量做时间序列分析就会很麻烦的了
我给你一个时间序列的编程模型 时间序列编程大概就是这样的 这个程序是拟合garch模型的
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; |
|
