|
|
||||||||
1 New Bolton Center, Center of Animal Health and Productivity School of Veterinary Medicine, University of Pennsylvania, Kennett Square 19348
2 Animal Breeding and Genetics Group, Wageningen University, The Netherlands
Corresponding author: H. Groenendaal; e-mail: huybert{at}risk-modelling.com.
| ABSTRACT |
|---|
|
|
|---|
Key Words: replacement economic model retention pay-off cost per extra day open
Abbreviation key: CI = calving interval, DO = days open, DP = dynamic programming, MNR = marginal net revenue, RPO = retention pay-off
| INTRODUCTION |
|---|
|
|
|---|
On dairy farms, the most observed reasons for culling cows are reproductive problems, low production, and mastitis (Morris and Marsh, 1985; Van Arendonk, 1988). Although culling decisions are of great economic importance for a dairy farm, they are often made in a non-programmed fashion and based partly on the intuition of the decision-maker (Lehenbauer and Oltjen, 1998). To improve expected future profits on the dairy farm, culling decisions should be based on economic principles rather than on biological considerations (Dijkhuizen, 1983; Lehenbauer and Oltjen, 1998). Economic analysis of the replacement decision should include the expectation of the cows future performance as well as that of the potential replacement (Dijkhuizen, 1983).
To evaluate decisions on dairy cattle breeding and replacement, 2 main techniques, marginal net revenue (MNR) and dynamic programming (DP), have been applied. Burt (1965) stated that the MNR approach is, in fact, a special case of DP. Both techniques rely on the production function approach in which the economic costs and revenues of a cow are modeled during her life span (Van Arendonk, 1985a). The 2 main differences between the MNR and the DP approach (and limitations of the MNR approach) are 1) the DP approach can take into account the variation in expected performances of both present and subsequent replacement, and 2) DP can take into account genetic improvement. Because DP can overcome both limitations of the MNR approach, many researchers (Giaever, 1966; Van Arendonk, 1985a, DeLorenzo et al., 1992; Jalvingh, 1993; Kristensen, 1993; Houben, 1995) have used DP techniques to provide guidelines for replacement and breeding decisions.
A problem with the DP technique, however, is that DP models can easily become very large and complicated depending on the number of states defined, incurring the risk of limited breadth of application because of intensive resources requirements (Smith et al., 1993). Although more efficient DP models have been developed (Kristensen, 1993), most of the DP models that have been developed so far are relatively complicated and need high computer skills to use. In addition, most DP models are compiled with many fixed parameters, thus limiting the number of parameters that could be changed by the user. Furthermore, most of the existing models are not very user friendly and have interfaces that are unfamiliar to the end users. The majority of effort on decision-supporting models in the dairy industry has been focused on constructing models (Van Arendonk, 1985a; Kristensen, 1993; Houben, 1995) rather than on using models as applied decision-making tools. This fact is illustrated by the observation that, to our knowledge, none of the existing models are directly available. As a consequence, little progress has been made at the farm level in making better culling decisions (Lehenbauer and Oltjen, 1998).
For use as a decision-supporting tool on the dairy farm, a model should be simplified as much as possible without compromising the accuracy of outputs. For that reason (Van Arendonk, 1985a) and for the possibility of structuring a replacement model in a spreadsheet program that is familiar and easily available to the end users, the MNR approach sometimes can be justified. An often-cited limitation of the MNR approach (Van Arendonk, 1985a; Kristensen, 1993) is its inability to easily account for genetic improvement. However, Van Arendonk (1985a) concluded that genetic improvement hardly affected the optimal breeding and replacement policy. A second limitation of the MNR approach is that it does not easily take into account the variation in the expected performance of present animals. However, as the expected performance of a replacement heifer, which is used in the MNR approach, is equal to the average of the probability distributions that are used in the DP model, both models are likely to give very similar results. In other words, although DP takes into account variation, the optimal decisions in both methods are based on the expected performances of the animals in the herd.
The goal of the current paper is to describe a spreadsheet dairy cattle replacement model. With the model, optimal replacement and breeding decisions can be supported for cows with different production characteristics. With the model, the costs per additional day open (DO) for cows with different production characteristics can also be calculated. The model is based on the MNR approach and is user friendly in that it allows users to easily change all input parameters under different production and economic situations on dairy farms.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
![]() | ([1]) |
where
| ANRj | = | annuity of the net revenue of a replacement animal per month;
| i | = | decision moment of replacement (1 I j) at the end of period i;
| j | = | period, at the end of which an animal can be replaced;
| r | = | discount rate per month;
| pi | = | probability of survival until the end of period i, calculated from the moment at which a young animal starts its first production (end of period 0);
| mi | = | length of period i (mo); and
| MNRi | = | MNR in period i, including a change in slaughter value and financial loss associated with involuntary disposal.
|
In the current application of the model in this paper, the replacement animal was assumed to have been bred to an average calving interval (CI) and kept until her maximal annuity value was reached.
The RPO value of a cow is equal to the total additional profits that a producer can expect from trying to keep the cow until her optimal age, taking into account the changes of involuntary premature removal compared with her immediate replacement (Huirne et al., 1997). In other words, the RPO value represents the maximum amount of money that could be spent to try to keep an animal in case of reproductive failure or health problems. The RPO value assumes that the only opportunity, other then keeping the cow, is a replacement heifer. Therefore, the maximization of net revenue per cow-place per year in the long run is the objective, and the opportunity costs have to be included in the calculation. The RPO value was calculated as follows (Huirne et al., 1997):
![]() | ([2]) |
where
| RPOi | = | RPO at decision moment i;
| d | = | optimal moment for replacement (when MNRj < ANRmax);
| r | = | discount rate per month;
| pj | = | probability of survival until the end of period j, calculated from decision moment i;
| j | = | period, at the end of which an animal can be replaced;
| mj | = | length of period j (mo);
| MNRj | = | MNR in period j; and
| ANRmax | = | expected maximum average net revenue per month.
|
In the model, the RPO values can be calculated for different production levels, which were defined relative to the herd average milk yield. The herd average milk yield in this application was set at 9072 kg (20,000 lb)/yr per cow. In the current models application, the RPO values of 5 production levels (76, 88, 100, 112, and 124% of average) were calculated, but any level can be used. To capture different reproductive efficiencies, the model calculates the RPO values for cows with CI of 11, 12, ...17 mo. The CI in the following lactations can take on a variety of levels, including the overall herd average (set as default). The herd average CI, without confounding by culling, was calculated by:
![]() | ([3]) |
where
| CIaver. | = | average CI (mo),
| RND | = | rounding function,
| VWP | = | voluntary waiting period (d),
| EDR | = | estrus detection rate (%),
| CR | = | conception rate (%), and
| LP | = | length of pregnancy (assumed to be 274 d).
|
In the current application, the average CI (CIaver.) was calculated at 15 mo (for input parameters, see Table 1
). This average CI was used as the expected CI of replacement heifers.
|
![]() | ([4]) |
where
| CDOCI | = | costs per DO ($),
| CI | = | CI (mo),
| RPOCI | = | RPO value in first month of a lactation with a CI measured in months, and
| RPOCI + 1 | = | RPO value in first month of a lactation with a CI of CI + 1 mo.
|
To calculate the effect of a difference in CI on the RPO values, the RPO values in the first month of lactation (between calving and the end of the voluntary waiting period) were used. One could also choose a later month during the lactation (as long as this month is before the month of conception of the shortest CI) to compare RPO values and calculate the costs per extra DO. To keep the model consistent for different CI, the first month of the lactation was chosen. To get more insight in the costs per extra DO both the costs with and without opportunity costs (of a replacement heifer) were calculated. The RPO values without opportunity costs represent a situation where the RPO value is equal to the total net present value of the current cow (representing a situation where no replacement is available), without confounding by voluntary replacement.
Input
General.
The spreadsheet model (http://cahpwww.vet.upenn.edu/software/econcow.html) was developed using Excel 2002 with Visual Basic for Applications (Microsoft, Redmond, WA). The user can customize the model by changing all input parameters and variables to calculate results for specific herd situations. In the model, the performances, revenues, and costs of individual dairy cows were calculated for each month (30.4 d) to allow replacement at regular intervals within the lactation period. All future costs and revenues were discounted at a 5% discount rate (Brealey and Myers, 2000). Previous work found that milk price, feed price, milk production, replacement costs, and carcass prices have the largest influence on optimal replacement decisions (Van Arendonk, 1985a). Therefore, to keep the number of input parameters low and the model simple, only these and a few other important parameters were included (Table 1
). The default parameter values for this papers application were chosen to represent a typical dairy farm in Pennsylvania.
Revenues.
Three lactation curves are currently available in the model (the user chooses the curve with a pop-down menu), but other equations can easily be included. The lactation equation used in this papers application was developed by Oltenacu et al. (1981), adapted by Marsh et al. (1988), and later modified by Skidmore (1990):
![]() | ([5]) |
where
| Y | = | daily milk yield (kg),
| A | = | ((GNRHA/100 a)/2.96),
| GNRHA | = | rolling lactation average (genetic rolling herd average) (kg/yr),
| DIM | = | days in lactation (milk),
| DP | = | days in gestation (days pregnant),
| e | = | base of natural logarithm, and
| a, b, c, g | = | constants that determine the shape of the lactation curves.
|
For lactation numbers 1, 2, and
3, different constants for the coefficients a, b, c, and g (Table 2
) were used (Skidmore, 1990). The constant g, which determines the effect of gestation stage on milk yield, was modified to reflect a 200- and 350-kg cumulative milk yield decrease in the 305-d production of first and second and higher lactation animals (Olori et al., 1997). Prediction of future milk production levels was done, assuming a repeatability of 0.55 for the next lactation and 0.50 for all lactations afterward (mean reversion), similar to the method used by Van Arendonk (1985b). Whereas the current application of the model used this formula for generating lactation curves, in other situations other curves might be more applicable.
|
![]() | ([6]) |
where
| CVi,j | = | carcass value of cow i in month j,
| LWi,j | = | live weight (kg) of cow i in month j,
| D%j | = | dressing percentage (%) in month j,
| P | = | average price per kilogram of carcass weight for a heifer 210 d in lactation ($/kg), and
| dpj | = | price in month j as a deviation from the average pj ($/kg).
|
The effect of lactation number and stage of lactation on dressing percentage and the price per kilogram of carcass weight is given in Table 3
. This price was expressed as a deviation from the price of a heifer at 7 mo in lactation, which was taken at $1.52/kg (Pennsylvanian Agricultural Statistic Service, 1997).
|
Costs.
In the model, the rearing of young stock was isolated from other activities under the assumption that pregnant heifers are purchased from the farms rearing enterprise or through the market. The costs of replacement heifers for this papers application were calculated by interpolation data from the Cornell Cattle System 4 model (Van Amburgh and Fox, 1996). These data imply that, between reasonable biological limits, the total costs per kilogram of weight gain are lower with a lower age at first calving. The total costs of a replacement heifer calving at an age of 26 mo were, therefore, calculated at $1132, which was consistent with the average costs to raise a replacement heifer of $1124 found by Gabler et al. (2000). However, the user can include costs of replacement animals at any price without using the interpolated data.
Feed costs were calculated by multiplying the DMI (kg) of each individual cow with the costs per kilogram of DMI. The DMI was calculated using the following formulas (Galligan et al., 1985).
For lactating cows:
![]() | ([7a]) |
For dry cows:
![]() | ([7b]) |
where
| DMIi,j | = | DMI (kg) for cow i in month j,
| BWi,j | = | mature BW (kg) for cow i in month j, and
| Pi,j | = | milk production (kg) for cow i in month j.
|
Although this model can account for a variety of feeding systems, a one-group TMR was used in this papers application. The costs of DMI in lactating cows were set at $0.20/kg (US $0.09/lb). The feed costs for growing heifers and dry cows were $0.15/kg of DMI ($0.07/lb) (Galligan et al., 1985).
To calculate the breeding costs for different CI, the estrus detection rate was set at 40%, and the conception rate was set at 40%, close to what has been observed in the field (Smith et al., 1988; Lucy, 2001). The voluntary waiting period was assumed to be 50 d. The breeding costs per month were calculated by multiplying the expected number of breedings per month with the insemination costs per breeding. The expected number of breedings per months was calculated as the potential number of breedings per month (1.45, which was calculated as 30.4 d ÷ estrous cycle of 21 d) multiplied by the estrus detection rate. The minimal number of total breedings per pregnancy was set at one.
Other costs included (Table 4
) were the costs associated with morbidity, disposal, and mortality. Typically, yearly veterinarian costs per cow were estimated at $50 for an average first lactation cow (Snow, 1993) and increased $5 each lactation. Van Arendonk (1985b) assigned 33% of these costs to the first month, 11% to the second and third months, and 5% to the later months of each lactation. The direct financial costs associated with mortality were set equal to the slaughter values.
|
|
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
For 5 different milk production levels, the calculated RPO values for typical Pennsylvania conditions are given in Table 6
, calculated for cows that have just become pregnant at 6 mo after calving (resulting in a 15-mo CI). The RPO values in Table 6
vary between $37 and +$1995; variability was mainly caused by the difference in milk production. A first lactation cow with a relative milk production of 76% has an RPO value of $37, which means that keeping her 1 additional mo instead of replacing her with an average replacement heifer (if available) would cost the producer $37. In contrast, if the producer would have to cull (involuntarily) a first lactation heifer that has a relative production level of 124% and replace her with an average heifer, he would have an economic loss of $1995.
|
|
124%. Finally, Houben (1995) found that not including within-lactation transitions results in an overestimation of high-producing cows and an underestimation of low-producing cows in the beginning of the lactation. As the current model allows for transitions between different levels of milk production only at the end of the lactation period, the models effect of milk production on RPO could be slightly overestimated.
Second, the RPO value of a cow is generally the highest just before calving. For cows with a 15-mo CI, the RPO is minimal around 7 to 9 mo after calving, depending on the milk production level. Thereafter, the RPO increases again (Figure 2
) because of the decreasing risk that she is involuntarily culled before calving, the increasing expected revenues of the next lactation, and because during the dry period the costs are becoming sunk cost.
Third, the maximum RPO values across lactations for high-producing cows gradually decline from lactation 1 to 12. An exception is that for lower-than-average-producing cows the maximum RPO values increase from lactation 1 to 2, but after lactation 2, the RPO values are also gradually decreased. The lower RPO value in the first lactation of low-producing animals was caused by the lower production in this lactation compared with later lactations (Equation 5
). For average- and higher-than-average-producing cows, lower production in the first lactation was offset by the high relative production level of the cow (
100%) compared with an average replacement heifer. After reaching the maximum, the RPO value of cows gradually declines with a higher lactation number for 3 reasons. First, the time until optimal culling becomes shorter, and, therefore, the total extra profits of keeping the cow until this optimal time of replacement (= RPO value) decrease. The 2 other reasons for the decreasing RPO values are the higher involuntary culling and higher morbidity rates in later lactations.
Fourth, the RPO of cows with a lower production level decreased <0 around 6 to 7 mo in lactation, depending on lactation number and reproductive efficiency. For example, during the third lactation, under a 15-mo CI, the RPO value of a cow at 76% production goes <0 at 7 mo, which means that this cow should be replaced after 7 mo in milk. If the RPO value of a cow is negative, replacing her is economically more attractive than keeping her. However, if the producer decides to keep the cow and successfully breeds her, her RPO value will increase again >0 around 4 mo before calving. Beyond that point, the cow should be kept again until the next optimal time of replacement when the RPO again goes <0.
The calculations of the opportunity costs (ANRmax; see Equation 1
) were based on the average performance of animals present in the herd, assuming this to be the best estimates for future net revenue of young replacement animals. Three underlying assumptions of the results shown above are 1) no genetic gain, 2) unlimited availability of identical replacement heifers, and 3) constant number of cows in the herd. Genetic gain would result in lower RPO values because the opportunity costs are higher (i.e., replacement heifers are on average better than the current average animals). However, Van Arendonk (1985a) showed that genetic gain had only a very small influence on the optimal breeding and replacement decisions. Secondly, opportunity costs are 0 if there are (temporarily) no replacement heifers available. Replacements on dairy farms are often dictated by the calving of new heifers (Kristensen, 1993), and then naturally the least profitable cows should be replaced. Such situations often arise in herds that only use homegrown heifers, which is a very common policy in Pennsylvania. Opportunity costs can also be 0 (temporarily) when circumstances of the dairy farm allow heifers to be added without requiring existing cows to be culled. This situation can occur because of planned long-term expansion or short-term fluctuation in cattle numbers (Lehenbauer and Oltjen, 1998). In both situations when opportunity costs are 0, current cows can be kept as long as her MNR is >0 (Huirne et al., 1997). However, in any situation, the ranking of cows according to their RPO values selects the least profitable cow, and if a heifer is available, this cow is replaced (Kristensen, 1993). Therefore, the ranking of cows is more important than the absolute RPO values (Kristensen, 1993). In addition, the ranking of RPO values is far less sensitive to changing opportunity costs and changing prices and interest rates than the absolute RPO values.
Costs per Extra DO
The model calculated the RPO value of cows that have different reproductive efficiencies as measured by differences in their CI. First, Figure 3
shows the RPO values of cows in the first month of the third lactation without and with opportunity costs. A comparison of the situation with and without opportunity costs reveals 2 main differences. First, the RPO values without opportunity costs are considerably higher than the RPO values with opportunity costs. This is caused by 2 factors. First, the RPO values with opportunity costs represent the total extra profits of keeping the cow compared with replacement. Therefore, the opportunity costs (expected maximal average revenues of the replacement heifer) are subtracted from the MNR of the present cow (see Equation 2
). In contrast, the RPO values without opportunity costs represent the extra profits of keeping the cow compared with an open place in the barn (zero opportunity costs subtracted). Second, the age until optimal replacement will be lower for situations where a replacement heifer is available than for situations where no replacement heifers are available. As a consequence, the total profits of keeping the cow were lower. The second difference between both situations is the slope. In the situation in which no opportunity costs were included, the graphs always had a negative slope (future value of the cow became lower); in the situation with opportunity costs, depending on the milk production, beyond a certain CI, the slope becomes zero. The reason is that if the CI increases, postponing breeding of the cow will not always decrease the RPO value of the cow because of the models assumption of economical optimal culling. Because the model calculates the total benefits of keeping the current cow until the optimal time of replacement, it does not take into account losses that occur after this optimal time of replacement. As a consequence, a longer CI (beyond the maximum interval that is allowed to still be economically attractive) will not result in a decreased RPO in the model and subsequently will not result in losses per additional DO (in other words, the cow should not be bred at all).
|
|
|
Costs per extra DO vary across the lactation for animals in different lactations (first, second, and third) or with a low (76%), average (100%), or high (124%) milk production levels (Figure 6
). Costs per extra DO are lower and increase slower for first lactation animals than for animals in the second and higher lactation. This effect is caused by higher persistency of milk production of first lactation animals (Skidmore, 1990). Also, for all 3 lactations, the costs per extra DO are higher for low-producing animals than for higher producing animals. An exception to this, are animals (for example first lactation animals that produce 76% of the average) that should not be bred because breeding would not increase their RPO value and, therefore, is not optimal. Hence, the costs per DO will be calculated as $0 (Equation 4
) for these animals (Figure 6
). Finally, the difference between an extra DO for average- and high-producing animals in the third lactation is small and smaller than the difference with first lactation animals (Figure 6
). This result is in agreement with Dijkhuizen (1983), Van Arendonk (1985b), and Strandberg and Oltenacu (1989).
|
Application on Dairy Herds
The RPO values of 5 dairy cows are shown in Table 7
. In addition to the herd data that are shown in Table 1
, the cow data that are needed to calculate individual RPO values include the cows lactation number, the current milk production per day, the number of DIM, and the number of days pregnant. To determine the cows milk production level as a percentage of the herd average mature equivalent milk, the inverse of the modified Oltenacu lactation curve (Equation 4
; Oltenacu et al., 1981) was used. With these data, the model determines the individual RPO values and, if applicable, costs per additional DO.
|
| CONCLUSIONS |
|---|
|
|
|---|
With the current spreadsheet model, based on the MNR approach, the dairy cow breeding and replacement problem was modeled accurately, and optimal replacement and breeding decisions can be supported. The results are in close agreement with former studies (Van Arendonk, 1985a; Jalvingh, 1993; Kristensen, 1993; Houben, 1995). The strength of the currently described model is the integral evaluation of age, production, reproductive efficiency, and survivability in a simple and user-friendly economic computer spreadsheet model to support replacement and insemination decisions on dairy farms. Strandberg and Oltenacu (1989) concluded that there are no magic numbers for the optimal breeding and replacement decisions nor for the losses per marginal DO that apply to all herds and cows. Rather, there is a need of more customized breeding decisions for each (type of) cow that are herd specific. The current model can help in this need by providing user-friendly input and output to customize the calculations for individual herds and cows. On-farm, the model can support optimal decisions regarding voluntary replacement and breeding decisions. In addition, it can determine the on-farm costs associated with involuntarily culling or with additional DO. Other potential users of the current model are researchers, economists, and governmental organizations that wish to calculate the (farm or cow-specific) losses of involuntarily culled dairy cows because of a particular disease (Van Schaik et al., 1996) or as part of a specific disease control program (Groenendaal et al., 2002). Finally, users of the model can obtain estimates on the farm-specific costs of an extra DO for cows with different production characteristics, which can be useful to calculate the economics of specific insemination policies.
In summary, the current user-friendly model determines the economic value of individual dairy cows under farm-specific circumstances and uses a new approach to calculate the costs per extra DO. The results of the model are very similar to results of more complex models that are more difficult to use. Therefore, we consider the model a valuable tool for dairy farms to support farm and cow-specific optimal breeding and replacement decisions. Although not shown in this paper, the model can also be useful to assess economic costs of involuntary culling under disease control programs.
Received for publication June 5, 2003. Accepted for publication February 12, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
U. S. Sorge, D. F. Kelton, K. D. Lissemore, W. Sears, and J. Fetrow Evaluation of the Dairy Comp 305 Module "Cow Value" in Two Ontario Dairy Herds J Dairy Sci, December 1, 2007; 90(12): 5784 - 5797. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. De Vries Economic value of pregnancy in dairy cattle. J Dairy Sci, October 1, 2006; 89(10): 3876 - 3885. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |