Article| Volume 83, ISSUE 11, P2709-2719, November 2000

Download started.


Comparison of Models for Describing the Lactation Curve of Latxa Sheep and an Analysis of Factors Affecting Milk Yield

      This paper is only available as a PDF. To read, Please Download here.


      The objectives of this work were 1) to compare the goodness-of-fit of empirical models of the lactation curve and 2) to analyze the factors affecting the shapes of the lactation curves, the parameters describing them, and the overall milk yield of Latxa dairy sheep. A total of 14,699 records from 2711 ewes, collected during three consecutive years (1995 to 1997) by the milk recording program of the Latxa ewe of the Basque Country (Spain), were used. Six mathematical models and three fitting procedures were compared. The estimation of model parameters by nonlinear fitting procedures was superior to that by linear regression methods. A nonlinear variable decay model fitted the data better than the other models, as judged by lower mean square prediction error, residual sums of squares, and a lack of first-order positive autocorrelation as assessed by the Durbin Watson coefficient. The effects of the flock, flock-year interaction, month of lambing, length of lactation nested within month of lambing, parity, and number of live lambs born had significant effects on the parameters of the model and the total milk yield (P < 0.01). The prediction of milk yield from the selected model was similar to the estimates obtained with the Fleischmann method currently used by the national breeding program for the Latxa breed.

      Key words

      Abbreviation key:

      CDM (centering date method for calculation of total milk yield per lactation), DMY (daily milk yield), DW (Durbin Watson statistic), ETMYL (estimated total milk yield per lactation), MSPE (mean square prediction error), RES (residuals), RSS (residual sums of squares), TMYL (total milk yield per lactation)


        • Aguilar J.A.
        • Aparicio J.J.
        • Fernández-Rivera S.
        • Tovar I.
        Lactation curves andmilk yield of Corriedale, Rambouillet and Suffolk ewes.
        J. Dairy Sci. 1989; 72 (Abstr.): 426-427
        • Arranz J.
        • Oregui L.M.
        • Bravo M.V.
        • Urarte E.
        • Lana M.P.
        • Torrano L.
        Estudio de la duración delamamantamiento en ovejas de raza Latxa.
        Inform. Téc. Econ. Agraria. 1995; 16: 714-716
      1. Barillet, F. 1985. Amélioration génétique de la composition du lait de brebis: l’exemple de la race Lacaune. Ph.D. Diss., Inst. Natl. Rech. Agron., Paris-Grignon, France.

        • Barillet F.
        Expression de la production laitiére à la traite des brebis en systéme allaitement × traite mécanique.
        in: Eitam M. Proc. 4th Int. Symp. Machine Milking Small Ruminants. Tel-Aviv, Israel1989: 463-495
        • Barillet F.
        • Astruc J.M.
        • de Brauwer P.
        • Casu S.
        • Fabbri G.
        • Feddersen E.
        • Frangos K.
        • Gabiña D.
        • Gama L.T.
        • Ruiz Tena J.L.
        • Sanna S.R.
        Int. Committee for Animal Recording Guidelines: International Regulations for Milk Recording in Sheep.
        Institut d’Elevage, Paris, France1992 (15 pp.+appendix)
        • Barillet F.
        • Boichard D.
        Use of first lactation test-day data for genetic evaluation of the Lacaune dairy sheep.
        in: Proc. 5th World Congr, Genet. Appl. Livest. Prod., Guelph, ON, Canada1994: 111-114 (18)
        • Beever D.E.
        • Rook A.J.
        • France J.
        • Dhanoa M.S.
        • Gill M.
        A review of empirical and mechanistic model of lactational performance by the dairy cow.
        Livest. Prod. Sci. 1991; 29: 115-130
        • Brody S.
        • Ragsdale A.C.
        • Turner C.W.
        The rate of decline of milk secretion with the advance of the period of lactation.
        J. Gen. Physiol. 1923; 5: 441-444
        • Camacho J.
        • Baaijen M.
        • Perez E.
        Comparison of different models for the lactation curve in a tropical country.
        in: van Arendonk J.A.M. Proc. 45th Annu. Mtg. Eur. Assoc. Anim. Prod. Wageningen Pers, Wageningen, the Netherlands1994: 250
        • Cappio-Borlino A.
        • Pulina G.
        • Rossi G.
        A non-linear modification of Wood's equation fitted to lactation curves of Sardinian dairy ewes.
        Small Ruminant Res. 1995; 18: 75-79
        • Carta A.
        • Sanna S.R.
        • Casu S.
        Estimating lactation curves and seasonal effects for milk, fat and protein in Sarda dairy sheep with a test day model.
        Livest. Prod. Sci. 1995; 44: 37-44
        • Cobby J.M.
        • Le Du Y.L.P.
        On fitting curves to lactation data.
        Anim. Prod. 1978; 26: 127-133
        • Durbin J.
        • Watson G.S.
        Testing for serial correlation in least squares regression.
        Biometrika. 1951; 38: 159-178
        • Gabiña D.
        • Arrese F.
        • Arranz J.
        • Beltrán de I.
        • Heredia
        Average milk yields and environmental effects on Latxa sheep.
        J. Dairy Sci. 1993; 76: 1191-1198
        • Gipson T.A.
        • Grossman M.
        Lactation curves in dairy goats: a review.
        Small Ruminant Res. 1990; 3: 383-396
        • Grossman M.
        • Koops W.J.
        Multiphasic analysis of lactation curves in dairy cattle.
        J. Dairy Sci. 1988; 71: 1598-1608
        • María G.
        • Gabiña D.
        Simplification of milk recording scheme in Latxa milking sheep.
        Livest. Prod. Sci. 1992; 31: 313-320
        • María G.
        • Gabiña D.
        Non-genetic effects on milk production of Latxa ewes.
        Small Ruminant Res. 1993; 12: 61-67
        • Molina M.P.
        • Sanna A.
        • Molle G.
        • Branca A.
        • Ruda G.
        • Casu S.
        Evolution de l’etat corporel des brebis Sardes pendant la lactation et relation avec la productivité en conditions d’élevage intensif.
        in: Purroy A. Options Méditerranéennes Ser. A: Semin. Centre International de Hautes Etudes Agronomiques Méditerranéennes, Zaragoza, Spain1992: 77-84 (Méditerraneens No. 13)
        • Morant S.V.
        • Gnanasakthy A.
        A new approach to the mathematical formulation of lactation curves.
        Anim. Prod. 1989; 49: 151-162
        • Nelder J.A.
        Inverse polynomials, a useful group of multifactor response functions.
        Biometrics. 1966; 22: 128-141
      2. Oregui, L. M. 1992. Estudio del manejo de la alimentación en los rebaños ovinos de raza Latxa y su influencia sobre los resultados reproductivos y de producción lechera. Ph.D.Diss. Servicio. Publ. Gobierno Vasco, Vitoria-Gasteiz, Spain.

        • Papajsik I.A.
        • Bodero J.
        Modelling lactation curves of Friesian cows in a subtropical climate.
        Anim. Prod. 1988; 47: 201-207
        • Portolano F.
        • Spatafora F.
        • Bono G.
        • Margiota S.
        • Todaro M.
        • Ortoleva V.
        • Leto G.
        Application of the Wood model to lactation curves of Comisana sheep.
        Small Ruminant Res. 1996; 24: 7-13
        • Ruiz R.
        • Oregui L.M.
        • Bravo M.V.
        Productive characteristics of dairy ewe flocks in the Basque Country.
        in: Sørensen J.T. Livestock Farming Systems: More than Food Production. Wageningen Pers, Wageningen, The Netherlands1997: 42-47 (EAAP Publ. No. 89)
      3. SAS® Procedures Guide, Version 6, Third Edition. 1990. SAS Inst., Inc. Cary, NC.

        • Sakul H.
        • Boylan W.J.
        Lactation curves for several US sheep breeds.
        Anim. Prod. 1992; 54: 229-233
        • Shanks R.D.
        • Berger P.J.
        • Freeman A.E.
        • Dickinson F.N.
        Genetic aspects of lactation curves.
        J. Dairy Sci. 1981; 64: 1852-1860
        • Stanton T.L.
        • Jones L.R.
        • Everett R.V.
        • Kachman S.D.
        Estimating milk, fat and protein lactation curves with a test day model.
        J. Dairy Sci. 1992; 75: 1691-1700
        • Torres-Hernández G.
        • Hohenboken W.D.
        Biometric properties of lactations in ewes raising single or twin lambs.
        Anim. Prod. 1980; 30: 431-436
        • Urarte E.
        • Gabiña D.
        • Arranz J.
        • Arrese F.
        • Gorostiza P.
        • Sierra I.
        Las razas ovinas Latxa y Carranzana I. Sistemas de producción.
        Inform. Téc. Econ. Agraria. 1992; 84: 7-25
        • Wilmink J.B.M.
        Adjustment of test-day milk, fat and protein yield for age, season and stage of lactation.
        Livest. Prod. Sci. 1987; 16: 335-348
        • Wilson L.L.
        • Varela-Alvarez H.
        • Hess C.E.
        • Rugh M.C.
        Influence of energy level, creep feeding and lactation stage on ewe milk and lamb growth characters.
        J. Anim. Sci. 1971; 33: 686-690
        • Wood P.D.P.
        Algebraic model of the lactation curve in cattle.
        Nature (Lond.). 1967; 216: 164-165
        • Wood P.D.P.
        Factors affecting the shape of the lactation curve in cattle.
        Anim Prod. 1969; 11: 307-316