用sas进行灰色预测..docVIP

data a1; input t year xt@@; yt+xt; /*生成一阶累加序列*/ index=1;zt=-(yt+lag(yt))/2; /*为数据矩阵B准备数据*/ datalines; 1 1990 24395 2 1991 25286 3 1992 26901 4 1993 27339 5 1994 27871 6 1995 28721 7 1996 29728 8 1997 30067 9 1998 30791 10 1999 31284 11 2000 33716 12 2001 34558 ; proc iml; use a1; read all var{zt index} into B where(zt^=.); read all var{xt} into yn where(zt^=.); ahat=inv(B`*B)*B`*yn; ahatt=ahat`; na={a u}; creat a2 from ahatt [colname=na]; append from ahatt; quit; data a3; set a2;index=1; data a4; set a1;if _n_=1;xt0=xt; keep xt0 index; data a5; merge a1 a3 a4; by index; if _n_=1 then xp=xt; e