圈养林麝生长曲线拟合

蔡永华， 杨营， 付文龙， 朱萍， 冯达勇， 雷毅， 张明

(1.四川省养麝研究所，四川都江堰611845； 2.四川农业大学旅游学院，四川都江堰611830；3.四川农业大学动物科技学院，四川雅安625014)



DOI:10.11984/j.issn.1000-7083.20150246

圈养林麝生长曲线拟合

蔡永华1， 杨营1， 付文龙1， 朱萍1， 冯达勇1， 雷毅2， 张明3*

(1.四川省养麝研究所，四川都江堰611845； 2.四川农业大学旅游学院，四川都江堰611830；3.四川农业大学动物科技学院，四川雅安625014)

为探索圈养林麝Moschusberezovskii体质量、体尺的生产发育规律, 选用320只林麝(雄麝160只，雌麝160只)，测定其1～400周龄的体质量和体尺(额宽、头长、耳长、脊柱长、胸围、肱骨长、大腿长和小腿长)，并用Logistic、Bertalanffy、Gompertz模型进行生长曲线拟合。3种模型均能很好地拟合周龄与体质量之间的回归关系，根据回归方程决定系数判断，Bertalanffy对雄麝(R2=0.966)、雌麝(R2=0.954)体质量的拟合优于Logistic和Gompertz模型。体尺与周龄的拟合结果显示：Bertalanffy模型优于Logistic和Gompertz模型；肢体体尺与周龄的拟合度较高，头部体尺与周龄拟合度较低。在所检测的体尺指标中，Bertalanffy对雄麝(R2=0.927)、雌麝(R2=0.933)小腿长的拟合度最高。林麝体质量、体尺的生产发育规律可用Bertalanffy模型拟合。

生长曲线；林麝；Logistic模型；Bertalanffy模型；Gompertz模型

动物生长曲线是动物个体生长发育规律的反映，通常采用非线性数学模型进行描述(Migonetal.，1990；盛志廉，吴常信，1995)。目前常用Logistic、Bertalanffy和Gompertz模型来探讨动物的生长发育规律(Woolson & Leeper，1980；von-Rosen，1991;Nunez-Antonetal.，2001)。生长曲线模型可以用于预测动物的体质量、体尺等生长发育指标，并进一步以此判断动物的饲养管理和营养水平(李娟等，2009)。

林麝Moschusberezovskii是国家Ⅰ级保护动物(Yangetal.，2003；Mengetal.，2006)，雄麝分泌的麝香具有极高的药用价值。随着林麝的人工驯养与繁殖的成功，已在林麝的基础生理学(姜海瑞等，2007；何迟酩等，2012)、行为学(Mengetal.，2011)、遗传学(Xiaetal.，2006；Zhangetal.，2007；Zhaoetal.，2008)和常见疾病(朱承嗣等，2012；周瑜等，2013)等方面进行了探索与研究。林麝体质量、体长和肩高(体高)的生长发育规律已有报道(竭航等，2014)，但已有的研究中，测量的个体数少，持续测定时间短，头长、耳长、胸围、大腿长和小腿长等体尺没有测定，而且使用的拟合模型不是经典模型。本研究利用Logistic、Bertalanffy和Gompertz模型对圈养林麝的体质量、额宽、头长、耳长、脊柱长、胸围、肱骨长、大腿长和小腿长的生长曲线进行拟合，为林麝生长发育研究提供更多的基础数据。

1 材料与方法

1.1 实验材料

四川省养麝研究所都江堰和马尔康养麝场的圈养林麝320只(雄麝160只，雌麝160只)。实验麝群自由采食、饮水。其他饲养管理措施参见程世国和邹真慧(1991)的研究。

1.2 测定指标

从出生到400周龄，每周测定实验林麝的空腹体质量，以及体尺指标：额宽、头长、耳长、脊柱长、胸围、肱骨长、大腿长和小腿长。

体尺指标定义如下：

额宽：两眼眶外侧突起点之间的距离；

头长：从两耳之间中部到鼻端的长度；

耳长：耳背根部到耳尖端的长度；

脊柱长：从头部枕骨后第一脊椎到尾根前第一尾椎之间，沿脊柱弯曲线的长度；

胸围：紧靠前肢腋下，垂直围绕胸部一周的长度；

肱骨长：前肢从肩关节凸出顶点到肘关节后端顶点间的直线距离；

小腿长：从后肢脚跟后端凸点到蹄冠后部间的直线距离；

大腿长：从后肢膝关节前段凸出顶点到脚跟后端突出点之间的直线距离。

1.3 生长曲线模型

采用3种最常见的动物生长曲线模型：Logistic、Bertalanffy和Gompertz(表1)。

表1 3种常用动物生长曲线模型

1.4 数据处理

用SPSS 20拟合模型，建立生长曲线模型，推算出拐点体质量和最大增重、增长，并且根据拟合度(R2)对生长曲线模型进行评价。

2 结果分析

2.1 林麝体质量与周龄拟合

雌雄麝体质量的生长拟合方程和拟合度(R2)见表2。结果显示：雄麝周龄与体质量的Logistic、Bertalanffy和Gompertz拟合方程决定系数分别为0.960、0.966和0.964，雌麝3种拟合方程的决定系数分别为0.938、0.954和0.950。结果表明：Logistic、Bertalanffy和Gompertz 3种模型均能很好拟合周龄与体质量的回归关系。但运用Bertalanffy模型计算雌麝的拐点周龄为0，因此相对较好的周龄-体质量拟合是Gompertz模型。

2.2 不同模型拟合的体质量生长曲线

利用Logistic、Bertalanffy和Gompertz回归方程绘制雌雄麝的体质量生长曲线(图1)。

在3种模型中，Bertalanffy模型在雌雄麝周龄-体质量拟合中具有最高的拟合度。在Gompertz模型中，雌雄麝的初生体质量(0周龄)分别是1.236 kg和1.339 kg。实测雌雄麝初生体质量为0.48 kg±0.13 kg和0.53 kg±0.23 kg。在Logistic和Bertalanffy模型中，雌雄麝的初生体质量分别为0.567 kg、0.785 kg和0.866 kg、0.762 kg。

2.3 头部体尺与周龄的拟合

表2 林麝体质量与周龄的3种拟合方程

Logistic模型 Bertalanffy模型 Gompertz模型

头部体尺指标额宽、头长和耳长与周龄的拟合方程和拟合度见表3。结果显示：Logistic、Bertalanffy和Gompertz模型中，周龄与额宽拟合方程的决定系数为0.717～0.804，周龄与头长拟合方程的决定系数为0.860～0.899，周龄与耳长拟合方程的决定系数为0.810～0.941。雄麝周龄与额宽、头长和耳长拟合方程的决定系数分别为0.799～0.804、0.895～0.899和0.810～0.821，雌麝周龄与额宽、头长和耳长拟合方程的决定系数分别为0.717～0.721、0.860～0.866和0.939～0.941。综上所述，雌麝耳长与周龄用3种模型均能很好地拟合，而额宽和头长的拟合程度相对较低。以拟合度判断周龄与额宽、耳长的回归关系中，Bertalanffy模型具有最好的拟合度；雄麝头长与周龄，Logistic模型具有最好的拟合度，雌麝头长与周龄，Bertalanffy模型具有最好的拟合度。

2.4 躯体体尺与周龄的拟合

躯体体尺指标：脊柱长和胸围与周龄的拟合方程和拟合度见表4。结果显示：Logistic、Bertalanffy和Gompertz中，雄麝脊柱长与周龄的拟合方程决定系数分别为0.940、0.943和0.942，而雌麝脊柱长与周龄的拟合方程决定系数分别为0.865、0.869和0.868。Logistic、Bertalanffy和Gompertz 3种模型中，雄麝胸围与周龄的拟合方程决定系数分别为0.913、0.925和0.918，而雌麝胸围与周龄的拟合方程决定系数分别为0.893、0.897和0.896。但Bertalanffy 和Gompertz模型中，均存在拐点周龄为0的情况。以拟合度判断，Bertalanffy模型在周龄与脊柱长和胸围回归关系中具有最好的拟合度。

2.5 肢体体尺与周龄的拟合

表5结果显示了林麝肢体体尺与周龄的3种拟合曲线。结果显示：周龄与肱骨长、大腿长和小腿长的拟合方程的决定系数分别为0.903～0.914、0.914～0.930和0.810～0.933。Logistic、Bertalanffy和Gompertz中，肢体体尺与周龄的拟合度分别为0.810～0.928、0.905～0.933和0.903～0.932。3种模型之间没有差异，但在肱骨和大腿与周龄的拟合中，Bertalanffy和Gompertz模型出现了拐点周龄为0的情况；在小腿与周龄拟合中，3种模型均出现了拐点周龄为0的情况。以拟合度判断，Bertalanffy模型在周龄与肱骨长、大腿长和小腿长回归关系中具有最好的拟合度。

表3 林麝头部体尺与周龄的3种拟合方程

表4 林麝躯体体尺与周龄的3种拟合方程

表5 林麝肢体体尺与周龄的3种拟合方程

3 讨论

生长曲线分析实质是把不同时期体质量、体尺资料提供的信息综合成少数几个参数。动物生长中常用Logistic、Bertalanffy和Gompertz模型来探讨动物的生长发育规律(Woolson & Leeper，1980；von-Rosen，1991；Nunez-Antonetal.，2001)。因此本研究中仅用了Logistic、Bertalanffy和Gompertz模型拟合林麝体质量与体尺的回归关系。增重是一个连续的过程，在动物的整个生产过程中，动物的生长主要表现为“S”曲线，但“S”曲线必须在非限制饲养条件下(竭航等，2014)。在本研究中，林麝体质量与周龄在Logistic模型中呈“S”曲线，在Bertalanffy和Gompertz模型中接近幂函数型，本研究中Bertalanffy模型具有最好的拟合度。研究结果反映出林麝前期增重非常迅速，这与已有的研究结果一致(张保良，1982；程世国，1984；竭航等，2014)。在已有的研究中，体质量与日龄在Gompertz、Logistic模型中的拟合度分别为0.988、0.980(竭航等，2014)，与本研究的拟合度相比更高，其原因主要是研究的个体数少，个体生长发育体质量变异较小；另一方面使用的Logistic模型与本研究中经典的Logistic模型有差异。

林麝生产发育过程中体尺生长模型的研究很少，竭航等(2014)研究了肩高(体高、体长)与日龄的拟合曲线。本研究尝试用Logistic、Bertalanffy和Gompertz模型来探讨林麝体尺发育规律。结果表明Bertalanffy模型的拟合度相对较高，而且肢体体尺(肱骨长、大腿长和小腿长)与周龄的拟合度较高，头部体尺(额宽、头长和耳长)与周龄的拟合度较低，躯体体尺(脊柱长和胸围)与周龄的拟合度介于二者之间。在动物的体尺测量上：体长和肩高与动物站姿有关，胸围测量比较麻烦，因此通过测量林麝肢体体尺和头部体尺来估测林麝周龄对于其管理非常方便。

程世国， 邹真慧. 1991. 麝的饲料和饲养[M]. 成都: 四川科技出版社: 4-32.

程世国. 1984. 麝的饲养[J]. 西南民族学院学报: 畜牧兽医版, 10(4): 45-49.

何迟酩, 林进生, 邓家波, 等. 2012. 林麝血液生理生化指标研究[J]. 四川动物, 31(3): 456-459.

姜海瑞, 徐宏发, 王晓黎, 等. 2007. 麝的生物学研究现状分析[J]. 生物学通报, 42(6): 4-6.

竭航, 封孝兰, 陈强, 等. 2014. 圈养林麝生长曲线拟合研究[J]. 四川农业大学学报, 32(1): 97-102.

李娟, 刘益平, 张明, 等. 2009. 大恒优质肉鸡早期体重发育规律[J]. 四川农业大学学报, 27(3): 345-349.

盛志廉, 吴常信. 1995. 数量遗传学[M]. 北京: 中国农业出版社: 114-115.

张保良. 1982. 麝的生长发育[J]. 野生动物, 4(4): 24-27.

周瑜, 赵克雷, 田永强, 等. 2013. 几种药物对圈养林麝化脓隐秘杆菌菌膜的影响[J]. 四川动物, 32(3): 330-333.

朱承嗣, 王方荣, 唐清山, 等. 2012. 家养林麝常见病诊断与防治研究[J]. 野生动物, 33(1): 3-4.

Meng XX, Zhao CJ, Hui CY,etal. 2011. Behavioral aspects of captive musk deer during non-mating season: gender differences and monthly pattern[J]. Asian-Australian Journal of Animal Science, 24(5): 707-712.

Meng XX, Zhou CQ, Hu JC,etal. 2006. Musk deer farming in China[J]. Animal Science, 82: 1-6.

Migon GS, Beaumont CL, Bihan DE,etal. 1990. Genetic parameters of growth curve parameter in male and female chicken[J]. British Poultry Science, 40(1): 44.

Nunez-Anton V, Zimmerman D, Gregoire TG,etal. 2001. Parametric modeling of growth curve data: an overview[J]. Test, 10(1): 1-73.

von-Rosen D. 1991. The growth curve model: a review[J]. Communications in Statistics-theory and Methods, 20(9): 2791-2822.

Woolson RF, leeper JD. 1980. Growth curve analysis of complete and incomplete logitudinal data[J]. Communications in Statistics-theory and Methods, 9(14): 1491-1513.

Xia S, Zou DF, Yue BS. 2006. Six microsatellite loci in forest musk deer,Moschusberezovskii[J]. Molecular Ecology Note, 6: 113-115.

Yang QS, Meng XX, Xia L,etal. 2003. Conservation status and causes of decline of musk deer (Moschusspp.) in China[J]. Biological Conservation, 109: 333-342.

Zhang SC, Yue BS, Zou FD. 2007. Isolation and characterization of microsatellite DNA marker from forest musk deer (Moschusberezovskii)[J]. Zoological Research, 18(1): 24-27.

Zhao SS, Chen X, Fang SG,etal. 2008. Development and characterization of 15 novel microsatellite markers from forest musk deer (Moschusberezovskii)[J]. Conservation Genetics, 9: 723-725.

Studies on Growth and Development Model of Captive Forest Musk Deer

CAI Youghua1, YANG Ying1, FU Wenlong1, ZHU Ping1, FENG Dayong1, LEI Yi2, ZHANG Ming3*

(1. Institute of Breeding and Reproduction of Forest Musk Deer in Sichuan, Dujiangyan, Sichuan Province 611845, China;2. College of Tourism, Sichuan Agricultural University, Dujiangyan, Sichuan Province 611830, China;3. College of Animal Science and Technology, Sichuan Agricultural University, Ya’an, Sichuan Province 625014, China)

In order to explore the rule of growth and development of body weight and size of caotive forest musk deer (Moschusberezovskii), body weight and size (forehead width, head length, ear length, spine length, chest circumference, humerus length, thigh length and leg length) of 320 forest musk deer (males: 160, females: 160) were measured from the birth to the 400thweek, and the regression of the measured values and week-age were fitted with Logistic, Bertalanffy and Gompertz model, respectively. The results showed that the three tested models, Logistic, Bertalanffy and Gompertz, could fit the regression of body weight and week-age of forest musk deer. The fitting degree (R2) of Bertalanffy model was better than that of Logistic and Gompertz for males (R2=0.966) and females (R2=0.954). The results of the regression body size and week-age revealed that: Bertalanffy model was better than Logistic and Gompertz models; The fitting degree of limb size and week-age was high, whereas that of head size and week-age was low. In all parameters of body size, the fitting degree of leg length and week-age was the highest in males (R2=0.927) or females (R2=0.933) with Bertalanffy model. The rule of growth and development of body weight and size of forest musk deer could be fitted with Bertalanffy model.

growth curve; forest musk deer; Logistic model; Bertalanffy model; Gompertz model

2015-08-11 接受日期：2015-10-29

四川省科技厅“‘林麝饲养和活体取香技术’成果产业化配套技术的研究及推广应用”项目

蔡永华(1962—)， 高级工程师

*通信作者Corresponding author， E-mail:zhm3000@163.com

Q959.8

A

1000-7083(2016)02-0238-05