S型曲线拟合PPT课件.ppt

第十章,曲线回归,1,本章介绍可以直线化的曲线回归的类型，以生长型曲线为例说明曲线的直线化配合，曲线回归方程的拟合度,2,第一节曲线回归的意义,3,直线回归的局限1、两变量之间的关系不完全是直线关系2、简单相关不显著并不表示两变量间无相关3、两变量间更普遍的关系是曲线关系4、直线回归仅是曲线回归的一种特殊形式5、直线回归是曲线回归中的一部分,4,曲线配合的一般步骤：1、确定回归关系的类型：线性非线性（曲线形状）2、确定回归关系的参数、相关指数、估计标准误3、对所得回归方程作显著性检验曲线方程可分为两种：可直线化的曲线方程不可直线化的曲线方程（多项式）因此，首先应确定两变量的曲线关系是哪一种,5,第二节曲线类型及其方程,6,本章仅讨论可以直线化的曲线方程函数型曲线方程（一）幂函数直线化：两边取对数：令：则有：对求A和b，并得即可得：a、b，建立方程（双对数转换，即对x、y均求对数后输入）,7,（二）指数函数或直线化：两边取对数：令：则有对求A并得即可得a、b，建立方程（单对数变换，即对y求对数后与x一起输入）,8,（三）双曲线函数令：则对x求X即可得中的a、b（倒数变换，即取x的倒数，与y一起输入）此外还有一些曲线方程：下面是几种可以转换为直线方程的曲线函数图形：,9,10,曲线回归的计算器计算方法：计算器将出现如下画面：,mode,3,LinLogExp123,2,3,12,（四）S型曲线陆生、水生动物的种群增长、微生物种群增长、细胞的生（增）长等都是这一模式因此，S型曲线又称为生长型曲线、logistic曲线，其变换形式有以下几种：,13,类似的生长型曲线还有Gompertz曲线：其变换形式：Bertalanffy曲线：,14,在这些曲线方程中，无一例外的都有3个需要计算的统计量：k、a、bK是当x趋向于+时y所能达到的最大值，往往是未知的，因此也是需要进行计算的这是生长曲线与其他可以直线化的曲线方程不同的地方这些曲线方程中的x往往是时间单位，因此一般可用t表示，而y往往是群体的增长量，或群体增长倍数，所以也可以用N表示我们这里仅对典型的S型曲线方程进行直线化，其他变换类型的方程直线化可以仿此进行,15,测得某微生物在一定温度下随时间变化的平均增长量数据如下：时间t123456789增长倍数N1.31.52.63.66.88.48.59.19.5从下面的散点图我们可以看出，可配合S型曲线：108642123456789,16,我们采用生长曲线的一般形式进行配合变换，两边取对数，得：并令：从数据表中取三个等距的点代入上式（一般总取始点、中点、末点）：（1，1.3）、（5，6.8）、（9，9.5）,17,解这一三元一次方程组，消去a、b，得：则这是一个通式，任何配置S型曲线的数据资料均可使用这一公式求得k值将上式中的代入式，得即为k的解将k=9.78代入可得和t相对应的各个Y值,18,将这些Y值写在数据表下方对应处，用最小二乘配置法配置直线时间t123456789增长倍数N1.31.52.63.66.88.48.59.19.51.881.711.020.54-0.82-1.81-1.89-2.59-3.52,19,得一级数据：或将时间t和Y值输入计算器直接进行计算,20,则将k、a、b代入方程，即得：或：,21,在这一类例子中，时间往往是有效单位时间，如一周、一月、一年、一个时间段等，如需换算成具体时间如天、小时、分等，则需将其换算值代入t值即可另外，在一般的通式中，我们往往以x、y作为自变量和依变量的符号，但在具体问题中，有时为了更形象、更直观地说明问题，可以用其他不同的字母（往往是相应的英文名词的首写字母）来代替,22,如长度用L、时间用t、增重倍数用N、体重用W等用统计软件进行计算时，可直接将原始数据输入数据库，调用相应的程序运算即可,23,第三节曲线配合的拟合度,24,曲线配合完成，其方程是否理想，同一批数据采用不同的曲线方程进行拟合，其效果如何，哪一种方程更好，可以用曲线方程的拟合度来衡量曲线方程的拟合度就是相关指数R2离回归平方和Q（实测值与预测值之差的平方和，即剩余回归平方和）在总平方和中所占的比例越小，说明方程的效果越好，因此可以用剩余回归平方和在总平方和中的比例来表示曲线配合的好坏：,25,在曲线回归方程中，我们必须实际求得每一个，然后求出，而不能象简单回归一样可以用有关公式求出在上例中：t123456789N1.31.52.63.66.88.48.59.19.50.94451.74813.00304.63866.33267.71678.64489.18749.47970.12640.06150.16241.07880.21840.46690.02100.00760.0004,26,R2的平方根R称为相关系数，为了和简单相关系数r有所区别，曲线回归方程和多元回归方程的相关系数称为复相关系数，写为R拟合度得到后，同样需要进行显著性检验，检验的方法还是查r表本例中，变量个数为m=2，自由度