全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
2926 3
2010-04-05
各位大虾:
  我是一位菜鸟,正在学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;

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-4-5 16:14:15
分析结果如上图
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-4-5 23:21:48
管家 发表于 2010-4-5 16:11
各位大虾:
  我是一位菜鸟,正在学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



DER. is used to get the partial derivatives.
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-4-7 10:32:13
谢谢指点

请问:改良高斯-牛顿法的具体算法(或提供学习资料)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群