nonlin3.sas

SAS for Windows统计分析系统教程,里面有源码

options nodate nonumber;
data nonlin3;
  input x y @@;
cards;
 7.5 121.2  8.5 125.7  9.5 130.7 10.5 135.3 11.5 139.9 12.5 145.2
13.5 151.8 14.5 158.3 15.5 163.8 16.5 167.0 17.5 168.6
;
proc nlin best=50 method=marquardt/* iteration method */;
  der.b0=1;
  der.b1=x;
  der.b2=x*x;
  der.b3=x*x*x;
  model y=b0+b1*x+b2*x*x+b3*x*x*x/* model y=b0+b2x+b2x*x+b3x*x*x */;
  output out=a predicted=yhat residual=yresid;
  parameters b0=150 to 170
             b1=-20 to -10
             b2=0 to 2
             b3=-1 to 0;
proc plot data=a;
  plot y*x='0' yhat*x='*'/overlay vpos=25;
run;