Prediction of random-regression coefficient for daily milk yield after 305 days in milk by using the regression-coefficient estimates from the first 305 days

  • Yamazaki, Takeshi (Dairy Cattle Group, Division of Dairy Production Research, Hokkaido Agricultural Research Centre, NARO) ;
  • Takeda, Hisato (Animal Breeding Unit, Division of Animal Breeding and Reproduction Research, Institute of Livestock and Grassland Science, NARO) ;
  • Hagiya, Koichi (Department of Life and Food Science, Obihiro University of Agriculture and Veterinary Medicine) ;
  • Yamaguchi, Satoshi (Computing Section, Milk Recording Division, Hokkaido Dairy Milk Recording and Testing Association) ;
  • Sasaki, Osamu (Animal Breeding Unit, Division of Animal Breeding and Reproduction Research, Institute of Livestock and Grassland Science, NARO)
  • Received : 2017.11.21
  • Accepted : 2018.03.06
  • Published : 2018.10.01


Objective: Because lactation periods in dairy cows lengthen with increasing total milk production, it is important to predict individual productivities after 305 days in milk (DIM) to determine the optimal lactation period. We therefore examined whether the random regression (RR) coefficient from 306 to 450 DIM (M2) can be predicted from those during the first 305 DIM (M1) by using a RR model. Methods: We analyzed test-day milk records from 85,690 Holstein cows in their first lactations and 131,727 cows in their later (second to fifth) lactations. Data in M1 and M2 were analyzed separately by using different single-trait RR animal models. We then performed a multiple regression analysis of the RR coefficients of M2 on those of M1 during the first and later lactations. Results: The first-order Legendre polynomials were practical covariates of RR for the milk yields of M2. All RR coefficients for the additive genetic (AG) effect and the intercept for the permanent environmental (PE) effect of M2 had moderate to strong correlations with the intercept for the AG effect of M1. The coefficients of determination for multiple regression of the combined intercepts for the AG and PE effects of M2 on the coefficients for the AG effect of M1 were moderate to high. The daily milk yields of M2 predicted by using the RR coefficients for the AG effect of M1 were highly correlated with those obtained by using the coefficients of M2. Conclusion: Milk production after 305 DIM can be predicted by using the RR coefficient estimates of the AG effect during the first 305 DIM.


Supported by : NARO


  1. ICAR. Milk recording surveys on cow, sheep and goats [Internet]. ICAR; 2017 [cited 2017 Oct 24]. Available from:
  2. Inchaisri C, Jorritsma R, Vos PLAM, van der Weijden GC, Hogeveen H. Analysis of the economically optimal voluntary waiting period for first insemination. J Dairy Sci 2011;94:3811-23.
  3. Knight C. Extended lactation: turning theory into reality. Adv Dairy Technol 2005;17:113-23.
  4. VanRaden PM, Dematawewa CMB, Pearson RE, Tooker ME. Productive life including all lactations and longer lactations with diminishing credits. J Dairy Sci 2006;89:3213-20.
  5. Bohmanova J, Miglior F, Jamrozik J. Use of test-day records beyond three hundred five days for estimation of three hundred five-day breeding values for production traits and somatic cell score of Canadian Holsteins. J Dairy Sci 2009;92:5314-25.
  6. Hagiya K, Terawaki Y, Yamazaki T, et al. Relationships between conception rate in Holstein heifers and cows and milk yield at various stages of lactation. Animal 2013;7:1423-8.
  7. Haile-Mariam M, Pryce JE. Variances and correlations of milk production, fertility, longevity, and type traits over time in Australian Holstein cattle. J Dairy Sci 2015;98:7364-79.
  8. Interbull. Interbull guidelines for national and international genetic evaluation systems in dairy cattle with focus on production traits. Interbull Bulletin 2001;28:1-27.
  9. ICAR. ICAR rules, standards and guidelines for dairy production recording. In: International Agreement of Recording Practices, Section 2. Roma, Italy: International Committee on Animal Recording; 2016. p. 25-94.
  10. Haile-Mariam M, Goddard ME. Genetic and phenotypic parameters of lactations longer than 305 days (extended lactations). Animal 2008;2:325-35.
  11. Jamrozik J, Schaeffer LR, Dekkers JCM. Genetic evaluation of dairy cattle using test day yields and random regression model. J Dairy Sci 1997;80:1217-26.
  12. Jensen J. Genetic evaluation of dairy cattle using test-day models. J Dairy Sci 2001;84:2803-12.
  13. Interbull. National genetic evaluation forms provided by countries [Internet]. Interbull; 2017 [cited 2017 Oct 24]. Available from:
  14. Kistemaker GJ. The Canadian test day model using Legendre polynomials. Interbull Bulletin 2003;31:202-4.
  15. Wilmink JBM. Comparison of different methods of predicting 305-day milk yield using means calculated from within-herd lactation curves. Livest Prod Sci 1987;17:1-17.
  16. Schaeffer LR, Jamrozik J, Kistemaker GJ, Van Doormaal BJ. Experience with a test-day model. J Dairy Sci 2000;83:1135-44.
  17. Gernand E, Wassmuth R, von Borstel UU, Konig S. Heterogeneity of variance components for production traits in largescale dairy farms. Livest Sci 2007;112:78-89.
  18. Yazgan K, Makulska J, Weglarz A, Ptak E, Gierdziewicz M. Genetic relationship between milk dry matter and other milk traits in extended lactations of Polish Holstein cows. Czech J Anim Sci 2010;55:91-104.
  19. National Livestock Breeding Center. Genetic evaluation for production traits [Internet]. National Livestock Breeding Center; 2018 [cited 2018 Feb 4]. Available from:
  20. Madsen P, Jensen J. A user's guide to DMU [Internet]. University of Aarhus Research Centre Foulum; 2013 [cited 2017 Oct 24]. Available from: dmuv6_guide.5.2.pdf
  21. SAS Institute Inc. SAS/STAT 14.1 User's guide. Cary, NC, USA: SAS Institute Inc.; 2015.