Research

Influence of milking frequency on genetic parameters associated with the milk production in the first and second lactations of Iranian Holstein dairy cows using random regression test day models

Moslem Moghbeli Damane1, Masood Asadi Fozi2, Ahmad Ayatollahi Mehrgardi2
Author Information & Copyright
1Department of Animal Science, Faculty of Agriculture, University of Jiroft, Jiroft, Iran
2Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, Iran

© Damane et al. 2016. Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Received: Nov 3, 2014 ; Accepted: Jan 6, 2016

Published Online: Feb 1, 2016

Abstract

Background

The milk yield can be affected by the frequency of milking per day, in dairy cows. Previous studies have shown that the milk yield is increased by 6–25 % per lactation when the milking frequency is increased from 2 to 3 times per day while the somatic cell count is decreased. To investigate the effect of milking frequency (3X vs. 4X) on milk yield and it’s genetic parameters in the first and second lactations of the Iranian Holstein dairy cows, a total of 142,604 test day (TD) records of milk yield were measured on 20,762 cows.

Results

Heritability estimates of milk yield were 0.25 and 0.19 for 3X milking frequency and 0.34 and 0.26 for 4X milking frequency throughout the first and second lactations, respectively. Repeatability estimates of milk yield were 0.70 and 0.71 for 3X milking frequency and 0.76 and 0.77 for 4X milking frequency, respectively. In comparison with 3X milking frequency, the milk yield of the first and second lactations was increased by 11.6 and 12.2 %, respectively when 4X was used (p < 0.01).

Conclusions

Results of this research demonstrated that increasing milking frequency led to an increase in heritability and repeatability of milk yield. The current investigation provided clear evidences for the benefits of using 4X milking frequency instead of 3X in Iranian Holstein dairy cows.

Keywords: Genetic parameters; Random regression; Milk production; Milking frequency; Iranian Holstein cows

Background

The purpose of animal breeding is to genetically enhance the livestock production where more efficient animals are produced to guard against future circumstances. An important method for maximizing response to selection program is the accurate prediction of breeding values of animal [34]. In dairy cows, to implement an efficient breeding program, the estimates of genetic parameters for the production traits are required. In addition, to predict expected selection response and to achieve the predicted breeding value by the mixed model procedures, the accurate heritability, repeatability and correlation estimates are required. Various test day (TD) models have been recommended for genetic evaluation of dairy cows such as multiple-trait models, covariance function models and random regression models [33]. However, in many countries, the random regression model (RRM) has been widely demonstrated to increase the accuracy of breeding value predictions [31]. The use of RRM makes it possible to study changes in TD records over time and a better understanding of lactation genetics [32].

It has been found that the milk yield could be affected by the frequency of milking cows per day [16]. Erdman [14] demonstrated that the increase in milk yield due to the increased milk frequency is fixed, and is not dependent to the level of milk production at the time of increase milking frequency. In high producing dairy cows, twice a day milking interval of 10–14 h and 12–14 h are suggested, whereas, once a day or skipping milking is not acceptable, since it may arise several problems such as an increased in the somatic cell counts as well as poor udder health [30]. By increasing milking frequency from 2 to 3 times per day, an increase in milk yield from 6 to 25 % per lactation is observed [11, 14, 22] and there is also a decrease in somatic cell counts [17]. Moreover, Armstrong [1] reported that the milk production is increased by 8–12 % when four times (4X) a day milking is applied instead of three times (3X) a day. Three times-a-day milking is currently the most frequently used milking schedule in Iranian Holstein dairy cows, however, four times a day are also rarely applied. The objectives of this study were to estimate variance components for test day (TD) milk yield using a random regression TD model to evaluate the effect of milking frequency (3X vs. 4X) on genetic parameters associated with milk yield in the first and second lactations of Iranian Holstein dairy cows and determines how much milk yield changes with increasing in milking frequency in different lactation yield. It is expected that increasing milking frequency from 3X to 4X causes augmentation of genetic parameters of milk production and also increase in milk production, thereby would be more profitable for farmers.

Methods

Data

A total of 142,604 TD records of milk yield were measured on 20,762 first and second lactations of Iranian Holstein dairy cows from 2009 to 2011. The cows were from 113 herds located in the different regions of Iran. In average, for each cow 6 records were available along the lactation period. The records of cows were deleted if they had fewer than 4 or more than 10 TD records or their ages at first calving were below 20 months or above 36 months, further edits excluded irregular data for daily milk yield (<1.0 and >70 kg). Number of times milked per day with 2 levels (3X or 4X). In both 3X and 4X, the cows were milked for 305-d period. Very small number of records were available before day 5 and after day 305 which were discarded. A summary of the pedigree information and the data used are presented in Tables 1 and 2, respectively. The residual variance was considered to be heterogeneous and modeled by grouping days in milk into 10 equal segments of 30 days containing similar variations.

Table 1. Summary of the pedigree information of milk yield data in 3X and 4X milking frequency
Number of
TraitMilking frequencyAnimalRecordsSireDamHerd
Milk yield318,764117,168100216,458112
428102543629621988
Download Excel Table
Table 2. Descriptive statistics of the test-day records of milk yield (kg) for some selected days in milk with 3X and 4X milking frequency
3X Milking4X Milking
Lact 1Lact 2Lact 1Lact 2
DIMNMeanSDNMeanSDNMeanSDNMeanSD
5–35639229.74.97691386.3114232.55131341.86.1
36–65702233.44.9821840.76.4129837.34.8137446.16.2
66–95709733.95.1824640.26.6129238.64.9137545.86.2
96–125743033.45877838.66.61336384.8146743.86.2
126–155681632.85.3781536.46.7133037.24.8144841.36.1
156–185629332.45.4705734.46.7137536.44.8139438.85.9
186–225523431.85.6545632.56.7133635.44.7126036.66.3
226–2454287315.6413830.96.8132434.44.6121734.26.3
246–2753105305.7270229.16.71318394.5104531.86.2
276–305193829.25.8145326.76.1104531.94.874729.86
Download Excel Table

Model

A single trait random regression model was applied to estimate genetic parameters of test day (TD) records of milk yield in the first and second lactations. The fixed effects of herd-test date, calving age and the number of days in milk were fitted in the model. Also the additive genetic effects and animal permanent environmental effects were added to the model analysis as the random effects. The data were analysed using ASReml software [15].

The genetic analysis of the milk yields were done using the following model:

yijnos=HTDi+n=1kbnageijnn+n=1kcndimijnon+n=1kα1αpnφndimijno+n=1kPE1γpnφndimijno+eijnos

Where: yijnos is the individual test-day records HTDi is ith Herd-Test date, bn is fixed regression coefficient of age at calving, ageijn is calving age, cn is fixed regression coefficient of days in milk, dimijno is days in milk, φn is nth Legender polynomial for days in milk, αpn is additive genetic effects, γpn is animal permanent environmental effects and eijnos is the residuals.

Models with different order of Legendre polynomials for the additive genetic effects and the animal permanent environmental effects were compared using Schwarz’s Bayesian information criteria (BIC) [29]:

BIC=2l+Klogn

Where l is the Log likelihood values, K is the number of estimated parameters and n is the number of observations.

Results and discussion

The best models for the genetic analyses were selected based on the BIC. Accordingly, in the selected models, the order of fit for genetic and permanent environmental effects, the number of estimated parameters, log likelihood values and BIC are presented in bold Tables 3 and 4.

Table 3. Number of estimated parameters, log likelihood values and Schwarz’s Bayesian information criteria (BIC) for 3X milking frequency in the first and second lactations
LactationKaaKpebNumber of parametersLog likelyhoodBIC
Lact 11116−106,018.02212,111.7
2119−105,796.05211,682
2222−105,772.28211,648.6
3225−105,735.65211,589.5
3328−105,703.71211,539.9
4228−105,680.26211,493
4331−105,656.98211,460.6
4434−105,645.56211,451.9
5231−105,661.18211,469
5334−105,648.43211,457.7
5437−105,635.52211,446c
5540−105,675.21211,539.6
Lact 21116−131,871.6263,819.6
2119−131,569.2263,229.1
2222−131,516.12263,137.3
3225−121,394.82242,909c
3328−141,350.47282,834.7
4228−131,347.03262,827.8
4331−131,303.31262,754.7
4434−131,275.85262,714.1
5231−131,134.43262,416.9
5334−131,086.81262,336
5437−131,050.94262,278.6
5540−131,041.94262,274.9

aka: order of fit for additive genetic effect

bke: order of fit for permanent environmental effect

cSelected order of legendre polynomials for the genetic and permanent environmental effects (best model)

Download Excel Table
Table 4. Number of estimated parameters, log likelihood values and Schwarz’s Bayesian information criteria (BIC) for 4X milking frequency in the first and second lactations
LactationKaaKpebNumber of parametersLog likelyhoodBIC
Lact 11116−24,329.1348,723.88
2119−24,234.3948,546.7
2222−24,226.5948,543.4
3225−24,174.0748,450.67
3328−24,171.0848,456.99
4228−24,173.3348,461.49
4331−24,169.1848,465.49
4434−24,167.3648,474.16
5231−24,112.2648,351.65c
5334−24,108.0948,355.62
5437−24,100.9148,353.56
5540−24,118.7148,401.46
Lact 21116−27,021.3654,108.25
2119−26,923.0953,923.99
2222−26,913.5753,917.24
3225−26,859.4153,821.2
3328−26,849.1553,812.97
4228−26,857.3753,829.41
4331−26,849.1753,825.3
4434−26,845.6753,830.58
5231−26,794.9353,716.82
5334−26,787.0953,713.42c
5437−26,784.1553,719.83
5540−26,791.1553,746.11

aka: order of fit for additive genetic effect

bke: order of fit for permanent environmental effect

cSelected order of legendre polynomials for the genetic and permanent environmental effects (best model)

Download Excel Table

Phenotypic analysis

The average milk yield in the first lactation were 31.8 and 35.5 kg for 3X and 4X milking frequency, respectively. For the second lactations, the values were 34.7 and 39 kg for 3X and 4X milking frequency, respectively. In comparison with 3X milking frequency, 11.6 and 12.2 % greater milk were obtained during the first and second lactations, respectively when 4X was applied (p < 0.01). This indicated that the 4X milking frequency is more effective during the second lactation than that of the first lactation. Figures 1 and 2 present the effect of milking frequency on milk yield for the first two lactations in different DIM, which is similar to the normal pattern of the milk production in dairy cows.

jast-58-0-5-g1
Fig. 1. Trend of milk yield by days in milk (DIM) for three and four times milking in first lactation
Download Original Figure
jast-58-0-5-g2
Fig. 2. Trend of milk yield by days in milk (DIM) for three and four times milking in second lactation
Download Original Figure

For both the first and second lactations, similar patterns were observed in which the greatest milk yield was obtained from day 65 to 155 and thereafter, gradually reduced until the end of lactation period. In addition, the greatest difference (14 %) between milk yield with 3X and 4X milking frequency was obtained at day 95 and 125, for the first and second lactations, respectively. The lowest difference (9 %) between the values for milk yield with 3X and 4X milking frequency was obtained at day 275 and 305, respectively. Similar results were reported by Armstrong [1].

Estimation of genetic parameters

A range of 0.09 to o.55 has been reported for the heritability of milk yield [5, 10, 27, 34]. Figures 3, 4 and Table 5 show milk yields heritability estimates by DIM for 3X and 4X milking frequency during the first two lactations. The overall heritability estimates for 3X and 4X during the first lactation were 0.25 and 0.19, respectively. The values were 0.34 and 0.26 for the second lactation. The results were in agreement with those reported by Muir et al. [25]. The heritability estimates were slightly higher than the results derived from a multivariate model using the same dataset [28]. The results show the estimated heritabilities of both milking frequency in first lactation were lower than those for the second lactation which is in consistent with the results reported by Miglior et al. [24]. During the lactation periods, the estimated heritability of 4X was higher than 3X (excepted d 305 in first lactation) but the differences between the heritability of 3X and 4X in the first lactation were higher than those in the second lactation (Figs. 3 and 4). The changes in heritabilities are because of the changes in the variance components which are explained in the following paragraphs.

jast-58-0-5-g3
Fig. 3. Trajectory of estimated heritability for milk yield by days in milk (DIM) for three and four times milking in first lactation
Download Original Figure
jast-58-0-5-g4
Fig. 4. Trajectory of estimated heritability for milk yield by days in milk (DIM) for three and four times milking in second lactation
Download Original Figure
Table 5. Additive genetic (σ2a), animal permanent environmental (σ2pe), residual (σ2e) variances and heritability (h2) of the test-day records of milk yield (Kg) for some selected days in milk and different milking frequency
VariancesDays in milk (DIM)
5–3536–6566–9596–125126–155156–185186–215216–245246–275276–305
Lact 1 (3X)
σ2a5.05.74.95.26.18.19.18.08.521.7
σ2pe9.211.011.711.311.913.415.217.820.018.1
σ2e11.410.09.18.59.48.58.99.27.57.7
h20.200.210.190.210.220.270.270.230.240.46
Lact 2 (3X)
σ2a4.04.16.36.87.39.914.016.414.712.3
σ2pe25.521.722.324.225.625.725.125.429.340.9
σ2e20.919.015.714.313.412.712.511.511.210.1
h20.080.090.140.150.160.200.270.310.270.19
Lact 1 (4X)
σ2a10.18.912.011.611.511.211.09.59.113.5
σ2pe20.314.311.610.710.610.911.412.715.621.6
σ2e8.99.58.68.38.57.37.86.06.56.0
h20.260.270.370.380.380.380.360.340.290.33
Lact 2 (4X)
σ2a16.612.010.410.08.710.214.416.619.215.5
σ2pe18.525.027.226.926.125.825.825.927.637.2
σ2e21.616.415.211.612.311.48.98.07.96.7
h20.290.230.200.210.180.220.290.330.350.26
Download Excel Table

Variances tended to be larger at the beginning and the end of lactation, which is probably because of the smaller number of records that corresponded to these time periods, and are possible artifacts of Legendre polynomials.

For the first lactation, additive genetic variances (VA) of 3X increased slowly and for and 4X decrease slowly during the lactation trajectory (Fig. 5). Additive genetic variance of 3X in the second lactation was decrease sharply from day 5 to 185 and then increased slowly up to 275 days. Additive genetic variance of 3X in the second lactation was increase sharply from day 5 to 245 and then decrease to end of lactation (Fig. 6).

jast-58-0-5-g5
Fig. 5. Trajectory of additive genetic variance for milk yield by days in milk (DIM) for three and four times milking in first lactation
Download Original Figure
jast-58-0-5-g6
Fig. 6. Trajectory of additive genetic variance for milk yield by days in milk (DIM) for three and four times milking in second lactation
Download Original Figure

As indicated in Table 5, for all traits, animal permanent environmental variances increased slowly during the lactation period (Figs. 7 and 8). Based on the results of Meyer et al. [23] and Zavandilova et al. [35] on Holstein cows, the highest additive genetic and permanent environment variances for dairy traits occurred in the first and last days of lactation.

jast-58-0-5-g7
Fig. 7. Trajectory of animal permanent environmental variance for milk yield by days in milk (DIM) for three and four times milking in first lactation
Download Original Figure
jast-58-0-5-g8
Fig. 8. Trajectory of animal permanent environmental variance for milk yield by days in milk (DIM) for three and four times milking in second lactation
Download Original Figure

In general, the trends in the VA and VEP throughout the lactations obtained in this study are comparable to trends found by Druet et al. [12] and Strabel et al. [35]. For all traits, except for 3X milking frequency in the first lactation, residual variance was smaller, compared with the total phenotypic variance. This indicates a good fit of the model. Since residual variance related to additive genetic variances was larger for 3X milking frequency in the first lactation, therefore, the model should be modified, so that the genetic parameters could be estimated with a greater accuracy (Figs. 9 and 10). Residual variances decreased slowly during the lactation for all traits; therefore, heritabilities increased slowly during the lactation. These results are similar to those observed by Cobuci et al. [8] and Biassus et al. [2].

jast-58-0-5-g9
Fig. 9. Trajectory of residual variance for milk yield by days in milk (DIM) for three and four times milking in first lactation
Download Original Figure
jast-58-0-5-g10
Fig. 10. Trajectory of residual variance for milk yield by days in milk (DIM) for three and four times milking in second lactation
Download Original Figure

Table 6 shows the genetic and phenotypic correlations for milk yield along DIM for the different milking frequency in first two lactations. The genetic correlations for 3X cow ranged from −0.18 to 0.97 for the first lactation, from 0.17 to 0.99 for the second lactation and for 4X cow from −0.71 to 0.97 and from −0.22 to 0.97 for first and second lactation respectively. In general, genetic correlations declined with an increasing interval in DIM for both lactations. Similar results have been reported by El Faro et al. [13] for Caracu cattle. Karacao et al. [20] reported genetic correlations for milk yield at calving to midlactation monotonically decreased to 0.40 and at the end of the lactation, decreased to −0.06 also genetic correlation between d 1 and 247 was estimated to be 0.03.

Table 6. Estimates of genetic (below diagonal) and phenotypic (above diagonal) correlations between milk yield on selected days in milk (DIM) for different milking frequency in the first and second lactations
MFa Lactation and DIMDIM
356595125155185215245275305
3X First
 35-0.340.290.270.230.190.170.180.190.10
 650.19-0.590.510.460.470.470.430.380.40
 950.140.87-0.620.570.550.520.500.460.41
 1250.170.670.94-0.640.600.560.530.500.42
 1550.070.740.950.97-0.660.620.560.510.49
 185−0.090.860.940.870.95-0.700.630.550.59
 225−0.180.890.940.810.880.97-0.690.630.63
 245−0.140.760.920.820.800.830.91-0.730.62
 2750.030.540.800.810.720.650.710.92-0.63
 3050.020.760.720.640.750.840.800.630.53-
3X Second
 35-0.510.430.380.350.330.310.290.270.22
 650.52-0.580.550.510.440.380.320.270.23
 950.290.96-0.650.610.550.460.390.320.26
 1250.320.900.97-0.680.630.560.480.400.30
 1550.450.730.790.92-0.700.660.590.500.36
 1850.530.470.500.700.92-0.730.690.600.44
 2250.540.270.280.500.800.97-0.750.690.52
 2450.510.180.180.410.730.930.99-0.760.63
 2750.430.200.230.440.730.910.960.98-0.75
 3050.170.370.460.600.720.740.710.730.85-
4X First
 35-0.260.190.240.250.200.140.140.170.04
 65−0.71-0.670.550.490.510.500.450.330.37
 95−0.570.90-0.680.620.610.580.540.450.43
 125−0.180.670.90-0.710.680.610.580.520.45
 1550.020.610.820.97-0.710.650.610.540.48
 185−0.010.750.830.900.95-0.730.680.580.55
 225−0.110.850.840.810.860.96-0.730.630.59
 245−0.090.750.800.790.800.860.94-0.740.63
 2750.000.470.680.790.760.720.750.91-0.67
 305−0.410.670.740.720.730.770.730.610.55-
4X Second
 35-0.490.470.450.350.240.200.260.320.17
 650.40-0.670.600.540.490.450.390.310.30
 950.360.87-0.710.650.570.500.440.380.32
 1250.310.610.91-0.730.640.560.520.470.37
 1550.140.570.860.95-0.710.660.610.540.45
 185−0.040.630.720.700.87-0.760.710.600.55
 2250.010.630.610.530.720.96-0.800.690.63
 2450.270.550.560.550.700.860.93-0.800.69
 2750.410.320.460.590.690.710.750.93-0.74
 305−0.290.170.290.420.730.970.970.880.73-

aMilking frequency

Download Excel Table

The phenotypic correlations were smaller than the corresponding genetic correlations, but they followed a pattern similar to the corresponding genetic correlations for all traits (Table 6). The least phenotypic correlations between 5 and 305 DIM, were 0.10 and 0.22 for 3X cow and 0.04 and 0.17 for 4X cow for the first and second lactations, respectively.

Negative phenotypic and genetic correlations were observed between initial and final test-days. After calving, the cow suffers from post-calving stress and also from an energy deficit. This can be caused the negative values. Negative genetic correlations estimated by RRM using different functions have also been reported for Holstein cattle by Jamrozik and Schaeffer [19], Olori et al. [26], Brotherstone et al. [6] and Kettunen et al. [21], and in Brazil by Cobuci et al. [7], Costa et al. [9] and Bignardi et al. [3, 4]. Jakobsen et al. [18] reported genetic correlations estimates higher than 0.40 for first lactation test-day milk yield of Holstein dairy cattle, therefore slowly lower than some estimates observed in this study. However, lower estimates, even close to zero, were obtained for genetic correlations between test-day milk yields in first lactation by Cobuci et al. [7] and Biassus et al. [2].

Conclusion

The results of this study show the milk yield in the first and the second lactations is significantly affected by milking frequency. The results also indicate that the milk yield is more affected by milking frequency in the middle of the lactation periods. In the first lactation, the Phenotypic and genetic correlations between the 3X milk yields in the studied DIM were higher than those for the 4X while they were similar in the second lactation. In both lactations, the heritability estimated for milk yield in 3X was smaller than those for 4X, therefore the 4X cows could be genetically more accurately evaluated. However, because of the cost of other inputs like labor and failure in the marketing, the application of 4 times milking per day in Iranian Holstein cows should be investigated economically.

Notes

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MMD performed the statistical analysis and provided a draft for the paper that was read and approved by MAF and AAM. All authors read and approved the final manuscript.

References

1.

Armstrong DV. Milking frequency. Proceedings of the 3rd Western Dairy Management Conference. 13–15 March 1997, Las Vegas, Nevada. 1997;79–84.

2.

Biassus IO, Cobuci JA, Costa CN, Roberto P, Rorato N, Neto JB, et al. Genetic parameters for production traits in primiparous Holstein cows estimated by random regression. Revista Brasileira De Zootecnia. 2011; 40:85-94.

3.

Bignardi AB, El Faro L, Cardoso VL, Machado PF. Random regression models to estimate test-day milk yield genetic parameters Holstein cows in southeastern Brazil. Livest Pro Sci. 2009; 123:1-7.

4.

Bignardi AB, El Faro L, Cardoso VL, Machado PF. Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models. J Dairy Sci. 2009; 92:4634-40.

5.

Bohlouli M, Alijani S. Estimation of genetic parameters of milk production traits of Iranian Holstein dairy cattle using multi-trait random regression model. Livest Res Rural Dev. 2012; 24:11.

6.

Brotherstone S, White IMS, Meyer K. Genetic modeling of daily milk yield using orthogonal polynomials and parametric curves. J Anim Sci. 2000; 70:407-15.

7.

Cobuci JA, Euclydes RF, Lopes PS, Costa CN. Estimation of genetic parameters for test-day milk yield in Holstein cows using a random regression models. Genet Mol Biol. 2005; 28:75-83.

8.

Cobuci JA, Costa CN, Neto JB, de Freitas AF. Genetic parameters for milk production by using random regression models with different alternatives of fixed regression modeling. Revista Brasileira De Zootecnia. 2011; 40:557-67.

9.

Costa CN, Melo CNR, Pacher IU, Freitas AF. Genetic parameters for test day milk yield of first lactation Holstein cows estimated by random regression using Legendre polynomials. Revista Brasileira De Zootecnia. 2008; 37:602-8.

10.

De Roos APW, Harbers AGF, de Jong G. Random regression test-day model in The Netherlands. Interbull Bull. 2001; 27:155-158.

11.

DePeters EJ, Smith NE, Acedo-Rico J. Three or two times daily milking of older cows and first- lactation cows for entire lactations. J Dairy Sci. 1985; 68:123-32.

12.

Druet T, Jaffrezic F, Ducrocq V. Estimation of genetic parameters for test-day records of dairy traits in the first three lactations. Genet Sel Evol. 2005; 37:257-71.

13.

El Faro L, Cardoso VL, Albuquerque LG. Variance component estimates applying random regression models for test-day milk yield in Caracu heifers (Bos Taurus Artiodactyla, Bovidae). Genet Mol Biol. 2008; 31:665-73.

14.

Erdman RA, Varner M. Fixed yield responses to increased milking frequency. J Dairy Sci. 1995; 78:1199-203.

15.

Gilmour AR, Gogel BJ, Cullis BR, Thompson R. ASReml User Guide Release 3.0 VSN International Ltd, Hemel Hempstead. HP1 1ES. UK. available at https://www.vsni.co.uk/downloads/asreml/release3/UserGuide.pdf 2009.

16.

Hale SAAV, Capuco TK, Erdman RA. Milk yield and mammary growth effects due to increased milking frequency during early lactation. J Dairy Sci. 2003; 86:2061-71.

17.

Hogeveen H, Ouweltjes W, de Koning CJAM, Stelwagen K. Milking interval, milk production and milk flow-rate in an automatic milking system. Livest Pro Sci. 2001; 72:157-67.

18.

Jakobsen JH, Madsen P, Jensen J, Pedersen J, Christensen LG, Sorensen DA. Genetic parameters for milk production and persistency for Danish Holsteins estimated in random regression models using REML. J Dairy Sci. 2002; 85:1607-16.

19.

Jamrozik J, Schaeffer LR. Estimates of genetic parameters for a test day model with random regressions for yield traits of first lactation Holsteins. J Dairy Sci. 1997; 80:762-70.

20.

Karacao B, Jaffrezic F, Kadarmideen HN. Genetic Parameters for Functional Traits in Dairy Cattle from Daily Random Regression Models. J Dairy Sci. 2006; 89:791-8.

21.

Kettunen A, Mantysaari EA, Poso J. Estimation of genetic parameters for daily milk yield of primiparous Ayrshire cows by random regression test-day models. Livest Pro Sci. 2000; 66:251-61.

22.

Klei LR, Lynch JM, Barbano PA, Lednor AJ, Bandler DK. Influence of milking three times a day on milk quality. J Dairy Sci. 1997; 80:427-36.

23.

Meyer K. Scope of random regression model in genetic evaluation of beef cattle for growth. Livest Pro Sci. 2004; 86:68-83.

24.

Miglior F, Gong W, Wang Y, Kistemaker GJ, Sewalem A, Jamrozik J. Genetic parameters of production traits in Chinese Holsteins using a random regression test-day model. J Dairy Sci. 2009; 92:4697-706.

25.

Muir BL, Kistemaker G, Jamrozik J, Canavesi F. Genetic parameters for a multiple-trait multiple- lactation random regression test-day model in Italian Holsteins. J Dairy Sci. 2007; 90:1564-74.

26.

Olori VE, Hill WG, McGuirk BJ, Brotherstone S. Estimating variance components for test-day milk records by restricted maximum likelihood with a random regression animal model. Livest Pro Sci. 1999; 61:53-63.

27.

Penasa M, Cecchinato A, Battagin M, De Marchi M, Pretto D, Cassandro M. Bayesian inference of genetic parameters for test-day milk yield, milk quality traits, and somatic cell score in Burlina cows. J Appl Genet. 2010; 51:489-95.

28.

Razmkabir M. Estimation of genetic trend for production traits in Iranian Holstein cattle. MS Thesis 2005, Agriculture Faculty, University of Tehran. Tehran. Iran.

29.

Schwarz G. Estimating the dimension of a model. Ann Stat. 1978; 6:461-4.

30.

Stelwagen K, Lacy-Hulbert SJ. Effect of milking frequency on milk somatic cell count characteristics and mammary cell damage in cows. Am J Vet Res. 1996; 57:902-5.

31.

Strabel T, Ptak E, Szyda J, Jamrozik J. Multiple-lactation random regression test-day model for Polish Black and White cattle. Interbull Bull. 2004; 32:133-6.

32.

Strabel T, Szyda J, Ptak E, Jamrozik J. Comparison of random regression test-day models for Polish black and white cattle. J Dairy Sci. 2005; 88:3688-99.

33.

Swalve HH, Guo Z. An illustration of lactation curves stratified by lactation yields within herd. Archiv 240 Tierzucht. 1999; 42:515-25.

34.

Yousefi Golverdi A, Hafezian H, Chashnide Y, Farhadi A. Genetic parameters and trends of production traits in Iranian Holstein population. Afr J Biotechnol. 2012; 11:2429-35.

35.

Zavandilova L, Jamrozik J, Schaeffer LR. Genetic parameters for test day model with random regression for production traits of Czech Holstein cattle. Czech J Anim Sci. 2005; 50:142-54.