THE ROBUST MIN-MAX NEWSVENDOR PROBLEM WITH BALKING UNDER A SERVICE LEVEL CONSTRAINT

In this paper we study an extension of a classic newsvendor model with balking under a service-level constraint. We also relax the assumption that the demand distribution is fully available to a decision-maker; we assume that only the mean and variance of the demand distribution are known, and discuss the procedure for determining the optimal order quantity for the concomitant model. We further extend our base model to two different cases: (i) when the fixed ordering cost is included, and (ii) when yield is uncertain. We illustrate the solution procedures for each case with numerical examples. Moreover, we discuss the performance and robustness of the approach through randomly generated test instances, and perform a numerical sensitivity analysis to evaluate the impact of the changes of a targeted fill-rate and variances of a demand distribution.

The newsvendor model has served as one of the fundamental models in stochastic inventory control problems [17], and has been used in numerous applications including production/inventory control, supply chain coordinations, option pricing and contract, and yield management.A comprehensive literature survey of the model can be found in Khouja [8] and Petruzzi & Dada [16].
It is well known that newsvendor-type problems can easily be solved when the full demand distribution is available, and that they have an elegant solution structure called the critical fractile solution [17].However, in practice it is often found that there may be weak or no justification that demand should follow a certain distribution.Thus a robust min-max (or distribution-free) approach for the newsvendor model, pioneered by Scarf [18], has been proposed, and various forms of extensions have been discussed in the extant literature [1,5,6,9,11,13].It allows us to make a decision about the worst case distribution when only the mean and variance information for demand can readily be captured or estimated from the data.This is indeed a risk-averse approach [15].
In this paper we specifically consider a robust min-max approach for an extension of the newsvendor model with customer balking and a service-level constraint.It is often observed that demand for a product can be influenced by product availability on the shelf.For example, customers may not purchase items if the inventory on the shelf is below a certain level.We often observe such customer behaviour in our daily life, especially for perishable products such as dairy products, fruits, or vegetables in grocery stores.If there are few items left on the shelf, customers are reluctant to buy them, possibly thinking that the products might not be in as good shape as they should to be, might no longer be fresh, or might be close to their expiry date.This is referred to as 'balking' [14].Thus, in this study, if the inventory level drops below a threshold , we assume that customers may balk with probability  − , and so demand is deflated by .
Moon & Choi [11] and Liao et al. [9] are the researchers who are most relevant to this study.Moon & Choi [11] provide the optimal order quantity that minimises the expected cost against the worst possible distribution of the demand with mean  and variance   with balking.However, they disregard the impact of lost sales, which is typically captured by the shortage cost in a newsvendor setting.Most of all, the phenomenon of customer balking exacerbates the magnitude of lost sales.It is thus reasonable and important to incorporate this aspect to reflect customer balking effects better.Liao et al. [9] examine a simple extension of the model in Moon & Choi [11] by also incorporating a linear lost-sales penalty cost in the profit function.However, the main concern about this extension is that, in many cases, the lost sales penalty is difficult to estimate in practice, since it includes intangible factors such as loss of goodwill or customers.In order to handle the lost sales effect in a tractable manner, we incorporate the notion of service level (or fill-rate) into the decision framework, and so our model is more viable in practice than the earlier models in the literature.
The main contribution of this work is that we incorporate the service level constraint for considering the lost-sales effect due to the understocking of items and customer balking, and we provide a practical guide to the interdependence of retail shelf management and associated inventory control in the broader sense.As Liao et al. [9] point out, the previous distribution-free newsvendor model with balking neglected this aspect, despite the significance of lost sales.They therefore attempt to consider this aspect with the explicit lost-sales penalty term in the profit function -although they admit that the penalty cost is hardly ever estimated in practice.In this paper we aim to discuss the optimal ordering policy for the robust min-max single-period model with balking under a fill-rate constraint, to improve the practicality of the concomitant model.Furthermore, we present two extensions of the base model: one including a fixed ordering cost, and the other assuming yield uncertainty.We note that, in the shipment of perishable goods, the amount that perishes during storage and transit in the United States is 10 to 15 per cent of perishable freight tonnage [10].The percentages are even higher in countries that do not have a sufficiently developed cold chain infrastructure [3].Thus the extension for yield uncertainty also readily captures the aspect in which the exact random shrinkage distribution during a shipment is almost never estimated.(Many items in a grocery store are perishable and are subject to freight shrinkage during a shipment.)We then evaluate the performance of the concomitant robust min-max model with a set of randomly generated instances, and discuss the conditions under which the robust min-max approach presented in this paper can be recommended for use by practitioners in the retail industry.
The remainder of this paper is organised as follows.Section 2 presents the base model with balking and a service level constraint, using cases where the exact demand distribution is either available or is not available.Sections 3 and 4 discuss the extensions of the base model in Section 2: the fixed ordering cost case and the random yield case respectively.Section 5 presents computational experiments for evaluating the performance and checking the robustness of distribution-free models.We conclude the discussion in Section 6.

BASE MODEL
We mainly adopt the notation presented in Moon and Choi [11] (listed in Table 1): Then the expected profit function π F (Q) can be expressed as which has explicit dependence on the demand distribution F(D) with mean µ and variance σ 2 .As presented in the equation above, there are three possible cases, depending on the demand realisation: a.
0 ≤ D ≤ Q − K : the resulting inventory level after the demand is met is still above the threshold K, and so no change occurs in the sales rate.Q − K + K/γ ≤ D: all the ordered items must be sold out; so a shortage of items occurs with the amount of D − (Q − K + K/γ).
As mentioned earlier, we introduce the service level constraint in our decision model.Service level or fill-rate β is defined as β = 1 − Expected lost sales Expected total demand ⁄ , which is the average fraction of demand that is satisfied immediately from stock on the shelf [12].Let ξ be the desired service level, and the decision-maker determines the optimal order quantity under the constraint β ≥ ξ.We note that, in the concomitant model, the expected lost sales equal ∫ , and so the fill-rate β(Q), which is a function of an order quantity Q, is given as Under the assumption that the demand distribution F(⋅) is fully available to the decisionmaker, the optimal order quantity Q F * is determined as follows: We note that π F (Q) can be rewritten as Note that maximising π F (Q) is equivalent to minimising   (), and it can easily be verified that   () is convex in  (by evaluating the secondorder derivative of   () with respect to ).
Similarly, the fill-rate constraint can be also expressed as Let Consequently, if the demand distribution information F(⋅) is available, the optimal order quantity Q F * is determined by solving the constrained convex optimisation problem min ∈Ω    ().We now discuss how to determine the optimal order quantity Q F * when the exact demand distribution information is known.Since   () is convex in , as mentioned earlier, we examine the first-order condition: † Let Otherwise the optimal order quantity must be a boundary point in Ω F .It can be shown that the fillrate , which is a function of , strictly increases in .Then the conditions that (i)   () is convex in  and (ii)   ′ ∉ Ω F imply that   * should satisfy the fill-rate constraint with equality (i.e.SL F (  * ) = 0).Algorithm 1 summarises the discussion above.
Algorithm 1: Procedure to determine the optimal order quantity   * when the exact demand distribution (⋅) is available So far we have discussed the newsvendor model with balking and a service level constraint under the full demand distribution information.We now present the corresponding robust † Let () where ℎ() is a function of .Then min-max model in which there is no assumption about the demand distribution apart from the mean  and variance  2 .Let Γ(, ) be the set of all possible demand distributions with mean  and variance  2 .In the robust min-max model, we derive an upper-bound function (Q) such that  ∈Γ(,)   () ≤ () for all  ≥ 0, and the following lemma provides the theoretical basis for the derivation of the upper bound function: Lemma 1.Let () be a random variable, which is a function of  and independent of , [()] be the expectation of (), and [()] be the variance of ().a.Then, (3) b.Suppose that () is a constant (i.e.[()] = () and [()] = 0).Then, for any (), there is always a demand random variable  whose distribution is   ∈ (, ) where the upper bound (right-hand side of Equation ( 3)) is tight.
Thus, (i) when , and (ii) when . Based on the results, The first and second derivatives of respectively by using the Leibniz rule.Since We next consider the fill-rate constraint when there is no assumption about the demand distribution.Due to the result from Lemma 1, and let () be the function of the right-hand-side of the inequality above.We identify the constraint () ≤ 0 as the distribution-free fill-rate constraint.We note that Lemma 1(ii) implies that there is a demand distribution  satisfying SL F () = ().Since it can be verified that () is also convex in , the set Ω = � ≥ 0: () ≤ 0� is convex, and moreover Ω ⊆ Ω  for all F ∈ Γ(µ, σ) due to the relation SL F () ≤ ().
Therefore, the robust min-max model is to determine Q that minimises the upper bound function C() over the set Ω, say, Q * all F ∈ Γ(µ, σ).We now discuss how Q * can be determined, and we can apply the same argument that we made for the case where the exact demand distribution is known because both () and () are convex in .Thus we present the simple solution procedure to determine Q * for this robust min-max model in Algorithm 2.
2. Find Q ′′ satisfying the distribution-free fill-rate constraint with equality -i.e. �Q ′′ � = 0 or Algorithm 2: Procedure to determine the optimal order quantity  *

for the robust minmax model
We note that Moon & Choi [11] suggested the notion of the expected value of additional information (EVAI) to measure the additional value of acquiring the exact demand distribution.In this case, the expected marginal value of knowing the exact demand distribution relative to the lack of distribution information can be expressed as ), and we introduce ) × 100 (%), which is the percentage loss due to the lack of demand distribution information, to evaluate the performance of the robust min-max approach.
Example 1.For the illustration, we assume that the mean and the standard deviation of the demand are  = 800 and  = 150 respectively.The demand deflation threshold and the corresponding deflation rate are  = 200 and  = 0.8 respectively.The cost parameters are as follows: the unit ordering cost is  = $35, the unit retail price is  = $60, and the unit salvage value is  = $15.Lastly, let the target fill-rate be  = 0.85.Suppose that the demand distribution is given as (540,1060), whose mean and standard deviation are 800 and 150 respectively.The optimal order quantity   * for the distribution is approximately 829; on the other hand, the optimal quantity Q * for the robust min-max model is close to 804.Since C F (  * ) ≈ $19319.764and   ( THE FIXED ORDERING COST CASE We now discuss an extension of the base model given in Section 2. Suppose (i) that whenever the retailer procures any positive number of items (i.e.Q > 0), a fixed cost A is charged; and (ii) that the retailer owns an inventory I ≥ 0 initially where possessing the initial inventory is irrelevant to any cost.It is well known that the classic newsvendor model with a fixed ordering cost follows the so-called (s, S) inventory replenishment policy [4].However, it will be shown that the extension we examine in this section does not follow the (s, S) policy because of the fill-rate constraint.
Let S be the inventory level after receiving an order (i.e. S = I + Q).Then the expected total profit function under the exact demand information π F (1) () is expressed as () + 1 > : () ≤ 0� where () = () + 1 > , -that is,  ′ satisfies the following first-order condition: ().In addition, there is a unique  ′ (<  ′ ) satisfying the equation The fact that () is strictly convex and is not bounded from above ensures the existence of such an  ′ .We also note that  ′ corresponds to the reorder point in a (, ) policy, and without the fill-rate constraint, the pair ( ′ ,  ′ ) forms such a policy, because the retailer places an order up to  ′ only when  <  ′ due to the fact that () > � ′ � +  for  <  ′ , and  ′ is the break-even point.
We now discuss how to determine the optimal order quantity  * for this robust min-max model.Let  ′′ be a  satisfying the distribution-free fill-rate constraint with equality, or We note that the fill-rate constraint � ′ � + .This suggests the following ordering policy: Algorithm 3 summarises the procedure for determining the optimal order quantity for the robust min-max model with a fixed ordering cost.
� ′ � + , and 2. Given the initial inventory level , the optimal order quantity  * is determined as follows: , the optimal policy in this case is (851 − ) if  < 851, and 0 otherwise.

THE RANDOM YIELD CASE
In this section we consider the case in which the ordered or produced quantity Q is not perfect.This type of problem is in the class of random yield models [19].We especially assume the multiplicative random yield model that appears frequently in the literature [2,7].We note that Jung & Lee [6] present a robust newsvendor model with balking and random yield; however, unlike ours, their model assumes that the yield distribution follows a binomial distribution.Let ρ be a random variable with support on (0,1] whose probability (cumulative) density function is y(ρ) �Y(ρ)�.In addition, we denote µ  and σ  2 as the expectation and variance of the random variable ρ respectively.We assume that ρ is independent of the order quantity Q and of the demand D. In this random yield case, it is assumed that out of Q items, only ρQ can be sellable or in a good shape.Then the expected profit function π F (2) () can be expressed as  ≥ .Under the assumption that the exact demand and yield distributions are available, the optimal order quantity   * is We now discuss the robust min-max approach where the exact distributions for both demand D and yield ρ are not available.Due to Lemma 1, by letting By applying the results above, the objective function for the corresponding robust min-max model, C � (2) (), is derived as Similarly, the distribution-free fill-rate constraint for the random yield case becomes () ≤ 0 where () = Therefore the corresponding robust min-max decision model with random yield reduces to Since it can be shown that  (2) () is convex in , we set the first derivative of () with respect to  equal to zero, which leads to And let   ′ be  satisfying the first-order condition above.If �  ′ � ≤ 0 holds, then   ′ is the optimal order quantity for the robust min-max model with random yield.Otherwise, in a similar way to the base model case   ′′ , satisfying the fill-rate constraint with equality (i.e. �  ′′ � = 0) is the optimal solution for the model.Thus we can summarise a simple procedure to seek the optimal order quantity   * for the robust min-max model with random yield as presented in Algorithm 4.

𝑆𝐿
(2)  �  ′ � ≤ 0, then, according to Algorithm 4, the optimal order quantity is   * ≈ 846.However, if  = 0.95, the fill-rate constraint is violated, and so, from Step 2 in Algorithm 4,   ′′ ≈ 986, which is the optimal order quantity in this case.We note that the order quantities under perfect yield assumptions in Example 1 are 804 if  = 0.85, and 851 if  = 0.95.Due to the the random yield effect of produced or ordered quantity, the order quantity increases, and its increasing rate seems higher as the targeted fill-rate increases.

COMPUTATIONAL RESULTS: PERFORMANCE OF THE ROBUST MIN-MAX APPROACH
We have performed extensive computational experiments to evaluate the performance and robustness of this distribution-free approach.For the experiments, 1,000 problem instances are randomly generated, and all the parameters for each instance are drawn from the uniform distributions specified in Table 2.We take some parameter values from Moon & Choi [11].For each instance, we evaluate three different distributions (i.e.normal, uniform, and triangle distributions), and compare them with results from the robust min-max approach.
We note that the information on the mode  of a triangle distribution is additionally required, and it is drawn from (750,900).In addition, for a given instance, we evaluate (%) for four different levels of the desired fill-rate  as {0.8, 0.85, 0.9, 0.95}.] for each distributional form, which are essential to calculate   (Q).We note that we replaced median with the sample average in the box plots.The results indicate (i) that the percentage gap (%) is less than 3 per cent on average for all the targeted service levels, and (ii) that when the targeted fill-rate  is less than or equal to 0.9 (90 per cent), the maximum value of (%) remains less than 3 per cent.However, we can observe that both the variability and the average of (%) increase significantly at  = 0.95, compared with the cases for  ≤ 0.9, and its maximum reaches to 10-12 per cent.
For  ≤ 0.95, it turns out that its application to the normal distribution results in better performance than in the other distributions.Furthermore, these computational results partly provide justification for the use of the robust min-max approach if (i) the desired fillrate is less than 90 per cent, and (ii) it is difficult to estimate the form of a demand distribution function, but the information of mean and variance is readily available or can be estimated.When the targeted fill-rate is 95 per cent, the performance of the robust min-max approach becomes degraded, and it seems that the degradation would be exacerbated for a fairly high fill-rate (close to 1).However, use of the robust min-max approach could still be recommended (at least for the cases we examined in this experiment), since the 75-percentile of (%) from the set of problem instances is still below 4 per cent.Therefore the robust min-max approach is likely to work well or to be acceptable in most cases.However, the increased mean and variance of (%EVAI) for the min-max model, compared with the base model, could mainly be because the optimal solutions for the robust min-max model are likely to lie on the constraint boundary of the distribution-free fill-rate constraint for a high ξ, and even more conservative decisions of the robust min-max model could be made when the service level is high We now investigate how the change in variance of the demand affects the performance of the robust min-max approach depending on the targeted fill-rate .Let  be the coefficient of variation (i.e.,  = /).We note that, in the previous numerical experiment, we restricted the value of  to between 0.1 and 0.5 (refer to Table 2).In order to see the effect of changes in the variance of demand for the uniform distribution, we choose 1,000 random instances in the same way as the previous experiments.For each value of  ∈ {0.8,0,85,0.9,0.95},we vary the value of  from 0.1 to 1.0 by a 0.1 increment, and calculate (%) as shown in Figure 2. As seen in Figure 2, the variance of demand has a significant impact on the performance of the robust min-max approach, and the impact increases as the desired service level increases.In general, the variability and average of (%) from a set of randomlygenerated problem instances increase as  increases.For  equal to either 0.8 or 0.85, the robust min-max approach could still be recommended on average, even for a higher variance in demand.The robust min-max approach performs moderately well in the cases when the demand variance is low ( ≤ 0.4) for  ≥ 0.9cases, and the performance degradation becomes significant when  is equal to 0.95.This suggests to practitioners that there may be a reason to explore the underlying true demand distribution if both the targeted fill-rate and the coefficient of variation are high enough.As mentioned earlier, decision-makers tend to make conservative decisions when they lack the exact demand distribution information for higher fill-rate cases.This often leads to a larger information gap and an increased variability of (%).In earlier works discussing robust min-max models [1,11], their numerical experiments used the value of  less than 0.3 to conclude that the distribution-free approach is robust.Although none of these studies attempted to investigate the impact of the variance of demand on the performance of the robust minmax approach, our results could indicate that robust the min-max models in the earlier literature might be robust even in the cases where demand distributions have high variance (see the cases when the desired service level  is not too high).As mentioned earlier, one reason that our model is more likely to be conservative as the targeted service level increases is due to the incorporation of the fill-rate constraint, which results in a performance degradation of robust min-max approaches.In harsh conditions, the fill-rate constraint could possibly raise the value of the order quantity  too high.However, their models have no such constraint, and we presume that the optimal order quantity obtained from the first-order condition may not be too conservative.
Lastly, we also performed numerical experiments to examine the effect of the threshold K, and found that its effect is insignificant.

CONCLUSION
In this paper we have investigated the robust min-max approach for the classic newsvendor model with customer balking under a service-level constraint.In the approach, we only require the first two moments of the demand distribution to determine an order quantity where the demand distributional form is hardly estimated, which in practice is usually the case.In order to control the lost-sales and customer balking effects in a tangible way, we incorporate the service rate constraint in our model, rather than include a lost-sales cost penalty, which is generally difficult to estimate in practice.We present the procedure to identify an optimal order quantity for the robust min-max model, and we extend the base model into two cases: one with a fixed ordering cost, and the other with imperfect yield.We also examine the effects of the variability of the demand and the targeted fill-rate on the information gap (%).Extensive numerical experiments indicate that the robust min-max approach generally performs well and efficiently; but we also observe that the quality of the solutions determined by the robust min-max approach is degraded for the very high targeted fill-rate, and the degradation accelerates as the fill-rate comes close to 1. Similarly, its performance becomes less robust as the coefficient of variation for the demand distribution increases for the high targeted fill-rate.Thus the distribution-free model is generally efficient and robust, especially for situations where the coefficient of variation  ≤ 0.4 or 0.5, and the targeted service level  ≤ 0.9.

Q
− K ≤ D ≤ Q − K + K/γ: the stock level becomes less than the threshold K, and so the amount of D � = (Q − K) + γ�D − (Q − K)� is the realised demand size, and the amount of Q − D � = K − γ(D − Q + K) will be salvaged at the end of the selling season.c.

where 1 𝛼
is an indicator function.The corresponding robust min-max optimisation model is min ≥ �

,
≤ 0 is equivalent to() = 1 − � 2 + ( −  + / − ) 2 − ( −  + / − ) 2 ≥ and it can be easily verified that the distribution-free fill-rate function () is increasing in .Thus we can interpret  ′′ as the minimum inventory level required for the targeted service level in the expected sense.Given  ′ , -rate constraint due to the property of (), and so the model becomes an unconstrained optimisation problem.The optimal order quantity is determined according to the ( ′ , the inventory level after receiving an order should be equal to  ′′ due to the convexity of the objective function  inventory level after receiving an order should be at least greater than or equal to  ′′ .If I ∈ ( ′ ,  ′′), the order needs to be placed up to  ′ 2) () = (1 − )( − ) , [( −  + ) + ] + ( − ) , �� − � −  +  �, and it is clear that maximising π F (2) () is equivalent to minimising C F (2) ().Like the previous cases, the fill-rate () is given as () = 1 −  , [�−(−+/)� + ]

Algorithm 4 :
Procedure to determine the optimal order quantity   * for the robust minmax model with random yield Example 3. In order to illustrate the random yield extension, we use the parameter values in Example 1.We additionally assume that   = 0.9 and   = 0.1.If  = 0.85,   ′ ≈ 846, and it satisfies

Figure 1 :
Figure 1: Effects of changing the targeted fill-rate Figure1presents the box plots of the ratio (%) from the 1,000 randomly-generated instances for each distributional form while varying the value of  in order to see the effect of the desired service level on the performance of the robust min-max approach.Appendix B provides the closed-form formulas for the cumulative density function ( ≤ ) and [( − ) + ] for each distributional form, which are essential to calculate   (Q).We note that we replaced median with the sample average in the box plots.The results indicate (i) that the percentage gap (%) is less than 3 per cent on average for all the targeted service levels, and (ii) that when the targeted fill-rate  is less than or equal to 0.9 (90 per cent), the maximum value of (%) remains less than 3 per cent.However, we can observe that both the variability and the average of (%) increase significantly at  = 0.95, compared with the cases for  ≤ 0.9, and its maximum reaches to 10-12 per cent.For  ≤ 0.95, it turns out that its application to the normal distribution results in better performance than in the other distributions.Furthermore, these computational results partly provide justification for the use of the robust min-max approach if (i) the desired fillrate is less than 90 per cent, and (ii) it is difficult to estimate the form of a demand distribution function, but the information of mean and variance is readily available or can be estimated.When the targeted fill-rate is 95 per cent, the performance of the robust min-max approach becomes degraded, and it seems that the degradation would be exacerbated for a fairly high fill-rate (close to 1).However, use of the robust min-max approach could still be recommended (at least for the cases we examined in this experiment), since the 75-percentile of (%) from the set of problem instances is still below 4 per cent.Therefore the robust min-max approach is likely to work well or to be

Figure 2 :
Figure 2: Effects of changing variance   of the demand when the demand follows uniform distribution