例3.9.sas

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

goptions vsize=7cm hsize=10cm;
data a;
input change_temp@@;
dif=dif(change_temp);
year=1880+_n_-1;
cards;
-.40    -.37    -.43    -.47    -.72    -.54    -.47    -.54    -.39    -.19
-.40    -.44    -.44    -.49    -.38    -.41    -.27    -.18    -.38    -.22
-.03    -.09    -.28    -.36    -.49    -.25    -.17    -.45    -.32    -.33
-.32    -.29    -.32    -.25    -.05    -.01    -.26    -.48    -.37    -.20
-.15    -.08    -.14    -.13    -.12    -.10     .13    -.01     .06    -.17
-.01     .09     .05    -.16     .05    -.02     .04     .17     .19     .05
 .15     .13     .09     .04     .11    -.03     .03     .15     .04    -.02
-.13     .02     .07     .20    -.03    -.07    -.19     .09     .11     .06
 .01     .08     .02     .02    -.27    -.18    -.09    -.02    -.13     .02
 .03    -.12    -.08     .17    -.09    -.04    -.24    -.16    -.09     .12
 .27     .42     .02     .30     .09     .05
;
proc print;
proc gplot;
plot dif*year;
symbol v=star i=join c=black;
proc arima data=a;
identify var=dif nlag=12;
estimate p=1 q=1;
run;