1、第八章第八章 参数估计方法参数估计方法第一节第一节 农业科学中的主要参数及其估计量的评选标准农业科学中的主要参数及其估计量的评选标准第二节第二节 矩法矩法第三节第三节 最小二乘法最小二乘法第四节第四节 极大似然法极大似然法第一节第一节 农业科学中的主要参数及其估计量的评选标准农业科学中的主要参数及其估计量的评选标准一、农业科学中的主要参数一、农业科学中的主要参数 (1)总体数量特征值参数，例如，用平均数来估计品种的产总体数量特征值参数，例如，用平均数来估计品种的产量，用平均数差数来估计施肥等处理的效应；量，用平均数差数来估计施肥等处理的效应；(2)在揭示变数间的相互关系方面，用相关系数来描述在
2、揭示变数间的相互关系方面，用相关系数来描述2个个变数间的线性关系；用回归系数、偏回归系数等来描述原因变数间的线性关系；用回归系数、偏回归系数等来描述原因变数变化所引起的结果变数的平均变化的数量，用通径系数变数变化所引起的结果变数的平均变化的数量，用通径系数来描述成分性状对目标性状的贡献程度等。来描述成分性状对目标性状的贡献程度等。农业科学研究中需要估计的参数是多种多样的，主要包括农业科学研究中需要估计的参数是多种多样的，主要包括:二、参数估计量的评选标准二、参数估计量的评选标准(一一)数学期望数学期望 样本平均数的平均数就是一种数学期望。样本平均数的平均数就是一种数学期望。例如，一个大豆品种的
3、含油量为例如，一个大豆品种的含油量为20%，测定一次可能，测定一次可能是大于是大于20%，再测定可能小于，再测定可能小于20%，大量反复测定后平均，大量反复测定后平均结果为结果为20%，这时，这时20%便可看作为该大豆品种含油量的数便可看作为该大豆品种含油量的数学期望，而每单独测定一次所获的值只是学期望，而每单独测定一次所获的值只是1个随机变量。个随机变量。抽象地，随机变量的数字特征是指随机变量的数学期抽象地，随机变量的数字特征是指随机变量的数学期望值。望值。对于离散型对于离散型(间断性间断性)随机变量随机变量y的分布列为：的分布列为：Py=yi=pi，其中，其中，i=1，2，那么随机变量，那
4、么随机变量y的数学期望的数学期望E(y)为：为：1iiipyyE)(81)这样可以求得总体平均值。这样可以求得总体平均值。对于连续型随机变数对于连续型随机变数y的数学期望的数学期望E(y)为：为：dyyyfyE)()(82)其中其中f(y)为随机变量为随机变量y的概率密度函数，这样可以求得总的概率密度函数，这样可以求得总体均值。体均值。用用D(y)表示方差，有表示方差，有 D(y)=E yE(y)2(83)这就是随机变量函数的数学期望。同理，离散型随机这就是随机变量函数的数学期望。同理，离散型随机变量方差的数学期望为：变量方差的数学期望为：12iiipyEyyD)()(84)连续型随机变量方差
5、的数学期望为：连续型随机变量方差的数学期望为：dyyfyEyyD)()()(2(85)数学期望有这样一些常用的性质：数学期望有这样一些常用的性质：(1)常数的数学期望为常数本身；常数的数学期望为常数本身；(2)随机变量与常数的乘积的数学期望是常数与随机随机变量与常数的乘积的数学期望是常数与随机变量的数学期望的乘积；变量的数学期望的乘积；(3)多个随机变量分别与常数的乘积的求和函数的数多个随机变量分别与常数的乘积的求和函数的数学期望是常数与多个随机变量的数学期望的乘积的和；学期望是常数与多个随机变量的数学期望的乘积的和；(4)多个相互独立的随机变量的乘积的数学期望是多多个相互独立的随机变量的乘积
6、的数学期望是多个随机变量的数学期望的乘积。个随机变量的数学期望的乘积。(二二)参数估计量的评选标准参数估计量的评选标准评价估计量优劣的标准主要有无偏性、有效性、相合性等评价估计量优劣的标准主要有无偏性、有效性、相合性等 (1)无偏性无偏性 参数估计量的期望值与参数真值是相等的，这参数估计量的期望值与参数真值是相等的，这种性质称为种性质称为无偏性无偏性，具有无偏性的估计量称为，具有无偏性的估计量称为无偏估计量无偏估计量。例如，在抽样分布中已经介绍了离均差平方和除以自由度得例如，在抽样分布中已经介绍了离均差平方和除以自由度得到的均方的平均数等于总体方差，即该均方的数学期望等于到的均方的平均数等于总
7、体方差，即该均方的数学期望等于相应总体参数方差，这就是说该均方估计量是无偏的。相应总体参数方差，这就是说该均方估计量是无偏的。估计量的数学期望值在样本容量趋近于无穷大时与参数估计量的数学期望值在样本容量趋近于无穷大时与参数的真值相等的性质称为的真值相等的性质称为渐进无偏性渐进无偏性，具有渐进无偏性的估计，具有渐进无偏性的估计量称为量称为渐进无偏估计量渐进无偏估计量。(2)有效性有效性 无偏性表示估计值是在真值周围波动的一个数无偏性表示估计值是在真值周围波动的一个数值，即无偏性表示估计值与真值间平均差异为值，即无偏性表示估计值与真值间平均差异为0，近似可以，近似可以用估计值作为真值的一个代表。用
8、估计值作为真值的一个代表。同一个参数可以有许多无偏估计量，但不同估计量的期望同一个参数可以有许多无偏估计量，但不同估计量的期望方差不同，也就是估计量在真值周围的波动大小不同。估计方差不同，也就是估计量在真值周围的波动大小不同。估计量的期望方差越大说明用其估计值代表相应真值的有效性越量的期望方差越大说明用其估计值代表相应真值的有效性越差；否则越好，越有效。不同的估计量具有不同的方差，方差；否则越好，越有效。不同的估计量具有不同的方差，方差最小说明最有效。差最小说明最有效。如果一个无偏估计量相对与其它所有可能无偏估计量，其如果一个无偏估计量相对与其它所有可能无偏估计量，其期望方差最小，那么称这种估
9、计量为期望方差最小，那么称这种估计量为一致最小方差无偏估计一致最小方差无偏估计量量。(3)相合性相合性 用估计量估计参数涉及一个样本容量大小问题，用估计量估计参数涉及一个样本容量大小问题，如果样本容量越大估计值越接近真值，那么这种估计量是如果样本容量越大估计值越接近真值，那么这种估计量是相相合估计量合估计量。除以上三方面标准外，还有除以上三方面标准外，还有充分性充分性与与完备性完备性也是常考虑也是常考虑的。的。充分性充分性指估计量应充分利用样本中每一变量的信息；指估计量应充分利用样本中每一变量的信息；完备性完备性指该估计量是充分的唯一的无偏估计量。指该估计量是充分的唯一的无偏估计量。第二节第二
10、节 矩法矩法一、矩的概念一、矩的概念 矩矩(moment)(moment)分为分为原点矩原点矩和和中心矩中心矩两种。两种。对于样本对于样本y1,y2,yn，各观测值的，各观测值的k次方的平均值，称为次方的平均值，称为样本的样本的k k阶原点矩阶原点矩，记为，记为 ，有，有 ,用观测值减去用观测值减去平均数得到的离均差的平均数得到的离均差的k次方的平均数称为次方的平均数称为样本的样本的k k阶中心矩阶中心矩,记为记为 或或 ，有，有 。kynikikyny11kyy)(knikikyynyy1)(1)(对于总体对于总体y1,y2,yN，各观测值的，各观测值的k次方的平均值，称为次方的平均值，称为
11、总体的总体的k k阶原点矩阶原点矩，记为，记为 ，有，有 ；用观测；用观测值减去平均数得到的离均差的值减去平均数得到的离均差的k次方的平均数称为次方的平均数称为总体的总体的k k阶阶中心矩中心矩，记为，记为 或或 ，有，有 )(kyENikikyNyE11)()(kyEkNikikyNyE1)(1)(二、矩法及矩估计量二、矩法及矩估计量 所谓所谓矩法矩法就是利用样本各阶原点矩来估计总体相应各阶就是利用样本各阶原点矩来估计总体相应各阶原点矩的方法，即原点矩的方法，即)(11knikikyEyny(86)也可以用样本各阶原点矩的函数来估计总体各阶原点矩也可以用样本各阶原点矩的函数来估计总体各阶原点
12、矩同一函数，即若同一函数，即若Q=f(E(y),E(y2),E(yk),则则)，(kyyyfQ2由此得到的估计量称为由此得到的估计量称为矩估计量矩估计量。例例8.1 现获得正态分布现获得正态分布 的随机样本的随机样本y1,y2,yn，要求正态分布要求正态分布 参数参数 和和 的矩估计量。的矩估计量。)，(2N)，(2N2 首先，求正态分布总体的首先，求正态分布总体的1阶原点矩和阶原点矩和2阶中心矩：阶中心矩：dyyydyyyfE(y)222)(exp21)(2222222)(exp21)()()()(dyyydyyyyE然后求样本的然后求样本的1阶原点矩和阶原点矩和2阶中心矩，为阶中心矩，为n
13、iiniiyynsyny12221)(1，1最后，利用矩法，获得总体平均数和方差的矩估计最后，利用矩法，获得总体平均数和方差的矩估计niiniiyynsyny12221)(1，1 故总体平均数和方差的矩估计值分别为样本平均数和样故总体平均数和方差的矩估计值分别为样本平均数和样本方差，方差的分母为本方差，方差的分母为n。单峰分布曲线还有二个特征数，即单峰分布曲线还有二个特征数，即偏度偏度(skewness)(skewness)与与峰峰度度(kurtosis)(kurtosis)，可分别用偏度系数和峰度系数作测度。，可分别用偏度系数和峰度系数作测度。偏度系数偏度系数(coefficient of
14、skewness)(coefficient of skewness)是指是指3阶中心矩与标准阶中心矩与标准差的差的3次方之比；次方之比；峰度系数峰度系数(coefficient of kurtosis)(coefficient of kurtosis)是指是指4阶中阶中心矩与标准差的心矩与标准差的4次方之比。次方之比。当偏度为正值时，分布向大于平均数方向偏斜；偏度当偏度为正值时，分布向大于平均数方向偏斜；偏度为负值时则向小于平均数方向偏斜；当偏度的绝对值大于为负值时则向小于平均数方向偏斜；当偏度的绝对值大于2时，分布的偏斜程度严重。当峰度大于时，分布的偏斜程度严重。当峰度大于3时，分布比较陡峭
15、，时，分布比较陡峭，峰态明显，即总体变数的分布比较集中。峰态明显，即总体变数的分布比较集中。由样本计算的偏度系数由样本计算的偏度系数 23121333)(1)(1niiniiyynyyncs(87)峰度系数峰度系数 24121444)(1)(1niiniiyynyynck(88)例例8.2 计算表计算表3.4数据资料数据资料(140行水稻产量行水稻产量)所属分布曲线所属分布曲线的偏度和峰度。的偏度和峰度。表3.4 140行水稻产量(单位：克)177215197 97123159245119119131149152167104161214125175219118192176175 9513619
16、9116165214 9515883137 80138151187126196134206137 98 97129143179174159165136108101141148168163176102194145173 75130149150161155111158131189 91142140154152163123205149155131209183 97119181149187131215111186118150155197116254239160172179151198124179135184168169173181188211197175122151171166175143190213
17、192231163159158159177147194227141169124159 首先，计算样本的首先，计算样本的2、3、4阶中心矩阶中心矩 ，以及标，以及标准差估计值：准差估计值：432，7351303)(1122.yynnii8913953)(1133.yynnii614410677294)(1.yynnii10736)(1122.yynni 然后，根据矩法原理，该分布的偏度与峰度估计值分别然后，根据矩法原理，该分布的偏度与峰度估计值分别为：为：0849033./cs752243./ck 因此，说明资料比较集中在平均数左右，分布曲线并不因此，说明资料比较集中在平均数左右，分布曲线并不是
18、特别陡峭。是特别陡峭。例例8.3 例例6.9为研究籼粳稻杂交为研究籼粳稻杂交F5代系间单株干草重的代系间单株干草重的遗传变异，随机抽取遗传变异，随机抽取76个系进行试验，每系随机取个系进行试验，每系随机取2个样品个样品测定干草重测定干草重(g/株株)。按单向分组方差分析进行分析，结果见。按单向分组方差分析进行分析，结果见表表6.9。此处用来说明由矩法估计误差、遗传方差和干草的。此处用来说明由矩法估计误差、遗传方差和干草的遗传力遗传力h2。因为因为76个系是随机抽取的，因而为随机模型。方差结个系是随机抽取的，因而为随机模型。方差结果说明系间差异显著，因而系间效应存在。果说明系间差异显著，因而系间
19、效应存在。根据矩法，首先应求出系间和误差变异来源的样本均根据矩法，首先应求出系间和误差变异来源的样本均方和总体期望均方方和总体期望均方(表表6.9)。然后，利用矩估计原理，令样本的均方与总体相应变然后，利用矩估计原理，令样本的均方与总体相应变异的期望均方相等，从而求出异的期望均方相等，从而求出 和和 的矩估计值。的矩估计值。22此处此处E(MS系统间系统间)=ETt-E(Tt)2，(Tt 为各个系统的总和数为各个系统的总和数)=22n E(MS误差误差)=E(e2)=，(e为误差为误差)2因而因而 77172.7972222.51.272)77.1779.72(2%76.6077.1751.2
20、751.27222222222eggpgh第三节第三节 最小二乘法最小二乘法 从总体中抽出的样本观察值与总体平均数是有差异的，从总体中抽出的样本观察值与总体平均数是有差异的，这种差异属于这种差异属于抽样误差抽样误差。因而，在总体平均数估计时要尽。因而，在总体平均数估计时要尽可能地降低这种误差，使总体平均数估计值尽可能好。可能地降低这种误差，使总体平均数估计值尽可能好。参数估计的最小二乘法就是基于这种考虑提出的。参数估计的最小二乘法就是基于这种考虑提出的。基本思想基本思想是使误差平方和最小，达到在误差之间建立一是使误差平方和最小，达到在误差之间建立一种平衡，以防止某一极端误差对决定参数的估计值起
21、支配种平衡，以防止某一极端误差对决定参数的估计值起支配地位。这有助于揭示更接近真实的状况。地位。这有助于揭示更接近真实的状况。具体方法具体方法是为使误差平方和是为使误差平方和Q为最小，可通过求为最小，可通过求Q对待对待估参数的偏导数，并令其等于估参数的偏导数，并令其等于0，以求得参数估计量。，以求得参数估计量。例例8.4 用最小二乘法求总体平均数用最小二乘法求总体平均数 的估计量。的估计量。若从平均数为的总体中抽得样本为若从平均数为的总体中抽得样本为y1、y2、y3、yn，则观察值可剖分为总体平均数与误差，则观察值可剖分为总体平均数与误差 ei 之和，之和，iiey 总体平均数的最小二乘估计量
22、就是使总体平均数的最小二乘估计量就是使 yi 与间的误差平方与间的误差平方和为最小，即和为最小，即niiiyeQ122)(为最小。为最小。为获得其最小值，求为获得其最小值，求Q对的导数，并令导数等于对的导数，并令导数等于0，可得：，可得：0)(2niiyQ1即总体平均数的估计量为：即总体平均数的估计量为：niiyn11 因此，算术平均数为总体平均数的最小二乘估计。这与因此，算术平均数为总体平均数的最小二乘估计。这与矩法估计是一致的。矩法估计是一致的。估计离均差平方和估计离均差平方和 的数学期望：的数学期望：2)(yyQi222222222)1()()()()(2)()()()(n/nnnyyE
23、yyyyEyyEyyEQEiiiii因而，因而，估计为：估计为：21)(1)(nyynQi22)(与矩法所得不同，而与常规以自由度为除数法一致。与矩法所得不同，而与常规以自由度为除数法一致。例例8.5 求例求例6.13的两向分组方差分析资料缺的两向分组方差分析资料缺1个小区个小区(表表8.1)的最小二乘估计量和估计值。的最小二乘估计量和估计值。iy表8.1 生长素处理豌豆的试验结果处 理(A)组(B)总和Ti 平均 对照(CK)6062616024360.8 赤霉素656568ye198+ye 动力精6361616024561.3 吲哚乙酸6467636125563.8 硫酸腺嘌吟626562
24、6425363.3 马来酸6162626525062.5总和 T j375382377310+yeT=1444+ye 从第从第6章可知，这种资料模式的线性模型为：章可知，这种资料模式的线性模型为：按照最小二乘法的估计原理，使按照最小二乘法的估计原理，使 ijjiijyaii10rjj10 该模型的约束条件为：该模型的约束条件为：，和误差项服和误差项服从正态分布。从正态分布。aijiijrjijyQ1212)(为最小时可以求出效应和缺失小区为最小时可以求出效应和缺失小区ye的估计量，即的估计量，即 airjjiijeaijiijjrjjiijiairjjiijyyQyQyQyQ1111110)2
25、(01)(2(01)(2(01)(2(从而，最小二乘估计量分别为：从而，最小二乘估计量分别为：airjijaiijrjijjieairjijaiijjjairjijrjijiiairjijyaryayryyaryayyyaryryyyary11111111111111111111因而表因而表8.1中，缺失小区的估计值可由下式求出：中，缺失小区的估计值可由下式求出：64144463104198eeeeyyyy解上述方程，最小二乘估计值为：解上述方程，最小二乘估计值为：ye=65.6。缺区估计是根据线性模型，以及最小二乘法的原理得缺区估计是根据线性模型，以及最小二乘法的原理得到的。不过，试验中尽可
26、能不要缺区，因为缺区估计尽管到的。不过，试验中尽可能不要缺区，因为缺区估计尽管可以估计缺区的值，但是误差的自由度将减少，本试验的可以估计缺区的值，但是误差的自由度将减少，本试验的误差自由度将减少误差自由度将减少1。一般地，若一般地，若m个自变数个自变数x1、x2、x3、xm与依变数与依变数y存在存在统计模型关系统计模型关系)，；，(kmxxxfy2121(89)k，21其中，其中，为待估参数。为待估参数。通过通过n次观测次观测(nk)得到得到n组含有组含有x1i,x2i,xmi,yi(i=1,2,n)的数据以估计的数据以估计 。其最小二乘估计。其最小二乘估计值为使值为使k，2122121112
27、)，；，(kmiiiininixxxfyQ(810)为最小的为最小的 。这种估计方法称为参数估计的。这种估计方法称为参数估计的最最小二乘法小二乘法(least squares)(least squares)，或，或最小平方法最小平方法。k21，第四节第四节 极大似然法极大似然法 所谓所谓极大似然法极大似然法(maximum likelihood method)(maximum likelihood method)是值选择是值选择使事件发生概率最大的可能情况的参数估计方法。使事件发生概率最大的可能情况的参数估计方法。极大似然法包括二个步骤：极大似然法包括二个步骤：(1)建立包括有该参数估计量的建
28、立包括有该参数估计量的似然函数似然函数(likelihood(likelihood function)function)(2)根据实验数据求出似然函数达极值时的参数估计根据实验数据求出似然函数达极值时的参数估计量或估计值。量或估计值。一、似然函数一、似然函数 对于离散型随机变量，似然函数是多个独立事件的对于离散型随机变量，似然函数是多个独立事件的概率函数的乘积，该乘积是概率函数值，它是关于总体概率函数的乘积，该乘积是概率函数值，它是关于总体参数的函数。参数的函数。例如，一只大口袋里有红、白、黑例如，一只大口袋里有红、白、黑3种球，采用复置抽种球，采用复置抽样样50次，得到红、白、黑次，得到红、
29、白、黑3种球的个数分别为种球的个数分别为12，24，14，那么根据多项式的理论，可以建立似然函数为：那么根据多项式的理论，可以建立似然函数为：143242121)()()(!14!24!12!50ppp 其中其中p1，p2，p3分别为口袋中红、白、黑分别为口袋中红、白、黑3种球的概种球的概率率(p3=1p1p2)，它们是需要估计的。，它们是需要估计的。对于连续型随机变量，似然函数是每个独立随机观测对于连续型随机变量，似然函数是每个独立随机观测值的概率密度函数的乘积，则似然函数为：值的概率密度函数的乘积，则似然函数为：)；()；()；()；，()(nnyfyfyfyyyLL2121(811)若若
30、yi 服从正态分布服从正态分布 ，则，则 ，上式可，上式可变为：变为：)，(2N)(,)()(212)(2)(221222221222nnyynyyeeeL)1(11)，(812)二、极大似然估计二、极大似然估计 所谓所谓极大似然估计极大似然估计就是指使似然函数为最大以获得总体就是指使似然函数为最大以获得总体参数估计的方法。其中，所获得的估计总体参数的表达式参数估计的方法。其中，所获得的估计总体参数的表达式称为称为极大似然估计量极大似然估计量，由该估计量获得的总体参数的估计，由该估计量获得的总体参数的估计值称为总体参数的值称为总体参数的极大似然估计值极大似然估计值。为了计算上的方便，一般将似然
31、函数取对数，称为为了计算上的方便，一般将似然函数取对数，称为对数对数似然函数似然函数，因为取对数后似然函数由乘积变为加式，其表，因为取对数后似然函数由乘积变为加式，其表达式为：达式为：niinyfyyyLL121,lnlnln)()；，()(813)求极大似然估计量可以通过令对数似然函数对总体参求极大似然估计量可以通过令对数似然函数对总体参数的偏导数等于数的偏导数等于0来获得，即当来获得，即当 ，有，有)，(l21)，；，(lnkyyyL2121ln0)，；(linikyf211（k=1，2，l）(814)由此获得总体参数的极大似然估计量。由此获得总体参数的极大似然估计量。例例8.6 设设y1
32、,y2,yn是正态总体是正态总体 的随机的随机样本，求正态分布样本，求正态分布 参数的极大似然估计量。参数的极大似然估计量。)，(2N)，(2N似然函数为：似然函数为：niinniiyyL12222122221exp212exp21)()，(取对数，得：取对数，得：niiynnL1222221ln22ln2ln)()()，(那么似然方程组为：那么似然方程组为：niniynLyL1242221220212ln01ln)()，()()，(解得：解得：niiniiynyyn122111)(因此，正态分布总体平均数的极大似然估计量为：因此，正态分布总体平均数的极大似然估计量为：当总体平均值为未知时，方
33、差估计量为当总体平均值为未知时，方差估计量为:niiyyn11当总体平均值为已知时，方差估计量为当总体平均值为已知时，方差估计量为:niiyyn1221)(niiyn1221)(例例8.7 求红、白、黑球事例中求红、白、黑球事例中p1，p2，p3的极大似然估计值。的极大似然估计值。由由 可获得对数似然函数可获得对数似然函数 14324212114241250)()()(！ppp)()(21213213211ln14ln24ln12 ln14ln24ln12,lnppppCpppCpppL其中，其中，C为常数。为常数。分别求分别求 对对p1，p2的偏导数，并的偏导数，并令为令为0，得似然方程组：
34、，得似然方程组：)1，(2121lnppppL0111424)1 ln01114121 ln2122121221121211)(，()()，(pppppppLppppppppLp联立求解，得：联立求解，得：257,2512256321/p/p/p,显然，极大似然估计值显然，极大似然估计值 等于其观测频率。等于其观测频率。321ppp，例例8.8 两个亲本的基因型分别为两个亲本的基因型分别为AABB和和aabb，这两，这两个亲本杂交后个亲本杂交后F2出现了出现了4种基因型，分别为种基因型，分别为A_B_、A_bb、aaB_和和aabb，得到四种基因型的个数分别为，得到四种基因型的个数分别为c、d
35、、e、f，已知已知AA和和BB两对基因间存在连锁关系，现欲估计重组率？两对基因间存在连锁关系，现欲估计重组率？设重组率为设重组率为r，根据遗传学推导，可以得到，根据遗传学推导，可以得到4种基因型种基因型的概率见表的概率见表8.2。表表8.2 F F2 2群体基因型的分离情况群体基因型的分离情况基基 因因 型型A_B_A_bbaaB_aabb总数总数观察得观察得到基因到基因型个数型个数c(289)d(26)e(29)f(76)n(420)概概 率率14122)(r4112)(r4112)(r412)(r 首先，通过表首先，通过表8.3介绍由两对连锁主基因控制的介绍由两对连锁主基因控制的F2群体群
36、体16种基因型的概率计算出种基因型的概率计算出4种表现型的概率种表现型的概率(表表8.2)。配子配子及概率及概率AB (1r)/2Abr/2aBr/2ab(1r)/2AB (1r)/2AABB(1r)2/4AABbr(1r)/4AaBBr(1r)/4AaBa(1r)2/4Abr/2AABbr(1r)/4AAbbr2/4AaBbr2/4Aabbr(1r)/4aBr/2AaBBr(1r)/4AaBbr2/4aaBBr2/4aaBbr(1r)/4ab(1r)/2AaBa(1r)2/4Aabbr(1r)/4aaBbr(1r)/4Aabb(1r)2/4表表8.3 F2群体的基因型及其概率群体的基因型及其
37、概率 按多项式分布，可以根据概率函数得到似然函数为：按多项式分布，可以根据概率函数得到似然函数为：fedcrrrrfedcnrL414114114122222!)(815)若以若以 代入上式，则似然函数和对数似然函数分别为：代入上式，则似然函数和对数似然函数分别为：21)(rfedcfedcnrL 4414142!)(4ln41ln)(42ln)(lnfedckL(k是常数项是常数项)(816)(817)对上式求导数，并令导数为对上式求导数，并令导数为0，可得方程：，可得方程：012fedc上式化解为一元二次方程上式化解为一元二次方程02222ffedcn)(nfnfedcfedc282222
38、2)()(818)在在 的两个解中取一个符合遗传规律的解，那么，的两个解中取一个符合遗传规律的解，那么，重组率的解为：重组率的解为：。1r重组率方差估计量为：重组率方差估计量为：)(1)()(2221nrD(819)对于本例，有对于本例，有6140.01226.04202420768)76292262289()76292262289(2取正根，取正根，=0.7366，736601.r0190)7366021(4202)736602)(736601()()(.rDrs由此，由此，=0.142。三、关于三种估计方法的讨论三、关于三种估计方法的讨论 上述上述3种参数估计方法对比：种参数估计方法对比：
39、(1)对于总体平均数的估计量，对于总体平均数的估计量，3种估计方法都具有无偏性、种估计方法都具有无偏性、有效性和相合性；有效性和相合性；(2)对于总体方差的估计量，由离均差平方和期望值所得的对于总体方差的估计量，由离均差平方和期望值所得的是无偏的，但由矩法和极大似然法所得两种估计量是有偏是无偏的，但由矩法和极大似然法所得两种估计量是有偏的，但都是相合的；最小二乘法无直接的总体方差估计量。的，但都是相合的；最小二乘法无直接的总体方差估计量。3种常用方法的不同要求种常用方法的不同要求:(1)极大似然法要求已知总体的分布，才能获得估计量极大似然法要求已知总体的分布，才能获得估计量。(2)另外两种方法
40、对分布没有严格的要求。另外两种方法对分布没有严格的要求。3种常用方法的应用范围：种常用方法的应用范围：(1)极大似然法估计结果大多具有无偏性、有效性和相合极大似然法估计结果大多具有无偏性、有效性和相合性等优良的估计量性质性等优良的估计量性质。(2)最小二乘法在估计线性回归模型参数时具有灵活方便最小二乘法在估计线性回归模型参数时具有灵活方便的特点。的特点。(3)矩估计方法由于不需要知道总体分布也是经常采用的矩估计方法由于不需要知道总体分布也是经常采用的方法，但该方法估计结果有时不具备优良的估计量性质，而方法，但该方法估计结果有时不具备优良的估计量性质，而且局限在与矩有关的估计量。且局限在与矩有关的估计量。
