proc iml 运行出错，大侠看看哦

shiyiming

[code:xsa3g8dd]data a1;
input t year xt@@;
        yt=xt;
        index=1;
        zt=-(yt+lag(yt))/2;
        jbi=lag(xt)/xt;
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
;[/code:xsa3g8dd]
[code:xsa3g8dd]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};
                create a2 from ahatt[colname=na];
                append from ahatt;
quit;[/code:xsa3g8dd]
[color=#FF0000:xsa3g8dd][size=150:xsa3g8dd]请教大侠，错在那里了，期望能够改正,不胜感激！！[/size:xsa3g8dd][/color:xsa3g8dd]

shiyiming

Are you willing to do linear regression? It sounds.
JingJu

[code:sz0n9hh3]...
ahat=inv(t(b)*b)*(t(b)*yn);

or

ahat=solve(b`*b, b`*yn);
...[/code:sz0n9hh3]

shiyiming

第二段程序中，
[code:3dyugrif] .....
na={a u};
......
[/code:3dyugrif]
是否表示这段需要计算的是a&u值，并且输出在数据集a2中。
若不是的话，表示什么意思呢？
怎么才能显示出a&u值？

shiyiming

to supply variable name(s).
it is equal to
[code:146uyw9e]create a2 from ahatt[colname='a'];[/code:146uyw9e]
since only 1 variable available in ahatt.

Jingju