Two-stage stochastic programming approach for limited medical reserves allocation under uncertainties

At the early stage of public health emergencies, when the conventional medical reserves prepared are insufficient, and productivity could temporarily not meet the surge in demand, donations can be used to cover excess demand for medical supplies to a large extent. This paper explicitly considers the allocation problem of limited medical reserves during a public health emergency, incorporating uncertainty in demand and donated supplies and the priorities of health care centers. The problem is formulated as a two-stage stochastic program that regards the donated supplies as an efficient recourse action, aiming to minimize the total losses. The optimal allocation strategy of limited medical reserves and donations is obtained by solving the model using Gurobi solver. Finally, the effectiveness of the proposed approach is verified by a series of computational results, which show that the solutions of our method not only benefit the emergency demand fulfill rate but reduce the total losses as well.


Introduction
The storage pool consists of a fixed quantity of medical reserves prepared for emergencies, including public health emergencies and natural disasters such as an epidemic [1,2], earthquake and flood [3][4][5]. Medical reserves are set for emergencies that occur occasionally and usually result in uncertain demand. A very high minimum stock level to ensure sufficient supplies will result in excessive space because of low mobility, as well as substantial waste due to expiration. Therefore, medical reserves in the storage pool for emergency supplies maintain a relatively low (limited) stock level [6]. This work explicitly focuses on how to allocate limited medical reserves when a public health emergency occurs. In the early stage of a large-scale public health emergency, it takes time to organize manufacturers to restore production to step up supplies, especially during unusual times. Consequently, productivity may temporarily not meet the surge in demand arising from emergencies. As medical reserves prepared for an emergency response might be insufficient, donated supplies can be served as a recourse action to cover some excess demand to a large extent.
Certainly, how to utilize scarce medical materials prepared for emergency response to maximize their value is a vital issue. In the aftermath of a large-scale public health emergency, to reduce mortality and raise admission and cure rates, health care centers are grouped into two classes, designated hospitals and general hospitals. Designated hospitals mainly receive and cure patients in severe conditions, while general hospitals converted from large buildings are for mild cases, suspected cases, and close contacts. That severe patients obtain medical treatment plays a critical role in preventing and controlling a large-scale public health emergency, indicating that the same reserves distributed to the designated hospital will yield more utility than general hospitals. Therefore, when the pre-positioned medical reserves are inadequate, the designated hospital for severe patients should be prioritized to maximize the value of limited reserves. Inspired by Shin and Lee [7], we introduce priority to describe the weight of importance. In addition, health 1 3 care centers' requirements for resources and donated supplies cause considerable variability in the decision-making process. We incorporate and analyze many important reallife issues: limited medical reserves, the priorities of health care centers, uncertainty in demand and donated supplies. This decision process is complex and involves the allocation of multiple supplies, which has a significant impact on both the quality and cost of service.
The rest of the paper is organized as follows: "Literature review" reviews the relevant literature in the area of the medical resources allocation. "Problem formulation" delineates the problem, notations, and the quantitative model addressed in this study. A modeling framework for the allocation of limited medical reserves is introduced, and it can be formulated as a two-stage stochastic program that involves uncertainty in demand and donations. "Computational experiments" reports computational experiment results performed by Gurobi for solving the model proposed in "Problem formulation". The paper ends with overall conclusions and some future research directions.

Literature review
The medical resources allocation problem has been studied in several recent studies, and it has considerable potential in the management of medical resources from emergency scenarios.
Our work is closely related to the resource allocation problem, which has attracted much attention because of its practical relevance. There are increasing and various studies concentrating on the medical resource allocation problem. Cao et al. [8] proposed a discrete event simulation model to determine the most reasonable scarce medical resources allocation principle from four principles and evaluated their effects by considering the trade-off between survival efficiency and moral issues in the emergency response. Su et al. [9] defined disaster response alliance as a form of teamwork and proposed a multi-constrained integer linear programming model, which was solved by an algorithm based on differential evolution and heuristic exploration to allocate multiple emergency resources to various events in a parallel way. To handle the problem of multi-resource and multi-location based on capability, Altay [10] established an integer programming model to allocate resources better. Yu et al. [11] paid more attention to resource allocation's efficiency, effectiveness and fairness. Two equivalent models were applied to address this problem which was discussed well through a rollout algorithm. Yu et al. [12] applied an improved dynamic programming method to measure critical resource allocation effectiveness and fairness.
There are several studies considering limited medical resources in the resource allocation problem. Yu et al. [13] developed a Timed Colored Petri net to describe and model multi-resource-multi-activity. The hospital's limited medical resources can be used more reasonably through the medical information integration system, and the activities become more efficient. Ho et al. [14] considered dynamic resource allocation and optimized sequential intervention policies under resource constraints with a multi-fidelity optimization model and method. In public healthcare, Lai et al. [15] modeled resource allocation by a team-DEA analysis which could be used to optimally allocate resources.
Several studies involve uncertain conditions. A possibilistic programming model was presented by Kargar et al. [16] for a liver transportation and allocation problem, which takes uncertainty into account and trade-off between quality indicators (i.e., efficiency and fairness). Their models have maximized the patients' survival probability and minimized transportation cost and time. Bastian et al. [17] presented a stochastic multi-objective automatic optimization model for allocating medical resources in fixed-input health systems. The model automatically determined the hospital level to reallocate the system's input resources for the optimization of the overall system performance. Huang et al. [18] studied the dynamic allocation and scheduling of stochastic medical resources in surgery. They analyzed the bottleneck resources' problem for optimizing medical resource scheduling through the drum-buffer-rope scheduling manner. To address the challenge of medical staff allocation and the scheduling in random environments, a two-stage algorithm based on goal planning has been proposed by Chen et al. [19] to identify the minimum number of medical staff needed and arrange the optimal schedule for them. Zonouzi and Kargari [20] constructed a new mathematical model based on the data mining methods to maximize ambulances' coverage to demand points under uncertain conditions of rescue stations. The data mining methods were also used to estimate allocation problems under uncertainty, Franco and Alfonso [21] used a random counterpart optimization approach and proposed two models considering different uncertainty sources. Zhang [22] studied a general resource planning and allocation problem under uncertain environment by two-stage stochastic programming.
Hu et al. [23] provided a bi-objective robust resource allocation model to maximize emergency resources' efficiency and fairness, considering scarce supplies and many types of uncertainty. Boostani et al. [24] investigated and simulated the uncertainty in disaster relief conditions by a multiobjective mixed-integer stochastic programming model. Yang et al. [25] formed a distributed robust optimization model for emergency supplies pre-positioning under high uncertainty and limited information. Gao [26] designed a multi-commodity rebalancing problem through a bi-level stochastic mixed-integer programming model. Zhang et al. [27] established a two-stage stochastic programming model and used the Lagrangian relaxation method to deal with mobile blood collecting facilities' location and allocation. Xiang et al. [28] formulated a two-objective robust model to study the berth allocation problem when considering uncertain factors.
Some researches specifically focus on the allocation under the demand uncertainty. To determine optimal policies for medical drug distribution, Nikzad et al. [29] used a twostage stochastic programming method to incorporate demand uncertainty into the optimization model and evaluated the effectiveness of the model by two different approaches based on chance constraints. Batista et al. [30] studied the tradeoff between the allocation of medical resources and the cost of service by a bi-objective stochastic programming model while considering the uncertainty in demand and capacity.
Furthermore, the problem of limited resource allocation in uncertain conditions has also been carefully explored by scholars. Feng et al. [31] constructed a stochastic bi-objective mathematical model in emergency departments, while their studies developed a multi-objective simulation method to deal with the difficulties in the multi-objective limited medical resources allocation. Tang and Wang [32] used an adjustable robust optimization model to explore surgical capacity allocation with uncertain demand and reduce the losses of income due to resource shortages.
In general, the works of literature above have conducted various researches on the allocation of medical resources. Nevertheless, most of the existing works assume that either medical resources for emergencies are sufficient so that demand arising from emergencies can be satisfied or the allocation is implemented, while several critical conditions are known, such as demand or supply. We observe that both uncertainties in demand and limitation of supplies are hardly avoidable, since emergencies are sudden and unpredictable, demand is not accurately known, and resources are not always sufficient.
Donated supplies can be used to cover some excess demand to a large extent when medical reserves prepared for emergency response are insufficient. Therefore, it is necessary to take donated supplies into account. In this research, we discuss donated supplies, limited medical reserves, and uncertainty in demand. This is an exploratory study on incorporating donated supplies in a limited resources allocation problem.

Problem description
We introduce priority which refers to a weight of importance. As medical reserves prepared for emergency response are insufficient, health care centers' demand with high priority, such as designated hospitals for severe patients, should be prioritized.
When both the demand of health care centers and donated supplies are uncertain, the allocation of limited medical reserves for each health care center should be performed. Since productivity may temporarily not meet the surge in demand caused by the emergencies at early stages, donations can be regarded as an efficient recourse action to cover some demand which medical reserves fail to satisfy, while the demand and donations become known.
There are several assumptions itemized: (1) When public health emergencies occur, the control and management of medical materials are performed by the government. Since the cost related to transportation distance is not included, the total of medical reserves pre-positioned in multiple locations (even the appropriate medical storage of each hospital) can be regarded as the inventory level stored in one storage pool. The government undertakes the task of distribution. (2) The probability distribution of the random variables is known and derived from historical data. (3) The lower and higher demand level based on a specific scenario is employed to describe each health care center's demand fluctuation. The number of supplies received by health care centers that is more than the higher demand level or less than the lower demand level will generate related costs, i.e., penalty cost and shortage cost. (4) Donated supplies are all distributed by charitable organizations without storage in the storage pool. (5) Inventory level of a single type of product in the storage pool for emergency supplies, that is, the amount of medical reserves is known.
The decisions to be made are: (1) how to allocate limited medical reserves in the storage pool; (2) how to allocate uncertain donated supplies.
The objective is to minimize the total cost, which has several terms: (1) the allocation costs involving distribution cost and labor cost; (2) recourse cost, i.e., to qualify "losses" related to the allocation of limited medical reserves, and other costs are incurred after the first-stage decision is made due to secondary allocation, shortage, and surplus.

A two-stage stochastic programming formulation
In this section, we propose a two-stage stochastic optimization model for the allocation problem of limited medical reserves. Uncertainty in demand and donations is taken into consideration. The problem has a two-stage characteristic, i.e., before the realization of uncertainty, there is a hereand-now decision to make. After the realization of uncertainty, recourse actions are possible. Before presenting the model, we introduce some notations used for the sets, indices, parameters, and decision variables of the formulation. (2) A two-stage stochastic program In our modeling framework, decisions are made in two stages. The first-stage variables are here-and-now decisions represented by the x i . These variables decide on allocating reserves in the storage pool for emergency medical supplies, which are taken before demand and donations are known. The second-stage decisions are taken after the realization of uncertainty, which can be regarded as recourse actions. In this problem, we have three possible such actions: allocation of donations, surplus and shortage of medical supplies. Note that the second-stage variables depend on the first-stage decision variables as well as the realization of the random variables denoted by the vector ( ) , whose realizations correspond to the various scenarios ∈ Ω . For simplicity, we use to replace ( ) . Accordingly, we denote the second-stage variables or recourse decisions as y i (x, ),b i (x, ) , and i (x, ) , respectively. When the random event occurs, the second-stage parameters d i ( ), h i ( ), e( ) , become known. Then, the second-stage variables y i (x, ),b i (x, ) , and i (x, ) must be taken. The recourse decisions thus involve the allocation of medical reserves.
The first-stage problem: The second-stage problem: The first-stage problem defines the allocation strategies of medical reserves for health care centers before demand and donations are known. The objective function (1) is to minimize the total cost, which consists of cost for allocation of medical reserves and recourse cost, where E Q(x, ) denotes the expected value of the second-stage problem; constraints (2) imply the total allocation of medical reserves must be less than or equal to the inventory level; constraints (3) represent that the first-stage variables are non-negative.
Note that in the second-stage problem, decision variables x i are fixed. After demand and donations are known, the allocation strategy of donated supplies for health care centers will be further made, and it is useful for supporting decisions on the first stage. The second-stage objective is to minimize the recourse cost related to the first-stage decision, which involves other costs incurred due to the allocation of donated supplies, shortage, and surplus.
Constraints (5) and (6) describe the relationship between the demand level and the total allocation, consisting of both medical reserves and donated supplies. Constraints (7) ensure the amount of donations allocation is no more than the maximum available donated supplies when scenario realized. Constraints (8), (9), and (10) state that the secondstage variables are non-negative.

Extensive form of linear equivalent
When the realizations of the scenario ∈ Ω are finite, i.e.,Ω contains a finite number of elements, the second-stage decision variables can be indexed and the two-stage stochastic program mentioned above can be formulated as an equivalent linear programming model. Uncertainties can be represented as a finite set of scenarios, which is denoted as S. Each scenario occurs with a known probability p s , where ∑ s∈S p s = 1 . We redefine the second-stage decision variables and parameters concerning demand level and donations as follows: d is : lower demand level of health care center i when scenario s occurs; h is : higher demand level of health care center i when scenario s occurs; e s : maximum available donated supplies when scenario s occurs; y is : amount of donations allocated to health care center i when scenario s occurs; b is : amount of surplus of health care center i when scenario s occurs; is : amount of shortage of health care center i when scenario s occurs; The extensive form of the equivalent linear programming model is formulated as follows: Note that the formulation above can be solved by commercial solvers directly, such as Gurobi.

Proposition 1 The constraint matrix of the equivalent linear programming model is totally unimodular (details of the proof are provided in Appendix).
According to Proposition 1, we can directly solve the linear equivalent to obtain the optimal solution of the corresponding integer programming model, i.e., the linear program has an optimal integral solution. Therefore, the proposed approach can be applied to large-scale integer programming of the limited medical reserves allocation problems.

Expected value problem and wait-and-see problem
The expected value problem replaces all random variables with their expected values to mitigate uncertainty. In our case, replacing the quantity of demand and donations with their expected values leads to a deterministic problem. The expected result of using the solution of the expected value problem is defined as EEV. The gap between EEV and the optimal value of a stochastic program (SP) is defined as the value of the stochastic solution or VSS [33].
The expected value of perfect information or EVPI measures the maximal amount the decision-maker would be ready ∑ i∈I x i ≤ r to pay to get the information in advance. To obtain EVPI, the optimization problem associated with each scenario , i.e., distribution problem, should be solved first. When the demand level of each health care center and donations are known before allocating the medical reserves, the decision process corresponds to the wait-and-see approach. Find the optimal objective values of the problem for all scenarios and compute the expected value, which is commonly denoted by WS. EVPI is the gap between the optimal value of a stochastic program (SP) and WS [33].
VSS and EVPI are two measures of uncertainty effects. We introduce alternative measures, VSS gap and EVPI gap , which are independent from the magnitude of the values involved. VSS gap and EVPI gap are calculated as follows:

Computational experiments
A series of computational experiments are conducted to evaluate the effectiveness of the proposed stochastic model. In "Data generation", we itemize the specifics of the instances used in the experiments. In "Computational results on different size instances", we concentrate on evaluating the stochastic model by comparing it with the deterministic model using the expected value of random variables. "Parameter sensitivity analysis" discusses the influence of both the inventory level of medical reserves and designated hospitals' priorities on the different costs. Finally, we perform a cost analysis to identify the factors that cause a more significant total loss in the deterministic problem in the last "Cost analysis".

Data generation
We generate 15 random instances, including different problem sizes with the number of health care centers ranging from 20 to 200, and the number of scenarios from 20 to 100. We define the set of designated hospitals as D, and the corresponding number is |D| . We assume that |D| accounts for 20% of |I|.
The priorities of designated hospitals are the same and are fixed to be v i = 2, i ∈ D and the priorities of other health care centers are set to be v i = 1, i ∈ I�D . The allocation cost per unit of supplies (including medical reserves and donated supplies) for each health care center is generated uniformly from (1,10). We simplify the experiments by assuming the unit penalty cost of each health care center is roughly identical and is set as 50, whereas the unit shortage cost is ten (20) VSS gap = (EEV − SP)∕SP (21) EVPI gap = (SP − WS)∕WS times greater than the penalty cost, since shortage may lead to severe consequences during a public health emergency. The lower demand level of each health care center is randomly (uniformly) generated between 0 and 100, while the higher is, between 100 and 200. Both of them are integers. We assume that the inventory level of medical reserves and the amount of donated supplies are related to the problem size, which corresponds to the severity and effects of a public health emergency.
The inventory level of medical reserves r is calculated as follows: In addition, the amount of donations under each scenario is drawn uniformly from where C UB is the upper bound of the integer uniform distribution followed by the higher demand level, and C LB is the lower bound of the integer uniform distribution followed by the lower demand level. They are 200 and 0, respectively.
Each scenario has the same probability of occurring:p s = 1∕|S|.

Computational results on different size instances
In this section, we present and discuss the results of the proposed model, the deterministic model using the expected value of random variables. We have performed the test instances and solved the equivalent linear programming model by the commercial solver Gurobi 9.0.2. All experiments have been conducted on a personal computer with a processor i5 at 2.3 GHz and 8 GB memory.
According to the construction of the test instance mentioned in "Data generation", we generate 10 instances separately for each data set and report on the expectation and variance of EVPI gap and VSS gap , and the average solution time of the stochastic program and the deterministic counterpart.
Detailed results for each data set (on average) are listed in Table 1. Columns |I| and |D| provide the number of health care centers and designated hospitals. Column |S| represents the number of scenarios in different sizes. Then, we report the averaged results corresponding to 10 instances, for example, the mean and variance of VSS gap and EVPI gap , respectively, in the next four columns. The "CPU Time" columns show the average computational time (in milliseconds) required to obtain the optimal solution of 10 independent invocations of two models.
The column " VSS gap Avg." and the column " EVPI gap Avg." contain the average relative VSS and EVPI in percentage and further details are referred to "Expected value problem and wait-and-see problem". From these columns, we notice that the expectation of VSS gap ranges from 3.392 to 6.192%, while the expectation of EVPI gap ranges from 1.095 to 5.074%. The variances of VSS gap and EVPI gap reported in columns 6 and 7 show that the differences among all cases are less significant and below 2.503(%). Analyzing these columns, one can observe that the relative values of VSS and EVPI increase as the number of health care centers increases. This is an indication that an increase in problem size leads to further improvements in the value of both the stochastic solution and the knowledge about the future. Moreover, concerning the limited medical reserves allocation problem, the stochastic program clearly outperforms the deterministic program that uses the expectation of random variables.
As can be observed from Table 1, the average running time, for solving the stochastic programming model of the largestsized instance with 200 health care centers and 100 scenarios, is less than 1.05 s. It indicates that the proposed method enables decision-makers to tackle large-scale allocation problems of limited medical reserves and obtain the optimal integer solution quickly due to its natural unimodularity.

Parameter sensitivity analysis
In this section, we analyze the impact of two parameters on the total cost (losses): the inventory level of medical reserves and the priority of designated hospitals.

Influence of inventory level on total cost and the first-stage cost
As mentioned in "Data generation", we generate a mediumsized instance with 50 health care centers and 10 scenarios.
The priorities of designated hospitals and other health care centers are set as v i = 2, i ∈ D and v i = 1, i ∈ I�D , respectively. Let the inventory level of medical reserves vary from 0 to 6000. The objective values of the stochastic program and the expected value problem are obtained by solving the corresponding models. Figure 1 describes the influence of different inventory levels on the stochastic program's objective values and the deterministic program, and we observe that the stochastic program improves more significantly in cost reduction in all cases.
It is evident that the cost for allocating medical reserves in the first stage of both models equals 0, as the inventory level of medical reserves is 0. It reflects that both models' objective values are similar and equal to recourse cost, which is the sum of costs incurred due to the allocation of donated supplies, shortage, and surplus. As shown in Fig. 1, the two models' optimal objective values (total cost) maintain a decreasing trend as the inventory level increased, since the shortage cost for unsatisfied demand of all health care centers decreases. In addition, the total cost corresponding to the deterministic program stabilizes when the amount of reserves reaches 4000. By contrast, the total cost of the stochastic program remains steady when it reaches 3000. We observe that the medical reserves are sufficient, and the inventory level can cover the lower demand levels of all health care centers in this case. The allocation of medical reserves obtained by solving the stochastic model is more reasonable, compared to the allocation from the deterministic model, since excess reserves received by some health care centers result in a higher total cost of the deterministic model, which includes more penalty cost. To summarize, the results suggest the benefit from the stochastic program in the allocation of limited medical reserves whether the amount of reserves is sufficient or not. Figure 2 displays the influence of different inventory levels on the first-stage allocation cost, indicating that the increase in inventory levels leads to an increase in the number of medical reserves allocated to hospitals in the first stage. In addition, the result obtained by the stochastic program pays more attention on the first-stage allocation, i.e., more medical reserves are employed to cover demand. Compared to the deterministic program, the stochastic program uses limited reserves materials as much as possible to meet demand instead of donated supplies.

Influence of priority on total cost and shortage cost of designated hospitals
We use the same medium-sized instance as described in "Influence of inventory level on total cost and the first-stage cost" and set the inventory level of medical reserves equal to 1000. Let the priority of designated hospitals vary from 1 to 10 and obtain the stochastic program's optimal objective values and the deterministic program by solving the corresponding models. Figure 3 clearly shows the influence of different priorities of designated hospitals on the total cost. It reveals that using a stochastic program approach for this problem instead of a deterministic approach allows for a considerable improvement both in cost (losses) reduction and in demand fulfill rate of designated hospitals, although with the increase of the priority of designated hospitals, the optimal values of two models also increase. The growth of the stochastic program's optimal value is at a noticeable slower rate than that of the deterministic program. When the priority of designated hospitals is set to be 1, which means there is no difference between designated hospitals and other health care centers, the unit shortage cost of all health care centers is identical and the total cost in this case is the lowest. However, the optimal value of the stochastic program is still smaller than that of the deterministic program.
Since designated hospitals for severe patients play a key role in the prevention and control of a large-scale public health emergency, their demand should be prioritized. For this reason, we investigate the impact of priority on the shortage cost of designated hospitals, which associates with demand fulfillment. Figure 4 plots different situations by changing the parameter. We observe that two curves show the same pattern: the shortage cost decreases sharply at first and then the rate of increase is slow. For designated Fig. 2 Influence of different inventory levels on the first-stage allocation cost: stochastic program vs. deterministic program Fig. 3 Influence of different priorities on the total cost: stochastic program vs. deterministic program hospitals, there is no shortage in the stochastic program, when the corresponding priority of equals to 6, whereas the shortage cost is still high in the deterministic counterpart. We exploit the parameter to maximize the value of limited reserves by prioritizing the designated hospitals.

Cost analysis
In this section, we perform a cost analysis to identify the factors that caused more significant losses in the deterministic problem and to further evaluate the effectiveness of the proposed stochastic program. We start setting the inventory level of medical reserves as 3000 and the priority of the designated hospital as 2. Other parameters are generated as described earlier. To obtain the proportion of various costs from stochastic program, deterministic program and waitand-see program, Gurobi is also used to solve these models to optimality. In this medium-sized case, two measures of uncertainty: VSS gap and EVPI gap are equal to 11.237 and 0.816, respectively. Figure 5 illustrates the proportion of various costs from three models. There is no penalty cost when the stochastic program and wait-and-see program are successfully solved to optimality.
Since it is too difficult to predict the relevant information of demand and donations accurately, random variables are often replaced with their expectation to tackle uncertainty by decision-makers, which corresponds to the deterministic program Fig. 5c. As the figure indicates, it is clear that 94.612% of the total cost comes from shortage cost, and penalty cost takes up approximately 0.041%. It suggests that although shortage costs are incurred after most of the scenarios occur, there still exist surplus supplies received by several health care centers under a few scenarios. Consequently, the optimal solution of the deterministic program using the expected value may cause severe losses in a large-scale public health emergency. By contrast, the proposed stochastic program improves more significantly in shortage cost reduction (7.466% vs. 94.612%), and there is no shortage cost when unknown information (both demand and donations) is revealed (Fig. 5a).
This suggests that we can employ the stochastic programming method to make more reasonable and practical decisions on the allocation problem of limited reserves under uncertainty. In addition, when both reserves and donations cannot cover excess demand, more donations should be appealed for help from the community. In addition, more supportive policies should be implemented to guarantee adequate emergency stockpiles.

Conclusions
This study focuses on the allocation of limited medical reserves in the storage pool, when a public health emergency occurs. What is more, accurate information such as the demand of health care centers and the donated supplies is unknown. To obtain the optimal allocation strategy of this research, a two-stage stochastic modeling framework is proposed and its equivalent linear programming model is written. Since the constraint matrix arising in the linear equivalent is totally unimodular, this research can be applied to large-scale integer programs of the allocation of limited medical reserves as well.
The experimental results demonstrate the proposed approach not only benefits the emergency demand fulfill rate but reduces the total losses as well. Compared to the deterministic program that uses the solutions of the expected value problem, the proposed approach results in a significant reduction in terms of losses, by explicitly considering the uncertainty in demand level and donated supplies, as well as the priorities of health care centers.
Future work is needed to handle some limitations in the proposed model. For example, incorporating the inventory level of medical reserves in the allocation problem is worth investigating to obtain the optimal strategy for both reserving and allocating. Various recourse actions to mitigate the uncertainty in demand, such as transshipments mechanism for improving demand fulfill rate, need to be further studied.