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;



雷达卡


京公网安备 11010802022788号







