reg_nlin6_11.sas

黄燕 吴平《SAS统计分析及应用》 书中一些例子的SAS程序代码

proc nlin data=data_prg.nlin6_11 method=newton;
parms a=1.7 to 2 by 0.05
      b=1.5 to 2.0 by 0.01
      c=0.5 to 0.8 by 0.01 
      d=0.1 to 0.2 by 0.005;
temp=(x/b)**c;
temp1=temp+1;
temp2=exp(x/b);
model  y= (a-d)/temp1+d  ;
der.a=1/temp1;
  der.b=((a-d)*c*temp**2)/(b*temp1**2);
  der.c=((d-a)*temp*temp2)/(temp1**2);
  der.d=1-1/temp1;
  OUTPUT OUT=data_prg.bb2 P=yhat R=resid SSE=scrs L95M=lm U95M=um;
run;
proc plot data=data_prg.bb2;
plot y*x='a' yhat*x='*'/overlay HAXIS=0.001 0.01 0.1 1 10 100 1000 10000;
run;
quit;