关于asp种畜blup育种值的程序设计

《关于asp种畜blup育种值的程序设计》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、关于ASP种畜BLUP育种值的程序设计关于ASP种畜BLUP育种值的程序设计1948年。美国学者Henderson提出了BLUP法，即最佳线性无偏预测(BestLinearUnbiasedPrediction)法。该方法具有估计值方差最小、可消除因选择或淘汰等原因造成的偏差等特性，获得的种畜育种值精确性较高。是进行种畜遗传评定的主要方法之一。该方法对于低遗传力性状和限性性状的效果明显，可广泛应用于奶牛各种性状的评定中。采用ASP语言，结合Excel文件开发了一个计算单性状BLUP的应用程序，实现种畜育种值的自动计算功能。l基本原理假设某一群体存在若干公畜。则任一公畜后代的观测值可用

2、以下模型表示：Yjk=U+hi+Sij+e(1)其中，Y是后代的成绩观察值，u是总体均值，h。是第i个畜群的固定效应，S.是第i个畜群第i个公畜的随机效应，而e是随机剩余效应。若(1)式采用矩阵形式表示，则为：Y=Xb十Zu+e(2)其中，Y为观察值向量，X为已知固定冈子的结构矩阵，b为未知固定向量，z为已知随机因子的结构矩阵，u为未知随机向量，e为随机效应向量。这里的就是要计算的家畜育种值估计量，即BLUP育种值。2ASP程序设计打开一个记事本文件，将文件名和后缀更改为index.asp，并保存。关于ASP种畜BLUP育种值的程序设计，其全部代码如下：0Ork(i，1)=k(j，

3、3)<>0ThenC(i,j)：O.5：c(i，i)=0.5Elsefk(i，2)=k(j，2)<>0Ork(i，3)=kO，3)<>0ThenⅡi=jThene(i，j)=0.25：cO，i)=0.25ElseIfk(i,2)：k0,2)<>0Andk(i,3)=k0,3)<>0ThenC(j，j)：O.5：c0，i)=0.5EndINext：NextFori=lToxIfk(i,2)+k(i,3)=0ThenC(i，i)=1ElseIfk(j，2)k0,3)=0ThenForj：1ToxIfk(i,2)=k0,1)Th

4、enp=jIfk(i，3)=k0,1)Thenp=jC(i，i)=1Form=lTOi一1C(i,m)=c(p，m)0.5：c(m，i)-c(i，m)Next：NextElseForj=1ToIfk(i，2)=k(j，1)Thenp=jIfk(i,3)=k0,1)Thenq=jForm=lToi-1C(i，m)=0.5(C(p，m)+c(q，m))：c(m，i)=c(i,m)Nexte(1一，1)=O。5e(P,q)+1NextEndIf：Next分子矩阵逆阵Dimbb，ee，kkFori=lTox：Fori-1Toxe(i,j)=c(1一，j)Ifj：iThene(i,j+x)=

5、1Ifj<>iThene(i,j+x)=0Next：NextFori：lTox：Fort=iToxIfe(t，i)<>0ThenFori-1To2xbb=e(i,j)：e(i,j)=e(t,j)：e(t,j)=bbNextcc=l/e(i，i)Forj=lTo2*xe(i'j)=eee(i,j)NextEndIf：：NextForj=1ToxIfj<>iThenkk=-e(j，i)Form-1To2*xe(i，m)-ej,m)+kke(i，m)NextEndIf：Next：NextFori=lTox：Fori-1ToxC(i,j)-e(i,x+j

6、)Next：Next输出矩阵Fori=lTox+z：For._lTox+zl(i，j)=0Ifi=jThenl(j，j)=b(i)Next：NextFori=lTox：Fori-1ToZl(z+i'j)=k(i,2j+2)：lO,z+i)=k(i，2j+2)Next：NextFori=z+lTox+z：Fori_z+1Tox+zl(i'j)=l(i'j)+he(i-z,j-z)Next：Next结构矩阵逆阵Dimbbb，CCC，kkkFori=lTox+z：Forj=1Tox+zf(i'j)=l)Ifj=iThenf(i'j+x+z)=lIfj<>iThenf(i,j+

7、x+z)=0Next：NextFori_1Tox+z：Fort=iTox+zIff(t，i)<>0ThenForj=1To2x+2Zbbb=f(i,j)：f(i,j)=f(t,j)：f(t,j)=bbbNextccc=l/f(i，i)Forj=1To2x+2Zf(i,j)=CCCf(i'j)NextEndIf：NextFori_lTox+zIfi<>iThenkkk=一f(j，i)Form=1To2*x+2*zfj,m)：fO,m)+kkkf(