1. Introduction
Kidney exchange programs (KEPs) offer a new opportunity for patients on a kidney transplantation list to find suitable kidneys. Because everyone is born with two kidneys and can survive with only one, kidney transplantations are possible. There is a higher chance of a match in family members, so kidney transplantations often occur among family members. However, problems arise when family members have kidneys that are incompatible with the patient. In this case, a kidney exchange program is a possible solution. Kidney exchanges can occur when there are two incompatible patient-donor pairs, but the donors are compatible with the patient in the other pairing. Kidney exchanges were introduced in 1986 by Felix Rappaport [
1] and first performed in South Korea in 1991 [
2,
3,
4]. Many countries, including Switzerland [
5], the UK [
6], and the USA [
7,
8,
9], have since performed kidney exchange transplantations. We here remark that the Korea Centers for Disease Control & Prevention (KCDC) has recently reviewed and discussed a plan for the development of an integrated IT healthcare platform for the KEP, and hospitals can participate in the provision of the essential information for their renal disease sufferers with their medical records to the platform. Thus, the KCDC believes that the nationwide KEP platform would possibly provide good matching plans to maximize social utility.
Operations research (OR) technique has been applied extensively across healthcare issues [
10], and KEP is one of the important issues where OR technique has been used. KEP is typically modeled as an integer linear programing network optimization problem. Depending on the nature of KEP, a lot of deterministic models for the optimization of these programs have been designed. However, considering that the risk of uncertainty in KEP has not been discussed much, it is still an ongoing issue. The risk of uncertainty in KEP can lead to failure in matching, which means that a surgery cannot be done. In this study, we use the Stochastic Programming (SP) [
11,
12] framework to consider failure scenarios before matching takes place. In particular, a two-stage stochastic programming is a promising method for this purpose, in that it incorporates the effect of data uncertainty in the second stage. We list possible exceptional events and classify them as two types of failure scenarios, i.e., node failure and arc failure. The former occurs within a single pair and the latter occurs between pairs. Node failure can occur when the condition of the patient in one of the pairs seriously deteriorates. Arc failure can occur when one of two pairs changes its mind, which means when a pair feels unsatisfied with its matching. Because a kidney exchange is a form of bartering, the patient-donor pairs play the role of customers. Depending on the nature of the barter, a single patient-donor pair may refuse to go through with a particular transplantation. We consider the fairness issue at this point and define fairness in different ways. We propose the notion of personal fairness, something which previous studies have not considered much. Personal fairness is the extent to which individual patient-donor pairs feel satisfied with the exchange. We propose a value indicator, unfairness, to represent the level of satisfaction. We give the penalty value to two representative examples of node failure and arc failure to reflect their effect in the second stage of the SP framework.
KEP is typically designed as follows.
Figure 1 illustrates how incompatible patient-donor pairs can be organized to allow successful surgery to take place by exchanging transplants.
Figure 1a illustrates the process by which two incompatible pairs exchange kidneys. If the donor of Pair 1 and the patient of Pair 2 are compatible and the donor of Pair 2 and the patient of Pair 1 are compatible, surgery can take place. Likewise, three incompatible patient-donor pairs can also be matched by exchanging kidneys (
Figure 1b).
The general KEP model can be expressed by directed graph
referred to as a
KEP graph—with set of vertices
and set of arcs
(
Figure 2). Vertex
represents an incompatible patient-donor pair and arc
represents compatibility between pairs
.
represents the weight of arc
, indicating the likelihood of surgical success.
Figure 2a presents possible links as dotted lines, while
Figure 2b presents a constructed cycle with vertices
and arcs
. A deterministic model for a KEP finds the set of cycles with the maximum weight using objective functions to maximize the sum of the weights.
The following are some of the studies that considered deterministic mathematical models. Roth et al. [
13] proposed both edge and cycle formulations for KEPs. Abraham et al. [
14] proposed a clearing algorithm that solves established models. Constantino et al. [
15] suggested two compact formulations, an edge-assignment (EA) formulation and an extended edge (EE) formulation. Yuh et al. [
16] applied the reformulation-linearization technique (RLT) to the EA and EE formulations for deriving a tighter and more compact formulation. The models in these previous studies use objective functions to increase the utility of the program as a whole. However, these models do not reflect equity as experienced by the individuals involved, nor did they attempt to incorporate the risk of failure.
Past research has identified several different perspectives on fairness in terms of KEPs. For example, Bertsimas et al. [
17] raised the issue of defining fairness in relation to kidney exchange transplants. They suggested that the concept of fairness can vary depending on the focused target, whether it be the individual or society as a whole. Indeed, most previous studies have defined fair exchanges as being those that increase the overall social benefit by prioritizing difficult-to-match pairs [
18]. This method achieves high levels of social fairness but fails to reflect the characteristics of bartering, in which individuals feel the exchange is fair when they are satisfied with it.
The risk of not considering personal fairness is that it can directly lead to the failure of surgery. A single patient-donor pair may refuse to go through with a particular transplantation when they feel unfair, even though the matched kidney is compatible, in order to wait for a better one. In fact, many kidneys that are medically compatible with some of the patients on transplant waiting lists are left unselected and unused [
19,
20]. Sabouri et al. [
21] suggested screening strategies for identifying patients who are ineligible for transplants. These studies imply that there are kidneys which are not suitable for transplant. The attractiveness of the exchange depends on the condition and the value of the kidney; based on these criteria, the patient-donor pair make a choice without considering the medical situation. Therefore, by considering the satisfaction of the individuals involved, it is possible to minimize the amount of unused kidneys and to reduce matching failure.
It is important to minimize matching failure in a KEP because a single failure can have a flow-on affect that disrupts other possible transplants as well. Previous studies have focused on failure in the context of the kidney exchange problem. For example, Zheng et al. [
22] were the first to consider arc failure, while Dickerson et al. [
23] mathematically outlined the impact of failure on exchange cycles and chains. Alvelos et al. [
24] proposed a compact integer programming model considering failure probabilities and they assumed that all arc and vertices have the same failure probabilities, while our study specified and considered individual pair’s characteristics and the effect of failure with the penalty value. Our work has uniqueness and advantage in its specified approach.
There have been some studies applying a stochastic and heuristic approach. Awasthi et al. [
25] proposed online stochastic optimization with trajectory-based algorithms. Dubey et al. [
26] utilized ant-lion optimization algorithm as a meta-heuristic approach to solve KEP.
Stochastic programming which is mentioned earlier can be a key to the limitation of uncertainty issues by utilizing different objective functions at each stage. A more robust model can be developed through using an objective function that increases total utility in the first stage, and an objective function that decreases the penalty for failure in the second stage. Because of the presence of the two objective functions, we expect two measurable improvements in this model. One is robustness to exceptional events, and the other is individual satisfaction while maintaining social utility.
In this paper, we propose a two-stage SP model that considers exceptional events that have not been included in previous deterministic models. We develop this two-stage SP model based on the compact EA model introduced by Constantino et al. [
15]. We redefine fairness as personal fairness in terms of an individual pair’s satisfaction and establish an indicator of unfairness. We classify failure scenarios and incorporate them into the second stage to obtain a more robust model. Computational results are provided to demonstrate the improvement of this model compared to the previous deterministic model. The contributions of our study can be summarized as follows.
We propose a stage-based SP model that reflects exceptional events before matching takes place. We consider failure scenarios in the second stage and classify failure scenarios as either arc failure or node failure. Experimental results illustrate that this stochastic model is more robust to failure than deterministic models when an exceptional situation occurs. This model thus proves that failure can be taken into account before matching, and this technique can be extended in many ways by adding failure scenarios.
We redefine and formalize fairness within a kidney exchange program as personal fairness, which is a measure of how satisfied individual pairs feel during an exchange. This is a unique approach compared to previous studies, which have defined fairness in terms of increasing social welfare. As such, this study proposes a new model that considers the nature of bartering. We propose an unfairness indicator and set a penalty function with this indicator as an input variable. The SP model that considers personal fairness produces lower total unfairness levels than that of existing deterministic models.
The rest of this paper is organized as follows.
Section 2 briefly introduces an overview of stochastic programming.
Section 3 then summarizes the two-stage stochastic programming approach to KEPs, illustrates how failure scenarios are classified and incorporated, and presents the stochastic formulation. Following this,
Section 4 presents computation results to assess the performance of our model. Finally, we conclude the study in
Section 5.
2. Stochastic Programming Overview
Stochastic programming is a framework used to model optimization problems that involve uncertainty. Unlike deterministic optimization problems, which are formulated with known parameters, stochastic programming incorporates unknown parameters. The general SP model can be presented as
Model (1) consists of two sections: first-stage and second-stage decisions. Here, is the first-stage decision variable, and is the second-stage variable. When the first-stage decision is made, after which a random event occurs affecting the outcome of the first-stage decision, a recourse decision can then be made in the second stage that compensates for any negative effects. is the cost of recourse action and is added up in the objective function. This second term in the objective function tries to repair the decisions which is made before the consideration of a random event. Through this procedure, SP can incorporate the first-stage decision’s effect as well as the second-stage decision’s effect. SP has a set of random data and random vector is used to represent random events with a certain probability. is the vector formed by components and and denotes mathematical expectations with respect to .
Figure 3 illustrates how scenarios are generated. In two-stage stochastic programming, random vector
can be represented with respect to scenario
. Each scenario
has the probability of occurrence
, and the sum of
with respect to
is 1. The second decision variable
appears in the objective function, and
and
, which are parameters for scenario
, appear in the mathematical formulation.
However,
in Model (1) can lead to a nonlinear model if
is continuous over
. Therefore, the model should be discretized for
as follows.
where
indicates the index of the scenarios, and
and
are the probability of scenario
, the parameter vector related to scenario
, and the second stage variable of scenario
, respectively. Model (2) is a compact formulation of a two-stage stochastic programming model. In this study, the scenarios are possible cases of failure in kidney exchanges. At the first stage, we consider an objective function that maximizes the sum of weights. At the second stage, we consider causes of failure and set the probability of each scenario
. We establish penalty value
with respect to
to incorporate the first stage’s objective function and the second stage’s failure scenarios. We derive a two-stage stochastic programming model in
Section 3 based on the well-known EA model by adding failure scenario information; we also demonstrate how the scenarios are classified and generated.
3. Two-Stage Stochastic Programming for Kidney Exchange Programs
3.1. Scenario Generation
In this section, we explain how we generate and classify failure scenarios. The KEP graph consists of two elements, nodes (vertices) and arcs (edges). Nodes represent an incompatible patient-donor pair and arcs represent compatibility between two pairs. Failure may occur between two pairs and within a pair. Therefore, failure scenarios for exceptional events are largely classified into two types of failure: node failure and arc failure. Node failure occurs if there is a problem in the patient-donor relationship, while arc failure occurs in the relationship between patient-donor pairs.
Figure 4 illustrates the two types of failure. Node
failure means that there is a problem within patient-donor pair
, which means that any compatible arcs related to node
become incompatible. Arc
failure means that there is a problem between pairs
and
. Therefore, if arc
failure occurs, the compatibility between donor
and patient
becomes incompatible.
In this study, two representative examples of these failure types are presented: a change of mind due to perceived unfairness for arc failure and the deterioration of a patient’s condition for node failure. However, there are various other causes of failure that could also be considered, as shown in
Figure 5.
A scenario
is a holistic and conceptual factor for possible cases of failure. It can be generated based on the classification that we propose and can be presented as below:
where
is a scenario of arc failure and
is a scenario of node failure on arc
. For notational simplicity, we simply denote
and
respectively in case that it is evident to understand the arc
associated with a scenario. In a mathematical model,
and
are realized based on the scenario
, as mentioned in
Section 2. The sum of the probabilities of these failures, that is
, should be 1.
is a penalty value given to each node and arc depending on the scenario
. In this study,
and
are proposed and calculated reflecting characteristics of their failure scenarios.
is the second-stage decision variable depending on the scenario
.
● Node failure
The deterioration of the health of the patient in a patient-donor pair is selected as a representative example of node failure. If the health of the patient is so poor that it is dangerous to perform the operation, the solution set may be broken. Other possible causes of node failure include the patient’s death and financial problems. We introduce a penalty value based on the patient’s health condition. This value is only affected by the health status of the patient in individual pair. The penalty value is given as follows.
a penalty value of arc depending on the health status of the patient in pair regardless of the health status of the donor in pair in scenario for .
● Arc failure
A change of mind due to perceived unfairness is selected as a representative example of arc failure. If a pair feels that the exchange is unfair, they may change their mind and thus the solution set may be broken. In this paper, we have assumed that all patients and donors have certain health conditions. Therefore, based on the health status of the patient and donor, the level of compatibility was calculated. Other possible causes of arc failure are shortages of operating rooms and a change of mind due to long distance between two pairs. We introduce unfairness indicator
for arc
as follows:
where
is the health status of the donor in pair
,
is the level of compatibility between pairs, and
is the value calculating unfair feelings between the relationship of pair
and pair
. From pair
’s perspective, it can be calculated by dividing the value that pair
receives into the value that pair
gives. Note that the larger the value of
and the smaller the value of
, the greater the value of
. It is natural to feel an exchange is unfair if the value of donating a kidney is greater than the value of receiving one. We introduce a penalty value based on this perceived unfairness. The following definition and equation is used to determine the penalty value.
: a penalty value of arc
depending on the unfairness indicator
in scenario
for
.
We remark that these parameters, and , have negative values imposed to penalize the effects of failure scenarios. When it comes to , the worse the patient’s condition, the higher the penalty. For , as unfairness increases, the absolute value of this negative number increases exponentially. This means that the effect of penalty increases when an individual pair j feels more unfair. Note that is divided by the constant to scale the penalty value in the range of −1 and 0.
3.2. Model
In this section, we present as a base model the reduced edge assignment (EA) formulation first developed by Constantino et al. [
15]. Our proposed stochastic programming model is then developed by adding new decision variables to this base model. These decision variables are related to the scenarios discussed earlier.
Let be an upper bound on the possible number of cycles in graph . The weight of edge indicates the level of compatibility between the kidney from the donor in patient-donor pair to the patient in pair for . Parameter represents the maximum cycle length. The decision variables in the EA formulation are defined as follows:
: 1 if the patient of pair receives a kidney from the donor of pair , and 0 otherwise,
: 1 if pair is included in cycle , and 0 otherwise.
With these decision variables, the reduced EA formulation can be presented as follows:
Deterministic model (EA formulation):
The objective function (3a) indicates that the total weights of the arcs involved in a kidney transplantation is maximized. This leads to a set of transplantations that maximizes the total compatibility of the given patient-donor pool to be selected by definition of weight
. Note that if every
is set to 1, the model then maximizes the number of transplantations in the pool. Constraint (3b) states that the number of kidneys that the donor in pair
provides must be equal to the number of kidneys that the patient in pair
receives. Constraint (3c) ensures that the maximum possible number of kidneys that a donor can donate is 1. Constraint (3d) states that the length of a cycle in the pool must be less than or equal to
. Constraint (3e) assigns a transplantation to a cycle. Constraint (3f) states that, if the donor in pair
gives the patient in pair
a kidney and pair
is included in cycle
, then pair
must also be included in cycle
. Constraints (3g) and (3h) are additional constraints that eliminate duplicate solutions, as discussed in detail by Constantino et al. [
15]. Constraints (3i) and (3j) set the decision variables
and
as binary variables.
Based on the above formulation, we now develop a stochastic programming formulation for a KEP to consider failure within the KEP graph. As defined in
Section 3.1, for each scenario
, we let
be the probability of the failure scenario occurrence, and
be the penalty value given to each node and arc depending on the scenario
We specify the additional decision variable
for the two-stage stochastic programming model which is given as follows:
With the above decision variable, the two-stage stochastic programming model can be expressed as
: 1 if the patient of pair
receives a kidney from the donor of pair
in scenario
, and
otherwise.
where
is the average
over scenarios
, and
is the optimal value of the following second-stage problem:
The objective function (5a) in the second-stage problem returns the total penalty for the possible failure of a solution in scenario . Constraint (5b) ensures that, in scenario , the pairs which were matched in the first stage are considered. The common constraint (5c) sets the decision variable as a binary variable.
We remark that we can easily derive the deterministic equivalent version of the two-stage stochastic problem with discretized scenarios where the corresponding probabilities are . (5a)–(5c) can be specified as below:
Finally, our proposed model, referred as a stochastic model, can be simply formulated as follows:
5. Conclusions
In this paper, stochastic programming is proposed as a method for incorporating uncertainty into KEP models. Based on the compact EA formulation, we develop a two-stage SP model by adding second-stage decision variables and parameters. In order to consider exceptional events in the model, we classify specific failure scenarios based on two types of failure: node failure and arc failure. The deterioration of the patient’s condition and a change of mind due to perceived unfairness are chosen as representative examples of node and arc failure, respectively. While previous studies have defined fairness as social fairness, we define it as personal fairness and consider how satisfied individual pairs feel. An indicator of unfairness is introduced; as this number increases, the penalty value imposed in the second stage increases. Computational simulations are run to compare an existing deterministic model with our proposed stochastic model. From these, it is shown that the stochastic model reduces overall perceived unfairness while maintaining total utility. We also demonstrate that, when the two failure scenarios occur, the stochastic model is more robust in terms of preserving the solutions. We believe that the proposed model could help to be utilized in the integrated KEP IT healthcare platform in South Korea to obtain optimized fair KEP exchange plans.
In future research, the two-stage SP model proposed in this paper can be tested with other types of failure, thus leading to a more robust model. Bertsimas et al. [
28] suggested the widely used Robust Optimization (RO) model, which made it possible to generate robust counterpart of Linear Programming (LP) problems. Büsing et al. [
29] also suggested robust counterpart of LP with uncertain coefficient matrix. These two approaches provide insight for future work because SP and RO have overlap in that they can be linearized. Robust models for KEPs have been presented in various perspectives by Dickerson [
30]. D’Andreagiovanni et al. [
31] considered data uncertainty in topology and applied RO. They applied RO to the Network Problem, in which our work belongs to. RO in KEP or Network Problem would be the topic for the future work.
Because a KEP is related to surgery, which is associated with significant time and money costs, robustness to failure is an important consideration. Our study only considers a static pool, but other studies have analyzed pools that move dynamically over time [
32,
33]. Therefore, even more promising results may be achieved if we incorporate patient-donor pairs that flow in time.