我的编程分为两部分，第一部分是回归：
proc reg data=BA;
model lnwage=exp exp2 exp3 exp4 schooling;
run;

然后我想把回归出来的结果转成矩阵，于是第二部分编程为：
data BA;
set BA;
keep exp exp2 exp3 exp4 schooling;
run;

proc iml;
use BA; read all into bu;

do j=1 to 42;
bi=bu[1,1]+bu[1,2]*j+bu[1,3]*j*j*0.01+bu[1,4]*j*j*j*0.001+bu[1,5]*j*j*j*j*0.0001+16*bu[1,6];

bu[1,2]对应exp,bu[1,3]对应exp2,bu[1,4]对应exp3,bu[1,5]对应exp4,bu[1,6]对应schooling,
但是keep命令部分出现了问题，我始终不能保留从回归得到的intercept部分，也就是bu[1,1]。
想请高人解答如何可以使用Keep命令保留Intercept部分，或者Drop命令也可以，谢谢。