|
|
||||||||


* Department of Agricultural Economics, University of Wisconsin-River Falls, River Falls 54022
Department of Agricultural Economics, Michigan State University, East Lansing 48824
1 Corresponding author: wolfch{at}msu.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: culling economics management
| INTRODUCTION |
|---|
|
|
|---|
High replacement costs and low milk prices tend to focus attention on culling decisions and optimal culling rates, as high replacement costs during periods of small profit margins make nonoptimal culling rates expensive. Past research has consistently estimated optimal herd-level culling rates ranging from 19 to 29% (Rogers et al., 1988; Bauer et al., 1993; Stott, 1994; Jones, 2001). Actual average annual culling rates are often above even this range. For example, the average annual culling rate for Upper Midwest DHIA herds in 2001 was 38% (Quaiffe, 2002). Because many dairy advisors, farm organizations, and agricultural lenders feel that the discrepancy between estimated optimal culling rates and observed culling rates is too large, dairy farmers are often advised to reduce culling rates.
Holding herd size constant, dairy cow replacement decision models instruct the manager to cull an incumbent cow when the challenger replacement cow is expected to generate more profit (Dijkhuizen et al., 1985). Without a challenger, the incumbent cow is appropriately culled when she is no longer profitable. Traditionally, an examination of culling patterns has included "voluntary" and "involuntary" culls. Culling due to low production, cow aggression, or when a cow is sold to another farm for dairy purposes is referred to as voluntary culling. Involuntary culling is all other factors including sales due to illness, injury, infertility, or death. If all culling decisions are voluntary, the manager can make decisions that maximize profit. Thus, it is important to understand how culling affects production, how many cattle are culled, why cattle are culled, and how management programs affect culling rates.
There are several examples of research on culling reasons, which are useful for identifying explanatory variables and comparing results. For example, Milian-Suazo et al. (1989) examined 34 New York dairy herds to explore the influence of production and disease on culling. They found that lower milk per calving-interval day, shorter previous lactation, and summer calving were associated with those cows sold for dairy purposes. Bascom and Young (1998) summarized reasons why farmers culled cows on 27 herds in New Hampshire, Vermont, and Massachusetts. They found that reproduction was the primary reason for culling, followed by production and mastitis. Higher producing herds were more likely to cull a cow for abortion and reproduction, but less likely to cull for mastitis. Smith et al. (2000) examined 11,250 herd summary records from states east of the Rocky Mountains to determine the effect of region, herd size, and milk production on cows leaving the herd. Low-producing herds reported a lower percentage of cows leaving than higher producing herds. High-producing herds reported more cows leaving for reproduction, mastitis, feet and legs, and disease than lower producing herds. Small herds reported more cows leaving for reproduction and mastitis than larger herds. Past research has largely focused on culling factors and patterns at the herd level. In contrast, we examine culling at a cow level.
We analyze DHI records for 5 Upper Midwest states and 5 Northeast states from 1993 through 1999 to determine the percentage of cattle culled each year, why cattle were culled, and when cattle were culled within a lactation. Probit regression models were used to determine how individual cow and herd characteristics affected the likelihood of an individual cow being culled.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Data
Data were supplied by Dairy Records Management Systems (DRMS, Raleigh, NC). The data contained production records of participating DHI dairy farms in the Northeast region [Maine, New Hampshire, Pennsylvania, and Vermont (1993 to 1999); and New York (1996 to 1999)] and the Upper Midwest region [Illinois, Iowa (19961999), Indiana (19931999), Michigan (19971999), and Wisconsin (19951999)] from 1993 through 1999 for individual cow records and 1995 through 1999 for herd level records in participating states. Data from 1995 through 1999 were used to capture the largest possible number of herd-related characteristics. The total numbers of herds with complete records for 1995 were 6,264 (about 9% of all operations with milk cows in those states; USDANational Agricultural Statistics Service, Milk Production Report). Other years and corresponding herd numbers with records were 14,028 herds in 1996; 18,290 herds in 1997; 20,132 herds in 1998; and 19,464 herds in 1999. These values represented approximately 22, 31, 33, and 35% of operations with milk cows in those states for 1996 through 1999, respectively (USDA-NASS, Milk Production Report). Herds could enter or exit at any time; thus, the data are an unbalanced panel from a herd perspective. We controlled for this in the estimation with a fixed-effects variable for herd. There were 7,087,699 individual cow lactation observations in the original data set.
Explaining Cow Cull Probability
For each cow and decision period, the decision-maker has a choice: the cow can be culled or retained. With a discrete, binomial choice, a nonlinear probability regression model was used for analysis. Specifically, a probit model was used to predict the probability that the cumulative effect of the explanatory variables had exceeded the threshold value for a cow to be culled. The probit model takes the form:
![]() |
where P(Cull = 1) refers to the predicted probability (0 or 1 for a given cow and lactation) that a culling event occurs in that lactation, ß0 is the intercept, and ß1, ..., ßn are coefficients relating to x1,..., xn, which are a set of independent variables that may influence whether a cull will occur. The marginal effect of the nth parameter in a probit model is calculated as ßn
(Yi), where ßn is the coefficient of the nth parameter,
is the standard normal density function, and Yi refers to the regression result at mean values (Gujarati, 1995). The SAS statistical software (SAS Institute, 2002) was used to calculate coefficients and resulting marginal effects. The probability threshold of a positive event was set at 0.6.
A variety of individual and herd-level characteristics contribute to a cow being culled. From the previous literature (Dijkhuizen et al., 1985; van Arendonk and Dijkhuizen, 1986; Rogers et al., 1988; Milian-Suazo et al., 1989; Delorenzo et al., 1992; Bascom and Young, 1998; Smith et al., 2000; Weigel et al., 2003), we identified a number of characteristics that were expected to influence the culling decision. Factors that influence culling decisions were categorized as relating to cow, herd, location (state or region), or time. Cow-related characteristics included calving season: spring (March, April, and May), summer (June, July, and August), fall (September, October, and November), or winter (December, January, and February), of the current lactation, lactation number (quadratic), milk production (kg), fat (%), protein level (%), persistency (%), as well as SCC score relative to the herd average, services per conception relative to the herd average, and breed. Persistency percentage is an index based upon each cows current projected 305-d mature equivalent (305ME) production compared with her projected 305ME production on the previous test day. A persistency value greater than 100 indicated that the lactation curve of the cow was more persistent than the average for her breed, age, and season of freshening. The formula to calculate persistency percentage is: (new projection/last months projection) x 100 (DRMS, 1999). The SCC difference measure was calculated by subtracting the herd average SCC in cells/mL from the individual cows SCC. Herd characteristics included whether the herd had breed-registered cows, herd size, whether the herd was in or recently had a herd expansion, and proportion of replacement heifers to milk cows. Two groups of 5 expansion dummy variables were designed to test the hypothesis that expansion was positively correlated with the likelihood of a cow being culled. The first expansion group, "small" expansions, included herds that increased herd size by between 20 and 300%. Dummy variables for this expansion group represented the first through fifth years during or following a herd expansion. The second group of expansions, "large" expansions, consisted of those farms that expanded by more than 300%. A dummy variable was also assigned to represent the first through fifth years during or following a large expansion. An explanatory variable, "heifer ratio current or lagged" representing the ratio of the number of heifers over 13 mo of age (corresponding to cattle that can be expected to calve within 10 to 13 mo) to the average annual milking and dry cow herd size (heifer ratio), was used to determine the effect that the number of available replacements had on the likelihood of a cow being replaced. One current period and a 1-yr lagged ratio were used. State level milk-to-feed and cull cow-to-replacement heifer price ratios for each month were included. The milk-to-feed price ratios were constructed from USDA reports for each state by month. This price ratio represents number of pounds of 16% mixed dairy feed equal in value to 1 pound of milk. All milk prices were from the Agricultural Prices Monthly publication (USDA-NASS). A representative feed price was constructed using prices of corn, alfalfa hay, and soybeans obtained from Agricultural Prices Monthly (USDA-NASS). In a few cases, prices were not available for a particular state (for example, a New Hampshire milk price). In these cases, the closest available neighboring state price was utilized (Vermont milk price for New Hampshire). Similarly, cull cow-to-replacement heifer price ratios for each month by state were calculated from Agricultural Prices Monthly data (USDA-NASS). Both of the price ratios are proxies for the contemporaneous dairy economic conditions under which a farmer made culling decisions. A state dummy variable captured other characteristics such as weather and geography unique to that state. A year variable captured potential time effects. Finally, a fixed-effect dummy variable was included for each farm observation to capture management and other unobservable characteristics. These characteristics are summarized in Table 1
and discussed in detail below.
|
![]() |
where s represents months prior to the current month for the variables that enter with a lag and r is the number of years following an expansion. Thus, the estimation is a cow model with explanatory variables that capture herd- and cow-specific characteristics. The farm fixed effect (ß2) captures characteristics unique to that operation that are not captured by other exogenous right-hand side variables. The estimated farm fixed effects are not reported below because there is one for each herd and individual values are only of relevance to the operation in question.
Individual cow data was available from 1995 through 1999. The estimation was performed for each region to examine how each parameter differed between the 2 regions rather than having the regional differences absorbed by each state-independent variable. Using data from Upper Midwest DHIA dairy farms in Illinois, Indiana, Iowa, Michigan, and Wisconsin, there were 1,126,610 useable observations (cow lactations) in which 335,711 culling events and 790,899 nonevents occurred. The estimation had an R2 of 0.28. The model was able to successfully classify 79.5% of the Upper Midwest culling events and nonevents.
To calculate the marginal effects of the independent variables on the probability of a cull, the model was designed to provide a culling probability for a baseline situation. The baseline situation for the Upper Midwest estimation was a first-lactation Holstein cow calving in Indiana in the spring of 1995. The farm from which this representative Upper Midwest cow came was equal in size to the sample average herd size of 193 cows and had not undergone an expansion before or during the 1995 through 1999 period. The representative farm had a current year heifer ratio that was equal to the sample average of 0.50. For the baseline Upper Midwest situation, there was a 26.1% chance that a first-lactation, average-producing Holstein cow would be culled during or at the termination of the first lactation on that particular farm in Indiana. The marginal effects that are reported below should be interpreted as the change in probability of a culling when that factor changes relative to the baseline cow all other factors held constant.
Using data from Northeast DHIA dairy farms in Maine, New Hampshire, New York, Pennsylvania, and Vermont, there were 225,987 culling events and 106,339 culling nonevents that took place. The probit model correctly predicted 79.9% of the cull decisions. The estimation had a R2 of 0.28. The baseline situation for the Northeast estimation was a first lactation Holstein cow calving on a particular farm in Vermont in the spring of 1995. The farm from which this representative cow came was equal in size to the sample average herd size of 200 cows and had not undergone an expansion before or during the 1995 through 1999 period. The representative farm had a current year heifer ratio that was equal to the sample average of 0.48. For the baseline Northeast farm situation, there was a 36.8% chance that a first-lactation, average-producing Holstein cow would be culled during or at the termination of the first lactation on a particular farm in Vermont.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Northeast state herds had average cull rates ranging from 29 to 31% (Table 2
). This average was lower than Upper Midwest region herds, which averaged 33 to 35%. Consistent with results of Smith et al. (2000), average annual culling rate increased with rolling herd average milk production. The average annual culling rate reached a plateau for herds producing more than 9,545 kg of milk at just over 32%. Average culling rate increased with herd size. Farms with a herd size of less than 150 cows had an average annual culling rate of 31.4%. The largest size group, herds with more than 600 cows, had the highest average culling rate of 36.8%.
|
|
|
|
Culling Probability Estimation Results
The data set in this study was large and the standard errors were small resulting in significant effects of many independent variables. Statistical significance does not necessarily translate to practical, economic significance and the results must be interpreted with care. Managers and advisers must consider what can be changed, what the costs and resource use will be, and what the potential benefit could be. A small percentage increase in risk of culling may be economically significant in a management sense, if a quick and inexpensive correction is available. The results presented here provide benchmarks and guidelines for producers.
Table 6
displays the coefficient estimate results for the culling probability model estimated with data from the Upper Midwest region. Table 7
shows the same estimation estimated with data from the Northeast region. As discussed above, the marginal effects should be interpreted as the change in probability of a culling when that factor changes relative to the baseline cow all other factors held constant.
|
|
Lactation number captured cow age, whereas lactation squared was included to capture the potential that the likelihood of being culled may have increased at a decreasing rate. In this study, each lactation subsequent to the first increased the probability of a culling more than 20% in the Northeast region and more than 30% in the Upper Midwest region. This increase in likelihood occurred at a decreasing rate as the lactation number was positive, whereas lactation squared was negative. Past research has shown that the incidence of culling increases with age. Dutch data used by van Arendonk and Dijkhuizen (1986) confirmed that the probability of an animal being culled due to health reasons (injury, disease, SCC, and reproduction) increased with age.
The herd production level of Dutch dairy farms (van Arendonk and Dijkhuizen, 1986) and United Kingdom dairy herds (Stott, 1994) did not affect the optimal cow longevity and culling rate. Rogers et al. (1988) found that larger milk yields supported only slightly higher culling rates for dairy farms in the United States. Weigel et al. (2003), however, reported a positive correlation between involuntary culling of high-producing cows, defined as cattle within the top 20% of their herd in terms of milk production. Based on these past models, the relative productivity of an individual cow compared with its herdmates should affect the probability of whether that cow is culled. Because US milk producers are paid based on milk components, differences in fluid milk production, butterfat production, and milk protein production were utilized to examine milk production effects. The 305ME milk, fat, and protein are production measurements that standardize production per cow so that the values represent the production as if all of the cattle were the same age, from the same location, milked twice a day, and had calved during the same season (DRMS, 1999).
As expected, the difference in milk production between a cow and its herdmates significantly affected the probability of a cow being culled. A cow that produced one additional hundredweight more than the average 305ME milk yield was 1.7% less likely to be culled than the average producing cow in the Upper Midwest region and about 0.5% less likely to be culled in the Northeast region. Cattle whose average difference in 305ME protein was 0.5 kg higher than the average producing cow were less likely to be culled in both regions. It is logical that cattle producing more milk and more protein than their herdmates are more profitable, all else being equal, and should remain in the herd longer than less profitable cattle. Counter to expectations, higher fat production increased the likelihood of a cow being culled, but the effect was small. This may result from managers trying to emphasize milk and protein production when making decisions. Prior to multiple component pricing, dairy farmers primarily were paid on fluid milk, butterfat, and SCC. Currently farmers are paid for fluid, butterfat, milk protein, SNF, and SCC. Thus, managers may now place more emphasis on components other than butterfat.
It was expected that more persistent cows have an increased probability of remaining in the herd as the less persistent cow will be less likely to cover variable costs in late lactation. Cows that were one persistency percentage higher than the average producing cow were about 2% less likely to be culled in both regions. Lactational persistency percentage is a DHI measure that compares the persistency of a cow with those of her herdmates. Dijkhuizen et al. (1985) determined that it became more advantageous to retain a long-open cow and continue breeding the animal as persistency increased.
Cattle with higher SCC score than the average were 1% more likely to be culled than the cow represented by the baseline situation. This is not surprising because high SCC scores decrease both the milk price per hundredweight as well as the amount of milk produced per cow (DRMS, 1999). In the United States, milk processors pay producers a premium for low SCC and discount prices for high SCC. As an individual cow SCC score increases relative to the rest of the herd, that cow becomes less profitable due to treatment costs, lost production, and quality discounts.
The number of services per conception in the previous lactation was associated with a decreased likelihood of a cull in the Upper Midwest (although the effect was small) region, but increased the likelihood of a cull in the Northeast region. Reproduction was the primary reason for culling cows (Bascom and Young, 1998). Cassell (2002) noted that cows that conceive more readily than their herdmates remain in the herd longer. Difference in services per conception from the prior lactation between a cow and her herdmates in the same lactation was used to represent individual cow reproductive efficiency. van Arendonk and Dijkhuizen (1985) hypothesized that a cow that had difficulty conceiving in the previous lactation may be viewed more critically during the present lactation.
The average first-lactation cow in Illinois, Iowa, and Wisconsin was less likely to be culled than similar cattle in Indiana, whereas a Michigan cow was more likely to be culled. The state differences were dramatic in the Upper Midwest region suggesting that infrastructure, economic, climatic, and geographic influences captured by the state variables were significant in explaining the likelihood of a cull. In contrast to the Upper Midwest results, the probability of a culling did not differ significantly across Northeast region states from the base Vermont case with the exception of Pennsylvania (where a culling was much more likely).
The effect of herd size on the likelihood of a cow being culled was significant, but small on a per-cow basis. Culling rates for herds with 600 or more cows were larger than the culling rates for herds with fewer than 100 cows, but that the difference was only 2% (Quaiffe, 2002). Weigel et al. (2003) found that high-producing cattle were more likely culled in larger herds due to herd health problems than in smaller herds.
The effect of expanding on the likelihood of culling was mixed. For the smaller category of expansions (an increase in herd size of less than or equal to 300%), the effect of the initial expansion year on the probability of a culling was negative, but insignificant for several years in the Upper Midwest. The time pattern for small expansions in the Northeast region was similar, but consistently more significant. For herds that had a large increase in herd size (a herd size increase of more than 300%), the first large expansion year was negatively correlated with culling likelihood, but often insignificant in both regions. Dairy farm expansion has been cited as a contributing factor increasing culling rates and high replacement heifer prices (Hoards Dairyman, 2003). Hadley (2001) found that dairy farms in Wisconsin and Michigan actually experienced a decrease in culling rates the first 2 yr of an expansion. Weigel et al. (2003) reported that expansion dairies were more likely to involuntarily cull high-producing cattle due to health problems, but less likely to cull low-producing cattle.
The heifer proportion from last year (i.e., lagged 1 yr) was included. For herds in both regions, the effect of the current year heifer ratio on the likelihood of a cow being culled was positive, but quite small. Thus, there was evidence suggesting that the number of heifers available had a small influence on culling rates. A manager may allow heifers to enter the herd through increased culling, selling the surplus heifers, or some combination of those actions. Radke and Lloyd (2000) asserted that many producers adjusted culling rates to accommodate the number of heifers available. If managers cull heavier to accommodate surplus heifers, this indicates that they believe it is more profitable to have the heifer enter the herd than to sell them as replacement heifers to other farms.
The milk-to-feed price ratio can serve as an important indicator of milk production profitability. As an increasing milk-to-feed price ratio represents an easier ability to cover variable costs, the marginal effect of each lactation months milk price ratio was expected to correlate negatively with the likelihood of a cow being culled. Renkema and Stelwagen (1979) concluded that the effects of milk price on optimal replacement decisions were small for Dutch dairy farms. Later work by Rogers et al. (1988) concluded that the same was true for US dairy farms. Bauer et al. (1993) found that the optimal terminal lactation for Alberta herds was not affected by changes in milk price. Jones (2001) found that the optimal longevity for Wisconsin dairy cows changed with milk price, but only in conjunction with opposite changes in replacement heifer prices. During periods of low milk price a producer may cull an animal without a replacement if the current financial return does not exceed the variable costs associated with production. This was the case in the majority of lagged price ratios by month in both regions.
It was expected that the marginal effect of the cull cow price to replacement heifer price ratio would be positively correlated with the likelihood of a cow being culled. As with the milk-to-feed prices, the cull-cow-to-replacement-heifer price ratio marginal effects were mixed being mostly positively correlated in the Northeast region with many of the parameters insignificant in both regions. Renkema and Stelwagen (1979) found that cull cow prices had a small effect on optimal culling policies for Dutch herds. Rogers et al. (1988) showed that replacement heifer prices had a large effect on culling decisions for US dairy farm managers. Bauer et al. (1993) estimated that changes in replacement heifer prices or a large decrease in cull cow prices should have large effects on the optimal terminal lactation for Alberta producers. Stott (1994) found that the optimal herd life for cows in the United Kingdom was sensitive to changes in replacement heifer prices. As the price farmers receive for selling cull cattle increases relative to replacement heifer price, it becomes less expensive from a capital investment perspective to cull and replace a given animal.
| CONCLUSIONS |
|---|
|
|
|---|
The largest marginal effects on likelihood of a cull were those related to lactation (cow age) and season. Summer and fall calving was associated with a smaller likelihood of a cull than spring calving. The milk market generally supplies a higher price in the fall months that encourages peak milk production in that season. Lower milk production and milk quality significantly increased the likelihood of culling the cow, although the effects were marginally small. Somatic cell count difference was the largest marginal effect and highlighted the continued importance of herd health programs. Expansion had a negative influence on the likelihood of a culling in the early years following the expansion. These results reflect the increased pressure to manage cash flow by keeping facilities at or near capacity following an expansion.
The study results are useful in describing patterns of culling and relating them to cow, herd, geographic, and time variables and can act as a benchmark for producers. Explaining discrepancies between farm management behavior and optimal culling decisions is an issue for future research.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication September 16, 2005. Accepted for publication January 11, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. W. Barkema, M. J. Green, A. J. Bradley, and R. N. Zadoks Invited review: The role of contagious disease in udder health J Dairy Sci, October 1, 2009; 92(10): 4717 - 4729. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. D. Dechow and R. C. Goodling Mortality, Culling by Sixty Days in Milk, and Production Profiles in High- and Low-Survival Pennsylvania Herds J Dairy Sci, December 1, 2008; 91(12): 4630 - 4639. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Bar, Y. T. Grohn, G. Bennett, R. N. Gonzalez, J. A. Hertl, H. F. Schulte, L. W. Tauer, F. L. Welcome, and Y. H. Schukken Effects of Repeated Episodes of Generic Clinical Mastitis on Mortality and Culling in Dairy Cows J Dairy Sci, June 1, 2008; 91(6): 2196 - 2204. [Abstract] [Full Text] [PDF] |
||||
![]() |
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] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |