nlin6_14.sas

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

title '负指数函数：Y=B0*(1-EXP(B1*X))' ;
DATA NLIN6_16;
  INPUT X Y @@;
  CARDS;
020  0.57  030  0.72  040  0.81  050  0.87  060  0.91  070  0.94  080  0.95
090  0.97  100  0.98  110  0.99  120  1.00  130  0.99  140  0.99  150  1.00
160	1.00  170  0.99  180  1.00  190  1.00  200  0.99  210  1.00 
 ;
PROC  NLIN  BEST=10  METHOD=MARQUARDT ;
  PARMS  B0=0  TO  2  BY  0.5     B1=0.01  TO  0.09  BY  0.01 ;
  MODEL Y=B0*(1-EXP(-B1*X)) ;
  DER.B0=1-EXP(-B1*X) ;
  DER.B1=B0*X*EXP(-B1*X) ;
OUTPUT  OUT=B  P=YHAT  R=YRESID ;
RUN;
PROC PLOT DATA=B ;
  PLOT   Y*X='A'  YHAT*X='P'  /OVERLAY  VPOS=25 ;
  PLOT  YRESID*X / VREF=0  VPOS=25 ;
RUN;