proc nlmixeddata=e.shanmuf;
parms a0=1a1=0.99 a2=0.38a3=0.49 a4=1.8s2u=200 s2e=3s3u=200 s3u1=200;
c=a0+u1;
b=a2+u2;
y=1.3+c*Ht**a1*(1-exp(-b*Dg**(-a3)*d))**a4;
model h~normal(y,s2e);
random u1 u2 ~ normal([0,0],[s2u, s3u,s3u1]) subject=class;
predict y out=e.shanmufm6;
run;



雷达卡






京公网安备 11010802022788号







