数学建模社区-数学中国

标题: SAS时间序列分析遇到参数检验不合格怎么办 [打印本页]

作者: hfut-liu    时间: 2009-9-1 20:08
标题: SAS时间序列分析遇到参数检验不合格怎么办
sas做时间序列预测比较好,有一些数据我先用差分进行处理然后去做预测,发现白噪声和残差自相关检验都通过,但是模拟参数你和不好,其中常数不显著,我就在程序p=3后加了noint,在运行,其中AR1,2通不过,该添加什么程序,或者用其他方法??+ O# e" t0 |1 Y% I
程序如下:; p3 S" L/ I6 P; u$ _) t
data example5_1;/ Z. M/ v0 e+ W2 e. s: n
input x@@;
9 l' h: t. N" W# O- B/ [; Bdifx=dif(x);
2 p) h9 a0 _  At=_n_;" ^9 V5 V' P3 ]9 z# V- t
cards;! g) K+ L3 H* Z0 m* I; N+ C
1601 5421 1890 4439 1703 3232 376 1167 1897 3737 1807 1628 1723
8 T' J3 Y) M$ H. J2 N& r2584 1551 2479 1199 4148 2449 2026 1690 3374 2015 2480 850 22492 h4 e7 _# q* e1 T
1674 3666 2029 1238
4 t- Q: B1 q4 n, p9 w0 f; C;
* `( G  T- s' h# Mproc gplot;1 v3 V; T! {* }% I
plot x*t difx*t;
' e5 e8 {* ]( [' ~, s9 j% E+ ?; Vsymbol v=star c=black i=join;; Z* b/ p. r% ^
proc arima;! L5 `2 P3 ~2 U0 G: m
identify var=x(1) nlag=10 minic p=(0:5) q=(0:5);
* ^1 j8 j8 W" X- W3 ^5 Qestimate p=3 noint;
( }! B- z: u& aforecast lead=7 id=t out=out;
$ z. I7 B' r3 I7 N" H. [6 P+ D/ Nproc gplot data=out;
, B- v7 D- }+ ]$ d2 p' d/ {+ c( Oplot x*t=1 forecast*t=2 l95*t=3 u95*t=3/overlay;7 d; K" S4 K; D" l/ Q( f4 |( w
symbol1 c=black i=none v=star;  m- w, I6 Z0 }' b3 D- q" o
symbol2 c=red i=join v=none;
. A) i9 L. U1 X/ b4 C  {9 }symbol3 c=green I=join v=none;  z( v2 y2 S5 Q) p& g3 ]0 L
run;
! U5 e7 i: }+ |( d. T9 h1 l3 ?: m, k! a% n* ~1 W0 s
部分结果分析如下(参数那一部分):/ K# P0 ]7 M% u& T- T$ w/ G
Conditional Least Squares Estimation
9 B/ y$ |/ i7 U! L9 s3 e2 u5 P. N( a* s, G& K$ V
                                                            Standard                 Approx
% u4 K9 E4 V9 ?                               Parameter      Estimate         Error    t Value    Pr > |t|     Lag1 M: m8 k% X$ X
5 P  u. t* f8 i3 n7 l4 U4 H) {
                               AR1,1          -0.69510       0.18063      -3.85      0.0007       10 f+ t, n0 B+ o7 d! ^
                               AR1,2          -0.35875       0.22142      -1.62      0.1173       2: ]8 z8 S" P* D' q
                               AR1,3          -0.49187       0.18485      -2.66      0.0132       3
' w; @6 A( d. Y, x3 a
- ^" G  l" E7 A2 @8 m; x2 T其中AR1,2是通不过检验的请假大家帮忙改怎么添加命令或者修改,谢谢
作者: hfut-liu    时间: 2009-9-2 14:20
给点意见吧
作者: C.J.    时间: 2011-8-11 15:31
estimate p=(1,3) noint;
作者: shuxuezaozhuang    时间: 2011-9-21 16:19
hao ma fan!!
作者: jeffreylhb    时间: 2011-9-30 11:16
努力回帖 努力作贡献
作者: rr_傻丫头    时间: 2012-8-11 16:48
这个软件怎么用哦……都没有头绪……




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5