汪聪勇 赵杨杨 陈江凌
摘要 利用EXCEL函数中的IF函数实现对荷斯坦牛20个体型线性性状(包括体高、胸宽、体深、腰强度、尻角度、尻宽、蹄角度、蹄踵深度、骨质地、后肢侧视、后肢后视、乳房深度、中央悬韧带、前乳房附着、前乳头位置、前乳头长度、后乳房附着高度、后乳房附着宽度、后乳头位置和棱角性)鉴定原始数据快速、简便地转化为对应功能分。根据各性状权重,计算个体体躯容量、尻部、肢蹄、泌乳系统和乳用特征5个部位得分以及个体总分。利用荷斯坦牛体型性状功能分,对牛场牛只选种选配提供技术支持。
关键词 荷斯坦牛;体型鉴定;功能分转化
中图分类号 S823 文献标识码 A
文章编号 0517-6611(2020)18-0098-03
Abstract IF function of EXCEL function was used to realize the rapid transformation of the original data of 20 linear type traits in Holstein cows.20 body shape linear traits were body height, chest width, body depth, waist strength, rump angle, rump width, foot angle, heel depth, bone structure, rear legs set, rear legs rear view, udder depth, central ligament, fore udder attachment, front teat placement, front teat length, rear udder height, rear udder width, rear teat placement, and angularity.According to the weight of each trait, the scores of five parts including truncus capacity, rump, feet and leg composite, udder composite, dairy form and the total scores of individuals were calculated.The function scores of type traits in Holstein cattle were used to provide technical support for cattle breeding.
Key words Holstein cow;Somatotype identification;Functional score transformation
奶牛体型外貌性状是奶牛育种工作中的一个重要的评估性状。我国荷斯坦牛体型线性鉴定评分方法先后采用了50分制[1]和9分制线性鉴定[2]。对奶牛工作者而言,如何高速、有效地将这些体型鉴定数据运用到奶牛育种工作实践中成为工作的重点和难点。徐夕水等[3]利用EXCEL函数中的VLOOKUP函数实现50分制线性评分数据转化百分制得分,并以软件方式呈现。易礼胜等[4]利用EXCEL函数中的IF函数三级嵌套实现了对体型性状的测定原始数据转化为50分制。胡伶俐等[5]利用Microsoft Access 2003软件开发了奶牛线性评定计算机管理系统。魏晓莉等[6]基于WEB采用ASP.NET技术开发了奶牛体型外貌在线评定系统。马彦男等[7]利用EXCEL函数功能和录制宏功能,自动检索奶牛线性鉴定原始数据转化为相应的功能分。刘和凤等[8]利用将9分制评分中单个性状表现分为高(7~9分)、中(4~6分)、低(1~3分)3个区域,最终实现功能分转化。
笔者根据多年奶牛工作经验和不断的探索,参考《中国荷斯坦牛体型鉴定技术规程》GB/T 35568—2017[2]标准,选取20个奶牛线性性状数据进行育种分析。通过EXCEL表格函数IF函数,实现在一张EXCEL表格中对这20个体型线性鉴定原始数据快速、准确地转化为对应的功能分。最终实现只要會简单使用EXCEL表格操作,就可以对数据量庞大、繁琐的体型鉴定数据实现准确、无误地快速转化,并准确掌握鉴定牛场牛只整体及个体体型性状优缺点,为牛场育种工作提供可靠的技术支持。
1 体型线性鉴定原始数据录入
将9分制原始数据导入2007版本的EXCEL表单SHEET1中,表单的表头参考中国奶牛数据中心网站公布的体型线性鉴定表格[9]。在A4~F4依次输入鉴定员、鉴定信息、牛号、牛场、管理号和胎次等信息。在G4~AC4依次输入体高、前段、大小、胸宽、体深、腰强度、尻角度、尻宽、蹄角度、蹄踵深度、骨质地、后肢侧视、后肢后视、乳房深度、乳房质地、中央悬韧带、前乳房附着、前乳头位置、前乳头长度、后附着高度、后附着宽度、后乳头位置和棱角性共23个性状。其中,前段、大小和乳房质地这3个性状在《中国荷斯坦牛体型鉴定技术规程》GB/T 35568—2017中删去了,可以空出不填[2]。将AD4~AZ4的表头依次设计成上述23个性状功能分转化后所在单元格。将BA4~BF4表头依次设计成体躯容量、尻部、肢蹄、泌乳系统、乳用特征和个体总分所在单元格。
2 EXCEL函数设定
2.1 各性状原始数据功能分的转化
参考《中国荷斯坦牛体型鉴定技术规程》GB/T 35568—2017中体型鉴定各性状线性分与功能分对照表设计EXCEL表格函数[2]。针对奶牛的20个线性性状(体高、胸宽、体深、腰强度、尻角度、尻宽、蹄角度、蹄踵深度、骨质地、后肢侧视、后肢后视、乳房深度、中央悬韧带、前乳房附着、前乳头位置、前乳头长度、后乳房附着高度、后乳房附着宽度、后乳头位置和棱角性),共依次设计了20个IF函数。将这20个函数依次录入到SHEET1表单的AD4~AZ4单元格中,按下电脑键盘Enter键,即可以得到对应性状的功能分。选中AD4~AZ4单元格,按下鼠标左键下拉,即可得到表格中所有牛只的功能分。20个函数公式如下:
体高=IF(G4=1,57,IF(G4=2,64,IF(G4=3,70,IF(G4=4,75,IF(G4=5,85,IF(G4=6,90,IF(G4=7,95,IF(G4=8,100,IF(G4=9,95)))))))))
胸寬=IF(J4=1,55,IF(J4=2,60,IF(J4=3,65,IF(J4=4,70,IF(J4=5,75,IF(J4=6,80,IF(J4=7,85,IF(J4=8,90,IF(J4=9,95)))))))))
体深=IF(K4=1,56,IF(K4=2,64,IF(K4=3,68,IF(K4=4,75,IF(K4=5,80,IF(K4=6,90,IF(K4=7,95,IF(K4=8,90,IF(K4=9,85)))))))))
腰强度=IF(L4=1,55,IF(L4=2,60,IF(L4=3,65,IF(L4=4,70,IF(L4=5,75,IF(L4=6,80,IF(L4=7,85,IF(L4=8,90,IF(L4=9,95)))))))))
尻角度=IF(M4=1,55,IF(M4=2,62,IF(M4=3,70,IF(M4=4,80,IF(M4=5,90,IF(M4=6,80,IF(M4=7,75,IF(M4=8,70,IF(M4=9,65)))))))))
尻宽=IF(N4=1,55,IF(N4=2,60,IF(N4=3,65,IF(N4=4,70,IF(N4=5,75,IF(N4=6,79,IF(N4=7,82,IF(N4=8,90,IF(N4=9,95)))))))))
蹄角度=IF(O4=1,56,IF(O4=2,64,IF(O4=3,70,IF(O4=4,76,IF(O4=5,81,IF(O4=6,90,IF(O4=7,100,IF(O4=8,95,IF(O4=9,85)))))))))
蹄踵深度=IF(P4=1,57,IF(P4=2,64,IF(P4=3,69,IF(P4=4,75,IF(P4=5,80,IF(P4=6,85,IF(P4=7,90,IF(P4=8,95,IF(P4=9,100)))))))))
骨质地=IF(Q4=1,57,IF(Q4=2,64,IF(Q4=3,69,IF(Q4=4,75,IF(Q4=5,80,IF(Q4=6,85,IF(Q4=7,90,IF(Q4=8,95,IF(Q4=9,100)))))))))
后肢侧视=IF(R4=1,55,IF(R4=2,64,IF(R4=3,75,IF(R4=4,80,IF(R4=5,95,IF(R4=6,80,IF(R4=7,75,IF(R4=8,65,IF(R4=9,55)))))))))
后肢后视=IF(S4=1,57,IF(S4=2,64,IF(S4=3,69,IF(S4=4,74,IF(S4=5,78,IF(S4=6,81,IF(S4=7,85,IF(S4=8,90,IF(S4=9,100)))))))))
乳房深度=IF(T4=1,55,IF(T4=2,65,IF(T4=3,75,IF(T4=4,85,IF(T4=5,95,IF(T4=6,85,IF(T4=7,75,IF(T4=8,65,IF(T4=9,55)))))))))
中央悬韧带=IF(V4=1,55,IF(V4=2,60,IF(V4=3,65,IF(V4=4,70,IF(V4=5,75,IF(V4=6,80,IF(V4=7,85,IF(V4=8,90,IF(V4=9,95)))))))))
前乳房附着=IF(W4=1,55,IF(W4=2,60,IF(W4=3,65,IF(W4=4,70,IF(W4=5,75,IF(W4=6,80,IF(W4=7,85,IF(W4=8,90,IF(W4=9,95)))))))))
前乳头位置=IF(X4=1,57,IF(X4=2,65,IF(X4=3,75,IF(X4=4,80,IF(X4=5,85,IF(X4=6,90,IF(X4=7,85,IF(X4=8,80,IF(X4=9,75)))))))))
前乳头长度=IF(Y4=1,50,IF(Y4=2,60,IF(Y4=3,70,IF(Y4=4,80,IF(Y4=5,90,IF(Y4=6,80,IF(Y4=7,70,IF(Y4=8,60,IF(Y4=9,50)))))))))
后乳房附着高度=IF(Z4=1,58,IF(Z4=2,65,IF(Z4=3,68,IF(Z4=4,70,IF(Z4=5,75,IF(Z4=6,80,IF(Z4=7,85,IF(Z4=8,90,IF(Z4=9,95)))))))))
后乳房附着宽度=IF(AA4=1,58,IF(AA4=2,65,IF(AA4=3,68,IF(AA4=4,70,IF(AA4=5,75,IF(AA4=6,80,IF(AA4=7,85,IF(AA4=8,90,IF(AA4=9,95)))))))))
后乳头位置=IF(AB4=1,57,IF(AB4=2,65,IF(AB4=3,75,IF(AB4=4,80,IF(AB4=5,85,IF(AB4=6,90,IF(AB4=7,85,IF(AB4=8,80,IF(AB4=9,75)))))))))
棱角性=IF(AC4=1,57,IF(AC4=2,64,IF(AC4=3,69,IF(AC4=4,74,IF(AC4=5,78,IF(AC4=6,81,IF(AC4=7,85,IF(AC4=8,90,IF(AC4=9,95)))))))))
2.2 五部位功能分和总分计算
根据中国荷斯坦牛体型鉴定技术规程》GB/T 35568—2017中各性状和5个部位性状权重,计算5个部位得分以及个体总分[2]。在SHEET1表单中,在BA4~BF4单元格依次输入个体体躯容量、尻部、肢蹄、泌乳系统、乳用特征和体型总分公式重复功能分转化操作即可得到所有牛只部位评分和体型总分。
体躯容量=AD4×0.25+AG4×0.35+AH4×0.25+AI4×0.15
尻部=AJ4×0.4+AK4×0.45+AI4×0.15
肢蹄=AL4×0.25+AM4×0.15+AN4×0.15+AO4×0.25+AP4×0.2
泌乳系統=(AQ4×0.55+AS4×0.45)×0.2+(AT4×0.45+AU4×0.25+AV4×0.18+AQ4×0.12)×0.35+(AW4×0.3+AX4×0.3+AY4×0.14+AQ4×0.12+AS4×0.14)×0.45
乳用特征=AZ4×0.8+AN4×0.2
体型总分=BA4×0.18+BB4×0.1+BC4×0.2+BD4×0.42+BE4×0.1
3 数据应用
在对牛只体型进行评价时,参照《中国荷斯坦牛体型鉴定技术规程》GB/T 35568—2017将中国荷斯坦牛母牛体型鉴定等级划分为6个等级:优EX(90~100)、很好VG(85~89)、好佳GP(80~84)、好G(75~79)、一般F(65~74)和差P(65分以下)[2]。以单个牛场的体型鉴定数据为例,转化过后的20个性状功能分、5个部位功能分和体型总分可以利用EXCEL的排序功能对它们分别进行由高到低或由低到高的排序,牛场鉴定牛只各个性状、部位和总分的优劣按排列顺序一目了然。依次选中各个等级范围内的数据,就可以依次得到各个等级牛只在整个鉴定牛群中的数量。同时,利用EXCEL函数中的AVERAGE函数,就可以得到各个性状的平均值,再参照等级鉴定标准,可以评价该牛场牛只整体的体型状况。
4 讨论
在EXCEL 2003软件中 IF函数可以嵌套7层[10],也就是单纯的IF函数最多可以有8个判断,与奶牛体型线性鉴定9分制中的9个判断不相符。但是,随着EXCEL软件的不断完善, EXCEL 2007版及以上的版本中函数可以嵌套到64层,可以实现奶牛体型线性鉴定9分制中9个可变量的判断,从而实现原始9分线性数据转化为百分制功能分,进而对奶牛体型性状的优缺点作出准确、全面的判断。
奶牛体型线性性状直接影响奶牛的功能寿命[11] 、305 d产奶量[12-13]、乳房炎发病率[14]等经济性状。近年来,申慧芳等[15]、刘松柏[16]、李培娟[17]、毛杰[18]和脱征军等[19]分别对山西、武汉、河北、上海和宁夏地区的荷斯坦牛体型线性性状进行了统计分析,结果发现各地区荷斯坦牛体型性状均未达到理想状态,还有很多的改良空间。
申慧芳等[15]、刘松柏[16]、李培娟[17]、毛杰[18]还对荷斯坦牛体型性状的遗传力进行了估计,估计值分别为0.10~0.45、0.03~0.37、0.07~0.61、0.09~0.64。上述研究表明,不同地区相同体型性状估计的遗传力不尽相同。按照遗传力划分标准(高遗传力h2<0.3;中等遗传力0.1
5 结论
EXCEL 2007及以上版本中的IF函数可以实现荷斯坦牛9分制体型鉴定原始数据简单、便捷地转化为百分制功能分。同时,可以利用这些功能分统计个体体躯容量、尻部、肢蹄、泌乳系统、乳用特征和体型总分。各个性状、部位功能分反应了牛只单性状和各个部位的情况,有利于育种工作者在选种选配时,及时改良个体劣势或缺陷性状,巩固并加强优良性状。评估一个场牛只体型鉴定等级情况,可以让管理者精准掌握牛场牛只整体情况,进而制定相应的育种计划。
参考文献
[1] 储明星,师守堃.奶牛体型线性评定及其应用[M].北京:中国农业出版社,1999:6-15.
[2] 张沅,石万海,张胜利,等.中国荷斯坦牛体型鉴定技术规程:GB/T 35568—2017[S].北京:中国标准出版社,2017.
[3] 徐夕水,陈惠华,邱雅洁.EXCEL在奶牛线性外貌评分中的应用[J].畜牧与兽医,2003(7):24.
[4] 易礼胜,付茂忠.Excel在奶牛线性评分中的应用[J].中国草食动物,2006,26(2):65.
[5] 胡伶俐,付海涵.基于Access的奶牛体型线性评定管理系统的设计与应用[C]∥中国奶业协会.第三届中国奶牛发展大会论文集.北京:《中国奶牛》编辑部,2008:255-257.
[6] 魏晓莉, 张长利, 兰翠.基于WEB的奶牛体型外貌线性评定系统的设计与开发[J].东北农业大学学报, 2008,39(12):93-95.
[7] 马彦男,吴建平,朱静.基于EXCEL的奶牛体型外貌线性评分系统[J].中国草食动物,2009,29(3):30-33.
[8] 刘和凤,汪湛.荷斯坦牛体型线性鉴定计算机系统的应用[J].中国奶牛,2012(13):34-36.
[9] 中国奶牛数据中心.tixjd.xls手机版体型评定数据导入模板[EB/OL].[2020-01-15].https://www.holstein.org.cn.
[10] 陈基容.用Excel中IF的嵌套函数进行数据的高级筛选[J].电脑编程技巧与维护,2015(3):60-61,75.
[11] 吴红超,樊新忠,徐凯勇,等.荷斯坦奶牛体型特征对功能寿命的影响[J].西南农业学报,2008,21(2):455-459.
[12] 王力生,蔡克周,程茂基,等.荷斯坦牛乳房诸性状与305d产奶量的相关[J].草食家畜,2007(4):12-13,18.
[13] 朱凯,刘光磊,黄黎明,等.上海地区荷斯坦牛体型外貌与产奶性状的相关分析[J].中国奶牛,2014(Z1):17-20.
[14] 刘贤侠,王建梅,张晓燕,等.荷斯坦牛乳房部分性状对乳房疾病的影响[J].中国奶牛,2011(2):34-37.
[15] 申慧芳,王钦德.山西省荷斯坦牛体型性状参数估计[J].山西农业大学学报(自然科学版),2006,26(2):189-191.
[16] 刘松柏.武汉地区荷斯坦牛体型线性性状的遗传分析[D].武汉:华中农业大学,2008.
[17] 李培娟.河北荷斯坦牛体型性状分析[D].邯郸:河北工程大学,2010.
[18] 毛杰.上海地区荷斯坦牛体型性状与产奶性状、SCS和寿命性状的遗传分析[D].南京:南京农业大学,2015.
[19] 脱征军,李委奇,田佳,等.宁夏规模化奶牛场荷斯坦牛体型线性评定浅析[J].中国奶牛,2017(10):20-24.
[20] 魏趁,赵俊金,黄锡霞,等.新疆地区西门塔尔牛核心群选择[J].中国农业科学,2019,52(5):921-929.