nlin6_14.sas

来自「黄燕 吴平《SAS统计分析及应用》 书中一些例子的SAS程序代码」· SAS 代码 · 共 20 行

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;