example4_3.sas

来自「第一章 时间序列分析简介 第二章 时间序列的预处理 第三章 平稳时间序列分析 第」· SAS 代码 · 共 29 行

data example4_3;                                                                                                                             
input x@@;                                                                                                                              
t=intnx('quarter','1jan1978'd,_n_-1);                                                                                                   
format t yyq4.;                                                                                                                         
cards;                                                                                                                                  
40777      41778      43160      45897                                                                                                  
41947      44061      44378      47237                                                                                                  
43315      43396      44843      46835                                                                                                  
42833      43548      44637      47107                                                                                                  
42552      43526      45039      47940                                                                                                  
43740      45007      46667      49325                                                                                                  
44878      46234      47055      50318                                                                                                  
46354      47260      48883      52605                                                                                                  
48527      50237      51592      55152                                                                                                  
50451      52294      54633      58802                                                                                                  
53990      55477      57850      61978                                                                                                  
;                                                                                                                                       
proc x11 data=example4_3;                                                                                                                               
quarterly date=t;                                                                                                                       
var x;                                                                                                                                  
output  out=out  b1=x d10=season d11=adjusted d12=trend d13=irr;                                                                                                      
proc gplot data=out;  
plot season*t=2 trend*t=2 irr*t=2;  
plot x*t=1 adjusted*t=2/overlay; 
symbol1 c=black i=join v=star;                                                                                                          
symbol2 c=red i=join v=none w=2;                                                                                                            
run;