家畜育種學(xué)07種畜的遺傳評(píng)估三blup育種值估計(jì)-展示頁(yè)
2025-05-27 02:06本頁(yè)面
【正文】 呈現(xiàn)連續(xù)性變異，通常是作為影響觀察值的協(xié)變量（回歸變量） ? 線性模型（ Linear model） ? 線性模型是指在模型中所包含的各個(gè)因子是以相加的形式影響觀察值，即它們與觀察值的關(guān)系為線性關(guān)系，但對(duì)于連續(xù)性的協(xié)變量也允許出現(xiàn)平方或立方項(xiàng)。 ? 隨機(jī)效應(yīng) 模型（ random model） ?若模型中除了總平均數(shù)外，其余的所有效應(yīng)均為隨機(jī)效應(yīng)，則稱此模型為隨機(jī)效應(yīng)模型或隨機(jī)模型 ? 混合 模型 （ mixed model） ?若模型中除了總平均數(shù)和隨機(jī)誤差之外，既含有固定效應(yīng)，也含有隨機(jī)效應(yīng)，則稱之為混合模型 ? 傳統(tǒng)的選擇指數(shù)法的基本假設(shè)是 不存在影響觀察值的系統(tǒng)環(huán)境效應(yīng) ，或者 在使用前剔除了系統(tǒng)環(huán)境效應(yīng) 。 ? 為克服以上缺陷， Henderson于 1948年提出了 BLUP方法，即 最佳線性無(wú)偏預(yù)測(cè) ，這個(gè)統(tǒng)計(jì)方法可同時(shí)估計(jì)固定效應(yīng)（例如系統(tǒng)環(huán)境效應(yīng)）和育種值。 ? 隨著計(jì)算機(jī)技術(shù)的高速發(fā)展，使這一方法的實(shí)際應(yīng)用成為可能，目前BLUP法已成為世界各國(guó)（尤其是發(fā)達(dá)國(guó)家）家畜遺傳評(píng)定的 規(guī)范方法 。 uMbK ???yL?yLuMbK ?????yL?uMbK ???b u ? BLUP 估計(jì)一般方程 ? ? ? BLUP法前提條件 1. 所用的表型信息必須真實(shí)可靠，系譜資料必須正確完整 2. 所用的模型是真實(shí)模型； 3. 模型中的隨機(jī)效應(yīng)的方差組分或方差組分的比值已知 yVXX)VX(b 11 ??? ????)bX(yVZGu 1 ?? ??? ? ?混合模型方程組的 一般形式 ?混合模型方程組的 簡(jiǎn)化形式 ???????????????????????????????????????yRZyRXubGZRZXRZZRXXRX1111111??????????????????????????????????? yZyXubAZZXZZXXX1 ??k22 uek ???2)( uV a r ?AGu ??2)( eV a r ?IRe ?? ?混合模型方程組的 度量 ??????????????????zzxxV arCG00Cub???????????????????? zzzxxzxxV arCCCCuub??kdduuC ovriiiiii uueauuuiiuu ????? 1/)()?,( 222?? ?????22 uek ??? 為 中與 個(gè)體對(duì)應(yīng)的對(duì)角線元素 iudzzC i 第三節(jié) BLUP的計(jì)算技術(shù) ?混合模型方程組的求解 ? 經(jīng)典解法 ?先求出方程組的系數(shù)矩陣和等式右邊的向量，建立方程組，然后迭代求解 ?缺點(diǎn)：混合模型方程組往往很大，容易受計(jì)算機(jī)內(nèi)存的限制，實(shí)際應(yīng)用范圍不廣 ? 間接解法 ?不需建立方程組，直接構(gòu)建觀測(cè)數(shù)據(jù)迭代公式，每次迭代讀入原始數(shù)據(jù)包括性狀觀測(cè)值和系譜記錄，并同時(shí)計(jì)算該次迭代的解 ?通用性不強(qiáng)，需要構(gòu)建特定的數(shù)據(jù)迭代公式 ?經(jīng)典的迭代方法 ? 高斯 賽德爾 迭代法 （ gaussseidel） ? 雅可比 迭代法（ jacobi） ? 松弛 迭代法（ relaxation） ?收斂標(biāo)準(zhǔn) （ convergence criteria) ? 一般標(biāo)準(zhǔn) ? 改進(jìn)標(biāo)準(zhǔn) iiijnijkjijkjijiki cxcxcrx )(11 1)1()()( ? ??? ??????)1(1)1()( )( ??? ??? ? kiiinjkjijiki xcxcrx)( )1()()()( ???? kikikiim p r ov e di xxxx ???? ? )1()(m a x titi xx??? ????? 21)(12)1()( )()( nitinititi xxx ? 列表法計(jì)算分子親緣矩陣 ，父母親號(hào)先于個(gè)體號(hào) ?個(gè)體 的父母未知時(shí)： ?個(gè)體 的父或母為 時(shí) ： ?個(gè)體 的父母已知為 或 時(shí)： 121 ?? ti ?、ttt pq1?tta ipitti aaa ?? 121 ?? ti ?、p1?tta 0?? itti aa pqtt aa ?? )( iqipitti aaaa ???121 ?? ti ?、 ? 分子親緣矩陣逆矩陣的計(jì)算 1. 構(gòu)造所有個(gè)體的系譜列表 ，父母親號(hào)先于個(gè)體號(hào) 2. 構(gòu)建三角矩陣 ? 個(gè)體 的父母未知時(shí)： ? 個(gè)體 的父或母為 時(shí)： ? 個(gè)體的父母已知為 或 ，假設(shè) ，這時(shí)： t p1?ttl 0?til 121 ?? ti ?、t?????????1210tppipill piti ??、ppititt fll 12 ???? ??L???????????????121021)(tqqiqppilpilll qiqipiti???、)(1 12qppiqitiqjpjtt ffllll ?????? ? ?? ?p q qp? ? 分子親緣矩陣逆矩陣的計(jì)算 3. 令 為 對(duì)角線元素組成的對(duì)角陣，讓 4. 按以下規(guī)則加入已知父母的個(gè)體的有關(guān)元素構(gòu)建 D L 2)( 11 DA ?? ?1A?? 如果雙親已知為 和 ： ? 如果個(gè)體父或母已知 為： p q1A?要加入的數(shù)值 中的位置 p1A?要加入的數(shù)值 中的位置 ? ),(),(),(),( qiiqpiip ),(),(),(),( qqpqqppp? ),(),(),(),( qiiqpiip ),(),(),(),( qqpqqppp ?如果是一個(gè) 非近交群體 ，則可直接構(gòu)建 ?如果雙親已知為 和 ： ?如果個(gè)體父或母已知為 ： p要加入的數(shù)值 中的位置 2 1 p q1A?),(),(),(),( qiiqpiip),(),(),(),( qqpqqppp),( ii要加入的數(shù)值 中的位置 3/4 2/3 1/3 1A?),( ii),(),( piip),(
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1