体长体重生长方程课件.ppt

1、渔业资源评估与管理渔业资源评估与管理第第2章章 生长方程与生长参数生长方程与生长参数1ppt课件学习目的1、生长方程的概念2、生长参数的生物学意义3、掌握用年龄-体长数据估计VBGF参数 的线性和非线性方法4、学会用FiSAT 软件估计生长参数5、样品随机取样的重要性6、Variation in length-at-age2ppt课件大量数据分析表明，鱼类（包括软体类、甲壳类）的长度大量数据分析表明，鱼类（包括软体类、甲壳类）的长度与重量呈一定函数关系，一般为幂指数：与重量呈一定函数关系，一般为幂指数：第第1节节 长度与体重关系长度与体重关系长度长度 L 重量重量 W全长全长体长体长 体重（净

2、体重）体重（净体重）叉长叉长胴长胴长头胸甲长头胸甲长),0(,2NeaLWb3ppt课件Biological implication),0(,2NeaLWbLbaWlnlnlnLog-normal error distributionImplication?)(BXAYnormal error distribution-linear regression4ppt课件如何估算参数如何估算参数 a,b跟踪测定一定数量个体的长度和体重跟踪测定一定数量个体的长度和体重标志标志放流；放流；用所有样品的长度、体重测定数据回归分析；用所有样品的长度、体重测定数据回归分析；根据长度频率，用长度组的平均长度与平

3、均体根据长度频率，用长度组的平均长度与平均体重进行回归；重进行回归；用年龄组的平均长度与平均体重回归分析。用年龄组的平均长度与平均体重回归分析。注意事项：取样的代表性注意事项：取样的代表性 5ppt课件线性回归方法变量之间呈线性函数关系（linear function）参数估算参数估算(拟合拟合)方法：最小二乘法方法：最小二乘法(least squares residual)6ppt课件2YX134LS minimizes$iin2112223242 objective function模型拟合图解7ppt课件线性回归方法的假设（Assumptions）误差的概率分布是正态分布；误差分布的方差

4、是常数；误差分布的均值是 0；误差是独立的；X 变量无误差。8ppt课件如何确定参数如何确定参数 a,b非线性方法直接估计非线性方法直接估计将将W=a Lb两边取对数后，两边取对数后，采用线性回归求取采用线性回归求取a,b：W=a+L Ricker(1975)认为，认为，a,b 可可以用来表示鱼类是否为匀以用来表示鱼类是否为匀速生长，速生长，a=3为匀速生长。为匀速生长。实际上许多种类实际上许多种类 a3，因此，因此有时近似地：有时近似地：W=a L3，那，那么：么：组数样品总重量长度组的组中值长度组的尾数:13NwlnlnwaiiNiii9ppt课件估算a、b,如何取样？随机取样，样品反映种

5、群的长度组成；随机取样难以做到时，可作选择性取样，数量并非越多越好，适量即可，关键是个体大小范围要较宽广，大个体和幼鱼都要有一定的数量。10ppt课件example11ppt课件概念：为了对鱼类个体的生长进行定量研究，人们用数学模型（或称方程）来描述鱼类体长和体重随时间（或年龄）增长的变化规律，该模型或方程称为生长方程生长方程。根据生长方程在坐标图上绘出的体长和体重随时间的变化曲线称为生长曲线生长曲线。第第2节节 生长方程生长方程(growth function/formula/model)12ppt课件key point：Average growth of an individual or

6、individuals over life spanVary with size and among individuals13ppt课件Typical growth modelVon Bertalanffy 生长方程生长方程(VBGF)(commonly used)指数生长方程指数生长方程Logistic 生长方程生长方程Gompertz 生长方程生长方程14ppt课件指数生长方程指数生长方程Ricker(1975)认为：把生长周期分成很短的间距，任何生长曲线可以作为指数生长来对待。tGiiieWW1平均体重增长率平均体重增长率平均体重增长率平均体重增长率15ppt课件Logistic 生长

7、方程生长方程Logistic 生长方程常用于研究人口的增长，但对许多鱼类和甲壳类的生长也符合。rtateLL1极限体长16ppt课件Gompertz 生长方程生长方程Ricker(1975)认为Gompertz曲线也可以很好的描述有些鱼类的生长情况。rtegteLL17ppt课件Von Bertalanffy growth function(VBGF)Von Bertalanffy(1934,1938)提出（推导过程略）。符合大多数鱼类，甚至甲壳类、软体类，运用非常广泛。可以简单的运用到B-H 资源评估模型中。目前极大多数综合性渔业资源评估模型采用的都是VBGF生长参数。18ppt课件Von

8、Bertalanffy growth function(VBGF)Key points:Variables in the model;Parameters in the model;Biological implications of the parameters;Measurement unit for parameters;Parameter estimation:nonlinear regression analysis（非线性回归方法）19ppt课件体长、体重生长方程：体长、体重生长方程：01ttKteLLbttKteWW01Von Bertalanffy growth functio

9、n(VBGF)20ppt课件东海黄姑鱼年龄与生长的初步研究（吴常文等，2005）21ppt课件东海带鱼年龄与生长的研究（严利平等，2005）22ppt课件如何估算如何估算Von Bertalanffy 生长参数？生长参数？1)线性回归方法(Linear Regression)2)非线性回归方法非线性回归方法(Non-linear Regression)01ttKteLL根据所利用的观测资料，生长参数估算有两类：1）Age-length based estimation2）Length frequency based estimation23ppt课件tKKtLeeLL11btKKbbtWeeW

10、W/11/1111）线性回归方法）线性回归方法估算估算VBGF参数参数(old method)线性模型回归：线性模型回归：BKBALeBeLAKKln,11则，设（同理）（同理）24ppt课件t0的估算，每一年龄的t0：tLLLKtttln1025ppt课件Von Bertalanffy 生长参数估算（例）：生长参数估算（例）：习题二，习题二，4，Posgay(1953)将海水扇贝进行标志方流，将海水扇贝进行标志方流，约在约在10个月的期间内获得个月的期间内获得和重捕的生长资料如下表：和重捕的生长资料如下表：贝壳直径贝壳直径(mm)标志放流时标志放流时尺寸尺寸64697194104105110

11、117117126重捕时尺寸重捕时尺寸9810293115120126125127136138试估算：试估算：(1)生长参数生长参数L和和K；(2)取取t00，计算，计算I至至龄各龄的壳径长度；龄各龄的壳径长度；(3)绘制长度生长曲线。绘制长度生长曲线。Excel demo26ppt课件2）非线性回归方法，直接求）非线性回归方法，直接求VBGF参数参数(non-linear least square method)：01ttKteLL渔业上渔业上,许多模型是非线性的！许多模型是非线性的！27ppt课件非线性回归参数估计非线性回归参数估计It uses the same objective fu

12、nction as in a linear regression analysis.The model parameters can not be estimated directly as in a linear regression analysis,and an approach called numerical iteration has to be used.The parameter estimation is commonly done on a computer with a statistical software that can conduct nonlinear par

13、ameter estimation.28ppt课件Excel demo using solver（规划求解）（规划求解）29ppt课件Length frequency based VBGF parameters estimation（体长频率法）方法产生并发展于 1980主要创立者 D.Pauly:Pauly,D.and N.David,1981.ELEFAN I,a BASIC program for the objective extraction of growth parameters from length-frequencies data.Meeresforsch.,28(4):2

14、05-211.Pauly,D.,1984.Length-converted catch curves:a powerful tool for fisheries research in the tropics(Part II).ICLARM Fishbyte,2(1):17-19.Pauly,D.and G.R.Morgan(eds).,1987.Length-based methods in fisheries research.ICLARM Conf.Proc.(13):468p.30ppt课件关于体长频率法长度频率法方法不仅适用于热带、亚热带鱼类，也适用于温带鱼类，只要长度频率分布合理，

15、样品中幼鱼和高龄鱼都有适当分布，长度频率数据进行种群参数估计可弥补通过年龄鉴定拟合的生长曲线无法描述1龄以内生长情况的缺陷。该方法的运用有3个假定条件：样品能代表所研究的种群；von Bertalanffy生长方程（VBGF）能描述鱼种生长；样品中所有同龄鱼具有相同长度，因而不同长度的鱼则年龄不同。31ppt课件体长频率法的实现运用FiSAT 软件的ELEFAN程序。体长频率法估计的原理，可参考文献：何宝全，李辉权.珠江河口棘头梅童鱼的资源评估J.水产学报，1988，12(6)：125134.SPARRE P,URSIN E,VENEMA S C.Introduction to tropica

16、l fish stock assessment Part 1 Manual R.FAO Fisheries Technical Paper 306/1,Rome,1989,109115.32ppt课件体长频率法估计的步骤(ELEFAN I)体长频率法拟合生长曲线主要分4个步骤：1、进行长度(体长或肛长)频率样品重组，突出频率的峰和谷。2、计算所有样品中各峰最大值的总和，这个值称为峰的最大和(ASP)，它表示单一生长曲线穿过样品累计所经长度组的频率值可能达到的最大值。3、适配生长曲线。选用某一参数初值确定的生长曲线，不断改变曲线的起点，选取一系列曲线中，穿过最多的峰而避开最多的谷的一条，计算此时

17、峰的解释和(ESP)。4、改变参数的初值，重复步骤三直到比率ESP/ASP(称为score，取值范围01)达到最大且曲线反映的生长在生物学上可接受，则相应的参数作为最后的估算结果。以上过程由软件执行。以上过程由软件执行。33ppt课件时间时间（月份）（月份）肛长组肛长组南海北部南海北部带鱼的体带鱼的体长频率长频率（19972001）34ppt课件Do an example in ELEFAN!35ppt课件文献阅读胡芬 李建生 刘勇 等，东海带鱼年龄与生长的研究，海洋渔业 2005年 27卷 2期，139-142吴常文 赵淑江 胡春春，东海黄姑鱼年龄与生长的初步研究，海洋渔业 2005年 27卷 3期，193-199水柏年，黄海南部、东海北部小黄鱼的年龄与生长研究,浙江海洋学院学报：自然科学版,2003年 22卷 1期 16-2036ppt课件课外作业复习课堂上的例子习题二：2，3阅读12篇 参考文献(见精品课程网站)37ppt课件