data a;
input year x y;
cards;
1987 357.4 12925
1988 364.79 14022
1989 371.39 15196
1990 380.88 16117
1991 386.1 16664
1992 392.16 21888
1993 397.08 23599
1994 402.71 16666
1995 411.11 17788
1996 416.47 18427
1997 422.31 18815
1998 425.98 18602
1999 429.95 18592
2000 438.18 19795
2001 442.16 20019
2002 448.1 20459
2003 456.6 19834
2004 444.68 20690
2005 455.7 20713
2006 469.85 23791
2007 478.9 24681
2008 486.73 25830
2009 491.43 27713
;
proc gplot data=a;
plot x*year=1 y*year=2;
symbol1 c=black i=join v=circle;
symbol2 c=black i=join v=star;
proc arima data=a;
identify var=x stationarity=(adf);
identify var=y stationarity=(adf);
proc arima data=a;
identify var=y crosscorr=x;
estimate input=x plot;
forecast lead=0 id=year out=result;
proc arima data=result;
identify var=residual stationarity=(adf=2);
run;
要求如题。我想探究y和x之间的关系,做了回归,检验了残差序列,如果有不妥当的地方,敬请指出!我会追加奖赏,谢谢大家!