摘要： 本项科研对南京市3头乳用种公牛（1059、1479、1391号）进行育种值的无偏估计。它们的女儿牛144头分布在江苏、安徽2个地区6个牧场，其女儿数不等。在地区差异较大，分布较广的情况下，作者应用混合线性模型（有交互作用）的数学模型，和最优线性无偏预测来估计公牛效应（Si）、地区效应（Dj）和（公牛×地区）交互作用效应。以公牛为随机效应、地区为固定效应，按照Dr. C. Henderson的方法在计算随机效应时，其系数矩阵的相应部分加k和λ值，并加约束条件ΣDj=0。通过Sharp PC 1501，计算结果为：(μ*+^)= 4385，(S*+^)#-1= -79，(S*+^)#-2= +47，(S*+^)#-3= 32，(D*+^)#-1=1153，(D*+^)#-2= -1153；（SD）(*+^)#-(11)= 2，（SD）(*+^)#-(12)= -42，（SD）(*+^)#-(21)= 24，（SD）(*+^)#-(22)= 0，（SD）(*+^)#-(31)= -26，（SD）(*+^)#-(32)= 42，经过数据分析认为：(1)种公牛的优劣名次为I#-(1479)＞I#-(1391)＞I#-(1059)；(2)地区效应差异很大，两地区相差2306公斤；(3)（公牛×地区）交互作用存在，其显著性和实用意义有待进一步研究 。
