数学建模社区-数学中国

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

作者: hfut-liu    时间: 2009-9-1 20:08
标题: SAS时间序列分析遇到参数检验不合格怎么办
sas做时间序列预测比较好,有一些数据我先用差分进行处理然后去做预测,发现白噪声和残差自相关检验都通过,但是模拟参数你和不好,其中常数不显著,我就在程序p=3后加了noint,在运行,其中AR1,2通不过,该添加什么程序,或者用其他方法??/ M, c' Q' I) w4 s( u
程序如下:6 d8 d) \3 _7 d! _# t9 g
data example5_1;
1 Y# g0 M: W- m. @4 [$ ?input x@@;  P/ `+ f+ Q, @  x# B( c
difx=dif(x);
* p" k% h1 Z7 @. Ut=_n_;
! ]! ]4 i2 k4 \& D2 dcards;, x8 ^% Y" A& p" e
1601 5421 1890 4439 1703 3232 376 1167 1897 3737 1807 1628 1723# p2 X! v( G4 B8 o7 [
2584 1551 2479 1199 4148 2449 2026 1690 3374 2015 2480 850 2249& w5 h& H" o5 \% I% j5 a" ]3 p* g8 V
1674 3666 2029 1238
) P9 f; S6 ?( }7 B% e; W, y;
5 p0 ?2 }3 B$ @2 `7 w5 t3 O) m' ~; [proc gplot;, T1 L- T5 P6 F
plot x*t difx*t;
  h* |2 k: x! C$ }8 G5 Msymbol v=star c=black i=join;
1 U, y6 T4 H* }' Z% Oproc arima;
9 G7 u( y: l6 G3 }5 d/ w5 ridentify var=x(1) nlag=10 minic p=(0:5) q=(0:5);
+ c( K3 ~; p5 w! B+ westimate p=3 noint;  O7 [# n1 w! r8 Y/ e
forecast lead=7 id=t out=out;0 \- X+ _2 K7 F4 N4 `5 P
proc gplot data=out;
; o' ~  D8 E  a( gplot x*t=1 forecast*t=2 l95*t=3 u95*t=3/overlay;' `# z' E" y( F) t! X5 l  p. D
symbol1 c=black i=none v=star;
; A* X! S3 e9 Q  ~; E3 Wsymbol2 c=red i=join v=none;
+ }5 e/ a, D& jsymbol3 c=green I=join v=none;
; ?. \2 N$ O# v) y9 W" W; _3 zrun;
: b0 z4 p) h! N$ J. B
( o$ o! @' F7 A4 c! e/ h& C部分结果分析如下(参数那一部分):1 U4 i3 q& \3 ^+ [4 y
Conditional Least Squares Estimation5 s2 }( [6 s1 U: h9 V' n1 Q

) _* x; ?1 c+ a9 a                                                            Standard                 Approx4 r% J. w' n9 @' x6 @, v
                               Parameter      Estimate         Error    t Value    Pr > |t|     Lag! T7 I# P. G5 H! V
' C- I/ b( L. L# K" H# K0 M
                               AR1,1          -0.69510       0.18063      -3.85      0.0007       1
1 L, }3 r; o' e6 H                               AR1,2          -0.35875       0.22142      -1.62      0.1173       24 d& j/ n4 Q+ f+ l4 ~: f( M
                               AR1,3          -0.49187       0.18485      -2.66      0.0132       33 ]2 c6 O" f' k/ ?( S6 j
) Z) w/ \& G8 J! `! K" ?
其中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