页数：48页文件大小：323.5 KB来源： 生物统计与田间试验ppt课件
本课件介绍农业科学中的主要参数及其估计量的评选标准，包括无偏性、有效性与一致性，并系统讲解矩法、最小二乘法和极大似然法三种常用参数估计方法的原理、计算过程与适用条件。
第八章参数估计方法
第一节农业科学中的主要参数及其估计量的评选标准
二、参数估计量的评选标准
(一) 数学期望样本平均数的平均数就是一种数学期望。例如，一个大豆品种的含油量为20%，测定一次可能是大于20%，再测定可能小于20%，大量反复测定后平均结果为20%，这时20%便可看作为该大豆品种含油量的数学期望，而每单独测定一次所获的值只是1个随机变量。抽象地，随机变量的数字特征是指随机变量的数学期望值
参数估计量的评选标准（续）
对于离散型(间断性)随机变量y的分布列为：P{y=yi}=pi，其中，i=1，2，…，那么随机变量y的数学期望E(y)为：(8·1)这样可以求得总体平均值。对于连续型随机变数y的数学期望E(y)为：(8·2)其中f(y)为随机变量y的概率密度函数，这样可以求得总体均值
参数估计量的评选标准（续）
用D(y)表示方差，有D(y)=E [y－E(y)]2(8·3) 这就是随机变量函数的数学期望。同理，离散型随机变量方差的数学期望为：(8·4) 连续型随机变量方差的数学期望为：(8·5)
数学期望有这样一些常用的性质
（二）参数估计量的评选标准评价估计量优劣的标准
主要有无偏性、有效性、相合性等(1) 无偏性参数估计量的期望值与参数真值是相等的，这种性质称为无偏性，具有无偏性的估计量称为无偏估计量。例如，在抽样分布中已经介绍了离均差平方和除以自由度得到的均方的平均数等于总体方差，即该均方的数学期望等于相应总体参数方差，这就是说该均方估计量是无偏的。估计量的数学期望值在样本容量趋近于无穷大时与参数的真值相等的性质称为渐进无偏性，具有渐进无偏性的估计量称为渐进无偏估计量
参数估计量的评选标准评价估计量优劣的标（续）
(2) 有效性无偏性表示估计值是在真值周围波动的一个数值，即无偏性表示估计值与真值间平均差异为0，近似可以用估计值作为真值的一个代表。同一个参数可以有许多无偏估计量，但不同估计量的期望方差不同，也就是估计量在真值周围的波动大小不同。估计量的期望方差越大说明用其估计值代表相应真值的有效性越差；否则越好，越有效。不同的估计量具有不同的方差，方差最小说明最有效。
如果一个无偏估计量相对与其它所有可能无偏估计量，其期望方差最小，那么称这种估计量为一致最小方差无偏估计量
参数估计量的评选标准评价估计量优劣的标（续）
(3) 相合性用估计量估计参数涉及一个样本容量大小问题，如果样本容量越大估计值越接近真值，那么这种估计量是相合估计量。除以上三方面标准外，还有充分性与完备性也是常考虑的。充分性指估计量应充分利用样本中每一变量的信息；完备性指该估计量是充分的唯一的无偏估计量
第二节矩法
一、矩的概念矩( moment )分为原点矩和中心矩两种。对于样本y1,y2,…yn，各观测值的k次方的平均值，称为样本的k阶原点矩，记为，有, 用观测值减去平均数得到的离均差的k次方的平均数称为样本的k阶中心矩, 记为或，有。对于总体y1,y2,…yN，各观测值的k次方的平均值，称为总体的k阶原点矩，记为，有；用观测值减去平均数得到的离均差的k次方的平均数称为总体的k阶中心矩，记为或，有
二、矩法及矩估计量
所谓矩法就是利用样本各阶原点矩来估计总体相应各阶原点矩的方法，即(8·6)也可以用样本各阶原点矩的函数来估计总体各阶原点矩同一函数，即若Q=f ( E(y),E(y2),…,E(yk) ), 则由此得到的估计量称为矩估计量
例8.1 现获得正态分布的随机样本y1,
y2 ,…yn，要求正态分布参数和的矩估计量。首先，求正态分布总体的1阶原点矩和2阶中心矩
然后求样本的1阶原点矩和2阶中心矩
然后求样本的1阶原点矩和2阶中心矩，为最后，利用矩法，获得总体平均数和方差的矩估计故总体平均数和方差的矩估计值分别为样本平均数和样本方差，方差的分母为n
单峰分布曲线还有二个特征数，即偏度
单峰分布曲线还有二个特征数，即偏度( skewness )与峰度( kurtosis )，可分别用偏度系数和峰度系数作测度。偏度系数( coefficient of skewness )是指3阶中心矩与标准差的3次方之比；峰度系数( coefficient of kurtosis )是指4阶中心矩与标准差的4次方之比。
当偏度为正值时，分布向大于平均数方向偏斜；偏度为负值时则向小于平均数方向偏斜；当偏度的绝对值大于2时，分布的偏斜程度严重。当峰度大于3时，分布比较陡峭，峰态明显，即总体变数的分布比较集中
例8.1 现获得正态分布的随机样本y1（续）
由样本计算的偏度系数(8·7)峰度系数(8·8)
例8.2 计算
本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。
表3.4数据资料(140行水稻产量)所属分布曲线的偏度和峰度。表3.4 140行水稻产量(单位：克)
首先，计算样本的2、3、4阶中心矩
首先，计算样本的2、3、4阶中心矩，以及标准差估计值
例8.2 计算（续）
然后，根据矩法原理，该分布的偏度与峰度估计值分别为：因此，说明资料比较集中在平均数左右，分布曲线并不是特别陡峭
例8.3 例6.9为研究籼粳稻杂交F5代
系间单株干草重的遗传变异，随机抽取76个系进行试验，每系随机取2个样品测定干草重(g/株)。按单向分组方差分析进行分析，结果见表6.9。此处用来说明由矩法估计误差、遗传方差和干草的遗传力h2。因为76个系是随机抽取的，因而为随机模型。方差结果说明系间差异显著，因而系间效应存在。根据矩法，首先应求出系间和误差变异来源的样本均方和总体期望均方(表6.9)。然后，利用矩估计原理，令样本的均方与总体相应变异的期望均方相等，从而求出和的矩估计值
例8.3 例6.9为研究籼粳稻杂交F5（续）
此处E(MS系统间)=E[Tt-E(Tt)]2，(Tt为各个系统的总和数)=E(MS误差)=E(e2)= ，(e为误差)因而
第三节最小二乘法
从总体中抽出的样本观察值与总体平均数是有差异的，这种差异属于抽样误差。因而，在总体平均数估计时要尽可能地降低这种误差，使总体平均数估计值尽可能好。参数估计的最小二乘法就是基于这种考虑提出的。基本思想是使误差平方和最小，达到在误差之间建立一种平衡，以防止某一极端误差对决定参数的估计值起支配地位。这有助于揭示更接近真实的状况。具体方法是为使误差平方和Q为最小，可通过求Q对待估参数的偏导数，并令其等于0，以求得参数估计量
例8.4 用最小二乘法求总体平均数的估计量
若从平均数为的总体中抽得样本为y1、y2、y3、…、yn，则观察值可剖分为总体平均数与误差ei 之和，总体平均数的最小二乘估计量就是使yi 与间的误差平方和为最小，即为最小
例8.4 用最小二乘法求总体平均数的估（续）
为获得其最小值，求Q对的导数，并令导数等于0，可得：即总体平均数的估计量为：因此，算术平均数为总体平均数的最小二乘估计。这与矩法估计是一致的
估计离均差平方和的数学期望
因而，估计为：与矩法所得不同，而与常规以自由度为除数法一致
例8.5 求例6.13的两向分组方差分析
本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。
资料缺1个小区(表8.1)的最小二乘估计量和估计值。表8.1 生长素处理豌豆的试验结果处理(A)组(B)总和Ti平均ⅠⅡⅢⅣ对照(CK) 赤霉素 +ye动力精 吲哚乙酸 硫酸腺嘌吟 马来酸 总和T j +yeT=1444+ye
从第6章可知，这种资料模式的线性模型为
从第6章可知，这种资料模式的线性模型为：按照最小二乘法的估计原理，使该模型的约束条件为：，和误差项服从正态分布。为最小时可以求出效应和缺失小区ye的估计量，即
从而，最小二乘估计量分别为
从而，最小二乘估计量分别为
例8.5 求例6.13的两向分组方差分（续）
因而表8.1中，缺失小区的估计值可由下式求出：解上述方程，最小二乘估计值为：ye=65.6。缺区估计是根据线性模型，以及最小二乘法的原理得到的。不过，试验中尽可能不要缺区，因为缺区估计尽管可以估计缺区的值，但是误差的自由度将减少，本试验的误差自由度将减少1
例8.5 求例6.13的两向分组方差分（续）
一般地，若m个自变数x1、x2、x3、…、xm与依变数y存在统计模型关系(8·9)其中，为待估参数。通过n次观测(n＞k)得到n组含有x1i , x2i ,…xmi , yi ( i=1,2,…,n )的数据以估计。其最小二乘估计值为使(8·10) 为最小的。这种估计方法称为参数估计的最小二乘法( least squares )，或最小平方法
第四节极大似然法
所谓极大似然法( maximum likelihood method )是值选择使事件发生概率最大的可能情况的参数估计方法。极大似然法包括二个步骤
一、似然函数
对于离散型随机变量，似然函数是多个独立事件的概率函数的乘积，该乘积是概率函数值，它是关于总体参数的函数。例如，一只大口袋里有红、白、黑3种球，采用复置抽样50次，得到红、白、黑3种球的个数分别为12，24，14，那么根据多项式的理论，可以建立似然函数为：其中p1，p2，p3分别为口袋中红、白、黑3种球的概率(p3= )，它们是需要估计的
似然函数（续）
对于连续型随机变量，似然函数是每个独立随机观测值的概率密度函数的乘积，则似然函数为：(8·11)若yi 服从正态分布，则，上式可变为：(8·12)
二、极大似然估计
所谓极大似然估计就是指使似然函数为最大以获得总体参数估计的方法。其中，所获得的估计总体参数的表达式称为极大似然估计量，由该估计量获得的总体参数的估计值称为总体参数的极大似然估计值。为了计算上的方便，一般将似然函数取对数，称为对数似然函数，因为取对数后似然函数由乘积变为加式，其表达式为：(8·13)
极大似然估计（续）
求极大似然估计量可以通过令对数似然函数对总体参数的偏导数等于0来获得，即当，有（k=1，2，…，l）(8·14)由此获得总体参数的极大似然估计量
例8.6 设y1 , y2 , … ,
yn是正态总体的随机样本，求正态分布参数的极大似然估计量。似然函数为：取对数，得
那么似然方程组为
解得
那么似然方程组（续）
因此，正态分布总体平均数的极大似然估计量为：当总体平均值为未知时，方差估计量为:当总体平均值为已知时，方差估计量为:
例8.7 求红、白、黑球事例中p1
p2，p3的极大似然估计值。由可获得对数似然函数其中，C为常数。分别求对p1，p2的偏导数，并令为0，得似然方程组
联立求解
联立求解，得：显然，极大似然估计值等于其观测频率
例8.8 两个亲本的基因型分别
为AABB和aabb，这两个亲本杂交后F2出现了4种基因型，分别为A_B_、A_bb、aaB_和aabb，得到四种基因型的个数分别为c、d、e、f，已知AA和BB两对基因间存在连锁关系，现欲估计重组率？设重组率为r，根据遗传学推导，可以得到4种基因型的概率见表8.2
表8.2 F2群体基因型的分离情况基因型A_B
本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。
表8.3 介绍由两对连锁主基因控制的F2群体1
本页含数据表，源文本中数字粘连、无法安全切分，表格未作自动还原，需对照课件原图人工补录。
配子及概率AB (1－r)/ (1－r)/2AB (1－r)/2AABB(1－r) (1－r)/4AaBBr(1－r)/4AaBa(1－r) (1－r)/ (1－r)/ (1－r)/ (1－r)/4ab(1－r)/2AaBa(1－r) (1－r)/4aaBbr(1－r)/4Aabb(1－r)2/4表8.3F2群体的基因型及其概率
按多项式分布
按多项式分布，可以根据概率函数得到似然函数为：(8·15) 若以代入上式，则似然函数和对数似然函数分别为：(k是常数项) (8·16) (8·17)
例8.8 两个亲本的基因型分别（续）
对上式求导数，并令导数为0，可得方程：上式化解为一元二次方程(8·18) 在的两个解中取一个符合遗传规律的解，那么，重组率的解为：。重组率方差估计量为：(8·19)
例8.8 两个亲本的基因型分别（续）
对于本例，有取正根，=0.7366，由此，=0.142
三、关于三种估计方法的讨论
上述3种参数估计方法对比
3种常用方法的不同要求
3种常用方法的不同要求:
下载地址：{{downurl}}
{{{a}}}
{{#b}}{{{b}}}
{{/b}} {{/title}} {{#nolog}}登陆注册
{{/nolog}} {{#djsx}}点我刷新
{{/djsx}} {{#gocz}} {{/gocz}} {{#ffqr}} {{/ffqr}}