我想用方程h=1.3+[(a1+a2*m+a3*n)exp(b/d)]  求出AIC和BIC以及-2LogLikelihood 的值 请高手指导一下我下面编的程序哪里不对 帮改一下 非常感谢！！！

proc import out=work.totalheight
            datafile="f:\sas kl.xls"
                  dbms=excel2000 replace;
            getnames=yes;
run;
data taper ;set totalheight;
run;
PROC nlmixed data=totalheight;
parms a1=0 to 10  a2=0 to 10 a3=0 to 10 sb1=0 sb2=0 sb3=0 s=0.05;
fm1=exp(b/d);
fm2=((a1+c1)+(a2+c2)*m+(a3+c3)*n;
mea=1.3+(fm1*fm2);
model h~ normal(mean,s);
random c1 c2  ~normal([0,0],[sb1,sb2,sb3])subject=plot;
predict mea out=p;
run;