各位大虾:
我是一位菜鸟,正在学sas,下面是工作中需用的一段分析代码。名称叫:线性加平台分析。它是一个肥料效应试验分析,施肥量在一个捌点值以下,产量呈线性增加,施肥量大于捌点值,产量不再增加(达到最佳产量)。下面代码计算出了在捌点值以下的线性的截距、线性系数和捌点值、平台产量。
关键语句俺不懂:
DER.A = 1;
DER.B = AA;
DER.K = 0;
请大虾明示!!!!
全部代码如下:
DATA A;
INPUT aa ry;
cards;
0
600
0
621
0
613
120
801
120
830
120
850
180
950
180
901
180
934
240
952
240
938
240
923
;
TITLE1
'Yield Response Data';
TITLE2
'LRP Model';
PROC
NLIN
DATA=A BEST=3;
PARMS A=600 to 1000 by 50 B=0 to 30 by 1 K=0 to 250 by 10;
IF AA<K THEN DO;
MODEL RY = A + B*AA;
DER.A = 1;
DER.B = AA;
DER.K = 0;
END;
ELSE DO;
MODEL RY = A + B*K;
DER.A = 1;
DER.B = K;
DER.K = B;
END;
FILE
PRINT;
IF _OBS_ = 1 AND _MODEL_= 0
THEN DO;
PLATEAU = A + B*K;
PUT PLATEAU=;
END;
PLATEAU = A + B*K;
ID PLATEAU;
OUTPUT
OUT=NEW P=PRY PARMS=A B K SSE=SSE;
RUN;
PROC
PLOT;
PLOT RY*AA='+' PRY*AA='*'/OVERLAY;
RUN;