|
|
||||||||
Department of Dairy Science Virginia Polytechnic Institute and State University Blacksburg 24061-0315
Corresponding author:
B. G. Cassell; e-mail:
bcassell{at}vt.edu.
| ABSTRACT |
|---|
|
|
|---|
Abbreviation key: DPL = days of productive life, HYS = herd-year-season, OC = opportunity cost, PL = productive life, RNI = relative net income, RNIOCF = RNI adjusted for opportunity cost in a fluid milk market, RNIOCM = RNI adjusted for opportunity cost in a manufactured milk market
Key Words: lifetime profit Net Merit sire evaluations
| INTRODUCTION |
|---|
|
|
|---|
Factors limiting use of AI include misconceptions of the true value of superior AI progeny, labor requirements for heat detection, quality of labor required for success, inability to consistently manage a herd for optimal performance through AI, perceptions of high cost for AI compared to natural service, and delay in receipt of benefits from use of genetically superior bulls. Dairy production is a competitive enterprise in the United States, with low per unit costs of production an essential ingredient for survival of the dairy business. Economic efficiencies of genetically superior cattle, if documented and presented to producers in a readily comprehensible manner, might be helpful in encouraging producers to more fully utilize the economic benefits of AI. Relationships between the many traits for which genetic evaluations are available and economic returns through daughter performance are complex and difficult to interpret, as shown by Weigel et al., 1997 and VanRaden, 2000. Long-term economic implications of alternative sire selection policies could be compared if useful prediction equations and estimates of genetic trend for a variety of traits were incorporated into computer graphics programs. Such an approach requires an estimate of lifetime economic merit. One function that has been used in several studies is lifetime relative net income, (Weigel et al., 1995a; Weigel et al., 1995b; Smith et al., 1998). Kulak, 1994, and Charagu, 1997, used similar but separately developed procedures to evaluate lifetime economic merit of dairy cows.
Many researchers and producers have embraced lifetime economic merit as the most useful breeding objective for genetic improvement programs in dairy cattle. VanRaden, 2000, described three indexes (Net, Fluid, and Cheese Merit) that included sire evaluations for production, fitness, and type traits, properly weighted to improve total lifetime economic merit of the dairy cow. References to Merit indexes throughout this manuscript refer to the indexes developed by VanRaden, 2000. Economic weights used in the Merit indexes (VanRaden, 2000) were determined from published research of expected costs and returns similar to Rogers, 1993, and by retrospective evaluations of lifetime economic merit from traits measured early in life similar to Weigel et al.(1995a) and Weigel et al. (1997). Documentation of the utility of the Merit indexes would encourage their utilization in sire selection decisions.
Our first objective in this study was to quantify differences in lifetime economic merit of daughters of bulls with no record of use in AI and bulls in AI service, using an appropriate lifetime economic merit function. Our second objective was to develop equations to predict an appropriate lifetime economic merit function using different combinations of bull proofs for production, fitness, and linear type. Our intent in this objective was to predict response to alternative sire selection practices, rather than to compare economic weights for different traits with the partial regression coefficients. Our third objective was to evaluate the utility of the Merit indexes using an alternative approach to estimate lifetime economic merit of cows. Here, our purpose was to compare the accuracy of prediction of lifetime economic merit by the Merit indexes to alternative weightings of the traits used in the Merit indexes. Ultimately, our purpose was to incorporate findings into a readily comprehensible, graphic representation of long term effects of alternative sire selection strategies that could be shown to producers using microcomputers.
| MATERIALS AND METHODS |
|---|
|
|
|---|
![]() |
+(number of lactations)(net value of a calf)
+net salvage value – rearing costs
–(total days in milk)(daily expenses per day in milk)
–(total days dry)(daily expenses per day dry)
Summation is over j components (milk or carrier, fat, and protein) and i lactations initiated before the end of the 84 mo herdlife opportunity. Further details of the economic assumptions summarized here are available (Weigel et al., 1995b). We used two milk prices, one appropriate for manufacturing milk markets and one appropriate for fluid milk markets. Milk prices are represented by "value" in the equation above. Feed costs for production of milk, fat, and protein were determined as suggested by Dado et al. (1994) and are represented by "cost" in the above equation. Protein is assigned a cost of production but has no value in the marketplace in fluid milk markets. Protein has the same cost of production, but also generates income in manufactured milk markets. The net value of the calf included costs to obtain a pregnancy in the dam, sex ratio, differential value of male and female calves, and average death losses. Sex of calf x lactation was not available. Net salvage value reflects the difference between the value of a fresh heifer and a cull cow, with cull cow value reduced for the average probability of death loss that was constant across cows and herds. Rearing costs were based on age at first freshening. Cows milked three times daily were charged $3.53/d compared with $2.80/d for cows milked twice daily for labor and facilities involved in milk harvesting. Daily expenses considered were labor costs for milking and handling of lactating cows when in milk, with lower costs during the dry period. Daily expenses also included feed costs for body maintenance for days milking and dry and fixed costs for facilities. Actual milk, fat, and protein production for up to the first 305 d of each lactation was used. All days beyond 305 d of production were considered to be dry days, even for cows that milked considerably longer than 305 d.
Lifetime Relative Net Income Adjusted for Opportunity Cost
We estimated opportunity cost (OC) from the average RNI on all heifers that freshened for the first time in each herd and year combination. Days of productive life (DPL) are necessary to calculate OC, and are the interval from first calving to (last calving + DIM for the last record). DPL is calculated differently from Productive Life (PL), the trait used for genetic evaluations of herd life in the United States. PL includes a maximum of 10 mo per lactation, whereas DPL includes the entire interval between first calf and culling or death. We used DPL to calculate OC, but PL for other references to herd life in this study. Opportunity cost was estimated by RNI/DPL, and is thus expressed as a cost per day in the herd. Opportunity cost can only be calculated directly for those herd and year of first freshening combinations that permit an 84-mo herdlife opportunity for first calf heifers. Production records initiated too recently for 84-mo herd life opportunity can contribute to calculation of RNI by the equation above, but cannot contribute directly to estimation of OC. We predicted RNI/DPL for each herd for those more recent years using linear regression from earlier OC for that herd. Details of estimation of OC for recent years are in Smith et al. (1998).
RNI adjusted for OC (RNIOCM for manufactured milk prices and RNIOCF for fluid milk prices) for each cow was the difference in her RNI and her total OC. Total OC was summed across lactations, multiplying the appropriate OC for the current lactation by the calving interval for that lactation on each individual cow in the study.
RNI increases with genetic trend and, since it is a function of RNI, so does OC. RNIOCM or RNIOCF, however, have almost no trend over time, because both cows born in later years and OC for those cows increase at a similar rate. RNIOCM or RNIOCF measure the ability of an individual cow to utilize herd resources relative to heifers freshening in the same year. A measure of lifetime economic merit without adjustment for OC (RNI in this study) measures the ability of a cow to generate return above costs, with no consideration for whether another cow would generate more or less return above costs from the same resources. We consider RNIOCM or RNIOCF to be the correct measure of lifetime economic merit to use to evaluate lifetime economic merit of cows from historical data when the objective is to determine the importance of individual traits to lifetime economic merit, or to compare past performance of cows produced by different breeding programs. However, prediction equations for RNI adjusted for OC by methods used in this study are inappropriate for RNIOCM or RNIOF, because the predicted values should sum to zero and display no trend. Prediction of RNI is appropriate and useful as long as it is understood that RNI does not consider alternative uses of the same resources (OC). We use RNI in this paper to predict future response to sire selection practices that remain consistent over time, since sire merit improves with genetic trend. The most appropriate use of such predictions would be to evaluate differences between alternative sire selection protocols, where one system would serve as "opportunity cost" for the other.
Data Used
Cows included in this study were from herds classified for type by the Holstein Association, Brattleboro, VT, from 1983 to 1990, described in detail by Smith et al. (1998). We did not use any classification scores on females or require that they be scored, but did use genetic evaluations for linear types of traits on sires of cows in those herds. Lifetime production records were obtained from the Animal Improvement Programs Laboratory (USDA, Beltsville, MD), with earliest first calvings in 1983 (birth year 1980), and most recent first calvings in 1991 (birth year about 1988). Calving dates for second and later records continued through 1995. Genetic evaluations for production and fitness traits PL and SCS on sires of cows in the study were provided by the Animal Improvement Programs Laboratory, while genetic evaluations for type on the same bulls were obtained from the Holstein Association. We used the November 1999 genetic evaluations. The age of the data in this study is dictated by herdlife opportunity requirements, a restriction that affects all lifetime performance studies based on field data.
RNI and RNIOC for fluid and manufacturing milk market prices (four variables per cow) were calculated by procedures described above. Sire evaluations for production, PL, and SCS were merged with lifetime economic merit of daughters, producing records for 2,493,585 cows in over 129,000 herd-years of first freshening. We used sire identification from cow records and USDA genetic evaluation history files to assign those cows to four groups based on characteristics of the bull that sired them. Sires of cows in Group 1 were bulls with no record of ever being in AI service. Sires of cows in Group 2 were in AI service, but were less than 4 years of age at the time the cow was conceived. Our interpretation is that these bulls would have been in an AI progeny test program when they were used. Sires of cows in Group 3 were in AI but were from 4 to 6 years of age at time of conception. Our interpretation is that these bulls were used without benefit of progeny test results, but were probably beyond their initial sampling period. The purpose of creating this group was to prevent such sires from being part of the other two groups of AI sires. Sires of cows in the fourth group were in AI and at least 6 years of age at time of conception of the daughter. Genetic evaluations based on progeny performance should have been available at the time of use for these bulls, and we interpret results from this group as indicative of response to breeding programs using proven AI bulls.
![]() | 1 |
Least squares contrasts between these groups were obtained from model [1] where herd-year of first freshening (hyi) was absorbed and sire group effect (Gi) was fit to various lifetime performance traits of cows. In model 1, Yijk represents RNIOCM, RNIOCF, or another lifetime performance trait of the kth cow calving for the first time in the ith herd and year (absorbed), and sired by a bull in the jth group as explained above with residual error eijk. The measure of lifetime economic merit used in this analysis was RNIOC because our interest was in comparing the use of herd resources between different kinds of sire selection policies, rather than in predicting future differences after sire merit changes through genetic trend. We used PROC GLM of SAS (1988) for analysis.
The second objective, predicting future profitability from long-term application of alternative sire selection strategies, used RNI rather than RNIOC as the response variable. The analysis also required a different model and additional restrictions on the data. For Model 2, we required that genetic evaluations for linear type traits be available on sires of cows in the data. This requirement left 2,136,141 cows in over 128,000 herd-years of first freshening. We no longer used the sire group effects described in Model 1 in this analysis. Data included cows sired by 6663 different sires. The model used was
![]() | 2 |
where yijk represented RNI (not RNIOC) from fluid or manufacturing milk markets on the kth daughter of the jth sire calving for the first time in herd-year of first freshening i. We used three different combinations of genetic evaluations (PTAs) on production, PL, SCS, and linear type traits to predict RNI in fluid or manufacturing milk markets of daughters. Subscript l in equation [2] represents the different traits, and n would vary depending on which of the three combinations of traits was included in an analysis. The effects of herd-year of first freshening were absorbed.
RNIOC from the 2,136,141 records evaluated by model [2] were also used to compare the accuracy of prediction of lifetime economic merit from the USDA Merit indexes to alternative weights applied to the same traits used in the Merit indexes. We divided the 2.1 million cows into two independent data sets of approximately equal size based on last digit (odd or even) of the DHI herd code number. Sire evaluations for the eight traits (or composites of linear type traits) used in the Merit indexes were used in Model 2 to develop two prediction equations, one from odd numbered herds and one from even numbered herds, for RNIOC using fluid or manufactured milk market prices. We then used the prediction equation developed from odd numbered herds to predict lifetime economic merit from genetic evaluations on sires of individual cows in even numbered herds and vice versa. This approach allowed us to predict lifetime economic merit on cows from prediction equations derived from an independent data source. Finally, we fit the linear regression of RNIOC for fluid or manufactured milk market prices on each cow on Net or Fluid Merit on each cows sire or on our prediction of lifetime economic merit. Herd-year of first freshening was absorbed, as in Model 2. The R2 values from these prediction equations were used to evaluate relative accuracy of the alternative indexes.
The Merit indexes were not available in the November 1999 genetic evaluations, which we used in this study, as USDA first released them in August 2000. The new formulas (VanRaden, 2000) were available to us and were used. Genetic evaluations from both USDA and HA, Inc. are included in the Merit indexes. One of the type composites used is the Feet and Legs Composite, and it includes the trait "Feet and Legs Score", a relatively recent addition to the type evaluation program of Holstein Association, Inc. Genetic evaluations for Feet and Legs Score is missing on 5099 of the 6663 sires involved, but daughters of those bulls were only 25% of the 2.1 million cows involved in this analysis. Feet and Legs Composite receives from 3 to 4% of total emphasis in the Merit indexes and feet and legs score is only half of the composite. Rather than discard a quarter of the data for an item of information that received little weight, we substituted the average Feet and Leg Score for the 1564 bulls with information available, 0.106, when Feet and Leg Score was missing on the sire of a cow. Other traits in the feet and leg composite were available on all bulls.
| RESULTS |
|---|
|
|
|---|
|
|
We fit genetic merit of sires of cows within groups in an analysis not shown, and obtained significant regressions of RNIOC on sire merit within all groups. Thus, use of carefully selected AI-proven bulls and non-AI bulls of poor pedigree merit would magnify the differences reported in Table 2
. Likewise, use of mediocre proven AI bulls and non-AI bulls carefully selected on pedigree merit would diminish or even reverse differences between groups of bulls. The decision to use AI rather than natural service has an impact on net income, but selection of better bulls remains critical to optimization of returns from a breeding program.
Daughters of AI-proven bulls had more months of productive life and produced more lifetime total milk, fat, and protein than daughters of bulls in any of the other three groups. Daughters of AI bulls in sampling at time of use (AI bulls <4 yr old) were the second most profitable of the four groups, but generated, on average $60 less lifetime-relative net-income adjusted for opportunity cost than daughters of AI bulls with proofs at time of use. This $60 can be viewed as income forfeited by foregoing the increased accuracy of selection resulting from progeny test results in sire selection. Our purpose in creating the 4 to 6-yr-old group was to prevent them from affecting comparisons of AI sample sires with proven AI bulls, rather than to evaluate the utility of such a group for genetic improvement. However, some producers do use AI bulls with only pedigree data available well after sufficient services have been performed for progeny test. Cassell et al. (1992) reported that heavily-used, AI-young sires were only marginally superior to all AI sample sires. Results of this study show that bulls used beyond the sampling period, but prior to progeny test results, produced slightly less profitable daughters than AI young sires used during their progeny testing period. Such bulls were not selected as accurately or intensely as proven AI bulls, and were genetically inferior to AI sample sires, presumably due to genetic trend.
Future Profit from Sire Selection Decisions
Most producers are interested in improving several traits at one time through sire selection, but may not be aware of the impact of selection for some of those traits on profitability of the resulting daughters. We wanted to show producers the impact of different sire selection approaches by predicting RNI using combinations of genetic evaluations for different traits on sires. RNI was used as the response variable because our intent was to predict future net income. RNI does not account for alternative uses of the same resources, as does RNIOC (Van Arendonk, 1990). However, RNI predicted from two different sire selection policies could be used as "opportunity cost" to each other. Care should be taken in interpreting predictions of future RNI without consideration of an alternative use of the same resources.
Table 3
includes three prediction equations each for RNI in fluid and manufactured milk markets based on genetic evaluations of the sire of the cow. Model 1 in Table 3
includes only those traits evaluated by AIPL from DHI data, milk, fat, protein, PL, and SCS. Model 2 includes those five traits and all 16 linear type traits evaluated by Holstein Association at the time of this study. Model 2 is affected by multiple co-linearity, as several of the linear type traits are highly correlated (Weigel et al., 1997). Model 3 removes body depth, fore udder attachment, and rear udder height from Model 2 as suggested by Weigel et al. (1997), because of multiple co-linearity problems with other traits and is the preferred prediction equation for future RNI with linear type data included.
|
Herd and year of first freshening alone produced R2 of .143 and .160 for fluid and manufactured milk market prices, respectively. Genetic evaluations of sires of these cows, while significant, did not explain much of the remaining variation between cows for RNI within herd and year of first freshening. There are several reasons for this result. Heritability of RNI in fluid and manufacturing milk markets (without adjustment for opportunity cost) is 15% (Smith et al., 1998). Much of the variation in profit functions for these cows is controlled by factors other than additive gene action. Genetic merit of the sire only partly controls half of the additive variation, as all gene samples transmitted to daughters are subject to mendelian segregation. Predictions of RNI are based on the traits listed for each model in Table 3
, and other traits such as disease resistance and reproduction are not included but also affect genetic control of RNI. Finally, the accuracy of PTAs for traits in Table 3
affect R2 values. Genes transmitted to each daughter are a random sample from the true breeding value of each sire, which is estimated with imperfect accuracy by PTA for each trait.
Partial regressions for production traits in Table 3
change considerably for RNI using fluid milk market or manufacturing milk market prices. Partial regressions for protein are negative for fluid markets, but positive for manufacturing milk markets. The RNI function includes a cost to produce protein under both markets, but protein generates revenue only in the manufactured milk market. Producers shipping milk to fluid markets face such economic conditions. AI bulls with high PTA for protein and premium semen prices may be more appropriate choices for producers who ship milk to manufactured milk markets than for producers shipping to fluid milk markets.
Models 1 and 2 in Table 3
show that linear type data added very little to the accuracy of prediction of RNI. R2 values increase only slightly under manufactured milk market prices, and not at all for fluid milk market prices. These results are consistent with Weigel et al. (1997), who reported that lower milk prices under manufacturing milk market conditions increased the relative value of functional type and fitness traits. Genetic evaluations for PL for these bulls were heavily influenced by direct evaluations of PL on progeny, because of the age of the data. This effect may reduce the relative value of linear type information compared to bulls with progeny not yet old enough for direct expression of herdlife. Models 2 and 3 reflect the effect of multiple co-linearity on partial regression coefficients for pairs of highly correlated traits such as stature, strength, and body depth or fore udder attachment, rear udder height and width, and udder depth. Elimination of strength, fore udder attachment and rear udder height changed size and sometime sign of partial regressions for remaining traits with no reduction in R2. For instance the partial regression coefficient for body depth predicting RNIM changed from –51.76 to 16.50 and stature changed from –20.41 to –10.57 with the elimination of strength. Partial regressions for SCS were also affected by elimination of linear traits, but SCS was significant for only one of the six models in Table 3
. SCS on individual cows was not used to calculate RNI in this study and would only improve accuracy of prediction through correlated effect on yields, length of productive life, reproduction and so forth. A more realistic evaluation of actual costs and returns from SCS would charge for treatment of mastitis and discarded milk and reward cows with low SCS for milk quality premiums.
Accuracy of Merit Indexes
Animal Improvement Programs Laboratory published three selection indexes in August 2000 (VanRaden, 2000) that were structured to improve lifetime economic merit of dairy cows. Net, Fluid, and Cheese Merit were intended for producers selling milk to markets that price milk according to use for a variety of milk products, predominantly for fluid consumption, or for the production of cheese, respectively. The lifetime economic merit functions we used did not mimic cheese market conditions, and we chose to study Net and Fluid Merit only. Those indexes had not been released when this study was conducted, but the formulas were available. We calculated Net and Fluid Merit indexes on the 6663 Holstein bulls represented in the data. Net and Fluid Merit averaged –$40 and –$29 with standard deviations of $171 and $179 respectively. By contrast, average Net Merit for active AI Holsteins in February 2002 averaged $338 with a standard deviation of $139. The difference in mean reflects genetic progress in the Holstein bull population between birth of sires of cows in this study and current active AI bulls, while smaller variation in Net Merit for current active AI bulls likely results from selection decisions by bull studs based on market potential.
Table 4
includes regression coefficients for prediction of RNIOC for fluid or manufactured milk prices using the eight traits or composites of linear type traits used in the Merit indexes. We used RNIOC in this analysis as our purpose was to predict lifetime economic merit calculated from existing data, and we did not intend to use the results to predict future lifetime economic merit. Table 4
contains two prediction equations each for fluid and manufactured milk market prices, one based on cows in herds with even DHI herd code numbers and the other from cows in odd numbered herds. Partial regression coefficients and R2 values were similar from odd and even herds within fluid or manufactured milk price assumptions. Partial regression coefficients for somatic cell score did change, but were never significantly different from 0.0 (P > 0.01) in any of the four prediction equations. SCS for individual cow records were not considered in calculation of lifetime economic merit due to missing SCS data (Smith et al., 1998). Partial regressions for PTA protein differed in sign for manufactured or fluid milk market price assumptions, as explained earlier. Partial regressions for PTA milk were larger with fluid than with manufactured milk prices, and partial regressions for PTA fat were not significantly different from 0.0 (P > 0.01) under fluid milk market prices. Partial regressions for production traits in Table 4
reflect the differences in market values and costs of production, and the genetic co-variances between all traits included as independent variables. Partial regressions for type composites were generally similar in size and sign for all analyses, but attempts at biological interpretation of the results are probably unjustified. Weigel et al. (1997) showed that linear type traits explained only small amounts of variation in profit functions calculated similarly to the ones used in this study.
|
Table 5
contains results of prediction of RNIOC on individual Holstein cows using different evaluations of lifetime economic merit from genetic evaluations on their sires. We used the least squares solutions from odd numbered herds to predict RNIOC on sires of cows in even numbered herds and then tested that prediction of lifetime economic merit on the sire by predicting actual RNIOC calculated for the cows. We also used the Net or Fluid Merit index of each cow for the same purpose. For instance, the regression of $1.03 in the first cell indicates that RNIOC calculated using manufactured milk prices on cows in odd numbered herds increased by $1.03 for each $1.00 increase in predicted RNIOC for manufactured milk prices when using the prediction equation developed using cows in even numbered herds. R2 for that prediction was 0.091. The regression coefficients in Table 5
were closer to $1.00 when RNIOC using manufactured milk market prices was predicted from an equation based on the same economic assumptions and vice versa. Regressions using Net or Fluid Merit were always less than $1.00. If these lifetime economic merit functions include and properly value all important sources of income and expense in the life of a dairy cow, then the Merit indexes overstate absolute profitability of progeny.
|
| CONCLUSIONS |
|---|
|
|
|---|
The lifetime profit functions we used are based on assumptions of costs of specified inputs and value of specified outputs. Other costs, reductions in costs, and sources of income, such as economic value attributable to lower SCC, that are widely measured on many cows could have been used with equal justification. All such assumptions affect results in this paper. Likewise, USDA weights for the Merit indexes are based on current research information that may change in the future. Results of this study demonstrate the robust nature of the assumptions on which the present Merit indexes are based. Results show that farmers can use the Merit indexes for sire selection decisions with confidence that the aggregate genotype favored is consistent with and effective in the improvement of a breeding objective of improved lifetime economic merit.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication April 10, 2002. Accepted for publication June 24, 2002.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. McDougall and C. Compton Reproductive Performance of Anestrous Dairy Cows Treated with Progesterone and Estradiol Benzoate J Dairy Sci, July 1, 2005; 88(7): 2388 - 2400. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. de Vries, C. Steenholdt, and C. A. Risco Pregnancy Rates and Milk Production in Natural Service and Artificially Inseminated Dairy Herds in Florida and Georgia J Dairy Sci, March 1, 2005; 88(3): 948 - 956. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |