استفاده از رگرسیون کوآنتایل در بررسی برخی خصوصیات منحنی شیردهی گاوهای شیری ایران

نوع مقاله : مقاله پژوهشی

نویسندگان

1 بخش علوم دام، دانشکده کشاورزی، دانشگاه بیرجند

2 مؤسسه تحقیقات علوم دامی، سازمان تحقیقات، آموزش و ترویج کشاورزی

چکیده

هدف از پژوهش کنونی، بررسی ویژگی‌های  منحنی شیردهی گاوهای شیری ایران با استفاده از رگرسیون کوآنتایل (چندکی) بود. داده­های مورد استفاده شامل 15647867 رکورد روز آزمون شیر بود که از 984690 رأس گاو در 579 گله موجود در 26 استان کشور گردآوری شدند. ابتدا، رکوردهای روز آزمون شیر در ترکیب­های مختلف نوبت زایش–فصل زایش (و مشتمل بر 20 ترکیب) با استفاده از یک مدل خطی تجزیه و تحلیل شدند. سپس، در هر یک از ترکیبات مزبور، میانگین حداقل مربعات رکوردهای روز آزمون شیر در فاصله­ روزهای 4 تا 305 شیردهی، در تابع خطی شده گامای وود و با روش رگرسیون کوآنتایل، مورد استفاده قرار گرفتند. نتایج نشان داد با بالا رفتن چندک­های رکوردهای روز آزمون شیر، مقدار تولید شیر در اوج شیردهی، افزایش می­یابد. زمان رسیدن به اوج شیردهی و میزان تداوم شیردهی در چندک­های پایینی، بزرگ­تر از چندک­های بالایی رکوردهای روز آزمون شیر بودند. در تمامی چندک­ها، کمترین و بیشترین میانگین زمان رسیدن به اوج شیردهی (92/50 و 36/93 روز) و میزان تداوم شیردهی (64/6 و 37/7) به­ترتیب مربوط به گاوهایی بود که در پاییز زایش داشتند و در دوره­های شیردهی سوم یا اول بودند (05/0P<). در رابطه با مقدار تولید شیر در اوج شیردهی، گاوهای زایمان کرده در بهار که دوره­ اول شیردهی را شروع کرده بودند، کمترین میانگین (84/31 کیلوگرم) و گاوهای زایمان کرده در پاییز که دوره­ چهارم شیردهی را شروع کرده بودند، بیشترین میانگین (43/40 کیلوگرم) را دارا بودند (05/0P<). با توجه به نتایج تحقیق حاضر، چنان­چه قرار باشد در امر ارزیابی ژنتیکی گاوهای شیری ایران برای صفت تولید شیر از اطلاعات دوره­های مختلف شیردهی استفاده شود، ضروری خواهد بود تا تغییرات یاد شده در صفات شیردهی، مدنظر قرار گیرند، به­ویژه هنگامی که از مدل­های روز آزمون جهت تجزیه ژنتیکی منحنی شیردهی استفاده می­ شود.

کلیدواژه‌ها

موضوعات


عنوان مقاله [English]

The use of quantile regression to evaluate some lactation curve characteristics of Iranian dairy cows

نویسندگان [English]

  • Z. Hesami Zostami 1
  • S. H. Farhangfar 1
  • M. A. Abbasi 2
  • H. Naeemipour Younesi 1
1 Animal Science Department, Agriculture Faculty, University of Birjand, Birjand, Iran
2 Animal Science Research Institute of IRAN (ASRI), Agricultural Research, Education & Extension Organization (AREEO), Karaj, Iran
چکیده [English]

Introduction: Over recent decades, numerous mathematical models have been developed to represent the lactation curve of dairy animals. These models are categorized into parametric and non-parametric types, with the three-parameter Wood's gamma incomplete function being extensively utilized for dairy cows. By using the parameters from Wood's model, key lactation characteristics such as peak time, peak yield, and persistency can be determined. These features are considered economic traits, suitable for genetic selection. Previous studies have primarily focused on employing mathematical models to calculate these characteristics based on the average of test day milk records. However, they can also be computed for individual quantiles of test day milk records (TDMR). To our knowledge, no research has yet explored the use of quantile regression for modeling the lactation curve of Iranian dairy cows. Therefore, the primary objective of this study was to assess certain lactation curve characteristics of Iranian dairy cows using the quantile regression technique.
Materials and methods: The data for this study were sourced from the Animal Breeding Centre of Iran. Initially, the dataset was merged and refined using SPSS software. The final dataset included TDMR data from 984,690 dairy cows across 579 herds in 26 provinces nationwide. A total of 15,647,867 TDMR records were collected from the first five parities of cows calving between 1995 and 2021. The TDMR was initially analyzed using a fixed linear model tailored to 20 different combinations of parity and calving season. For each combination, the least squares means of TDMR, measured from 4 to 305 days in milk, were applied to Wood’s incomplete gamma function (y=atbe-ct) in a quantile-based approach. This parametric three-parameter model (a, b, c) is extensively used to describe lactation curves. The linearized form of Wood’s model was fitted to these means using the Quantreg procedure in SAS software, estimating parameters across quantiles from 5 to 95 in increments of 5. Based on these parameters, lactation characteristics such as peak time, peak yield, and persistency were calculated across all parity-calving season combinations. Statistical comparisons of the combination means were conducted using SPSS software.
Results and discussion: The results of this research indicate that peak time, peak yield, and the persistency of the lactation curve in Iranian dairy cows vary across different quantiles of TDMR, parities, and calving seasons. The findings reveal that peak yield increases with higher TDMR quantiles. Conversely, peak time and persistency are greater in the lower TDMR quantiles compared to the upper ones. Across the TDMR quantiles, the minimum (50.92 days) and maximum (93.36 days) peak times, as well as persistency values (6.64 and 7.37), were observed in cows that calved in autumn and began their first or third lactation period (P<0.05). Regarding peak yield, cows that calved in spring and started their first lactation had a lower mean yield (31.84 kg), while the highest mean yield (40.43 kg) was noted in cows that calved in autumn and began their fourth lactation period (P<0.05). As anticipated, cows in their first parity exhibited greater lactation persistency compared to others. Peak yield was significantly correlated with peak time (-0.9520) and persistency (-0.9516) (P<0.0001), suggesting that cows with lower milk yield at peak time tend to reach their peak later and maintain higher persistency post-peak.
Conclusions: The research findings reveal that the peak time, peak yield, and persistency of the lactation curve in Iranian dairy cows differ across the quantiles of test-day milk records. Specifically, peak time and persistency are higher in the lower quantiles compared to the upper ones. Additionally, these lactation characteristics vary across different parities and calving seasons. Based on these results, it is recommended that these variations be considered when using records from various parities for the genetic evaluation of Iranian dairy cows. This consideration is crucial when employing test-day models for genetic analysis of the lactation curve.

کلیدواژه‌ها [English]

  • Lactation characteristics
  • Milk production
  • Quantile regression
  • Wood function
  • Dairy cow
Angeles-Hernandez, J. V., Aranda-Aguirre, E., Muñoz-Benítez, A. L., Chay-Canul, A. J., Albarran-Portillo, B., Pollott, G. E., & Gonzalez-Ronquillo, M. (2021). Physiology of milk production and modelling of the lactation curve. CAB Reviews, 60(056), 1-22. doi: 10.1079/PAVSNNR202116056
Atashi, H., Zamiri, M. J., Akhlaghi, A., Dadpasand, M., Sayyadnejad, M. B., & Abdolmohammadi, A. R. (2013). Association between the lactation curve shape and calving interval in Holstein dairy cows of Iran. Iranian Journal of Veterinary Research, 14(2), 88-93. doi: 10.22099/IJVR.2013.1580
Biswas, J., Kulkarni, H., & Das, K. (2017). Quantile regression in biostatistics. Biostatistics and Biometrics, 2(5), 102-105. doi: 10.19080/BBOAJ.2017.02.555596
Bouallegue, M., & M’Hamdi, N. (2020). Mathematical modeling of lactation curves: a review of parametric models. In: Lactation in Farm Animals - Biology, Physiological Basis, Nutritional Requirements, and Modelization, Edited by Naceur M’Hamdi. Part of IntechOpen Book Series: Veterinary Medicine and Science, 3, 95-114. doi: 10.5772/intechopen.78900
Boujenane, I., & Hilal, B. (2012). Genetic and non genetic effects for lactation curve traits in Holstein-Friesian cows. Archiv Tierzucht, 55(5), 450-457. doi: 10.5194/aab-55-450-2012
Capuco, A. V., Wood, D. L., Baldwin, R., Mcleod, K., & Paape, M. J. (2001). Mammary cell number, proliferation, and apoptosis during a bovine lactation: relation to milk production and effect of bST. Journal of Dairy Science, 84(10), 2177-2187. doi: 10.3168/jds.S0022-0302(01)74664-4
Castillo-Gallegos, E. (2018). Characterization of the lactation curve of f1 Holstein-Zebu cows. Ecosistemas y Recursos Agropecuarios, 5(14), 335-343. doi: 10.19136/era.a5n14.1524
Chen, Y. (2024). Lactation curve modelling in dairy production: applications at cow and herd level. Ph.D. thesis, Faculty of Veterinary Medicine, Utrecht University, the Netherlands. doi: 10.33540/2148
Cobby, J., & Le Du, Y. L. P. (1978). On fitting curves to lactation data. Animal Production, 26(2), 127-133. doi: 10.1017/S0003356100039532
Dado, R. G., & Allen, M. S. (1994). Variation in and relationships among feeding, chewing, and drinking variables for lactating dairy cows. Journal of Dairy Science, 77(1), 132-144. doi: 10.3168/jds.S0022-0302(94)76936-8
Dekkers, J. C. M., Ten Hag, J. H., & Weersink, A. (1998). Economic aspects of persistency of lactation in dairy cattle. Livestock Production Science, 53(3), 237-252. doi: 10.1016/S0301-6226(97)00124-3
Dijkstra, J., López, S., Bannink, A., Dhanoa, M. S., Kebreab, E., Odongo, N. E., Fathi Nasri, M. H., Behera, U. K., Hernandez-Ferrer, D., & France, J. (2010). Evaluation of a mechanistic lactation model using cow, goat and sheep data. Journal of Agricultural Science, 148(3), 249-262. doi: 10.1017/S0021859609990578
Elahi Torshizi, M. (2016). Effects of season and age at first calving on genetic and phenotypic characteristics of lactation curve parameters in Holstein cows. Journal of Animal Science and Technology, 58(8), 1-14.doi: 10.1186/s40781-016-0089-1
Farhangfar, S. H., Jafari, M., Montazar Torbati, M. B., & Sayyadnezhad, M. B. (2021). Genetic analysis for the trait of first lactation milk production and evaluation of the impact of selection for it on the phenotypic value of age at first calving of Iranian Holsteins using quantile regression. Animal Production Research, 10(4), 61-72. doi: 10.22124/AR.2022.18580.1585 [In Persian]
Farhangfar, S. H., Namjou, M., & Eghbal, A. R. (2022). Estimation of the impact of somatic cell count on daily milk yield in different lactation stages of Iranian dairy cows using quantile regression. Breeding and Improvement of Livestock Journal, 2(2), 43-52. doi: 10.22034/BILJ.2022.153230 [In Persian]
Farhangfar, S. H., Rashidi Toghroljerdi, M. S., Montazer Torbati, M. B., & Sayyad Nezhad, M. B. (2023). A study on the lactation curve characteristics of grade and Iranian purebred Holstein cows with the use of raw, fat corrected, and energy-corrected milk test day records. Animal Production Research, 12(2), 71-84. doi: 10.22124/AR.2023.22771.1718 [In Persian]
George, D., & Mallery, P. (2020). IBM SPSS Statistics 26 Step by Step: A Simple Guide and Reference. 16th Edition. Taylor & Francis.
Ghavi Hossein-Zadeh, N. (2014). Comparison of non-linear models to describe the lactation curves of milk yield and composition in Iranian Holsteins. Journal of Agricultural Science, 152(2), 309-324. doi: 10.1017/S0021859613000415
Ghavi Hossein-Zadeh, N. (2017). Application of growth models to describe the lactation curves for test-day milk production in Holstein cows. Journal of Applied Animal Research, 45(1), 145-151. doi: 10.1080/09712119.2015.1124336
Ghavi Hossein-Zadeh, N. (2019). Comparison of the parameters of the lactation curve between normal and difficult calvings in Iranian Holstein cows. Spanish Journal of Agricultural Research, 17(1), 1-13. doi: 10.5424/sjar/2019171-13673
Grossman, M., & Koops, W. J. (1988). Multiphasic analysis of lactation curves in dairy cattle. Journal of Dairy Science, 71(6), 1598-1608. doi: 10.3168/jds.S0022-0302(88)79723-4
Ihsanullah, I., Qureshi, M. S., Akhtar, S., & Suhail, S. M. (2020). Seasonal stress affects reproductive and lactation traits in dairy cattle with various levels of exotic blood and parities under subtropical condition. Pakistan Journal of Zoology, 52(1), 147-155. doi: 10.17582/journal.pjz/2020.52.1.147.155
Jianga, H., Hicksona, R. E., Woodsa, O. T., Morandeaub, M., Burkea, J. L., Correa-Lunaa, M., Donaghya, D. J., & Lopez-Villalobos, N. (2020). Persistency and lactation curves modelled using nonlinear random regression in dairy cows milked once a day. New Zealand Journal of Animal Science and Production, 80, 131-136.
Koenker, R., & Bassett, G. (1978). Regression quantiles. Econometrica, 46(1), 33–50. doi: 10.2307/1913643
Madsen, O. (1975). A comparison of some suggested measures of persistency of milk yield in dairy cows. Animal Production, 20(2), 191-197. doi: 10.1017/S0003356100035182
Koloi, S., Pathak, K., Behera, R., Mandal, D. K., Karunakaran, M., Dutta, T. K., & Mandal, A. (2018). Factors affecting the persistency of milk production in Jersey crossbred cattle. Journal of Dairy, Veterinary and Animal Research, 7(6), 268-271. doi: 10.15406/jdvar.2018.07.00225
Kopec, T., Chládek, G., Kučera, J., Falta, D., Hanuš, O., & Roubal, P. (2013). The effect of the calving season on the Wood’s model parameters and characteristics of the lactation curve in Czech Fleckvieh cows. Archiv Tierzucht, 56(80), 808-815. doi: 10.7482/0003-9438-56-080
Li, M., Rosa, G. J. M., Reed, K. F., & Cabrera, V. E. (2022). Investigating the effect of temporal, geographic, and management factors on US Holstein lactation curve parameters. Journal of Dairy Science, 105(9), 7525-7538. doi: 10.3168/jds.2022-21882
Marumo, J. L., Lusseau, D., Speakman, J. R., Mackie, M., & Hambly, C. (2022). Influence of environmental factors and parity on milk yield dynamics in barn-housed dairy cattle. Journal of Dairy Science, 105(2), 1225-1241. doi: 10.3168/jds.2021-20698
Mellado, M., Antonio-Chirino, E., Meza-Herrera, C., Veliz, F. G., Arevalo, J. R., Mellado, J., & de Santiago, A. (2011). Effect of lactation number, year, and season of initiation of lactation on milk yield of cows hormonally induced into lactation and treated with recombinant bovine somatotropin. Journal of Dairy Science, 94(9), 4524-4530. doi: 10.3168/jds.2011-4152
Montiel-Olguín, L. J., Jesús Ruiz-López, F. D., Gómez-Rosales, S., Cantó-Alarcón, G. J., Estrada-Cortés, E., & Vera-Ávila, H. R. (2025). Comparison of different models to describe the lactation curve of Holstein cows in a small-scale dairy system. Emirates Journal of Food and Agriculture, 36, 1-7. doi: 10.3897/ejfa.2025.129407
Naeemipour Younesi, H., Shariati, M. M., Zerehdaran, S., Jabbari Nooghabi, M., & Løvendahl, P. (2019). Using quantile regression for fitting lactation curve in dairy cows. Journal of Dairy Research, 86(1), 19-24. doi: 10.1017/S0022029919000013
Namjou, M., & Farhangfar, S. H. (2022). A study on the curve of far percentage, protein percentage and fat to protein ratio in test day milk records of Iranian dairy cows using quantile regression statistical technique. Journal of Animal Science, 32(4), 15-29. doi: 10.22034/AS.2022.29102.1449 [In Persian]
Pedrosa, V. B., Schenkel, F. S., Chen, S. Y., Oliveira, H. R., Casey, T. M., Melka, M. G., & Brito, L. F. (2021). Genome wide association analyses of lactation persistency and milk production traits in Holstein cattle based on imputed whole-genome sequence data. Genes, 12(11), 1-28. doi: 10.3390/genes12111830
Radjabalizadeh, K., Alijani, S., Gorbani, A., & Farahvash, T. (2022). Estimation of genetic parameters of Wood’s lactation curve parameters using Bayesian and REML methods for milk production trait of Holstein dairy cattle. Journal of Applied Animal Research, 50(1), 363-368. doi: 10.1080/09712119.2022.2080211
Rezaei, R, Wu, Z, Hou, Y, Bazer, F. W., & Wu, G. (2016). Amino acids and mammary gland development: nutritional implications for milk production and neonatal growth. Journal of Animal Science and Biotechnology, 7(20), 1-22. doi: 10.0.4.162/ s40104-016-0078-8
Ruban, S., Danshyn, V., Matvieiev, M., Borshch, O. O., Borshch, O. V., & Korol-Bezpala, L. (2022). Characteristics of lactation curve and reproduction in dairy cattle. Acta Universitatis Agriculturae et Silviculturae Mendelianae Brunensisdoi, 70(6), 373-381. doi: 10.11118/actaun.2022.028
SAS Institute. (2013). Base SAS® 9.4 Procedures Guide Statistical Procedures. Second Edition. SAS Institute Inc., Cary, NC, USA.
Sauvant, D. (1988). Modelling of lactation and nutrition. In: Modelling of livestock production systems. Eds: Korver S., & J. A. M. Van Arendonk. Kluwer Academic Publishers, Dordrecht, The Netherlands. Pp. 149-155.
Steri, R. (2009). The mathematical description of the lactation curve of ruminants: issues and perspectives. University of Sassari, Sassari, Italy.
Tekerli, M., Akinci, Z., Dogan, I., & Akcan, A. (2000). Factors affecting the shape of lactation curves of Holstein cows from the Balikesir province of Turkey. Journal of Dairy Science, 83(6), 1381-1386. doi: 10.3168/jds.S0022-0302(00)75006-5
Vahedi Darmian, R., Farhangfar, S. H., & Sayyadnezhad, M. B. (2020). Evaluation of the impact of inbreeding coefficient on milk production, lactation length and first calving age of Iranian dairy cows using quantile regression. Animal Science Journal, 30(2), 25-39. doi: 10.22034/AS.2020.11452 [In Persian]
Vijayakumar, M., Park, J. H., Ki, K. S., Lim, D. H., Kim, S. B., Park, S. M., Jeong, H. Y., Park, B. Y., & Kim, T. I. (2017). The effect of lactation number, stage, length, and milking frequency on milk yield in Korean Holstein dairy cows using automatic milking system. Asian-Australasian Journal of Animal Sciences, 30(8), 1093-1098. doi: 10.5713/ajas.16.0882
Wahinya, P. K., Jeyaruban, M. G., Swan, A. A., Gilmour, A. R., & Magothe, T. M. (2020). Genetic parameters for test-day milk yield, lactation persistency, and fertility in low-, medium-, and high-production systems in Kenya. Journal of Dairy Science, 103(11), 10399-10413. doi: 10.3168/jds.2020-18350
Walsh, S. W., Williams, E. J., & Evans, A. C. O. (2011). A review of the causes of poor fertility in high milk producing dairy cows. Animal Reproduction Science, 123(3-4), 127-138. doi: 10.1016/j.anireprosci.2010.12.001
Wathes, D. C., Cheng, Z., Bourne, N., Taylor, V. J., Coffey, M. P., & Brotherstone, S. (2007). Differences between primiparous and multiparous dairy cows in the inter-relationships between metabolic traits, milk yield and body condition score in the periparturient period. Domestic Animal Endocrinology, 33(2), 203-225. doi: 10.1016/j.domaniend.2006.05.004
Webster, H. H., Lengi, A. J., & Corl, B. A. (2024). Mammary epithelial cell exfoliation increases as milk yield declines, lactation progresses, and parity increases. Journal of Dairy Science Communications, 5(2), 707-712. doi: 10.3168/jdsc.2023-0534
Wood, P. D. P. (1967). Algebraic model of the lactation curve in cattle. Nature, 216, 164-165. doi: 10.1038/216164A0