数学建模社区-数学中国

标题: SAS时间序列分析遇到参数检验不合格怎么办 [打印本页]
作者: hfut-liu    时间: 2009-9-1 20:08
标题: SAS时间序列分析遇到参数检验不合格怎么办
sas做时间序列预测比较好，有一些数据我先用差分进行处理然后去做预测，发现白噪声和残差自相关检验都通过，但是模拟参数你和不好，其中常数不显著，我就在程序p=3后加了noint,在运行，其中AR1,2通不过，该添加什么程序，或者用其他方法？？
- b. _1 F* ]+ q% d程序如下：2 a2 g& N0 ?5 x3 l
data example5_1;! m2 v7 M6 r4 j; d7 A" X4 W
input x@@;# s  G3 h! J# V$ ?1 C
difx=dif(x);
  S. [5 u- l0 G( Mt=_n_;
7 g2 y  G/ i2 h/ [$ @cards;" L6 b; c: K. q/ @& f) k
1601 5421 1890 4439 1703 3232 376 1167 1897 3737 1807 1628 17232 H5 t% I5 U  t; R
2584 1551 2479 1199 4148 2449 2026 1690 3374 2015 2480 850 2249
8 f& i$ p3 s) s9 D1674 3666 2029 1238
6 p0 A1 K, [# K4 Q6 p  ~;
/ h% }- q- f1 E/ W$ gproc gplot;8 g, U+ M: U3 l7 Y# p1 H, B
plot x*t difx*t;
  {6 A" w1 b0 u- x" I+ vsymbol v=star c=black i=join;
/ t; @  a. [5 g6 t7 wproc arima;
  F' l+ G0 B% f/ K; G7 z9 ?identify var=x(1) nlag=10 minic p=(0:5) q=(0:5);# S/ j( R9 H* ]) X. h" E# Y
estimate p=3 noint;! G: F: {+ {" e9 q) F6 \
forecast lead=7 id=t out=out;
4 F1 J  n. e' i; {0 c4 o. }* Mproc gplot data=out;
5 o2 f5 @7 g2 `3 b) s4 yplot x*t=1 forecast*t=2 l95*t=3 u95*t=3/overlay;
* ]2 w/ p( W& J. Esymbol1 c=black i=none v=star;
, g1 e7 _- h6 g3 Fsymbol2 c=red i=join v=none;
  Q& z5 \, R/ Q, i3 N* r# @symbol3 c=green I=join v=none;$ H9 n  d( |* h1 s( \! l- I
run;
- O" B$ b9 @+ {! ~" V/ b0 H+ u8 u
- d4 o- F. m  M2 H部分结果分析如下（参数那一部分）：
7 k5 r5 E3 a5 p) V4 m1 M& X Conditional Least Squares Estimation* k6 j: I: t. Q7 U, Y2 _4 g
6 ^' T" |" a/ |' x! L  d1 z; v
                                                            Standard                 Approx
, O2 n9 o* c, U; T' C# A                               Parameter      Estimate         Error    t Value    Pr > |t|     Lag
$ Q; b" D* V' N9 h) ]5 Z# u4 \
9 E$ L1 R) o: q: c                               AR1,1          -0.69510       0.18063      -3.85      0.0007       1
' D" \4 v& c; s/ `                               AR1,2          -0.35875       0.22142      -1.62      0.1173       21 u& c) Y2 X& d* j9 a
                               AR1,3          -0.49187       0.18485      -2.66      0.0132       3
, a; a, [' a" ]6 Q( S* y0 B+ c- z+ @# |, @( a1 I; l
其中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