家畜育種學07種畜的遺傳評估三blup育種值估計(存儲版)
【正文】 0 0 .0 0 6 81 .3 8 8 91 .3 8 8 90 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 .7 7 7 80 0 0 0 0 0 0 0 0 0 0 0 .3 8 8 91 .3 8 8 91 .3 8 8 90 0 0 0 0 0 0 0 0 0 0 0 0 .1 6 6 70 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 .0 0 6 8 0 .0 0 6 8 0 0 0 0 0 0 0 0 0 0 0 .0 1 3 60 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 .0 0 6 8 0 .0 0 6 8 0 .0 0 6 8 0 0 0 0 0 0 0 .0 2 0 42524232221151413121122211211aaaaaaaaaahhhh ?方程組的解為： ?綜合育種值為： ? ? ? ?1 2 . 9 3 6 91 3 . 0 8 5 91 5 0 . 9 6 7 61 4 3 . 0 8 0 7???? 22211211 ?hhhh ? ? ? ?2 . 1 5 4 41 . 0 8 9 53 . 7 4 8 04 . 3 7 1 72 . 8 6 5 7????? 1514131211 ?aaaaa? ? ? ?0 . 0 6 3 90 . 0 6 2 30 . 5 4 0 20 . 5 1 0 20 . 2 2 7 8????? 2524232221 ?aaaaa? ? ? ?1 . 3 4 3 80 . 6 0 3 92 . 6 8 1 0 3 . 0 3 1 21 . 9 0 1 7????? 54321 ?aaaaa ?如果說模型是 BLUP法的關鍵，那么，計算問題則是 BLUP法的難點。 uMbK ???yL?yLuMbK ?????yL?uMbK ???b u ? BLUP 估計一般方程 ? ? ? BLUP法前提條件 1. 所用的表型信息必須真實可靠，系譜資料必須正確完整 2. 所用的模型是真實模型； 3. 模型中的隨機效應的方差組分或方差組分的比值已知 yVXX)VX(b 11 ??? ????)bX(yVZGu 1 ?? ??? ? ?混合模型方程組的 一般形式 ?混合模型方程組的 簡化形式 ???????????????????????????????????????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 ??? 為 中與 個體對應的對角線元素 iudzzC i 第三節(jié) BLUP的計算技術 ?混合模型方程組的求解 ? 經典解法 ?先求出方程組的系數矩陣和等式右邊的向量，建立方程組，然后迭代求解 ?缺點：混合模型方程組往往很大，容易受計算機內存的限制，實際應用范圍不廣 ? 間接解法 ?不需建立方程組，直接構建觀測數據迭代公式，每次迭代讀入原始數據包括性狀觀測值和系譜記錄，并同時計算該次迭代的解 ?通用性不強，需要構建特定的數據迭代公式 ?經典的迭代方法 ? 高斯 賽德爾 迭代法 （ gaussseidel） ? 雅可比 迭代法（ jacobi） ? 松弛 迭代法（ relaxation） ?收斂標準 （ convergence criteria) ? 一般標準 ? 改進標準 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 ? 列表法計算分子親緣矩陣 ，父母親號先于個體號 ?個體 的父母未知時： ?個體 的父或母為 時 ： ?個體 的父母已知為 或 時： 121 ?? ti ?、ttt pq1?tta ipitti aaa ?? 121 ?? ti ?、p1?tta 0?? itti aa pqtt aa ?? )( iqipitti aaaa ???121 ?? ti ?、 ? 分子親緣矩陣逆矩陣的計算 1. 構造所有個體的系譜列表 ，父母親號先于個體號 2. 構建三角矩陣 ? 個體 的父母未知時： ? 個體 的父或母為 時： ? 個體的父母已知為 或 ，假設 ，這時： t p1?ttl 0?til 121 ?? ti ?、t?????????1210tppipill piti ??、ppititt fll 12 ???? ??L???????????????121021)(tqqiqppilpilll qiqipiti???、)(1 12qppiqitiqjpjtt ffllll ?????? ? ?? ?p q qp? ? 分子親緣矩陣逆矩陣的計算 3. 令 為 對角線元素組成的對角陣，讓 4. 按以下規(guī)則加入已知父母的個體的有關元素構建 D L 2)( 11 DA ?? ?1A?? 如果雙親已知為 和 ： ? 如果個體父或母已知 為： p
教學課件相關推薦
鄂ICP備17016276號-1