《生物统计学课件--16曲线拟合(回归)讲解》由会员分享，可在线阅读，更多相关《生物统计学课件--16曲线拟合(回归)讲解（35页珍藏版）》请在金锄头文库上搜索。
1、第五节 曲线拟合（非线性回归分析） （可化为直线的曲线回归方程） 一元线性回归方程是一种最简单的回归，但在实际工作中，遇到 的却也很多。因为任何一种曲线回归，在一个很小的区间内，都可以认 为是直线回归。 尽管如此，直线回归并不能代替曲线回归。 生物学中呈曲线关系的例子很多： 细菌生长的数量与时间的关系 年龄与身高的关系 作物种植密度与作物产量的关系 辐射强度与或药物的浓度与致死率的关系 我们将数据做适当的变换，将曲线化为直线，再按直线回归分析的 方法处理，就可以求出相应的曲线回归方程，这一分析过程，就称为曲 线拟合。 为了将回归曲线直线化，首先要了解两个变量间所呈的函 数关系。 两个变量间的函
2、数关系，可以根据专业知识去判定。 例如： 生长曲线一般呈S形曲线关系。 细菌的生长量与时间的关系为指数函数关系等。 此外还可以由散点图来判定两个变量的关系。 一、对数函数曲线的拟合 1、对数方程的一般表达式： 2、对数曲线 的图象 3、 直线化方法： 若令 ，则有 4、求 a 和 b 的值： 例：在水稻塑料薄膜青苗床内空气最高温度和室外空气最高温 度资料列于下表，试求它们之间的相关关系。 序号12345678910 X7.27.911.816.912.018.718.920.221.822.7 Y13.821.424.933.632.339.540.136.940.242.6 X0.8573
3、0.897 6 1.0719 1.227 9 1.079 2 1.271 8 1.2765 1.305 4 1.3385 1.356 0 序号11121314151617181920 X22.923.123.323.623.827.027.628.630.731.4 Y44.636.635.144.444.143.948.348.546.350.4 X 1.359 8 1.363 6 1.367 4 1.372 9 1.376 5 1.431 4 1.440 9 1.456 4 1.487 1 1.496 9 X：室外气温 Y：苗床内气温 首先,根据上表的数据(x，y)绘制散点图。 由散点图看
4、x与y的关系符合曲线类型，很象b0时的对数函数曲线 y=a+blgx，所以做x=lgx的对数转换，并求x=lgx。以（x，y）作散 点图，散点近似于线性分布。 求一级和二级数据： 所以有 ： 回归关系的检验可以进行回归系数或相关系数的检验 。 苗床内最高气温y与空气最高气温x（lgx）之间的线性关系 二、指数曲线的拟合 1、一般化方程： 2、描述的现象： 细菌的繁殖 土壤中某些杀虫剂的分解曲线 放射性同位素的衰变 水果、蔬菜的腐败和温度的关系 3、直线化的方法： 方程的两边取常用对数，有 lg y = lga + xlgb 另lgy = y，lga = a，lgb = b，则有：y=a+bx
5、4、指数曲线 的图象 例：棉花红铃虫的产卵数与温度有关，试根据下表数据，建立棉花红铃 虫产卵数与温度的回归方程。 x 温度21232527293235 y 产产卵数711212466115325 y=lgy0.84511.04141.32221.38021.81922.06072.5119 由(x，y)的散点图可见，x 与 y 有曲线回归关系,散点分布近于指数函数曲 线的分布，将y=lgy与x做散点图(x，lgy)，x与lgy的分布近似于线性分布 。 X 与 lgy 的分布近似于线性分布，做线性回归： 解： 回归关系的检验：可以利用 b 或者 r 进行检验，主要是对线 性关系的检验，线性回归或
6、相关显著，则指数回归关系的拟 合就显著。 三、幂函数曲线的拟合 1、一般化方程： 2、幂函数 曲线的图象 3、直线化的方法： 两边取常用对数，使之线性化，得： 令：lgy = y，lgx = x，lga = a，则：y= a+bx 上式符合直线回归方程的一般形式。 4、怎样知道利用对数转换 最好的方法依然是绘制两个变量的散点图，如果散 点不符合直线或狭长形的椭圆，则可在双边对数纸上， 再将散点绘制出来，若散点的分布符合线性或狭长椭圆 ，则可以进行x与y的双对数转换。 例题：在突变实验中，用不同剂量的射线，照射植物的种子，结果苗期高 度与成活株之间有一定的关系。用X射线照射大麦的种子，记处理株第
7、一 叶平均高度占对照株高度的百分数为X，存活百分数为Y，得到以下结果 ，试对两者的关系进行分析。 X283240506072808085 Y81218283055618580 1、画两个变量的散点图： 由不同数据变换的散点图可见，X与Y的关系是幂函数的关系 。 按照幂函数的直线化方法，求二级数据： Y 对 X 的回归方程为： 因为：a = lga，所以：a = 10-1.9577 = 0.0110 所以： 五、曲线的检验 有时将同一组数据，我们将其做指数函数或幂函数形式的变 换，都能得到X与Y的拟合曲线，并且可能在做线性回归关 系检验的时候，线性关系都显著。 那么，究竟哪一条拟合曲线是最好的呢
8、？ 一般情况下，以剩余平方和或称之为误差平方和的大小来判 断，即SSe最小时的拟合曲线为最好的曲线。 可以对一组数据尝试进行多种拟合，找出SSe最小，其为最好 的拟合曲线。 必须用 来计算， 不可以用 来计算。 四、概率对数变换 在辐射育种和药理学实验中，经常遇到寻找半致死剂量的 问题。 致死率（y）与剂量（x）间的关系曲线往往呈S形。在半致 死处，曲线的斜率最大，与死亡率的交点最清楚，因此在 实际工作中，常常采用半致死剂量这一标准。 确定半致死剂量，最常用的方法就是概率对数转换，当我们 把死亡率做概率转换，辐射剂量做对数转换后，两者会呈很 好的线性关系。 1、概率转换： P（u up）= p
9、， 其中， p 为死亡率，即为 y 值，up为y的值。 2、对数转换： x= lgx，x 为辐射剂量。 则：y=a+bx，y与x间有很好的线性关系。 3、半致死剂量： 所谓的半致死剂量，即为Y=50%时的辐射或药物剂量X，用 LD50来表示，当p = y = 0.5时，P（uup）= 0.5，则 up= 0。 将up= y= 0 代入 y= a + bx， 则有 ：0 = a + bx， 则有：x= -a/b， 因为 x = lgx，所以 此时的x即为半致死剂量，即LD50。 例题：用不同剂量的 射线照射小麦品种库斑克，调查死苗 率，得到以下结果： 剂剂量（Kr）x 14 16 18 20 2
10、2 24 26 死苗率（%）y 6 10 40 70 80 93 95 x=lgx1.146 1.204 1.255 1.301 1.342 1.380 1.415 y=up-1.555 -0.842 -0.253 0.524 0.842 1.476 1.645 利用原始数据绘制散点图，x与y为S形曲线形状，将y与x 做概率转换和对数转换，则两者的散点图呈线性关系。 因此，对转换后的数据（褐色）进行线性分析： 求二级数据： 五、曲线的检验 有时将同一组数据，我们将其做指数函数或幂函数形式的变 换，都能得到X与Y的拟合曲线，并且可能在做线性回归关 系检验的时候，线性关系都显著，那么，究竟哪一条拟合曲 线是最好的呢？ 一般情况下，以剩余平方和或称之为误差平方和的大小来判 断，即SSe最小时的拟合曲线为最好的曲线。 可以对一组数据尝试进行多种拟合，找出SSe最小，其为最好的拟合曲线 。 必须用 来计算， 不可以用 来计算。
