• Asian-Australasian Journal of Animal Sciences
• /
• v.29 no.6
• /
• pp.775-781
• /
• 2016

A single trait linear mixed random regression test-day model was applied for the first time for analyzing the first lactation monthly test-day milk yield records in Karan Fries cattle. The test-day milk yield data was modeled using a random regression model (RRM) considering different order of Legendre polynomial for the additive genetic effect (4th order) and the permanent environmental effect (5th order). Data pertaining to 1,583 lactation records spread over a period of 30 years were recorded and analyzed in the study. The variance component, heritability and genetic correlations among test-day milk yields were estimated using RRM. RRM heritability estimates of test-day milk yield varied from 0.11 to 0.22 in different test-day records. The estimates of genetic correlations between different test-day milk yields ranged 0.01 (test-day 1 [TD-1] and TD-11) to 0.99 (TD-4 and TD-5). The magnitudes of genetic correlations between test-day milk yields decreased as the interval between test-days increased and adjacent test-day had higher correlations. Additive genetic and permanent environment variances were higher for test-day milk yields at both ends of lactation. The residual variance was observed to be lower than the permanent environment variance for all the test-day milk yields.

• Asian-Australasian Journal of Animal Sciences
• /
• v.19 no.12
• /
• pp.1696-1701
• /
• 2006

A random regression model was applied for the first time for the analysis of test day records and to study the genetic persistency of first lactation milk yield of Indian Murrah buffaloes. Wilmink's Function was chosen to describe the shape of lactation curves. Heritabilities of test day milk yield varied from 0.33 to 0.58 in different test days. The highest heritability was found in the initial test day ($5^{th}$ day) milk yield. Genetic correlations among test day milk yields were higher in the initial test day milk yield and decreased when the test day interval was increased. The magnitude of genetic correlations between test day and 305 day milk yield varied from 0.25 to 0.99. The genetic persistencies of first lactation milk yield were estimated based on daily breeding values using two methods. $P_1$ is the genetic persistency estimated as a summation of the deviation of estimated daily breeding value on days to attain peak yield from each day after days to attain peak yield to different lactation days. $P_2$ is the genetic persistency estimated as the additional genetic yield (gained or lost) from days to attain peak yield to estimated breeding value on different lactation days relative to an average buffalo having the same yield on days to attain peak yield. The mean genetic persistency on 90, 120, 180, 240, 278 and 305 days in milk was estimated as -4.23, -21.67, -101.67, -229.57, -330.06 and -388.64, respectively by $P_1$, whereas by $P_2$ on same days in milk were estimated as -3.96 (-0.32 kg), -23.94 (-0.87 kg), -112.81 (-1.96 kg), -245.83 (-2.81 kg), -350.04 (-3.28 kg) and -407.58 (-3.40 kg) respectively. Higher magnitude of rank correlations indicated that the ranking of buffaloes based on their genetic persistency in both methods were similar for evaluation of genetic persistency of buffaloes. Based on the estimated range of genetic persistency three types of genetic persistency were identified. Genetic correlations among genetic persistency in different days in milk and between genetic persistencies on the same day in milk were very high. The genetic correlations between genetic persistency for different days in milk and estimated breeding value for 305 DIM was increased from 90 DIM to 180 DIM, and highest around 240 DIM which indicates a minimum of 240 days as an optimum first lactation length might be required for genetic evaluation of Indian Murrah buffaloes.

• Asian-Australasian Journal of Animal Sciences
• /
• v.15 no.3
• /
• pp.315-320
• /
• 2002

Test-day records (n=13677) sampled from 896 ewes in 5-9 (${\mu}$=7.5) monthly test-days were used to estimate genetic and phenotypic parameters of test-day yields, lactation milk yield (TMY), length of the milking period (DAYS) and three measures of persistency of milk yield in Boutsico dairy sheep. Τhe measures of persistency were the slope of the regression line (${\beta}$), the coefficient of variation (CV) of the test-day milk yields and the maximum to average daily milk yield ratio (MA). The estimates of variance components were obtained under a linear mixed model by restricted maximum likelihood. The heritability of test-day yields ranged from 0.15 to 0.24. DAYS were found to be heritable ($h^2$=0.11). Heritability estimates of ${\beta}$, CV and MA were 0.15, 0.13, 0.10, respectively. Selection for maximum lactation yields is expected to result in prolonged milking periods, high rates of decline of yields after peak production, variable test-day yields and higher litter sizes. Selection for flatter lactation curves would reduce lactation yields, increase slightly the length of the milking period and decrease yield variation as well as litter size. The most accurate prediction of TMY was obtained with a linear regression model with the first five test-day records.

• Asian-Australasian Journal of Animal Sciences
• /
• v.16 no.9
• /
• pp.1285-1290
• /
• 2003

The present study was carried out to investigate the effect of test day milk yield, test day evening milk yield, parity, stage of lactation and body weight on milk urea and milk protein concentration. A total of 319 milk samples was collected from buffaloes over four month's period and subjected to urea and protein analysis. Milk urea concentration (mg/dl) was significantly (p<0.01) increased with increasing test day milk yield. The lowest value ($57.03{\pm}1.13$) was observed in the milk yield group ${\leq}4.5kg/day$ and the highest value ($64.15{\pm}1.13$) in the group 7.7-10.7 kg/day. However, test day evening milk yield had no significant effect on milk urea concentration. Milk protein did not vary significantly with the test day milk yield as well as test day evening milk yield. A clear decreasing trend of milk urea concentration (mg/dl) was found with the increasing parity. The highest MU concentration ($64.03{\pm}1.14$) was found in the first parity and the lowest ($55.67{\pm}1.22$) was found in the sixth and above parity. Whereas, stage of lactation had no effect on milk urea concentration. Moreover, parity and stage of lactation did not have any significant effect on milk protein concentration. Body weight (kg) was also found negatively (p<0.05) related with urea content (mg/dl) in milk. The highest mean MU concentration ($64.34{\pm}0.88$) was found when body weight was between 532 and 598 kg and lower mean values ($59.24{\pm}0.94$ and $59.33{\pm}1.23$) were observed in 599 to 665 kg and ${\geq}666kg$ group. Body weight also had significant (p<0.05) effect on milk protein content. The highest milk protein content (%) was found in ${\geq}666kg$ group and the lowest in <531 kg group. In conclusion, for proper interpretation of milk urea values to monitor protein nutrition status of the buffaloes parity, milk yield and body weight should be considered.

• Animal Bioscience
• /
• v.35 no.10
• /
• pp.1499-1511
• /
• 2022

Objective: This study was aimed to estimate the genetic parameters, including genetic and phenotypic correlations, of milk yield, lactation curve traits and milk composition of Thai dairy cattle from three government research farms. Methods: The data of 25,789 test-day milk yield and milk composition records of 1,468 cattle from lactation 1 to 3 of Holstein Friesian (HF) and crossbred HF dairy cattle calved between 1990 and 2015 from three government research farms in Thailand were analysed. 305-day milk yield was estimated by the Wood model and a test interval method. The Wood model was used for estimating cumulative 305-day milk yield, peak milk yield, days to peak milk yield and persistency. Genetic parameters were estimated using linear mixed models with herd, breed group, year and season of calving as fixed effects, and animals linked to a pedigree as random effects, together with a residual error. Univariate models were used to estimate variance components, heritability, estimated breeding values (EBVs) and repeatability of each trait, while pairwise bivariate models were used to estimate covariance components and correlations between traits in the same lactation and in the same trait across lactations. Results: The heritability of 305-day milk yield, peak milk yield and protein percentage have moderate to high estimates ranging from 0.19 to 0.45 while days to peak milk yield, persistency and fat percentage have low heritability ranging from 0.08 to 0.14 in lactation 1 cows. Further, heritability of most traits considered was higher in lactation 1 compared with lactations 2 and 3. For cows in lactation 1, high genetic correlations were found between 305-day milk yield and peak milk yield (0.86±0.07) and days to peak milk yield and persistency (0.99±0.02) while estimates of genetic correlations between the remaining traits were imprecise due to the high standard errors. The genetic correlations within the traits across lactation were high. There was no consistent trend of EBVs for most traits in the first lactation over the study period. Conclusion: Both the Wood model and test interval method can be used for milk yield estimates in these herds. However, the Wood model has advantages over the test interval method as it can be fitted using fewer test-day records and the estimated model parameters can be used to derive estimates of other lactation curve parameters. Milk yield, peak milk yield and protein percentage can be improved by a selection and mating program while days to peak milk yield, persistency and fat percentage can be improved by including into a selection index.

• Asian-Australasian Journal of Animal Sciences
• /
• v.19 no.7
• /
• pp.915-921
• /
• 2006

Various random regression models with different order of Legendre polynomials for permanent environmental and genetic effects were constructed to predict future milk yield of Holstein cows in Korea. A total of 257,908 test-day (TD) milk yield records from a total of 28,135 cows belonging to 1,090 herds were considered for estimating (co)variance of the random covariate coefficients using an expectation-maximization REML algorithm in an animal mixed model. The variances did not change much between the models, having different order of Legendre polynomial, but a decreasing trend was observed with increase in the order of Legendre polynomial in the model. The R-squared value of the model increased and the residual variance reduced with the increase in order of Legendre polynomial in the model. Therefore, a model with $5^{th}$ order of Legendre polynomial was considered for predicting future milk yield. For predicting the future milk yield of cows, 132,771 TD records from 28,135 cows were randomly selected from the above data by way of preceding partial TD record, and then future milk yields were estimated using incomplete records from each cow randomly retained. Results suggested that we could predict the next four months milk yield with an error deviation of 4 kg. The correlation of more than 70% between predicted and observed values was estimated for the next four months milk yield. Even using only 3 TD records of some cows, the average milk yield of Korean Holstein cows would be predicted with high accuracy if compared with observed milk yield. Persistency of each cow was estimated which might be useful for selecting the cows with higher persistency. The results of the present study suggested the use of a $5^{th}$ order Legendre polynomial to predict the future milk yield of each cow.

• Asian-Australasian Journal of Animal Sciences
• /
• v.18 no.9
• /
• pp.1231-1236
• /
• 2005

Genetic parameters for first lactation milk production based on test day (TD) records of 56319 Iranian Holstein cows from 655 herds that first calved between 1991 and 2001 were estimated with restricted maximum likelihood method under an Animal model. Traits analyzed were milk yield and milk fat percentage. Heritability for TD records were highest in second half of the lactation, ranging from 0.11 to 0.19 for milk yield and 0.038 to 0.094 for milk fat percentage respectively. Estimates for lactation records for these traits were 0.24 and 0.26 respectively. Genetic correlations between individual TD records were high for consecutive TD records (>0.9) and decreased as the interval between tests increased. Estimates of genetic correlations of TD yield with corresponding lactation yield were highest (0.78 to 0.86) for mid-lactation (TD3 to TD8). Phenotypic correlations were lower than corresponding genetic correlations, but both followed the same pattern. For milk fat percentage no clear pattern was found. Results of this study suggested that TD yields especially in mid-lactation may be used for genetic evaluation instead of 305-day yield.

• Asian-Australasian Journal of Animal Sciences
• /
• v.21 no.8
• /
• pp.1088-1096
• /
• 2008

Data from milk performance testing were used to analyze genetic and environmental trends for purebred Tsigai, Improved Valachian and Lacaune sheep. 103,715 (Tsigai), 212,962 (Improved Valachian) and 2,196 (Lacaune) test-day records gathered by the State Breeding Institute of the Slovak Republic entered the analyses. The respective pedigree data comprised 23,724 (Tsigai), 51,401 (Improved Valachian) and 438 (Lacaune) records. The multiple-trait, mixed model methodology was used to predict the breeding values for daily milk yield, fat and protein content and to estimate the fixed and remaining random effects assumed to affect the above mentioned traits, separately for each breed. The breeding values for daily milk yield were adjusted for 150-day standardized lactation length by multiplying with the constant 150, as the breeding goal of the selection scheme in Slovakian sheep is to increase 150-day milk production and constant heritability throughout the whole lactation is assumed. The genetic trends were expressed as changes in averages of breeding values across birth years of animals. For Tsigai and Lacaune breeds, cumulative genetic changes over the analyzed period were 3.8 and 5.1 kg for 150-day milk, 0 and -0.16% for fat content and 0 and -0.12% for protein content. For Improved Valachian breed, either a low (1.6 kg for 150-day milk yield) or zero (fat and protein content) cumulative genetic change was found. The environmental trends were calculated as averages of solutions for flock-test day effect across years and months in which measurements were taken. A distinctive cyclical pattern which reflected short-time variation in milk production traits was found. Possible explanations for this phenomenon are given and discussed.

• Journal of Animal Science and Technology
• /
• v.58 no.2
• /
• pp.5.1-5.9
• /
• 2016

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.

• Asian-Australasian Journal of Animal Sciences
• /
• v.29 no.5
• /
• pp.607-614
• /
• 2016

The objectives of the study were to estimate genetic parameters for milk production traits of Holstein cattle using random regression models (RRMs), and to compare the goodness of fit of various RRMs with homogeneous and heterogeneous residual variances. A total of 126,980 test-day milk production records of the first parity Holstein cows between 2007 and 2014 from the Dairy Cattle Improvement Center of National Agricultural Cooperative Federation in South Korea were used. These records included milk yield (MILK), fat yield (FAT), protein yield (PROT), and solids-not-fat yield (SNF). The statistical models included random effects of genetic and permanent environments using Legendre polynomials (LP) of the third to fifth order (L3-L5), fixed effects of herd-test day, year-season at calving, and a fixed regression for the test-day record (third to fifth order). The residual variances in the models were either homogeneous (HOM) or heterogeneous (15 classes, HET15; 60 classes, HET60). A total of nine models (3 orders of $polynomials{\times}3$ types of residual variance) including L3-HOM, L3-HET15, L3-HET60, L4-HOM, L4-HET15, L4-HET60, L5-HOM, L5-HET15, and L5-HET60 were compared using Akaike information criteria (AIC) and/or Schwarz Bayesian information criteria (BIC) statistics to identify the model(s) of best fit for their respective traits. The lowest BIC value was observed for the models L5-HET15 (MILK; PROT; SNF) and L4-HET15 (FAT), which fit the best. In general, the BIC values of HET15 models for a particular polynomial order was lower than that of the HET60 model in most cases. This implies that the orders of LP and types of residual variances affect the goodness of models. Also, the heterogeneity of residual variances should be considered for the test-day analysis. The heritability estimates of from the best fitted models ranged from 0.08 to 0.15 for MILK, 0.06 to 0.14 for FAT, 0.08 to 0.12 for PROT, and 0.07 to 0.13 for SNF according to days in milk of first lactation. Genetic variances for studied traits tended to decrease during the earlier stages of lactation, which were followed by increases in the middle and decreases further at the end of lactation. With regards to the fitness of the models and the differential genetic parameters across the lactation stages, we could estimate genetic parameters more accurately from RRMs than from lactation models. Therefore, we suggest using RRMs in place of lactation models to make national dairy cattle genetic evaluations for milk production traits in Korea.