System.Exception: 字符串 'jiangyouxinxiwanghaiboerzhongzhuyuzhongyouxiangongsilaibinxinxiwangliuhenongmukejiyouxiangongsichengdouxinjinxinhaonongmuyouxiangongsishandongxinxiwangliuhejituanyouxiangongsixinxiwangliuheguf' 后的引号不完整。 “jiangyouxinxiwanghaiboerzhongzhuyuzhongyouxiangongsilaibinxinxiwangliuhenongmukejiyouxiangongsichengdouxinjinxinhaonongmuyouxiang”附近有语法错误。 在 DBUtility.DbHelperSQL.RunProcedure(String storedProcName, IDataParameter[] parameters, Int32& rowsAffected) 在 DAL.APT_Common.Update(String tbName, String filter, String where) 在 zhuanliShow.pingyinfenlei(String id)
【技术实现步骤摘要】
多性状育种值的测定方法、装置、设备及存储介质
[0001]本申请涉及生物育种
,尤其涉及一种多性状育种值的测定方法、装置、设备及存储介质。
技术介绍
[0002]目前,育种值测定算法主要包括BLUP(Best Linear Unbiased Prediction)算法和Bayesian算法,BLUP算法和Bayesian算法均是基于无穷小位点模型(infinitesimal model)的假设。其中,算法采用最小二乘法对线性模型进行回归,得到育种值的最佳线性无偏估计值;Bayesian算法采用贝叶斯参数估计方法进行育种值估计。
[0003]然而,BLUP算法和Bayesian算法在针对多性状基因组的育种值测定时,都会出现严重的性能瓶颈。具体地,在多性状基因组育种值计算中,内存消耗与性状数量n之间的关系正比于n2,即每增加1个性状,则内存消耗量增加((n+1)2‑
n2)/n2倍的。而当前国际育种的性状计算规模通常达到30+,则其内存消耗量比单性状的内存消耗量至少多900倍,导致育种值测定成本增大。
技术实现思路
[0004]本申请提供了一种多性状育种值的测定方法、装置、设备及存储介质,以解决当前育种值测定算法针对多性状基因组存在测定性能瓶颈的技术问题。
[0005]为了解决上述技术问题,第一方面,本申请提供了一种多性状育种值的测定方法,包括:
[0006]获取多性状育种参数,所述多性状育种参数包括多种形状对应的表型文件参数、系谱文件参数和基因组文件参数;
...
【技术保护点】
【技术特征摘要】
1.一种多性状育种值的测定方法,其特征在于,包括:获取多性状育种参数,所述多性状育种参数包括多种形状对应的表型文件参数、系谱文件参数和基因组文件参数;利用克罗内克定理,根据BLUP算法的预设等价式模型,对所述多性状育种参数进行分解,以提取每种性状对应的克罗内克分量;对每个所述克罗内克分量进行单性状测定,得到每个所述克罗内克分量对应的第一育种值;根据多个所述第一育种值,计算所述多性状育种参数对应的第二育种值。2.如权利要求1所述的多性状育种值的测定方法,其特征在于,所述利用克罗内克定理,根据BLUP算法的预设等价式模型,对所述多性状育种参数进行分解,以提取每种性状对应的克罗内克分量,包括:根据所述多性状育种参数,建立第一设计矩阵和第二设计矩阵,以及残差效应矩阵与单位矩阵之间的关系矩阵;利用所述克罗内克定理,根据所述第一设计矩阵、所述第二设计矩阵和残差效应,对所述预设等价式模型进行矩阵转换,得到每种性状对应的所述克罗内克分量。3.如权利要求2所述的多性状育种值的测定方法,其特征在于,所述根据所述多性状育种参数,建立第一设计矩阵和第二设计矩阵,以及残差效应与单位矩阵之间的关系矩阵,包括:根据所述多性状育种参数,建立所述第一设计矩阵、所述第二设计矩阵和残差效应矩阵;基于目标行数与性状数之间的比值,确定所述单位矩阵,所述目标行数为所述第一设计矩阵或所述第二设计矩阵的行数;对所述残差效应矩阵与所述单位矩阵进行克罗内克乘积,得到所述关系矩阵。4.如权利要求1至3任一项所述的多性状育种值的测定方法,其特征在于,所述预设等价式模型为:其中,X为第一设计矩阵,X
′
为第一设计矩阵的转置矩阵,R为残差效应矩阵,I为单位矩阵,Z为第二设计矩阵,Z
′
为第二设计矩阵的转置矩阵,G0为加性效应的多性状方差组分,G为加性效应的方差,为固定效应的最佳无偏线性估计值,为育种值的最佳无偏线性估计值,y为表型值,表示克罗内克乘积。5.如权...
【专利技术属性】
技术研发人员:张兴,丁偌楠,杨雨婷,牛安然,荆晓燕,张锡飞,杨昌雨,张丽萍,马康,梁志刚,龚华忠,闫之春,
申请(专利权)人:江油新希望海波尔种猪育种有限公司来宾新希望六和农牧科技有限公司成都新津新好农牧有限公司山东新希望六和集团有限公司新希望六和股份有限公司新希望集团有限公司新希望六和股份有限公司北京分公司,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。