The study was conducted to identify a suitable milk-recording protocol for small-scale dairy production. Data from Holstein-FriesianxLocal and SahiwalxLocal crossbred cows on 2 selected farms in Chittagong metropolitan area of Bangladesh form 2010 to 2011 was used to compare and evaluate the co-efficients of 3 different mathematical models of lactation curve for daily and test-day milk yield. The estimated model co-efficients and the 210 days predicted lactation milk yields varied significantly between crossbreds. The models were evaluated with 3 fit statistic values (Coefficient of Determination (R2), Root Mean Square Error (RMSE) and Coefficient of Variation (COV)). The highest lactation milk yield was found for FriesianxLocal on Farm 1. The average daily milk yield of FriesianxLocal was higher than Sahiwal crossbreds in both the farms. The Wood model was found to be suitable to transform test-day milk yield into a 210-days predicted milk yield for all crossbreds based on higher R2 and COV and lower RMSE values. The test-day milk yield of 1 and 2 weeks interval was found to be best fitted with the actual 210-days milk yield by using Wood model for both crossbreds.
PDF Abstract XML References Citation
How to cite this article
In order to increase the productivity of the small-scale farmers, proper management practice is needed and the most important management tool is proper record keeping (Silver, 2006). The purpose of a recording system is to improve the level of herd performance by running the farm more competently, or reducing the possibility of poor future performances. Studies have shown that, the efficiency of milk production can be increased through the use of a simple, precise, understandable and easy to maintain recording system (Chagunda et al., 2006).
In order to monitor and improve the production of cows, milk recording is necessary: in a way, the persistence of satisfactory milk production in low-to medium-input environments is a sign of acceptable sustainability (Duclos et al., 2008).
Mathematical models of lactation curve can be used to predict the total amount of milk yield for the entire milking period using daily or test-day recorded yields. Three parametric mathematical models of lactation yields were developed in the 1960s by Wood (1967) which defines the shape of the lactation curve. However, before Wood, other workers (Ganies, 1927; Nelder, 1966) estimated lactation yields with simple mathematical models but those did not define the shape of the lactation curve accurately.
For establishing the sustainable dairy industry it is necessary to develop and adapt a suitable milk-recording protocol/system for small-dairy herd to undertake future nutritional, management and genetic improvement program. Therefore, the current study was undertaken for small-scale dairying with the objectives: (1) To estimate the actual milk yield of different genotypes (2) To simplify the milk-recording protocol through comparing actual and predicted milk yield by using different lactation models and (3) To chose a suitable milk recording interval for small-holder dairy production.
MATERIALS AND METHODS
Data used in this study was 210 days daily milk yield records collected from November 2010 to June 2011 on 137 cows from two selected dairy farms in Chittagong metropolitan area of Bangladesh. The available genotypes were Holstein-FriesianxLocal (86) and SahiwalxLocal (51) crossbreds and the cows were in third lactation.
Milk recording was performed by the appointed data collector. Milk recording data was collected on daily (24 h) basis in a standard record sheet. The daily milk yield was estimated by adding the morning (AM) and evening (PM) milk yield.
Actual milk yield: The actual 210 days milk yield per cow was calculated by adding all days recorded milk yield data. The mean and standard error for 210 days milk yield of Holstein-FriesianxLocal and SahiwalxLocal crossbred cows under two farming conditions was analyzed using PROC GLM of SAS (2000). The mean differences were estimated using Least Significant Differences (LSD) test (Steel et al., 1997).
Simplification of records by using different models: For the simplification of milk recording, the recording was simulated by putting the intervals between the milk recording (test-days (TDs)). The interval between 2 test-days (TDs) milk recording was gradually increased from 4 up to 7, 14 or 30 days, respectively. The schematic diagram for this sampling is presented in Fig. 1.
Models for the lactation curve: Following three lactation curve models were used to predict the 210 days milk yield from the different milk recording protocol:
|•||Quadratic polynomial model, Y= a+ bx+cx2|
|•||Wood (1967) model, Y = axbe-cx|
|•||Wilmink (1987) model, Y = a+bx+ce-dx|
where, Y is milk yield (kg day-1), x is the days in milk and a, b, c and d are co-efficient that define the shape of the lactation curve.
Statistical analyses: The Wood model was log linear transformed; then for the Wood model and quadratic polynomial regression equation, the test day milk yield was set as dependent variable and days in milk was set as independent variable for obtained the model co-efficient. The models were analyzed by using the SAS (2000). The Wilmink model was reduced to a 3-co-efficient linear model by setting the d exponent to a suitable fixed value. In the present study, d was assumed equal to 0.05 according to Khan (2009). Along with model co-efficient the fit statistics, R2 (coefficient of determinant), COV (coefficient of variation) and RMSE (Root mean square error) were also obtained by SAS analysis. The predicted milk yields were obtained by using the model co-efficient with days in milk as independent variables. To obtain the mean value of different components (model co-efficient, fit statistics and predicted 210 days milk yields), the linear mixed procedure of SAS (2000) was used in the following model:
Yijk = μ + Bi + Fj + eijk
where, Yijk represents the value of predicted milk yield, a, b, c, d, k and fit statistics:
|μ||=||Mean of the data|
|Bi||=||Effect of the ith breed (FriesianxLocal and SahiwalxLocal)|
|Fj||=||Effect of jth farm (Farm 1 and Farm 2)|
|eijk||=||Random residual effect distributed as N (0, σ2)|
|Fig. 1:||Illustration of the recording procedure to simulate the simplification of milk-recording protocols|
The differences of model co-efficient, predicted milk yield and fit statistics were tested with the probability value of p≤0.5.
Actual and predicted lactation milk yield: The mean value along with their Standard Error (SE) of actual and predicted lactation milk yield of FriesianxLocal and SahiwalxLocal in two different farms are presented in Table 1. From Table 1 it is seen that the 210 days actual milk yield was found to be higher (1689±11.74) for Holstein-FriesianxLocal in farm 1 and lower (1553.08±20.75) for SahiwalxLocal genotype in farm 2.
The Predicted Lactation Milk Yield (PLMY) for different models was varied between breeds within the model (Table 1). However, for Wood model, the actual lactation yields calculated from the raw data was close to that of the PLMY. The PLMY for various intervals fitting after Wood model were different between breed groups within the farm and also between breed groups within the intervals (Table 2). The PLMY for 4 and 30 day intervals were different with the actual milk yield.
|Table 1:||Estimated model co-efficient, fit statistics and the predicted milk yield of various models for different genotypes in two farms|
|a: Intercept, b: Shape of the curve, c: Placement of the curve, R2: Coefficient of determinant, RMSE: Root Mean Square Error, COV: Coefficient Of Variation, HFxL: FriesianxLocal, SahxL: SahiwalxLocal, PMY: Predicted milk yield. Mean with different superscripts are different at p<0.05 between breeds within model. NS : Not significant, *,**Significant at 0.05 and 0.01 level, respectively|
But the PLMY for 7 and 14-days interval were similar with the actual milk yield except for SahiwalxLocal (1481.1±110.32) under farm 1.
Use of 3 different lactation models to the base data to predict the 210 days milk yield: The mean value of the estimated model co-efficients (a, intercept; b, shape of the curves and c, placement of the curve; fit statistics (R2, RMSE and COV) and predicted lactation milk yield for three different models along with their standard error for different genotypes (Holstein-FriesianxLocal and SahiwalxLocal) under two different farming conditions are shown in Table 1.
The model co-efficient a, b and c for different models were differed between breed groups within and between the models. In the Polynomial model, the a was different between breed groups and farm but the b and c was not differed between breed groups. The difference in a, b and c was observed between breed groups and farms using Wilmink model. But in case of Wood model, the a was varied between breed groups within the farm. However, the b and c was different between breed groups and within farm.
|Table 2:||Estimated model co-efficients, fit statistics and the predicted milk yield of various intervals/test-day for different genotypes under two farming condition by using wood model|
|a: Intercept, b: Shape of the curve, c: Placement of the curve, R2: Coefficient of determinant, RMSE: Root Mean Square Error, COV: Coefficient Of Variation, HFxL: FriesianxLocal, SahxL: SahiwalxLocal, PMY: Predicted milk yield. Mean with different superscripts are different at p<0.05 between breeds within model. NS : Not significant, *,**significant at 0.05 and 0.01 level, respectively|
The values of three fit statistics differed between breed groups and also in between the farms. The R2, COV and RMSE values for both Polynomial and Wilmink models indicated its inconvenience as fitness. In the case of Wood model, the smaller value of RMSE (0.117±0.02 and 0.118±0.02), for R2 (0.49±0.08 and 0.79±0.08) and COV (5.75±0.529 and 5.91±0.555), the bigger values were observed which indicated the superiority of model in the Holstein-FriesianxLocal crossbred in both farming conditions.
Use of wood models to the different test-day data to predict the milk yield: The Wood model was observed to be good fitted in all genotypes under both the farms therefore the Wood model was used to predict various test-days milk yield from the different milk recording protocols. The estimated model co-efficients (a, b and c), fit statistics (R2, RMSE and COV) and predicted lactation milk yield for Wood model along with their standard error of Holstein-FriesianxLocal and SahiwalxLocal crossbreds under two different farming conditions are shown in Table 2. Only 4, 7 and 14 days intervals values are shown in Table 2.
For 4 days intervals, the Predicted Milk Yield (PMY) varies between the genotypes and within the farms. On the other hand, for 7 days intervals, the PMY differed only in farm but not in breed groups. For 14 days intervals, the PMY differed slightly between the crossbreds and within the farms. While for 30 days interval, the PMY varied highly in between breed groups and also in between the farms.
The model co-efficients (a, b and c) for different interval (4, 7, 14 and 30 days) were differed between genotypes within the model and also between the farms. Among various intervals, values the 14 days intervals were uniformly fitted with the actual milk yield value. The fit statistics (R2, COV and RMSPE) values were also varied between the breeds of different farms and within the breeds in same farm for various intervals.
From the study it has been observed that the highest 210 days milk yield was found in FriesianxLocal and lowest in SahiwalxLocal genotypes in both the farms. Al-Amin and Nahar (2007) found average lactation yield of FriesianxLocal genotype was higher than other breeds groups. The lactation milk yield for both the genotypes was found to be higher in Farm 1 than in Farm 2. The variation is due to the effect of feeding and also due to the effect of different environment surrounding both the farms. The higher production could be due to the effect of GenotypexEnvironment interactions in that particular farm. Rehman et al. (2008) observed that total milk yield per lactation affected by the difference of herd. Similarly, milk production between breed groups, seasons and management systems were reported by other researchers (Val-Arreola et al., 2004; Perochon et al., 1996).
The model co-efficients a, b and c for different models were varied between breed groups within the model and also between the models. The differences of model co-efficients were differed due to the differences in breeds were previously reported by other workers (Khan et al., 2012; Alam et al., 2009; Perochon et al., 1996).
The variation of fit of the models between breeds may have arisen from the differences in breeds, mathematical functions of the models and amount of records. The Polynomial and Wilmink models were not considered further due to their unsuitable model co-efficients and poor fit statistics values. In case of Wood model, the smaller value of RMSE (Val-Arreola et al., 2004) was considered to be superior but for R2 and COV, the bigger values indicated the superiority of FriesianxLocal in both farming conditions. A similar finding was observed by Alam et al. (2009) for the Wood model. The differences of lactation curve traits can be differed with genetic group (Perochon et al., 1996; Khan et al., 2012) and also their fitting ability (Alam et al., 2009).
There was a variation in predicted lactation milk yield among the different genotypes under both farming conditions. For Polynomial model, the lactation yield varied with in breed in between the two farms. In case of Wilmink model, lactation milk yield also varied between the breeds within the farms. However, for the Wood model, the predicted lactation yield was slightly varied between breed but a considerable variation occurred within the two farms. The total lactation milk yields calculated from the raw data for different genotypes under two farms were similar with the total lactation yields predicted by Wood model than the lactation yields predicted by Polynomial and Wilmink model. Finding of Alam et al. (2009) and Khan et al. (2012) were similar to the results of the study.
For 4 days intervals, the predicted lactation milk yield (PLMY) varies between the crossbreds and within the farm. On the other hand, for 7 days intervals, the PLMY differed only in farm but not in breeds. For 14 day intervals, the PLMY differed slightly between the breed and within the farms. While for 30 day interval, the PLMY varied highly in between breeds and also in between the farms. However it was seen that the total lactation yields predicted at 14 day interval by Wood model was similar for both genotypes, except for the SahiwalxLocal under farm 2. Duclos et al. (2008) investigated on simplified milk recording protocol at various test-day intervals and found similar result.
The model co-efficient a, b and c for different interval were varied between breed groups within the models and also between the farms. The b (shape of the curve) was significantly of SahiwalxLocal at 4 and 30 days intervals in both farms and FriesianxLocal also higher in farm 2 at 14 days interval than other genotypes. Three fit statistics (R2, COV and RMSPE) values were used to evaluate the model performance at various intervals. The values of fit statistics varied between the breeds of different farms and within the breeds in same farm for various intervals. The smaller value of RMSE was considered to be superior but for R2 and COV, the bigger values. The fitness of the 3 fit-statistics value indicated that the 7-days interval was good but best for the interval at 14 days with the fitting of Wood model (Khan et al., 2012) repeated the similar result.
The variation of fit of models between breeds may have arisen from differences of test-day yield, the amount of data, the number of test day records and the intervals between tests. The effects of the number of test day records and the intervals between tests on the estimation of fit statistics and also their fitting ability, was reported by Tekerli et al. (2000) and Wiggans et al. (2002).
This study indicated that the Wood model was the most suitable to transform test-day milk yields into a 210 days predicted milk yield for all crossbreds based on higher R2 and COV and lower RMSE values. It was seen that one week (7 days) interval test-day milk recording was good fit with actual 210 days yield. However, two weeks (14 days) interval test-day milk recording was best fit with that of actual 210 days milk yield by using Wood model for both crossbreds based on the different fit statistics value. The Holstein-FriesianxLocal genotypes fitted better than SahiwalxLocal for all models. Such study could be useful for the scientist to undertake further research on genetic improvement programme of dairy cattle under different management system where herd size is small. However, research is needed with more order polynomials and other non-linear and logistic models for the estimation of predicted total yield to confirm the results.
The study was undertaken with the financial support from DANIDA through Poultry Research and Training Center (PRTC) and are hereby acknowledged to the Director PRTC. The authors are grateful to the authority of Chittagong Veterinary and Animal Sciences University, for the permission to conduct this research.
- Alam, M.R., M.K.I. Khan and J. Khanom, 2009. Prediction of lactation milk yield from test day records using wood model. Wayumba J. Anim. Sci., 1: 18-19.
- Al-Amin, M. and A. Nahar, 2007. Productive and reproductive performance of non-descript (local) and crossbred dairy cows in costal area of Bangladesh. Asian J. Anim. Vet. Adv., 2: 46-49.
- Chagunda, M.G.G., A.C.M. Msiska, C.B.A. Wollny, H. Tchale and J.W. Banda, 2006. An analysis of smallholder farmers willingness to adopt dairy performance recording in Malawi. Livestock Res. Rural Dev., Vol. 18.
- Duclos, D., S. Gokhale, R. Bacilieri and V. Ducrocq, 2008. Simplified milk-recording protocols adapted to low-input environments with very small herd size. Animal, 2: 160-166.
- Khan, M.K.I., H.T. Blair and N. Lopez-Villalobos, 2012. Lactation curves of different cattle breeds under cooperative dairying conditions in Bangladesh. J. Applied Anim. Res., 40: 179-185.
- Nelder, J.A., 1966. Inverse polynomials, a useful group of multi-factor response functions. Biometrics, 22: 128-141.
- Perochon, L., J.B. Coulon and F. Lescourret, 1996. Modelling lactation curves of dairy cows with emphasis on individual variability. Anim. Sci., 63: 189-200.
- Rehman, Z.U., M.S. Khan, S.A. Bhatti, J. Iqbal and A. Iqbal, 2008. Factors affecting first lactation performance of Sahiwal cattle in Pakistan. Arch. Tierz. Dummerstorf, 51: 305-317.
- Steel, R.G.D., J.H. Torrie and D.A. Dickey, 1997. Principles and Procedure of Statistics-A Biomedical Approach. 3rd Edn., McGraw-Hill Companies, New York, London, IBN: 9780070610286, Pages: 666.
- Tekerli, M., Z. Akinci, I. Dogan and A. Akcan, 2000. Factors affecting the shape of lactation curves of Holstein cows from the Balikesir province of Turkey. J. Dairy Sci., 83: 1381-1386.
- Val-Arreola, D., E. Kebreab, J. Dijkstra and J. France, 2004. Study of the lactation curve in dairy cattle on farms in central Mexico. J. Dairy Sci., 87: 3789-3799.
- Wiggans, G.R., P.M. Van Raden, J. Bormann, J.C. Philpot, T. Druet and N. Gengler, 2002. Deriving lactation yields from test-day yields adjusted for lactation stage, age, pregnancy and herd test date. J. Dairy Sci., 85: 264-275.
- Wilmink, J.B.M., 1987. Adjustment of test-day milk, fat and protein yield for age, season and stage of lactation. Livestock Prod. Sci., 16: 335-348.
- Wood, P.D.P., 1967. Algebraic model of the lactation curve in cattle. Nature, 216: 164-165.