可直线化曲线回归分析

1、非线性回归分析与SAS的智能化实现实验报告实验名称 第 3章班 级 数学与应用数学可直线化的曲线回归分析一、上机操作与练习1、曲线直线化分析、研究免疫球蛋白A 做火箭电泳，测得的火箭高度，x 为免疫球蛋白 A 的浓度， y 为火箭高度。（1）程序：datanr3_1;/*命名将要建立的数据集为nr3_1*/inputx y ;/*要输入变量为x、 y, 并且连续输入*/x1=log10(x);x2=1/x;x3=log(x);y1=1/y;y2=log(y);cards ;/*数据块开始 */0.27.6 0.4 12.30.6 15.7 0.8 18.21.0 18.7 1.2 21.41.

2、4 22.6 1.6 23.8;/*数据块结束*/run ;odshtml ;procgplot;/*绘制散点图 */ploty*x/haxis=0 to1.6/*散点图中点表示在以y 为纵轴 ,x 为横坐标的直角坐标系中by0.4vaxis=0 to25 by5 ;symbolvalue=dot;run ;procreg ;/*标志 reg 过程的开始, 进行含有截距项的直线回归分析*/model y=x1;model y1=x2;model y2=x3;run ;odshtmlclose;*/（2）运行结果及其解释：该图是根据原始资料绘制的 IgA 浓度和火箭高度之间的数据对用 sas程序

3、运行得到的散点图，可发现其趋势与对数曲线、双曲线或幂函数曲线很相似，故对两变量进行相应变换。上图是根据原始资料绘制的IgA 浓度和火箭高度之间的数据对用曲线专家画出的散点图，可发现其趋势与对数曲线很相似，故对两变量进行相应变换。以上是结果变量火箭高度y 和 IgA 浓度 x 经对数变换后产生的新变量x1 之间进行直线回归分析的结果，给出了拟合的直线回归方程中的截距项和斜率的估计值，以及对它们与零的差别是否具有统计学意义进行假设检验的结果。可以看出，截距和斜率与零的差别均有统计学意义。因此，直线回归方程为：y=19.47272+17.15637x。因 x1=log(x)，故y=19.47272+

4、17.15637 log(x)以下 是结果变量火箭高度y 和 IgA 浓度 x 经倒数变换后产生的新变量y1 和 x2 之间进行直线回归分析的结果。可以看出，截距和斜率与零的差别均有统计学意义。因此，直线回归方程为： y1=0.03122+0.02006x2。因 x2=1/x ，y1=1/y ，故 y=x/(0.03122+0.02006x)以上是结果变量火箭高度y 和 IgA浓度 x 经倒数变换后产生的新变量y2 和 x3 之间进行直线回归分析的结果。可以看出，截距和斜率与零的差别均有统计学意义。因此，直线回归方程为： y2=2.95526+0.55655x3。因 x3=Inx，y2= In

5、x y，故y=exp(2.95526+0.55655Inx)2、比较上例三条回归曲线的拟合效果免疫球蛋白浓度x 和火箭高度y 的关系（1）程序：datanr3_2;/*命名将要建立的数据集为nr3_2*/inputx y ;/*要输入变量为x、 y,并且连续输入*/cards;/* 直接输入数据，数据块开始*/y1= 19.74512+17.90734*log10(x);y2=x/(0.03029*x+ 0.02029 );y3=19.32481*x*0.53667;residual1=y1-y;residual2=y2-y;residual3=y3-y;ssrs1=residual1*2;s

6、srs2=residual2*2;ssrs3=residual3*2;cards;0.27.60.4 12.30.6 15.70.8 18.21.0 18.71.2 21.41.422.61.6 23.8;run ;symbol1color=greenvalue =dot;symbol2color=redvalue=pointline=1 interpol=spline;symbol3color=bluevalue=pointline= 2interpol=spline;symbol4color=greenvalue =pointline=3 interpol=spline;legend1ac

7、ross=1 cborder=blacklabel=none value =( ' 实际值 ''对数曲线 ''双曲线'' 幂曲线 ' )position=(bottom right inside)mode=protect;procgplotdata=nr3_2;ploty*xy1*xy2*xy3*x/haxis =0 to1.6by0.4 vaxis =0 to 25by5 overlaylegend =legend1;run ;procsql;createtabletest1asselectcss(y)astotal,sum(s

8、srs1)assum_resi1,sum(ssrs2)assum_resi2,sum(ssrs3)as sum_resi3fromnr3_2;quit;datatest;settest1;s1=sqrt(sum_resi1/(&n-2);s2=sqrt(sum_resi2/(&n-2);s3=sqrt(sum_resi3/(&n-2);r_square1=1-sum_resi1/total;r_square2=1-sum_resi2/total;r_square3=1-sum_resi3/total;run ;procprintdata =nr3_2;run ;proc

9、printdata =testnoobs ;run ;（2 ）运行结果及其解释：下图是对数曲线、双曲线和幂函数曲线对上例免疫球蛋白浓度x 和火箭高度 y 的关系拟合效果的展示图。以上是 3 条曲线对结果变量的拟合情况。其中， x，y 分别代表原始资料中的原因变量与结果变量， residual1，residual,2，residual3 分别表示对数曲线、双曲线和幂函数曲线拟合资料的残差，即估计值与真实值之差， ssrs1，ssrs2， ssrs3为对应的残差平方。以上是衡量曲线拟合情况的几个统计量的值。其中，total 为变量 y 的离均差平方和，sum_resi1， sum_resi2， s

10、um_resi3 依次表示对数曲线、双曲线和幂函数曲线拟合资料的残差平方和，s1，s2，s3 为对应的剩余标准差， r_square1，r_square,2，r_square3 分别代表对应的相关系数。（3）CurveExpert的拟合)stinu(轴YS = 0.74528715r = 0.9920756125.4222.1818.9415.7012.469.225.980.10.30.60.91.21.51.7X 轴 (units)用户模型 : y=a*xbCoefficient Data:a =1.92550057763E+001b =4.94059120119E-001（4）总结：三个

11、曲线回归方程所含参数个数相同，其拟合资料的残差平方和分别为1.64615,1.60966,4.50771。比较最大残差平方和4.50771 与最小残差平方和1.60966，由于每个曲线回归方程中均有两个待估计参数，故其自由度均为6。F=4.50771/1.60966=2.80，由 F 临界值表可知， F=2.80<5.82，所以残差平方和最大和最小的曲线回归方程拟合此资料的效果没有统计学差异，因而三个曲线回归方程拟合此资料的效果没有统计学差异。但是，由于双曲线回归方程对应的残差平方和1.61 比其他两个曲线方程的小，故此资料采用双曲线来拟合效果更好。3、指数曲线直线化分析红铃虫产卵数y

12、和温度 x 的关系（1 ）程序：datanr3_3;inputx y;y1=log10(y);cards;21 7231125 2127 24 29 66 32 11535 325;run;odshtml;procgplot;plot y*x/haxis =21 to 35 by 2 vaxis =0 to 350 by 50;symbolvalue =dot;run;procreg ;modely1=x;run;（2）运行结果及其解释：上图是根据原始资料绘制的产卵数和温度之间数据对用SAS 程序运行得到的散点图，可发现其趋势与指数曲线很相似，故对产卵数进行以10 为底的对数变换。Logist

13、ic 模型 : y=a/(1+b*exp(-cx)Coefficient Data:a = -4.44359029217E+002b = -7.26278315976E+003c = 2.29376490358E-001上图是根据原始资料绘制的产卵数和温度之间数据对用曲线专家画出的散点图，可发现其趋势与指数曲线很相似，故对产卵数进行以10 为底的对数变换， 该模型为 logistic模型。由运行结果可以看出， 结果变量产卵数y 经对数变换后产生的新变量y1 和温度 x 之间进行直线回归分析，截距和斜率与零的差别均有统计学意义。因此，直线回归方程为：y1=-1.67168+0.11814x4、比

14、较上例两条条回归曲线的拟合效果红铃虫产卵数y 和温度x 的关系（1）程序：%letn=7;datanr3_4;inputx y;y1=0.0213 * 10 *(0.11814*x);residual1=y1-y;ssrs1=residual1*2;cards ;21 7231125 2127 24 29 66 32 11535 325;run ;symbol1color=green value =dot;symbol2color=redvalue=pointline =1interpol=spline;legend1across=1cborder=blacklabel=none value

15、=( ' ê ? ê ? ''? ê y? ú ?')position=(bottom right inside)mode=protect;procgplotdata=nr3_4;ploty*x y1*x/haxis=21 to35 by2 vaxis = 0 to 350 by50 overlaylegend =legend1;run;procsql;createtabletest1asselectcss(y)as total,sum(ssrs1)as sum_resi1from nr3_4;quit;datatest;

16、settest1;s1=sqrt(sum_resi1/(&n-2);r_square1=1-sum_resi1/total;run;procprintdata=nr3_4;run;procprintdata=testnoobs ;run;（2）运行结果及其解释：该图是指数曲线y=0.0213100.11814x对上例资料的拟合效果展示图。 以下是所得指数曲线对结果变量的拟合情况。其中， x，y 分别代表原始资料中的原因变量与结果变量， residual1 表示指数曲线拟合资料的残差，即估计值与真实值之差，ssrs1为对应的残差平方。从其拟合资料的残差平方和1534.72、剩余标准差 1

17、7.5198 及相关指数 0.98036 这三个统计量的值可以看出，所得曲线对资料的拟合效果比较理想。5、指数曲线直线化分析某县年月累计麻疹的发病情况（1）程序：datanr3_5;inputx y;y1=log(19.34-y)/y);cards;10.27920.92032.07844.39358.144612.174717.306818.435918.7281018.9651119.2301219.328;run;odshtml;procgplot;ploty*x/haxis=0 to12 by 3 vaxis =0 to 20 by 2;symbolvalue=dot;run;proc

18、reg ;modely1=x;run;odshtmlclose;（2）运行结果及其解释： 下图是根据原始资料绘制的麻疹累计发病率和月份之间数据对用sas 程序运行得到的散点图，可发现其趋势近似S 型，故对麻疹累计发病率进行logistic 变换。观察散点图，可发现下渐近线约为y=0，,上渐近线约为y=19.34，两条渐近线之间的距离约为19.34，所以取 L=0，K=19.34，对其进行变量变换。下图是根据原始资料绘制的麻疹累计发病率和月份之间数据对用曲线专家画出的散点图，可发现其趋势近似S 型，故对麻疹累计发病率进行logistic 变换，该模型为Richards模型。)stinu(轴YS

19、= 0.36166389r = 0.9992515721.2317.7014.1610.637.103.560.030.12.34.46.68.810.913.1X 轴 (units)Richards 模型 : y=a/(1+exp(b-cx)(1 /d)Coefficient Data:a =1.91496505259E+001b =9.00929738629E+000c =1.44722142155E+000d =2+000由运行结果可以看出， 结果变量麻疹累计发病率y 经 logistic 变换后产生的新变量y1和月份 x 之间进行直线回归分析，截距和斜率与零的

20、差别均有统计学意义。因此，直线回归方程为：y1=5.09934-0.97293x。6、比较上例 Logistic 曲线拟合效果某县年月累计麻疹的发病情况（1）程序：%letn=12;datanr3_6;inputx y;y1=19.34 /( 1+163.91369*exp(-0.97293*x);residual1=y1-y;ssrs1=residual1*2;cards ;10.27920.92032.07844.39358.144612.174717.306818.435918.7281018.9651119.2301219.328;run;symbol1color=greenvalue

21、=dot;symbol2color=redvalue=pointline=1interpol=spline;legend1across=1 cborder=blacklabel=none value =( ' ê ? ê ? ''logistic?ú ?' )position=(bottom right inside)mode=protect;procgplotdata =nr3_6;ploty*x y1*x/haxis =0 to12 by3 vaxis =0 to20 by 2overlaylegend=legend1;run;procsql;createtabletest1asselectcss(y)as total,sum(ssrs1)as sum