JDS
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by McConnel, M. B.
Right arrow Articles by Galligan, D. T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by McConnel, M. B.
Right arrow Articles by Galligan, D. T.
J. Dairy Sci. 87:3542-3549
© American Dairy Science Association, 2004.

The Use of Integer Programming to Select Bulls Across Breeding Companies with Volume Price Discounts

M. B. McConnel and D. T. Galligan

University of Pennsylvania, School of Veterinary Medicine Center for Animal Health and Productivity 382 West Street Road, Kennett Square, PA 19348

Corresponding author: M. B. McConnel; e-mail: mmcconne{at}wharton.upenn.edu.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Optimization programs are currently used to aid in the selection of bulls to be used in herd breeding programs. While these programs offer a systematic approach to the problem of semen selection, they ignore the impact of volume discounts. Volume discounts are discounts that vary depending on the number of straws purchased. The dynamic nature of volume discounts means that, in order to be adequately accounted for, they must be considered in the optimization routine. Failing to do this creates a missed economic opportunity because the potential benefits of optimally selecting and combining breeding company discount opportunities are not captured.

To address these issues, an integer program was created which used binary decision variables to incorporate the effects of quantity discounts into the optimization program. A consistent set of trait criteria was used to select a group of bulls from 3 sample breeding companies. Three different selection programs were used to select the bulls, 2 traditional methods and the integer method. After the discounts were applied using each method, the integer program resulted in the lowest cost portfolio of bulls. A sensitivity analysis showed that the integer program also resulted in a low cost portfolio when the genetic trait goals were changed to be more or less stringent. In the sample application, a net benefit of the new approach over the traditional approaches was a 12.3 to 20.0% savings in semen cost.

Key Words: integer programming • semen selection • optimization • dairy cattle

Abbreviation key: LP = linear program.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Artificial insemination is an important technology used by the dairy industry. Surveys show that over 80% of dairy herds in the United States use at least some AI today (Plourd, 2000) in contrast to only 70% of herds in 1992 (Norman and Powell, 1992). Not only is AI being used by more herds, but it is also being used more within each herd. This change is probably due to the widespread adoption of systematic breeding programs where timed insemination is a component (Stevenson, 2001).

Like any important decision, semen selection (which ultimately defines the production levels and characteristics of future cows) must start with concrete goals. Dairy producers must decide which traits are particularly important to them and they must specify target values for each of these traits (Dematawea et al., 1998; Tozer and Stokes, 2002). Once these goals have been defined, a systematic method can be used for efficiently choosing semen.

Many computer programs have been developed to aid in the systematic process of semen selection. Most of these models use simple linear programming as their basis (McGilliard and Clay, 1983a, 1983b; Galligan and Ferguson, 1995; Weigel and Lin, 2000), but others use multi-objective linear programming (Tozer and Stokes, 2001), integer programming (Erba et al., 1991; Nash and Rogers, 1995) or quadratic programming (Schneeberger and Freeman, 1982). The objectives of semen selection models also vary, from minimizing the average cost of semen (Galligan and Ferguson, 1995; Eicker, 1995), while meeting or exceeding target trait values to maximizing the target trait values while remaining within a specified price range (McGilliard and Clay, 1983b). Still other models have been proposed that minimize risk (Schneeberger and Freeman, 1982) or that minimize inbreeding (Weigel and Lin, 2000).

Despite the differences between these selection programs, they share a common characteristic; all of these approaches ignore the impact of volume price discounts. This is not to say that all of these programs ignore price discounts completely. For example, the well-known program developed by McGilliard and Clay (1983a), MAX-BULL, allows the user to specify a constant percentage discount across all straws of semen. However, a constant percentage discount is not the same as a volume discount. Volume discounts are discounts that vary depending on the number of straws purchased. Thus, there may be no discount if 200 straws are purchased, but a 10% discount may apply if 300 straws are purchased.

The approach taken by MAXBULL is sufficient if the breeding company offers a constant discount on all straw purchases regardless of quantity. It can also be sufficient if the breeding company offers volume discounts, but the producer plans to purchase semen from only one company. In this case, the producer can estimate how many straws he will need for the year and given that purchase quantity, he can project his volume discount rate and manually input the appropriate discount into the model, thereby creating a constant discount model.

However, when the producer would like to consider the impact of volume discounts on his semen and breeding company choices, an optimization program which incorporates volume discounts (not just a constant discount) is required. If the producer is free to purchase semen from 2 or more breeding companies and at least one of these companies offers volume discounts, it may be optimal to purchase all of the needed semen at one breeding company in order for example, to receive the highest discount level. Or, it may be optimal to split the semen purchase between 3 different breeding companies to receive the first discount level from each company. Thus, when multiple breeding companies are available and these companies offer quantity discounts, it is not adequate to simply indicate a single discount rate for each breeding company. Volume discounts are dynamic since the discount rate changes as the quantity of straws purchased from that breeding company changes. Similarly, the optimal quantity changes as the discount rate changes. This interdependence means that to find an optimal, low-cost portfolio of bulls, the effect of volume discounts must be imbedded in the optimization program.

Recently, breeding companies have used volume discounts as marketing incentives to attract potential customers and to ensure the complete use of their bulls on a given herd. This allows breeding companies to compete for herds rather than for the sale of each straw of semen. Since volume discounts have become such an important part of the semen industry, these discount opportunities should be considered as part of the semen selection decision. Yet, current programs treat volume discounts as afterthoughts, and these types of discounts are only applied if by chance sufficient straws have been purchased from a given breeding company. Under these traditional programs, when a volume discount is received, it is the result of good fortune not of systematic optimization. Traditional methods fail to capture the potential economic benefits of optimally selecting and combining breeding company discount opportunities in the decision process itself.

The purpose of this project was to look at the optimal selection of bulls across multiple breeding companies given unique price/volume discounts. This was accomplished by using integer programming. In the past, integer programming has been used to ensure that the semen optimization program solved for only integer values of straws of semen (Erba et al., 1991; Nash and Rogers, 1995) rather than for fractions of straws. In the current application, integer programming was not used in this manner. Instead, it was used as an on/off switch. Binary integer programming was used to turn different discount rates ‘on’ when certain minimum quantities were purchased. This allowed the program to model decisions regarding which bulls should be used, from which breeding companies, and at which discount prices. Data from a local Holstein herd of 625 cows was used to demonstrate the application of this methodology and to estimate its benefits over traditional approaches for this herd.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The Programming Approach
A linear programming model (LP) was used to define the objective of minimizing total semen cost subject to a set of trait constraints defined by the producer. These constraints were identified in a manner similar to previously published selection models (Galligan and Ferguson, 1995) with the general format being:

Decision variables:


Objective Function:


Constraints:




where

Bulln,s=is the number of straws of bull n in breeding company s (for studs 1, 2, and 3) used in the breeding program (rounded to integer units),

Pricen,s=is the cost of a straw of semen from bull n in breeding company s,

N=the number of straws required by the farm,

Traitn,s=reported trait value for bull n in breeding company s,

Traitgoal=trait level goal for the particular dairy herd, and

t=the number of bulls in each breeding company.

The first constraint ensures that the total number of straws of semen required is purchased. The second constraint limits each bull to a maximum of 15% of the portfolio. This reduces the risk of inbreeding and propagation of unknown deleterious genes (Galligan and Ferguson, 1995). The final constraint ensures that specific trait objectives are met on average by the portfolio of bulls. All of these constraints can be changed to reflect the owners’ individual preferences.

Integer programming techniques can be used to augment the general linear program above and to capture the opportunity to identify which combination of breeding companies should be used at which price discount level. Furthermore, integer programming can help the producer to decide whether more straws than are needed should be purchased in order to capture an additional discount tier. For example, a given breeding company may offer a 30% discount off all straws purchased once a certain volume (Qs) is reached. Under this scenario, producers might find it advantageous to purchase additional semen, which would not be used in the breeding program but would allow for the capture of a volume price discount. A further complication arises when breeding companies offer different volume discounts for different quantities of semen purchased. To capture these types of questions, binary decision variables (integer 1,0 variables) must be integrated into the LP model above.

For example, to select the optimal breeding company to use when one company offers a flat 30% discount on all straws purchased once a certain volume (Qs) is purchased, the following decision variables and constraints may be added to the general linear programming model:

Additional decision variables are defined:



New Objective Function:


Additional constraints:


([1])


([2])

Original constraints:




where

Xs=a binary variable that is 1 if the discount is used (d = 1) and 0 if it is not used (d = 2).

Bulln,s,d=the number of straws of bull n from breeding company s purchased at price point d,

where

d=1, at the first price point (no discount is used),

d=2, at the second price point (the discount is used),

r=the total number of price points available at each breeding company,

Qs=the quantity of semen required before price point 2 (d = 2) is available,

Bullexcessn,s,d=number of straws of each bull selected to increase the "volume" to achieve a discount. Thus, these straws do not contribute to the breeding target values, but do contribute to the total cost of the breeding program,

M=the total number of straws needed.

Bullexcess is the number of excess straws that are purchased in order to receive the discount. As such, they contribute to the total price (and thus are added to the objective function), but do not contribute to the trait values (and thus are not added to the original goal trait constraints). Constraining the decision variable, Xs, to be binary forces this variable to be either one (if the discount is used, d = 2) or zero (if the discount is not used, d = 1). If more than Qs straws are purchased, the 30% discount is used and the binary variable has the following value: Xs = 1.

Equations 1 and 2 are called linking constraints (Ragsdale, 2001). These constraints link Bulln,s,d and Xs, ensuring that the required relationship between these variables is enforced. The first linking constraint (equation 1) sets a minimum order quantity (Qs) before sprice point 2 can be used. It guarantees that when the discount rate from breeding company s is used (i.e., when d = 2 and Xs = 1), the sum of bulls from that breeding company is greater then the semen quantity (Qs) required for the discount to apply. The second of these constraints (equation 2) uses the variable M (a number equal to the total number of straws required) to force the total number of discounted semen straws purchased at the discounted price ({Sigma} Bulln,s,d=2) to equal zero if Xs equals 0. Together these constraints ensure that the discount is applied to all straws if more than Qs straws are purchased and to zero straws if less than Qs straws are purchased.

The use of these linking constraints ensures that the decision program considers quantity discounts during the selection decision. Moreover, linking constraints maintain the linearity of the program.

The simple discounting structure above describes a situation in which a certain number of straws are purchased (Qs), resulting in all of the straws receiving a discount. Other purchasing incentives are currently offered in the industry. These involve some form of marginal discounting where, for example, the first 100 straws may be available at full price while the next 200 straws may be available at a 30% discount. In effect, if 250 straws are purchased, a herd can pay 2 different prices for the same straw of semen. The linking constraints for this problem are shown below and are similar to the previous constraints but are modified to reflect the marginal issues:


([3])


([4])

where

Bulln,s,d+1=the number of straws of bull n from breeding company s purchased at the next marginal price point.

The first linking constraint (equation 3) ensures that the undiscounted price point (d = 1) is fully utilized before the discounted price point (d = 2) is permitted. The last linking constraint (equation 4) links the first price point to the second price point. It forces Xs to equal one if semen is purchased at the second price point. Together, these 2 linking constraints do not allow any straws of semen at price point 2 to be purchased unless the required number of straws at price point 1 has been purchased. This captures the marginal dimension of the decision problem. In other words, for the discounted price d = 2 to apply (where Xs = 1), the sum of bulls and bullexcess at price point 1 must be greater than or equal to Qs. The objective function to minimize semen cost will ensure that the higher discounts are applied to the more expensive semen.

Application Data Set
To demonstrate the application of these principles, a local 625-cow Holstein dairy provided information about the opportunities presented to them from 3 large commercial breeding companies (studs A, B, and C). The trait and semen price data were obtained from the National Association of Animal Breeders May 2003 list of active sires and are summarized in Table 1Go. Eligible bulls were required to have a minimum reliability index of 95%, resulting in a total of 52 bulls available for selection across the 3 breeding companies. Four traits (average PTA milk, PTA fat, PTA protein, and PTA calving ease) were identified as important to the producer, and the goal values for these traits were stipulated (Table 1Go).


View this table:
[in this window]
[in a new window]
 
Table 1. Summary of available bulls.
 
The producer’s objective was to obtain a portfolio of semen that, on average, met or exceeded the trait goals while minimizing cost (LP model). This herd required 1500 straws of semen from proven bulls each year to meet its expected breeding needs.

To understand the value of missed opportunities (due to a failure to consider breeding company selection and volume discounts within the optimization program), 3 selection methods were evaluated. The first method (method 1) reflected a traditional approach, where the LP format presented above was used to select bulls across all 3 breeding companies at their published (non-discounted) prices. If by chance a sufficient quantity of semen was purchased to qualify for a volume discount within a breeding company, then the discount was applied to determine the total discounted semen cost for the herd. Note that the discount was applied after the selection of bulls.

The second method (method 2) used the LP program to select bulls from a single breeding company. This method reflects the fact that producers sometimes choose to use just one breeding company and they often choose the company before selecting individual bulls. The LP format was used to create a separate optimization program for each of the 3 breeding companies. In this case, due to the ability to predict the volume discount (because it was known that 1500 straws would be purchased from one of the 3 companies), the producer could have anticipated the discount level before the bulls were selected or he could have applied the discount after the bulls were selected; both approaches would have led to the same result. In this model, the applicable discounts were applied after the bulls were selected.

The third approach (INTEGER) incorporated the volume discount opportunities into the decision process. The 3 breeding companies offered this farm different volume discounts. Stud A offered a flat 50% discount on all straws if more than 250 straws were purchased. Thus, if 300 straws were purchased, every single straw received a 50% discount. Stud B offered a tiered, fixed dollar discount, where a $6 discount applied to every straw if more than 500 straws were purchased and a $7 discount applied to all straws (even the first 500) if more than 1000 straws were purchased. Stud B also stated that the minimum price a straw could cost was $5. Stud C offered a marginal discount where the first 750 straws were full priced, the next 250 straws received a 25% discount, and any additional straws received an additional 30% discount.

INTEGER solved for an objective function that minimized the total cost of discounted semen (in contrast to method 1 and 2 which minimized the total cost of nondiscounted semen). To adequately account for these various price discounts, the following integer program was developed:

Decision variables:




Objective Function:


Constraints:





Stud A constraints


([5])


([6])

Stud B constraints


([7])


([8])

Stud C constraints


([9])


([10])


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The combinations of bulls that met or exceeded the producer’s breeding goals, while minimizing the cost of semen were calculated for each of the strategies. Method 1 selected bulls from across all breeding companies, resulting in a low cost portfolio total of $11,040 once discounts were appropriately applied, retroactively (Table 2Go). When the optimization program was constrained to select bulls from within breeding companies (method 2), stud A had a final total semen cost of $12,127 and stud B had a final total semen cost of $11,038, once discounts were applied. Stud C did not have a combination of bulls that met the producer’s breeding goals. This failing was probably due to the fact that only 8 bulls were available for use from stud C due to the 95% reliability requirement. This requirement, in conjunction with the previously defined maximum of 15% of the total semen coming from any one bull, made it impossible for stud C to provide a feasible solution. Hence, the optimal breeding company selection from method 2 was to use stud B at $11,038.


View this table:
[in this window]
[in a new window]
 
Table 2. Selection goals and optimization results for three optimization methods.
 
The integer program strategy (INTEGER) incorporated breeding company choice into the optimization routine and resulted in a total discounted semen price of $9,677. This represents a 12.3% decrease in the cost of semen compared with the discounted method 1 portfolio and a 20% and 12.3% decrease in the cost of semen compared with the discounted portfolio strategies for just stud A and stud B, respectively. A sensitivity analysis showed that the INTEGER model produced cost savings regardless of whether the producer set higher or lower genetic trait goals (Table 3Go). The INTEGER program consistently resulted in the low cost solution since it optimized the use of quantity discounts.


View this table:
[in this window]
[in a new window]
 
Table 3. Sensitivity of price to changes in trait goals.
 
The INTEGER program also provided the option of purchasing extra bull semen (Bullexcessn,s,d), which would not be used in the breeding program, but would allow the producer to reach the quantity threshold for discounts. This semen was assumed to be thrown away and thus, it did not contribute to the breeding goals, but was incorporated into the cost per straw of the semen required by the producer. While Bullexcess was offered as an opportunity, it was not used in this particular application. A sensitivity analysis was performed to identify when Bullexcess might be incorporated into the decision. The INTEGER program was solved again (under the original trait goals) using different values for the total straws of semen required (Table 4Go).


View this table:
[in this window]
[in a new window]
 
Table 4. Sensitivity of Bullexcess to changes in number of straws required.
 
As was expected, the cost per straw of semen increased when fewer straws of semen were required. This was due to the decreased ability to take advantage of quantity discounts. Furthermore, when fewer straws of semen were required, the option to purchase extra straws that would not be used in the breeding program became more important. This opportunity arises because straws of semen are priced differently; thereby creating an opportunity to purchase inexpensive straws of semen to drive up the volume and achieve a quantity discount, which is applied to the entire semen purchase. Often the first discount is the most substantial, thus it may be economical to purchase excess semen to reach the first discount level.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Current optimization programs fail to address the decisions that producers face regarding volume price discounts. These discounts can be significant determinants in deciding which bulls to use and from which breeding companies. This application used binary variables to incorporate breeding company selection and quantity discounts into the semen selection program and, it resulted in a more economical portfolio of bulls with a lower total semen cost. In the applied example, a 12 to 20% savings was observed over traditional optimization methods.

While the problems to be optimized may vary and thus may require different configurations of linking constraints and decision variables, binary programming (which is now readily available on common software; Solver in Excel) can be a powerful tool for making the model more realistic. Incorporating binary variables, instead of ignoring these nuances, will result in better solutions that more closely emulate real life decisions.

Received for publication August 31, 2003. Accepted for publication April 19, 2004.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 


Dematawewa, M. B., P. J. Berger, and B. E. Melton. 1998. Optimization of sire selection based on maximization of guaranteed income and risk associated with sire merit. J. Dairy Sci. 81:807–816.[Abstract]

Eicker, S. W. 1995. Using linear programming to select a least-cost portfolio of services for dairy herds. American Dairy Science Association, 78. (Abstr.)

Erba, E. M., J. G. Fadel, and T. R. Famula. 1991. The use of integer programming in dairy sire selection. J. Dairy Sci. 74:3552–3560.[Abstract]

Galligan, D. T., and J. D. Ferguson. 1995. Application of linear programming in bull selection for a dairy herd. JAVMA 206:173–176.

McGilliard, M. L., and J. S. Clay. 1983a. Breeding programs of dairymen selecting Holstein sires by computer. J. Dairy Sci. 66:654–659.

McGilliard, M. L., and J. S. Clay. 1983b. Selecting groups of sires by computer to maximize herd breeding goals. J. Dairy Sci. 66:647–653.

Nash, D. L., and G. W. Rogers. 1995. Herd sire portfolio selection: A comparison of rounded linear and integer programming. J. Dairy Sci. 78:2486–2495.[Abstract]

Norman, H. D., and R. L. Powell. 1992. Genetic Change Attained and Possible. Pages 59–66 in Large Dairy Herd Management. American Dairy Science Association, Savoy, IL.

Plourd, R. 2000. The drive for efficiency. Dairy Herd Management. Online. Available: http://www.dairyherd.com.

Ragsdale, C. T. 2001. Spreadsheet Modeling and Decision Analysis, 3E. South-Western College Publishing, Cincinnati, OH.

Schneeberger, M., and A. E. Freeman. 1982. Application of portfolio theory to dairy sire selection. J Dairy Sci. 65:404–409.[Abstract/Free Full Text]

Stevenson, J. S. 2001. Reproductive management of dairy cows in high producing herds. J. Dairy Sci. 84(E. Suppl.):E124–E143.

Tozer, P. R., and J. R. Stokes. 2001. Using multiple objective programming in a dairy cow breeding program. J. Dairy Sci. 84:2782–2788.[Abstract]

Tozer, P. R., and J. R. Stokes. 2002. Producer breeding objectives and optimal sire selection. J. Dairy Sci. 85:3518–3525.[Abstract/Free Full Text]

Weigel, K. A., and S. W. Lin. 2000. Use of computerized mate selection programs to control inbreeding of Holstein and Jersey cattle in the next generation. J. Dairy Sci. 88:822–828.



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by McConnel, M. B.
Right arrow Articles by Galligan, D. T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by McConnel, M. B.
Right arrow Articles by Galligan, D. T.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS