6 生物数学学报 第15举 d=41-h42+43-44+451-52+53-45, h2=4（-41-42+3+h+451+452-53-4s4) h=4-41+42-4g++1-452+53-45 同理可获得其它模型参数的极大似然估计.关于二阶遗传参数的估计结合1和2.1有关内容 获得.关于分离世代个体的主基因型归类可参考文献2,4,8刚 3多基因存在的鉴定 从1可知，，号和可剖分为多基因方差组分和环境方差组分两部分.在只有B, 和B2群体或只有F2群体的情况下，不能将多基因方差组分和环境方差组分分开.因此，在 多基因鉴定时，尚需利用亲本和下1群体以获得环境方差的估计，从而鉴定多基因的存在 记x1,22:和3，以及n1,n2和n3分别为P1,F1,P2群体观测值和样本容量.由假定可知： 工~N(4,)(行=1,2,3).因此，由B:和B2(或F2)群体鉴定多基因存在的样本似然函数分 别为2,4 e)(). 1 (6a) La()=II()If()IIf) i∑maf(iHa,) (6b) 通过构造H0:不存在多基因(=和=或=，多基因效应的平均数为0)： H。:存在多基因（σ>0和>0或σ>0，多基因效应的平均数不为0）.计算两种假设 下的最大似然函数值Lo和L。,构造出似然比统计量入=2血L。-nLo)~x以鉴定多基因 是否存在，其中，山为两种假设下相差的遗传叁数个数.采用IECM算法获得分布参数的极 大似然估计值，E步骤和迭代CM步骤的CM1与S2.2相似，下面给出利用P1,F1,P2和F2 鉴定多基因存在的迭代CM2和迭代CM步骤.迭代CM2是固定分布平均数与误差方差条件 下用选代公式求多基因方差组分，选代CM:是固定分布平均数和多基因方差组分条件下用选 代公式求误差方差.。和？的选代公式如下 3 ne (8) =1j=1 11 =1=1 19-09 China Academie JouaEleetronic Pablishing House.All rights reserved.hutp:/www.enki.net生 物 数 学 学 报 第 卷 砚一 。 一 。 。 一 。 。 一 。 。 一 。 , 一 壹 一 。 一 。 。 。 。 。 一 。 一 。 , · 一 。 十 。 一 。 。 。 一 。 。 一 。 八 同理可获得其它模型参数的极大似然估计 关于二 阶遗传参数的估计结合 互 和 互 有关内容 获得 关于分离世代个体的主基因型归类可参考文献 , , 多基因存在的鉴定 从 盯 可知 , 叹 , 峭 和 峭 可剖分为多基 因方差组分和 环境方差组分两部分 在只有 和 群体或只有 群体的情况下 , 不能将多基因方差组分和环境方差组分分开 因此 , 在 多基 因鉴定时 , 尚需 利用 亲本和 群体以 获得 环 境方差 的估计 , 从而 鉴定多基因的存在 记 、, ‘ 和 、 , 以及 , 和 分别为 , , 群体观测值和样本容量 由假定可知 勺 一 州户 , 嵘 , , , 因此 , 由 和 或 群体鉴定多基因存在 的样本似然函数分 别为 【 , 〕 了 ‘ 。, , 。置 了 、 。 , 。了 了 ‘ 。 , , 署 感 一 天 几 无 艺 二 。 , ‘ 。 。, 孑 艺 二 、 。 亡, 叮爹 , 茗 茗 几 , ,‘ 。 , 。营 了 、 。 , , 置 。 , , 奢 ‘ 葱 百 性 人 艺 二 亡 ‘ 。 。, 口言 百 通过构造 。 不存在多基因 峭 心 和 峭 嵘 或 峭 二 心 , 多基因效应的平均数为 。 存在多基因 二军 。 和 二 。 或 嵘 。 , 多基因效应的平均数不为 计算两种假设 下的最大似然函数值 和 。 , 构造 出似然 比统计量 入 。 一 一 二 夯以鉴定多基因 是否存在 , 其中 , 为两种假设下相差 的遗传参数个数 采用 算法获得分布参数的极 大似然估计值 , 步骤和迭代 步骤 的 与 壮 相似 , 下面给出利用 , , 和 鉴定多基因存在的迭代 和迭代 屿 步骤 迭代 是固定分布平均数与误差方差条件 下用迭代公式求多基因方差组分 , 迭代 屿 是固定分布平均数和多基因方差组分条件下用迭 代公式求误差方差 峭 。 和 心 的迭代公式如下 二 , 、 , 一 尸 ‘ 。 一 , 了 艺、 , 概艺阁一 ‘, 口 , 人 介 心 艺公 夕 艺闪 一 。、 ’ 艺 ”矛艺 切 , 、 一 。 ‘ ’ 艺 ‘ 艺 。‘ 艺 二 , ‘ , 玄 夕 二 了