7.2一元非线性回归与相关

《7.2一元非线性回归与相关》由会员分享，可在线阅读，更多相关《7.2一元非线性回归与相关（64页珍藏版）》请在金锄头文库上搜索。

1、7.2 一元非线性回归与相关一元非线性回归与相关 在实际应用中在实际应用中, 当当 x 的取值范围较广泛时的取值范围较广泛时, x与与y的关系并非都是直线的关系并非都是直线, 散点图往往呈现非线散点图往往呈现非线性关系性关系, 如施肥量如施肥量 x 和产量和产量 y, 光照强度光照强度 x 和光合和光合作用作用 y, 药剂浓度药剂浓度x 和害虫死亡率和害虫死亡率y 等等, 就呈现出就呈现出曲线关系。非线性的形式是多种多样的曲线关系。非线性的形式是多种多样的*, 本节本节只讨论生物学常遇到的非线性回归问题。只讨论生物学常遇到的非线性回归问题。* * 注注: : 一个恰当的非线性回归模型一个恰当的

2、非线性回归模型, ,需经过多次配合检验需经过多次配合检验, ,并要借助并要借助 于一定的专业知识才能确定。于一定的专业知识才能确定。7.282021/6/162021/6/161 1对幂函数直化的主要手段是两边取常用对数。对幂函数直化的主要手段是两边取常用对数。(或或自然对数自然对数), 取常用对数如下：取常用对数如下：1. 幂函数直化幂函数直化7.2.1 可化为直线的非线性回归可化为直线的非线性回归图图7.5 幂函数幂函数-1-1by = a x0b1b1 b=1y由对数性质得由对数性质得0 这种类型的回归方程一般只确定两个参数。这种类型的回归方程一般只确定两个参数。7.29幂函数式为幂函数

3、式为:(见图见图7.5)2021/6/162021/6/162 2然后再转换然后再转换 a, 即：即：lg a = a回归方程回归方程 - -1令令:可得幂函数的直化回归方程为可得幂函数的直化回归方程为 例例7.7 已测得云杉平均胸径已测得云杉平均胸径x(cm)与平均树高与平均树高 y(m)的资料的资料, ,试配合幂函数的回归方程试配合幂函数的回归方程 ;再试配合二次抛物线方程再试配合二次抛物线方程* 。7.30* * 注注: : 一元二次抛物线方程的配合过程见本章一元二次抛物线方程的配合过程见本章7.2.2所述。所述。所以可得曲线所以可得曲线2021/6/162021/6/163 3 x y

4、 lg(x) lg(y) 1 15 13.9 1.176 1.143 2 20 17.1 1.301 1.233 3 25 20.0 1.398 1.301 4 30 22.1 1.477 1.344 5 35 24.0 1.544 1.380 6 40 25.6 1.602 1.408 7 45 27.0 1.653 1.431 8 50 28.3 1.699 1.452 9 55 29.4 1.740 1.46810 60 30.2 1.778 1.48011 65 31.4 1.813 1.497表表7.4 云杉平均胸径与平均树高数据云杉平均胸径与平均树高数据图图7.6 云杉胸径树高散点

5、图云杉胸径树高散点图7.3115 20 25 30 35 40 45 50 55 60 65 3028 26 24 22 20 18 1614 .xy 对表对表7.4 x,y 画出散点图画出散点图初步判定为幂函数初步判定为幂函数, 分分别对别对 x, y 取常用对数得取常用对数得表表7.4 第第4, 5列列2021/6/162021/6/164 4 解：解：利用计算器统计功能按下面步骤计算利用计算器统计功能按下面步骤计算得直化后的回归方程得直化后的回归方程: :7.322021/6/162021/6/165 5所以可得幂函数的回归方程所以可得幂函数的回归方程 对曲线相关的研究只研究其相关指数。

6、曲线相对曲线相关的研究只研究其相关指数。曲线相关指数可定义为关指数可定义为对对a进行反对数转换得进行反对数转换得7.332021/6/162021/6/166 6 计算例计算例7.7相关相关指数（见表指数（见表7.5）15 13.9 14.63 - -0.73 0.5320 17.1 17.12 - -0.02 0.00 25 20.0 19.33 0.67 0.4530 22.1 21.36 0.74 0.5535 24.0 23.23 0.77 0.59 40 25.6 24.99 0.61 0.3745 27.0 26.65 0.35 0.1250 28.3 28.23 0.07 0.0

7、055 29.4 29.73 - -0.33 0.1160 30.2 31.18 - -0.98 0.96 65 31.4 32.57 - -1.17 1.37 xi yi yi (yi - - yi) (yi - - yi) 2 440 269 Q =5.06表表7.5 剩余平方和剩余平方和Q Q的计算过程的计算过程 注注: : 对曲线回归的检验较复杂对曲线回归的检验较复杂, ,主要是主要是Q的值较难求。的值较难求。7.342021/6/162021/6/167 7(1) 曲线关系检验曲线关系检验表表7.7 曲线回归的方差分析表曲线回归的方差分析表差异来源差异来源 平方和平方和 自由度自由度

8、 均方均方 均方比均方比 F0.01 回回 归归 314.307 1 314.307 559.27* 10.6 剩剩 余余 5.06 9 0.562 总总 和和 319.367 10Lyy =319.367 , Q = 5.06 U = 319.367 - - 5.06=314.3077.35 对回归关系的检验对回归关系的检验: :2021/6/162021/6/168 8 (2) 直化关系检验直化关系检验表表7.6 直化后的回归方差分析表直化后的回归方差分析表差异来源差异来源 平方和平方和 自由度自由度 均方均方 均方比均方比 F0.01 回回 归归 0.12478 1 0.12478 36

9、7 * 10.6 剩剩 余余 0.00308 9 0.00034 总总 和和 0.12786 107.36注：注：两种检验效果是一样的，如果曲线回归不能直化那么就只能用两种检验效果是一样的，如果曲线回归不能直化那么就只能用 曲线方法检验了。曲线方法检验了。2021/6/162021/6/169 92. 指数函数直化指数函数直化指数函数指数函数 , 两边取两边取自然对数自然对数可得可得(注意注意: 这里这里 x 变量并未转换）变量并未转换）令令:b0图图7.7 指数函数图示指数函数图示x y则直化方程为则直化方程为:7.372021/6/162021/6/161010 例例7.8 在栽培试验中在

10、栽培试验中, 测得某夏季绿肥播测得某夏季绿肥播种种15天后天后,每每5天天 x为一期的生长量为一期的生长量y(kg), 所得数所得数据如下据如下, 试对该资料进行回归与相关分析。根试对该资料进行回归与相关分析。根据散点图据散点图7.8选择指数函数选择指数函数模型模型 y = ae 来研究。来研究。bx图图7.8 绿肥绿肥5 5天一期生长量曲线天一期生长量曲线两边取自然对数两边取自然对数:直化方程直化方程:153045500100200400300.7.382021/6/162021/6/161111表表7.8 绿肥绿肥5 5天一期生长量数据天一期生长量数据 x y y= lny x y1 15

11、 58 4.060 60.902 20 67 4.205 84.103 25 79 4.369 109.23 4 30 140 4.941 148.235 35 200 5.298 185.436 40 320 5.768 230.727 45 480 6.173 277.78 210 1344 34.814 1096.367.392021/6/162021/6/161212回归方程的配合回归方程的配合:直化后的回归方程直化后的回归方程:7.402021/6/162021/6/161313取反对数得取反对数得:转曲线回归方程：转曲线回归方程：表表7.9 剩余平方和剩余平方和Q的计算过程的计算过

12、程 xi yi yi (yi - - yi) (yi- - yi ) 15 58 47.43 10.57 111.72 20 67 68.77 - -1.77 3.13 25 79 99.81 - -20.81 433.06 30 140 144.58 - -4.58 20.98 35 200 209.41 - -9.41 88.55 40 320 303.32 16.68 278.22 45 480 439.35 40.65 1652.42210 1344 2588.08相关指数计算：相关指数计算：27.412021/6/162021/6/1614143. 双曲线函数直化双曲线函数直化a 0

13、，b0(c)(b)yyxxbxa 0，b0b00b(a).y令令: 可直化为下式可直化为下式7.422021/6/162021/6/161515相关指数的意义和直线回归决定系的意义是一相关指数的意义和直线回归决定系的意义是一样的。这里不能用原始数据求直线相关系数样的。这里不能用原始数据求直线相关系数r, r 只能用于直线只能用于直线, 不能用于曲线。不能用于曲线。可得：可得：回归方城配合回归方城配合:对相关指数的说明对相关指数的说明:7.432021/6/162021/6/161616 例例7.9 研究研究“岱字岱字”棉自播至齐苗（以棉自播至齐苗（以80%出出苗为准）的天数（苗为准）的天数（x

14、）与日平均土温）与日平均土温(y)的关系的关系, 得得下表下表7.10, 试作回归分析。试作回归分析。表表7.10 “岱字岱字”棉至齐苗棉至齐苗 天数与土温的关系天数与土温的关系xy3330272421183.0034.005.008.0015.0yx4 8 12 16 40 30 20 10 0.图图7.10 “岱字棉岱字棉”至齐苗至齐苗 天数土温的关系图天数土温的关系图7.44322021/6/162021/6/161717表表7.11 “岱字岱字”棉至齐苗天数与棉至齐苗天数与 土温的关系数据变换土温的关系数据变换xyx= 1/x3330272421183.0034.005.008.001

15、5.00.33330.27270.25000.20000.12500.066610.99898.18106.75004.80002.62501.1988xy所选模型：所选模型：23直化式如下直化式如下7.45xy22 0.1111 0.0744 0.0625 0.0400 0.0156 0.0044 1089 900 725 576 441 324153 1.2436 34.5537 0.308 40592021/6/162021/6/161818直化回归方程：直化回归方程：曲线回归方程：曲线回归方程：相关指数：相关指数：7.46注注: : 本例还可用本例还可用 y=a+bx 转换转换, 其中

16、其中 y= xy , 其结果相同。其结果相同。 2021/6/162021/6/161919 在农业试验中经常遇到一些曲线回归配合问在农业试验中经常遇到一些曲线回归配合问题题, 现列出几种常用的可直化的函数模型现列出几种常用的可直化的函数模型*如如: :密度与有效穗密度与有效穗数间的关系。数间的关系。如如: :毛细管水的上升毛细管水的上升高度与时间的关系。高度与时间的关系。如如: :叶绿素与透叶绿素与透光度间的关系。光度间的关系。如如: :棉纤维长分布棉纤维长分布与次数间关系。与次数间关系。如如: :大豆密度与青荚大豆密度与青荚产量间关系。产量间关系。* * 注：注：农业试验统计农业试验统计莫

17、惠栋莫惠栋 编著编著 上海科技出版社上海科技出版社 1984如如: :洋葱鳞茎的直径洋葱鳞茎的直径与重量间的关系。与重量间的关系。7.472021/6/162021/6/1620204. S 形曲线直化形曲线直化 在生物界中在生物界中, 反应生物的累积发生量反应生物的累积发生量(率率)y与与时间、浓度、药品用量时间、浓度、药品用量 x 等关系时等关系时, 可以用可以用 S 形形曲线来刻画。如曲线来刻画。如: 动植物的累积发病率、累积死动植物的累积发病率、累积死亡率、累积发生进度、繁殖过程与时间、浓度与亡率、累积发生进度、繁殖过程与时间、浓度与药品用量的关系等。药品用量的关系等。图图7.11 对

18、称的对称的S S形曲线形曲线图图7.12 不对称的不对称的S S形曲线形曲线时间时间 浓度浓度 用量用量累累积积率率yx时间时间 浓度浓度 用量用量x累累积积率率y7.482021/6/162021/6/162121 (1) 对称的对称的S 形曲线形曲线: 一些作物的日生长量或一些作物的日生长量或发生量等往往为正态或近似正态分布发生量等往往为正态或近似正态分布(下图左下图左), 其累积量其累积量(率率) y%与时间与时间 x 则呈对称的则呈对称的S 形形曲线曲线(下图右下图右)。xy日日发发生生量量日日生生长长量量日期日期图图7.13 正态分布曲线正态分布曲线图图7.14 对称对称 s 形曲线

19、形曲线xy%日期日期累累积积率率累累积积量量7.492021/6/162021/6/162222 可将资料中的累积率可将资料中的累积率y%视为正态分布的分视为正态分布的分布函数布函数(u)=PUu, 若将若将y% 转化为正态分布的转化为正态分布的u值值, 则则u与时间与时间 x 就成为直线关系了就成为直线关系了, 但是但是u变变量的取值有时为正或为负计算起来很不方便量的取值有时为正或为负计算起来很不方便, 可可采用采用 P = u+ +5 来解决。则来解决。则 P =a+ +bx 。P 称概率单称概率单位值位值(不要和概率符号不要和概率符号P 相混淆相混淆), 在一些统计书在一些统计书中可查到

20、百分率与概率单位换算表中可查到百分率与概率单位换算表。7.50对称的对称的S 形曲线直化过程：形曲线直化过程：2021/6/162021/6/162323 例例7.10 研究某地越冬代棉花红铃虫的化蛹研究某地越冬代棉花红铃虫的化蛹进度进度 y%与时间与时间 x 的关系的关系, 试进行回归分析。试进行回归分析。 解解: 首先将观察日期首先将观察日期转换为以转换为以5月月31日为日为0的数的数值值(x), 化蛹进度化蛹进度 y%与与 x 成成 S 形曲线形曲线, 再将再将 y%查查表转换为概率单位值表转换为概率单位值P, 而而P与与x即为直线了。即为直线了。 5 10 5 10 15 20 25

21、30 35 40 45 5015 20 25 30 35 40 45 50 y%y%505025257575100100. . . . . . . . . . .图图图图7.17.15 5 棉红铃虫棉红铃虫棉红铃虫棉红铃虫s s形曲线图形曲线图形曲线图形曲线图x x 5 10 15 20 25 30 35 40 45 505 10 15 20 25 30 35 40 45 50 y%y%505025257575100100. . . . . . . .x x7.512021/6/162021/6/162424表表表表7.127.12 棉花红铃虫化蛹进度表棉花红铃虫化蛹进度表棉花红铃虫化蛹进度表

22、棉花红铃虫化蛹进度表 月月月月/ / / /日日日日 x x 天天天天 y%y% P P 概率单位概率单位概率单位概率单位 P P6 6/ / / /5 5 5 5 3.5 3.5 3.18813.1881 3.1264 3.1 3.1264 3.16 6/ / / /10 10 10 10 6.4 6.4 3.47803.4780 3.5627 7.5 3.5627 7.5 6 6/ / / /15 15 15 15 14.6 14.6 3.94633.9463 3.9991 15.8 3.9991 15.86 6/ / / /20 20 2020 31.4 31.4 4.51554.515

23、5 4.4354 28.6 4.4354 28.66 6/ / / /25 25 25 25 45.6 45.6 4.88954.8895 4.8718 44.9 4.8718 44.96 6/ / / /30 30 3030 60.4 60.4 5.26375.2637 5.3082 62.1 5.3082 62.17 7/ / / /5 5 3535 75.2 75.2 5.68085.6808 5.7445 77.2 5.7445 77.27 7/ / / /10 10 4040 90.2 90.2 6.29306.2930 6.1809 88.1 6.1809 88.17 7/ / /

24、 /15 15 4545 95.4 95.4 6.68496.6849 6.6172 94.7 6.6172 94.77 7/ / / /20 20 50 50 97.5 97.5 6.96006.9600 7.5360 98.0 7.5360 98.0 275275 50.899850.8998y%y%7.522021/6/162021/6/162525得回归方程得回归方程得回归方程得回归方程: : : :7.532021/6/162021/6/162626例如例如: 测报测报 7月月2日时即日时即 x=32, 代入回归方程代入回归方程得得P =5.4719 , 查表得查表得P对应的累积率为

25、对应的累积率为68.15%P概概率率单单位位值值P=2.756 + 0.0848x 1020304050x. .73456日期日期图图7.16 概率单位与日期直线关系图概率单位与日期直线关系图可计算直线回归的相关系数可计算直线回归的相关系数7.542021/6/162021/6/162727 (2) 不对称的不对称的 S 形曲线形曲线: 生物生长量或发生物生长量或发生量生量 y 与时间与时间 x 的曲线并不都是正态的曲线并不都是正态(图图7.17), 其累积率也呈现出不对称的其累积率也呈现出不对称的 S 形曲线形曲线(图图7.18), xyxy%发发生生量量日期日期日期日期生生长长量量图图7.

26、18 不对称不对称S S形曲线形曲线图图7.17 偏态曲线偏态曲线累累积积率率7.552021/6/162021/6/162828 例例7.11 用氯氰菊酯对菜粉蝶五龄幼虫进行触用氯氰菊酯对菜粉蝶五龄幼虫进行触杀毒力测定杀毒力测定(LD50), 累积死亡率累积死亡率y%与单位体重药与单位体重药量量 x ( / /L)之间的图形呈不对称的之间的图形呈不对称的 S 形曲线形曲线(见下见下图图7.19), 试对其直化后进行回归分析。试对其直化后进行回归分析。g如果将累积率如果将累积率 y% 转换成转换成 P, 则则P 和时间和时间 x 的关的关系仍然是曲线系仍然是曲线, 达不到直化的目的达不到直化的

27、目的, 前面的方法前面的方法不能直接搬用不能直接搬用, 但是只要将但是只要将 x 取对数转换取对数转换 (lg x), 则则P与与 lgx 就为直线了就为直线了, P = a + + blgx 。7.562021/6/162021/6/162929 解解: 将累积更正死亡率查百分率与概率单位将累积更正死亡率查百分率与概率单位值转换表得概率单位值值转换表得概率单位值P, 再将单位体重药量再将单位体重药量 转换转换为为lgx 建立直化回归方程建立直化回归方程 0.007999 19.8 0.01559 32.2 0.0376 57.8 0.0684 81.7 0.131 92.3 x(g/ /L)

28、 死亡率死亡率y% 表表7.13 菜粉蝶五龄幼虫毒杀力测定菜粉蝶五龄幼虫毒杀力测定12345图图7.19 氯氰菊酯触杀菜氯氰菊酯触杀菜 粉蝶毒力测定图粉蝶毒力测定图单位体重药量单位体重药量y%累累积积死死亡亡率率.10000.040.120.08806040200.16x7.572021/6/162021/6/163030 0.007999 - -2.097 19.8 4.1512 0.01559 - -1.807 32.2 4.5379 0.0376 - -1.425 57.8 5.1968 0.0684 - -1.165 81.7 5.9040 0.131 - -0.883 92.3 6.

29、4255 - -7.377 26.2154 平均平均 - -1.4754 5.2431 x( / /L) x= = lgx y y% P P 概率单位概率单位概率单位概率单位 g表表7.14 菜粉蝶五龄幼虫毒杀力测定数据菜粉蝶五龄幼虫毒杀力测定数据123457.582021/6/162021/6/1631317.592021/6/162021/6/163232即即: lgx概率单位值概率单位值 P- -0.8- -1.6- -2.0- -2.4.76543- -1.4图图7.20 直线回归方程直线回归方程 直线回归方程直线回归方程7.602021/6/162021/6/1633331 0.00

30、7999 - -2.097 19.8 4.1512 4.050 17.1 7.292 0.01559 - -1.807 32.2 4.5379 4.607 34.7 6.253 0.0376 - -1.425 57.8 5.1968 5.339 63.3 30.34 0.0684 - -1.165 81.7 5.9040 5.839 80.9 0.645 0.131 - -0.883 92.3 6.4255 6.379 91.6 0.49 - -7.377 26.215 44.92 平平 均均 - -1.475 5.2431 x( / /L) x=lgx y% Pg表表7.15 菜粉蝶五龄幼虫

31、毒杀力测定预测及相关指数计算表菜粉蝶五龄幼虫毒杀力测定预测及相关指数计算表Py%(yi yi)2ii7.612021/6/162021/6/163434预测预测 x = 0.04( / /L)时时, 菜粉蝶五龄幼虫死亡率菜粉蝶五龄幼虫死亡率?概率单位值概率单位值 P = 8.0746 + + 1.9191lg 0.04 = 5.392, 将概率单位值转换为百分率将概率单位值转换为百分率: 5.392 转换为百分率转换为百分率65.25% , (查附表查附表) 相关指数：相关指数：r = 0.98 g g 2 27.62 5. Logistic 生长曲线生长曲线( (S 形曲线形曲线) ) 动植

32、物的生长或繁殖动植物的生长或繁殖, 假定在无限空间、无假定在无限空间、无限营养来源等无约束条件下进行限营养来源等无约束条件下进行, 则其生长或繁则其生长或繁2021/6/162021/6/163535殖量与时间的关系开始时为指数函数殖量与时间的关系开始时为指数函数: 这种关系通常只能在一开始的一段时间内成立这种关系通常只能在一开始的一段时间内成立, 但实际上由于环境条件的恶化、营养的不足、衰但实际上由于环境条件的恶化、营养的不足、衰老的加快、繁殖的减慢、死亡的增多等约束条件老的加快、繁殖的减慢、死亡的增多等约束条件伴随而来伴随而来, 不能保持稳定不能保持稳定, 它们将改变这这种关它们将改变这这

33、种关系而呈现出系而呈现出 Logistic 曲线形式。曲线形式。7.632021/6/162021/6/163636 称上式为自然生长方程。其中称上式为自然生长方程。其中: y是生长量是生长量, 繁殖量繁殖量(或累积率或累积率), K 是时间无限时是时间无限时 y 的极限量的极限量(率率), x 为演变时间。该曲线与前两种为演变时间。该曲线与前两种 s 形曲线有所不同形曲线有所不同, 它没有正态分布或近似正态分布假定它没有正态分布或近似正态分布假定, 所以具有更所以具有更广泛的适应性。广泛的适应性。(1) Logistic 函数表达式函数表达式或或7.642021/6/162021/6/163

34、737 (2) Logistic 函数的性质函数的性质1. 当当x = 0 时可算得起始量时可算得起始量2. 当当 x =- ln(1/ /a)/ /b, y = K/ /2 时为拐点时为拐点, 即变即变化率从逐增过度到逐减的分界点。化率从逐增过度到逐减的分界点。1.3. 当当 x 在在0 , - - ln(1/ /a)/ /b 区区2.间时变化率逐增间时变化率逐增; 3.在在 - -ln(1/ /a)/ /b , 变化率变化率4.逐减。见图逐减。见图7.22xy 时间时间累累积积量量累累积积率率 % %Kk2x = - - ln(1/ /a)/ /b0图图7.22 Logistic曲线曲线.

35、7.65拐点拐点2021/6/162021/6/163838未知参数的确定：未知参数的确定：如果如果 y 是累积率则极限值是累积率则极限值 则则K=100。如果如果 y 是生长量或繁殖量是生长量或繁殖量, 其极限其极限 值值 K 可由按下式近似求得可由按下式近似求得在实际应用时在实际应用时 (x1, y1), (x2, y2), (x3, y3), 是数据中是数据中的三对观察值的三对观察值, x1, x2, x3 是等间距的量是等间距的量, 间距不间距不应太小应太小, 间距小了算出的间距小了算出的K值误差较大。值误差较大。7.662021/6/162021/6/163939 (3) Logis

36、tic 生长曲线直化过程生长曲线直化过程将函数将函数 变形得变形得等式两边同除以等式两边同除以 y 并取自然对数得并取自然对数得令：令：则有则有: 再求解直线回归系数再求解直线回归系数：a,- - b 7.672021/6/162021/6/164040 再对再对 a取反自然对数得曲线回归系数取反自然对数得曲线回归系数 a , 可得方程式可得方程式 。 注意注意: 如果对如果对式直化则有式直化则有 (对对回归系数回归系数a 不再取反对数不再取反对数)。即。即: , 式的式的 a 是不同的是不同的, ,这里必须指出用同一资料所配合的这里必须指出用同一资料所配合的 曲线方程曲线方程 和和 是等价的

37、。是等价的。注注: : Logistic 曲线方程是由比利时数学家曲线方程是由比利时数学家P.F.Verhulst t于于1838 年导出年导出, ,但长期湮没但长期湮没, ,至至20世纪世纪20年代才又被统计学家年代才又被统计学家R.Pearl 和和L.J.Reed重新发现。该方程在动植物的饲养、栽培、资源、重新发现。该方程在动植物的饲养、栽培、资源、 生态、环境保护等方面的模拟研究中已有较广泛的应用。生态、环境保护等方面的模拟研究中已有较广泛的应用。7.682021/6/162021/6/164141 例例7.13 试以生长曲线描述例试以生长曲线描述例7.10越冬代棉红越冬代棉红铃虫化蛹进

38、度和时间的关系。利用铃虫化蛹进度和时间的关系。利用 Logistic 函数函数的数学性质的数学性质, 求出拐点并解释其实际意义。求出拐点并解释其实际意义。 解：解：因为因为 y 变量是化蛹进度为百分数变量是化蛹进度为百分数(%), 所所以当发生时间以当发生时间 x 无限时无限时y 的极限值的极限值 K = =100 分析计算过程：分析计算过程：7.692021/6/162021/6/164242表表7.17 棉红铃虫化蛹进度和时间表棉红铃虫化蛹进度和时间表x(天天) y% 100 - -yy ln( )100 - -yy 5 3.5 27.571 3.316810 6.4 14.625 2.6

39、82715 14.6 5.8493 1.766320 31.4 2.1847 0.781525 45.6 1.1930 0.176530 60.4 0.6556 - -0.422235 75.2 0.3298 - -1.109340 90.2 0.1068 - -2.219645 95.4 0.0482 - -3.032050 97.5 0.0256 - -3.6636iiii7.702021/6/162021/6/164343反对数求反对数求 a 可得可得:得曲线回归方程得曲线回归方程:或或7.712021/6/162021/6/164444当当 x = 0 时化蛹进度时化蛹进度 y% =

40、=1.58%, 当当 x = =- -ln(1/ /62.08864)/ /0.156395 = 26.4为拐点为拐点, 见图见图7.23。说明在说明在6月月27日时化蛹日时化蛹进度达进度达 y% = 50%, 6月月27日前化蛹进度是日前化蛹进度是逐增的逐增的, 6月月27日以后日以后化蛹进度就逐渐减慢。化蛹进度就逐渐减慢。 图图7.23 棉红铃虫化蛹进度棉红铃虫化蛹进度xy%化化蛹蛹进进度度日期日期10 20 3040 5025507510050%100%6月月27日日. . 07.722021/6/162021/6/164545 例例例例7.147.14 测定某种肉用鸡在良好饲养条件下测

41、定某种肉用鸡在良好饲养条件下的生长过程的生长过程, 每两周测定一次每两周测定一次, 得数据如下得数据如下, 试试以以 Logistic方程描述之。方程描述之。(x 周次周次, y 重量重量kg)先求先求极限生长量极限生长量K值：由资料可选等间距值：由资料可选等间距 x 对应的对应的 y 值为值为 7.73解：解：由下式可求得由下式可求得K值值2021/6/162021/6/164646表表7.18 肉用鸡周次体重肉用鸡周次体重( (ykgkg) )数据表数据表2.827- -yy ln( )2.827- -yy 1 2 0.30 8.4233 2.1310 4.2620 2 4 0.86 2.

42、2872 0.8273 3.3092 3 6 1.73 0.6341 - -0.4556 - -2.7336 4 8 2.20 0.2850 - -1.2553 - -10.0424 5 10 2.47 0.1445 - -1.9345 - -19.345 6 12 2.67 0.0588 - -2.8336 - -34.0032 7 14 2.80 0.0096 - -4.6460 - -65.0224 56 - -8.1667 - -123.597平均平均 8 - -1.1667 周次周次x 重量重量 yiiiii7.742021/6/162021/6/164747配合回归方程：配合回归方

43、程：7.752021/6/162021/6/164848 当当: 周周 (为拐点为拐点 ) 是生长过程的速率由快转慢的界点是生长过程的速率由快转慢的界点, 当当 x5.76或或周时速率逐增周时速率逐增, x5.76时速率逐减。时速率逐减。x拐点拐点时时, 可用指数函数描述可用指数函数描述, x拐点时拐点时, 则需用则需用 Logistic 方程才能描述全过程。方程才能描述全过程。7.762021/6/162021/6/164949一元二次抛物线回归一元二次抛物线回归(一元二次多项式回归一元二次多项式回归) 7.2.2 7.2.2 不可直化的一元非线性回归不可直化的一元非线性回归不可直化的一元非

44、线性回归不可直化的一元非线性回归样本回归方程样本回归方程:1. 一元二次抛物线回归的配合一元二次抛物线回归的配合数学模型：数学模型：上式不能直接用取对数的方进行直化上式不能直接用取对数的方进行直化, ,只能用解只能用解方程的方程的方法求解回归系数方法求解回归系数 a, b, c 。7.772021/6/162021/6/165050图图7.24 抛物线回归方程与图形抛物线回归方程与图形y =a+bx+cx2yx 在研究动植物发育速度与温度间的关系、肥在研究动植物发育速度与温度间的关系、肥料或饲料与产量或重量的关系有时可用抛物线回料或饲料与产量或重量的关系有时可用抛物线回归来描述。如：苜蓿盲蝽卵

45、不同温度与卵的发育归来描述。如：苜蓿盲蝽卵不同温度与卵的发育速度间的关系就为抛物线。速度间的关系就为抛物线。2. 二次抛物线方程系数的最小二乘估计：二次抛物线方程系数的最小二乘估计： 令令: 求求函数函数 (因为因为Q是二次函数是二次函数, 所以一定有极值且有所以一定有极值且有极小极小7.782021/6/162021/6/165151值值), 对函数对函数Q(a,b,c)分别求一阶偏导分别求一阶偏导, 令其为令其为0：即即:7.792021/6/162021/6/165252三元一次正规方程组三元一次正规方程组对上面正规方程组求解对上面正规方程组求解 a, b, c (行列式解法行列式解法)

46、将方程两边各同除以将方程两边各同除以2, 进行整理之进行整理之, 可得下面可得下面7.802021/6/162021/6/165353得：得：回归方程回归方程:7.812021/6/162021/6/165454 例例例例7.157.15 对苜蓿盲蝽卵进行研究对苜蓿盲蝽卵进行研究, 发现不同温发现不同温度度x与卵的发育速度与卵的发育速度 y 呈二次曲线型呈二次曲线型, 见表见表7.19 图图7.25 试对其进行回归分析。试对其进行回归分析。 解：解：盲蝽象卵的发育速度盲蝽象卵的发育速度 y 是按照发育天是按照发育天数的倒数乘以数的倒数乘以100转换来的。转换来的。图图7.25 苜蓿盲蝽卵温度与

47、苜蓿盲蝽卵温度与发育速度关系图发育速度关系图 发育速度发育速度y温度温度 x.16481215 20 25 30 35 40表表7.19 苜蓿盲蝽卵温度与苜蓿盲蝽卵温度与 发育速度关系发育速度关系 温温 度度 x 发育速度发育速度 y 15 3.85 20 7.58 25 9.43 30 11.35 35 14.16 37 12.14 40 10.357.822021/6/162021/6/165555 7.832021/6/162021/6/165656分别求行列式分别求行列式 、 a、 b、 c7.842021/6/162021/6/165757回归方程：回归方程：分别求回归方程参数的估计

48、值分别求回归方程参数的估计值a、b、c7.85* * 注注: :本例计算量较大本例计算量较大, ,也可直接上计算机用程序计算。也可直接上计算机用程序计算。2021/6/162021/6/1658583. 抛物线回归的检验与预测抛物线回归的检验与预测表表7.20 方差分系表方差分系表 F0.01(1, 3) = 4.70差异来源差异来源 平方和平方和 自由度自由度 均方均方 均方比均方比 回回 归归 61.735 1 61.735 52.00* 剩剩 余余 5.933 5 1.187 总总 和和 67.668 6 检验：检验：7.862021/6/162021/6/165959 预测：预测：将将

49、 x 变量值代入回归方程变量值代入回归方程, 可求可求 y 变量的点预测值变量的点预测值, 例例7.15中当温度为中当温度为32 C度盲蝽度盲蝽象卵的发育速度为何？象卵的发育速度为何？ t(5, 0.05)=2.571 在对二次抛物线回归方程的研究中在对二次抛物线回归方程的研究中, 回归值回归值 y 的方差可近似等于剩余方差的方差可近似等于剩余方差 Se, 所以区间式为所以区间式为:207.872021/6/162021/6/1660607.88 例例例例7.17.16 6 现在配合例现在配合例7.7资料的二次抛物线方资料的二次抛物线方程并给于检验程并给于检验 , 再与幂函数进行比较。再与幂函

50、数进行比较。 x y x x x y xy x y 1 15 13.9 2 20 17.1 3 25 20.0 4 30 22.1 5 35 24.0 6 40 25.6 7 45 27.0 8 50 28.3 9 55 29.410 60 30.211 65 31.4表表7.21 云杉胸径与树高云杉胸径与树高配合二次抛物线方程数据计算表配合二次抛物线方程数据计算表42322225400625900122516002025250030253600422533758000156252700042875640009112512500016637521600027462550625160000390

51、625810000150062525600004100625625000091506251296000017850625193.21193.21292.41292.41400400488.41488.41576576655.36655.36729729800.89800.89864.36864.36912.04912.04985.96985.96208.5208.53423425005006636638408401024102412151215141514151617161718121812204120413127.53127.568406840125001250019890198902940

52、0294004096040960546755467570750707508893588935108720108720132665132665 440 269 20350 1034000 55783750 6897.64 11677.5 568462.52021/6/162021/6/166161将下面所得数值代入正则方程组并解之：将下面所得数值代入正则方程组并解之：7.892021/6/162021/6/166262通过上机运算解方程组得二次抛物线方程为通过上机运算解方程组得二次抛物线方程为3. 对对抛物线方程的检验抛物线方程的检验表表7.20 方差分系表方差分系表 F0.01(1, 3) = 4.70差异来源差异来源 平方和平方和 自由度自由度 均均 方方 均方比均方比 回回 归归 318.5963 1 159.3 1659.38* 剩剩 余余 0.77077 9 0.096 总总 和和 319.3627 10* * 注注: :该例用二次抛物线配合优于幂函数方程配合。该例用二次抛物线配合优于幂函数方程配合。7.902021/6/162021/6/166363 结束语结束语若有不当之处，请指正，谢谢！若有不当之处，请指正，谢谢！