Next Article in Journal
Generalized Data–Driven Predictive Control: Merging Subspace and Hankel Predictors
Next Article in Special Issue
Approximate Subdifferential of the Difference of Two Vector Convex Mappings
Previous Article in Journal
Improving Mechanical Oscillator Cooling in a Double-Coupled Cavity Optomechanical System with an Optical Parametric Amplifier
Previous Article in Special Issue
Application of an Intensive Longitudinal Functional Model with Multiple Time Scales in Objectively Measured Children’s Physical Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Data-Driven Decision-Making Model for Configuring Surgical Trays Based on the Likelihood of Instrument Usages

1
Stetson-Hatcher School of Business, Mercer University, Atlanta, GA 30341, USA
2
Department of Engineering Education, College of Engineering, The Ohio State University, Columbus, OH 43210, USA
3
Enterprise Supply Chain Services, Geisinger Health, Danville, PA 17822, USA
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(9), 2219; https://doi.org/10.3390/math11092219
Submission received: 19 March 2023 / Revised: 20 April 2023 / Accepted: 2 May 2023 / Published: 8 May 2023
(This article belongs to the Special Issue New Advance in Operations Research and Analytics)

Abstract

:
In order to perform a surgical procedure, substantial numbers of sterile instruments should be readily available to surgeons through the containers referred to as surgical trays and peel packs. After the procedure, all instruments in the opened containers, regardless of whether they have been used or not, must go through the labor-intensive re-sterilization process. Empirical studies have shown that the utilization rate of instruments within trays is very low due to not having optimized tray configurations. Additionally, surgical trays often include instruments that are not likely to be used but are included “just in case”, which imposes an additional cost on hospitals through unnecessary instrument re-sterilization. This study is the first analytical attempt to address the issue of configuring surgical trays based on the likelihood of instrument usage through formulating and solving a probabilistic tray optimization problem (PTOP). The PTOP model can serve as a decision support for surgeons by providing them with the tray’s probability of being used for optimally configured trays and its associated reprocessing costs. The PTOP is constructed upon an integer non-linear programming model. A decomposition-based heuristic and metaheuristic method coupled with two novel local search algorithms are developed to solve the PTOP. The application of this model can be illustrated through a case study. We discuss how hospitals could benefit from our model in reducing the costs associated with opening trays unnecessarily before a procedure. Additionally, we conducted a risk analysis to estimate the level of confidence for the recommended solution.

1. Introduction

In the United States, health expenditures are projected to grow by 5.4 percent annually to reach $6.2 trillion by 2028 [1]. Hospitals account for about one-third of health expenditure [1], and operating rooms are the main cost drivers within hospitals [2,3]. An improvement in managing sterile instruments and disposable surgical supplies is one area where operating rooms could reduce costs. Surgeons request these items on their preference cards for performing procedures. The surgeons’ preference cards also indicate which surgical supplies must be opened before a procedure. Several recent studies have introduced methodologies for managing disposable surgical supplies based on historical data and determining the quantity to include in a surgeon’s preference card to be opened before procedures [4,5]. For sterile instruments, the ones that surgeons require to perform a procedure must have already been packed into surgical trays and peel packs. However, the trays are not surgeon specific, and a tray type may contain instruments requested by multiple surgeons, which leads to some instruments remaining unused during a procedure. After the procedure, if a tray is opened, all instruments in the tray, regardless of whether the instruments are used or not, must proceed to the Sterilization Processing Department (SPD) for cleaning, inspections and grouping into trays. Our research focuses on managing sterile instruments.
There is evidence that shows a high proportion of instruments in trays are not used during the procedures, which causes millions of dollars of cost burden on hospitals. Stockert and Langerman [6] observed that about 82.9% of instruments had not been used across four busy surgical services (i.e., Bariatric Surgery, Plastic Surgery, Otolaryngology, and Neurosurgery). In another study, Mhlaba et al. [7] reported that the utilization rate for the plastic soft tissue tray and the major laparotomy tray was 14% and 29%, respectively. Koyle et al. [8] have also indicated a utilization rate of 42% for the pediatric inguinal hernia repair tray.
A direct observation of tray usage has been conducted by the authors of this manuscript at the SPD of Geisinger Health System for the Laparoscopic and Plastic major trays. The Geisinger Health System is an integrated healthcare organization located in Pennsylvania, USA that has 12 hospitals serving more than three million people. For a single observation of Laparoscopic and Plastic trays, the utilization rates were 25% and 18%, respectively. Figure 1 illustrates a Plastic tray with 14 of the 79 instruments used during a procedure. Since the operations of the SPD are very labor-intensive, the cost associated with the unnecessary processing of instruments is significant. Farrokhi et al. [9] reported an annual cost saving of $60,000 by reducing 70% of instruments in the trays for two procedures. Then, they projected that the institution could achieve a $2.8 M annual saving by similar instrument reductions in all trays. Harvey et al. [10] also reported that two procedures of the female pelvic medicine and the minimally invasive gynecologic could annually save $151,691 by safely removing unused instruments from their surgical trays. Nast and Swords [11] indicated that improvements in the management of surgical trays in a way that more than 50% of their instruments could be reduced could lead to at least 20% of cost savings.
Better management of surgical trays would not only result in immediate cost savings but also allow for indirect cost reductions. The frequent handling of the trays in the SPD and operating rooms is the main reason for the development of Work-Related Musculoskeletal Disorders (WRMSDs) among Perioperative Nurses and Technicians (PNT) [12,13,14,15]. The occurrence of WRMSDs diminishes the productivity of the PNTs, which has financial implications for hospitals through the absenteeism of PNTs from work. In addition, improved surgical trays could allow for cost savings through a reduced setup time required for counting the instruments before and after the procedures [8,11], as the average cost of running an operating room is estimated to be $63.64 per minute [16].
For improving the management of surgical trays and their sterile instrument components, a hospital has to make three types of decisions: (1) determine the configuration of the trays regarding the type and the number of instruments; (2) assign the configured trays to surgeons and procedures; and (3) determine the quantity of each tray type [2,17,18]. The two former decisions depend on the surgeon’s preferences in the type and quantity of requested instruments. The last decision is driven by the frequency and scheduling of the cases. This problem of configuring surgical trays is known as the tray optimization problem (TOP).
Despite the importance and significant impact of managing surgical trays on hospital costs, there have been relatively few studies devoted to addressing it. Van De Klundert et al. [19] developed an integer linear programming model for the TOP and proved that the TOP is an NP-hard problem. Florijn [20] focused on developing an optimization model for the TOP with the objectives of minimizing the number of trays, the handling of trays, and the number of instruments. Reymondon et al. [21] developed a non-linear mathematical model to solve the TOP. Dobson et al. [22] designed a heuristic for solving the TOP and showed that the optimal solution of the TOP depends on both the surgeon’s preferences and the operation rooms’ schedules. Dollevoet et al. [17] also formulated the TOP and defined capacity constraints to incorporate the capacity of trays into modeling. They proposed three solution methodologies: row and column generation, greedy heuristics, and genetic algorithm. Ahmadi et al. [2] designed a bi-objective optimization model for the TOP to concurrently enhance the utilization rate of the instruments and decrease the tray types required for each procedure. They took the weight of the trays into account to avoid composing overweight trays to decrease the risk of WRMSDs. They also developed an iterated local search heuristic for solving the designed model. For a comprehensive review of sterile instrument management papers, readers are referred to Ahmadi et al. [18] and Dos Santos et al. [23]. Harris and Claudio [24] formulated the TOP using a goal programming approach. Most recently, Deshpande et al. [25] developed an integer programming model and used real instrument usage data to address the TOP.
All aforementioned papers developed deterministic models for the TOP based on the assumption that the instances of instruments required for performing a procedure are precisely known in advance. However, this is not true in practice. Even if an instrument is requested by a surgeon for a procedure, it may be unneeded for a specific patient. Instruments are included to ensure the quality of care, which leads to having to reprocess instruments that were not used during a procedure, which imposes a high cost on the system.
A recently published paper utilized the empirical usage distributions of disposable surgical supplies to determine the optimal quantity to be included on the surgeon’s preference cards and to be opened before the procedure [4]. The authors used a newsvendor approach to formulate the problem of designing surgeons’ preference cards and evaluated their model under different service levels. However, this research does not address the configuration of surgical trays. In our research, we considered the fact that the true usage of instruments highly depends on the surgeon’s preferences, as well as the condition of patients during the procedure. Therefore, we extended the conventional deterministic TOP to the probabilistic tray optimization problem (PTOP), which enabled us to configure trays based on the likelihood of instrument usage. To this aim, the cost components (i.e., the cost associated with inventory, sterilization, and processing of instruments), along with the historical data about instrument usage, could be incorporated to configure trays. To solve the PTOP, we constructed a heuristic algorithm based on the well-known p-median problem. In addition, we developed an empowered genetic algorithm with two local searches to solve the PTOP. The performances of different solution approaches were investigated through a case study.
The designed probabilistic optimization method captures instruments that are never (or rarely) used and has them grouped separately from the instruments with a high chance of usage. It provides surgeons with the tray’s probability of being used and the associated costs. This information, along with case consideration, enables surgeons to better decide whether to open a tray before a procedure or to only open the tray if it is needed during a procedure. Since hospitals incur high costs for the supplies and instruments to function properly, and a great quantity of these items are left unused, the implementation of the proposed model is expected to have a broad financial impact by eliminating unnecessary healthcare costs.
The rest of this paper is organized as follows. In Section 2, the PTOP is described in detail. In Section 3, the mathematical model describing the PTOP is formulated, followed by solution methods presented in Section 4 and Section 5. In Section 6, the incorporated case study, along with the computational results, are reported. In Section 7, managerial insights are discussed. Finally, the conclusion and discussion about implementing our model are covered in Section 8.

2. Problem Description

In order to perform a procedure in an operating room (OR), surgeons request a list of sterile instruments, which have already been grouped into trays and peel packs, to be available in the OR. The trays and peel packs associated with these sterile instruments are included on the surgeons’ preference cards. A collection of surgeons’ preference cards can be shown using a matrix whose rows and columns represent surgeon procedures and instruments, respectively. The elements of this matrix indicate the number of each instrument a surgeon has requested to perform a given procedure. A sample of the matrix constructed based on two surgeons and three procedures (i.e., six surgeon-procedure combinations) and five instruments has been shown in Table 1. For example, surgeon-procedure 5 requested one instance of instrument 3, while surgeon-procedure 6 requested two instances of instrument 3. Therefore, if two instances of instrument 3 are placed in the same tray type, and the tray type is assigned to surgeon-procedure 5 and surgeon-procedure 6, every time surgeon-procedure 5 performs that particular procedure, one instance of instrument 3 that was not requested would be sent in the tray. These unrequested instruments cause unnecessary reprocessing such as inspection, re-sterilization, packaging, and handling, all of which can be reduced through an improved tray configuration.
In addition, not all of the instruments requested by surgeons are used. Whether an instrument is used for a procedure depends on the patient’s condition during the procedure, so surgeons often request extra instruments to be included in the trays just in case. By grouping instruments that are unlikely to be used in separate trays from those that are frequently used, surgeons can decide whether to open the trays before or during the procedure. Not opening the trays that contain less frequently used instruments before the procedure can lead to the decreased reprocessing of unused instruments.
The above-discussed issues can be addressed by solving the probabilistic tray optimization problem (PTOP), which consists of two decisions: configuring surgical trays based on the likelihood of instrument usage and then assigning the configured trays to the surgeon procedures. In order to formally describe the PTOP, we can define the usage of the j th copy of instrument i for surgeon-procedure k as a random variable with P i j k probability of being used. In other words, there is a P i j k chance that the j th copy of instrument i will be used during surgeon-procedure k , and 1 P i j k chance that the instrument will not be used. Throughout the rest of this paper, for the sake of simplicity, the procedure is used instead of the tuple surgeon-procedure so that each surgeon performing a given procedure is treated as a unique procedure.
To show the probability of usage associated with each instance of instruments, the matrix describing the PTOP (i.e., Table 1) is converted to a zero-one matrix, namely a procedure-instrument incidence matrix (PIIM), with the same number of rows as the original PTOP matrix. Each column i in the original PTOP matrix is converted to N i columns in the PIIM, where N i = max k { A i k } and A i k is the requested quantity of item i for procedure k . The entries in the PIIM, a i j k , are equal to “1” if the j th copy of instrument i is requested for procedure k and 0 otherwise, where a i j k a i j k if j < j . The PIIM corresponding to the PTOP example shown in Table 1 is provided in Table 2.
Each column of the PIIM stands for an instance of an instrument, which can be treated as a single instrument. For example, the first column of PIIM stands for the first instance of instrument 1, the second column stands for the second instance of instrument 1, the third column stands for the first instance of instrument 2, and so on. One of our assumptions in this study is that each surgeon performed a procedure only one time in a given day, and therefore, instruments are not reused within a day. Given this assumption, it is obvious that for each instrument i , at least N i instances needed to be stocked in an inventory to ensure that what surgeons requested on their preference cards would be available within a tray in the OR. For the instruments and procedures in the PIIM, we could show the probability of each instrument’s usage, denoted by P i j k , as presented in Table 3, where P i j k P i j k if j < j . For example, Table 3 shows that the first instance of instrument 1 had an 80% chance of being used during procedure 3, while its second instance had a 15% chance of usage during the same procedure.
The specific question to be answered in this research is how to cluster instruments into containers (i.e., trays and peel packs) to decrease the expected number of instruments that were reprocessed without being used. A container is defined as a “peel pack” if it holds a single instrument by itself; otherwise, the container can be referred to as a “tray type”. The next section describes the formulated mathematical model followed by a solution approach for answering this question.

3. Mathematical Formulation

We assumed that instrument usages were mutually independent and did not account for correlated instruments that are often used together. In the absence of data, we were not able to directly verify this assumption.
  • Indices
i : Index for instrument, i { 1 , 2 , 3 , , I } .
i j : Index for j th copy of instrument i , j { 1 , 2 , 3 , , N i } where N i is the maximum number of instruments of type i that can be required in one procedure.
k : Index for procedure, k { 1 , 2 , 3 , , K } .
t : Index for container (i.e., tray or peel pack), t { 1 , 2 , 3 , , T } .
  • Parameters
a i j k : The elements in the PIIM.
P i j k : The probability that jth copy of instrument i will be used during procedure k .
W i : Unit weight of instrument i .
F k : Frequency of procedure k performed in a year.
L : Weight limit for each tray type established by AAMI and AORN.
C 1 : The unit cost of reprocessing an instrument in a tray.
C 2 : The unit cost of reprocessing an instrument in a peel pack.
C 3 : The unit cost associated with administration and handling of a tray type.
C 4 : The unit cost associated with administration and handling of a peel pack.
M : A big number.
  • Objective function
M i n   f = f 1 + f 2 + f 3 + f 4
f 1 :   C 1 [ t = 1 T ( 1 e t ) k = 1 K y t k F k ( i = 1 I j = 1 N i z i j t ) ( 1 i = 1 I j = 1 N i ( 1 z i j t P i j k ) ) ]
f 2 : C 2 [ t = 1 T e t k = 1 K y t k F k i = 1 I j = 1 N i z i j t P i j k ]
f 3 :   C 3 t = 1 T k = 1 K F k x t y t k ( 1 e t )
f 4 :   C 4 t = 1 T k = 1 K F k x t y t k e t
This is subject to:
i = 1 I j = 1 N i z i j t W i x t L ,       t { 1 , , T } ,
t = 1 T y t k z i j t a i j k ,     i { 1 , , I } ,   j { 1 , , J } , k { 1 , , K } ,
i = 1 I j = 1 N i z i j t ( 1 e t ) M + 1 ,       t { 1 , , T } ,
i = 1 I j = 1 N i z i j t 2 2 e t ,       t { 1 , , T } ,
t = 1 T z i j t = 1 ,       i { 1 , , I } , j { 1 , , J } ,
z i j t { 0 , 1 } ,   x t { 0 , 1 } ,   y t k { 0 , 1 } ,   k { 1 , , K } , i { 1 , , I } , e t { 0 , 1 } , t { 1 , , T } ,   j { 1 , , N i } .  
In the first component of the objective function, the term i = 1 I j = 1 N i ( 1 z i j t P i j k ) indicates the probability that none of the instruments in tray type t are used for procedure k , i.e., the probability that tray type t remains unused. Therefore, ( 1 i = 1 I j = 1 N i ( 1 z i j t P i j k ) ) describes the likelihood of tray t being used during a procedure k . Multiplying this expression by the total number of instruments in the tray (i.e., i = 1 I j = 1 N i z i j t ) yields the expected number of instruments that are exposed during the procedure. Likewise, i = 1 I j = 1 N i z i j t P i j k describes the expected number of instruments in peel packs that are exposed during the procedure. Thus, component (1a) is the total expected reprocessing cost of the instruments in trays, including the cost of re-sterilization, the cost of inspection, and the cost of packing instruments into trays. Component (1b) is the total expected reprocessing cost of instruments in peel packs. Component (1c) is the total administrative cost associated with handling and stocking tray types. Creating more tray types as a result of separating the existing trays into additional smaller trays induces higher logistics costs that are related to stocking and transporting them. Finally, component (1d) is the total administrative cost of peel packs. It should also be noted that the model assumes that the reprocessing costs only depend on whether an instrument is in a tray or a peel pack, and the costs associated with instruments that are in the same tray are the same. Even though the sterilization time varies with the type of instrument, it is very difficult to quantify the time and, therefore, the cost of sterilizing every single instrument.
Constraint set (2) imposes a weight limit for trays to not exceed a user-specified value. In order to prevent or reduce ergonomic, sterilization, and drying issues, the Association for the Advancement of Medical Instrumentation (AAMI) recommends that the weight of surgical trays not exceed 25 pounds [26]. Constraint set (3) respects the surgeons’ preferences in the type and quantity of instruments. Constraint sets (4) and (5) can determine whether an instrument is in a peel pack or placed in a tray. It can be recalled that the peel packs and tray types are distinguished by the number of instruments in the container. If two or more instruments are grouped together to be placed in a container, this is referred to as a tray, but if a single instrument is placed in a container by itself, this is defined as a peel pack. Constraint set (6) ensures that each instrument has only been assigned to one tray. Constraint set (7) imposes binary conditions of decision variables. Since the objective function (1) and constraint set (3) involve non-linear terms, this model is an integer non-linear programming model (INLP).
Lemma 1. 
The PTOP described by (1)–(7) is NP-hard.
Proof of Lemma 1. 
Karp [27] proved that the set-covering problem is an NP-hard problem. Here, we showed that the set-covering problem was not more complex than the deterministic version of the PTOP. This meant that the set-covering problem was reducible to PTOP in polynomial time complexity.
We defined an instance of a set-covering problem by a universe U = { 1 , , n } and a set S containing subsets of U . The objective was to find the cover C , where C S , had the smallest cardinality such that its union was U . Then, we defined U as the set of instruments, S as the collection of all possible containers, and a single procedure K = 1 .
For simplicity, instruments of the same type were assumed to be indexed sequentially. For example, the first copy of instrument 1 was indexed 1, the second copy was indexed 2, and so on. We also set C 3 = C 4 = 0 , C 1 = C 2 = 1 , and a i 1 = 1 for i U . For a given container t S , z i t = 1 if i t   . Thus, for a given container t S , z i t , e t   , and x t become a parameter instead of a variable.
We also assumed that P i 1 = 1 for all instruments i U (i.e., deterministic version of the PTOP), and F 1 = 1 . The PTOP can be written as:
M i n f = t S ( 1 e t ) y t 1 i U z i t + t S e t y t 1 i U z i t
For a given t S if | t | = i U z i t = 1 then e t = 1 (i.e., the container is a peel pack), and if | t | = i U z i t 2 then e t = 0 (i.e., the container is a tray). Thus, the objective function becomes:
M i n f = t S y t 1 i U z i t
This is subject to:
t S y t 1 z i t 1 , i U ,
y t 1 { 0 , 1 } , t S .
This matches the ILP formulation of the set-covering problem. □
Lemma 1 implies that the PTOP cannot be solved in a reasonable time with any INLP solver, even for small sizes. Given this challenge, a heuristic and a metaheuristic algorithm were developed to solve the formulated PTOP. To this aim, the PTOP can be decomposed into two sub-problems, namely the probabilistic tray configuration problem (PTCP) and the probabilistic tray assignment problem (PTAP).
The PTCP deals with the problem of grouping instruments into trays and/or peel packs. More specifically, the PTCP determines which instruments and in what quantities should be grouped into the same tray, as well as which instruments should be wrapped as a peel pack. In the next step, the PTAP utilizes the solution of the PTCP and determines which tray types and peel packs should be assigned to each surgeon and procedure to provide all the instruments listed on the preference card. In the following sections, we design and discuss a heuristic and a metaheuristic for dealing with PTCP. These approaches determine the values for three out of four decision variables of the PTOP (i.e., z i j t , x t , and e t ). In order to cope with PTAP, which is the problem of determining values of the y t k decision variables, the solution generated by solving the PTCP can be fixed to the PTOP so that the PTAP becomes the trivial problem of solving an integer linear programming with objective function (1) subject to the constraint set (3).

4. Heuristic for Solving PTCP

The clustering heuristic provides an approximate solution to the PTCP by formulating the PTCP as a p -median model. The p -median model has a long history of use in several problems, such as location-allocation problems [28,29] and cell formation problems [30,31]. In this approach, the PTCP is formulated as an undirected graph, where each node represents an instrument, and each arc represents the expected number of exposed instruments when the respective two instruments in an arc are placed in the same tray. We wish to select some of the instruments as median vertices (to represent containers) and assign other instruments to these containers in such a way that can minimize the total distance between the selected containers and other instruments. We can define an undirected complete graph G = ( V , E ) , where the set of vertices V is the set of instances of the instruments ( | V | = i = 1 I max k { A i k } ). Additionally, there are | V | potential containers to which instruments can be assigned. Two vertices i j and m n (i.e., jth copy of instrument i and n th copy of instrument m , respectively) are connected by an arc ( i j , m n ) E . Associated with each arc ( i j , m n ) , there is a distance d i j m n . The d i j m n = k = 1 K F k 2 ( 1 ( 1 p i j k ) ( 1 p m n k ) ) which indicates the expected number of instruments i j and m n that are exposed in the required procedures−and therefore need to be re-sterilized−when two instruments i j and m n are placed in the same container and are assigned to all procedures. It should be noted that if the container holding two instruments i j and m n is assigned to a procedure k , for which these instruments are not requested (i.e., p i j k = 0 and p m n k = 0 ), then it contributes a value of 0 in the d i j m n calculation. The PTCP decides which instruments are stored in each of the predetermined P containers, which might be trays or peel packs, and assigns instruments to these P containers. The formulated integer programming model for solving PTCP in the p -median context is provided below.
M i n i = 1 I j = 1 N i m = 1 I n = 1 N m d i j m n h i j m n
m = 1 I n = 1 N m h i j m n = 1 ,   i { 1 ,   ,   I } ,   j { 1 ,   ,   N i } ,
m = 1 I n = 1 N m h m n m n = P
i = 1 I j = 1 N i h i j m n W i h m n m n L ,   m { 1 ,   ,   I } , n { 1 ,   ,   N i } .
where h i j m n is a binary decision variable, which is equal to one if the jth copy of instrument i and the n th copy of instrument m are placed in the same tray type and is equal to 0. The set of decision variables { h i j m n } m = i , n = j indicates the candidate containers that need to be created. The d i j m n values can be collected into a symmetric distance matrix whose diagonal values { d i j m n } m = i , n = j indicate the expected number of instruments to be re-sterilized when the instrument is placed in a container by itself (i.e., peel packed). The objective function (8) minimizes the sum of the expected number of instruments to be re-sterilized. The constraint set (9) ensures that each instrument is assigned to exactly one potential container. The constraint set (10) ensures that P containers are assembled, and constraint set (11) prevents tray configurations that violate the weight limit.
Solving the formulated p -median problem provides an approximate solution to the PTCP, and in the next step, by fixing the decision variables corresponding to the PTCP in the model described in Section 3, we can solve the PTAP. However, the p -median model only accounts for a pre-specified number of containers. Therefore, the p -median model and then PTAP can be solved iteratively by choosing the extreme value for P = i = 1 I max k { A i k } , which implies that all instruments are peel-packed and decreases the value of P by one in each iteration until no feasible solution can be generated for the selected P . A summary of the proposed p -median based heuristic is presented in Algorithm A1 of Appendix A.
In order to show how the p-median-based heuristic works, let us consider the PTOP example presented in Table 1. According to the PIIM (see Table 2), there exist 13 instances of instruments, each of which can be a candidate for a container. Therefore, the highest possible value for P is 13. Using the probabilities of instrument usage presented in Table 3, the values for d i j m n by assuming F k = 1 for all procedures k are shown in Table A1 of Appendix A. In the following, the results of solving the p -median model are reported, where the instruments have a unit weight ( W = 1 pound) and the capacity of the tray is restricted to not exceed 5 pounds. Therefore, each container can hold a maximum of five instruments, which implies that at least three containers are needed to be able to generate a feasible solution. It should be noted that the values of the weight of the instruments and the capacity of the trays for this small example were chosen arbitrarily just to show the applicability of the p-median model when solving the PTCP.
Figure 2 illustrates the obtained solution for the PTCP when four containers are considered (i.e., P = 4 ). As Figure 2 shows, the instruments are clustered into three tray types and one peel pack. The capacities of tray type 1 and tray type 2 are fully utilized as each contains five instruments. Tray type 3 contains two instruments of {1, 6}, and instrument 12 is placed in a peel pack.

5. Metaheuristic for Solving PTCP

The Genetic Algorithm (GA) is a kind of evolutionary algorithm that can solve intractable optimization problems based on stochastic search. Researchers have reported the successful application of GA to many different large space optimization problems, such as the shortest path [32,33] and inventory management [34,35,36]. The first step in implementing a GA is to encode a solution to the problem, denoted either as a chromosome or individual. Each chromosome is composed of genes, and each gene carries some information. The GA starts by generating a group of initial chromosomes. This group constitutes the population of the first generation. Each individual in the population is evaluated through the objective function and receives a fitness value. Individuals with a better fitness value have a higher chance of surviving for the next generation. Individuals are selected for removal to the mating pool for reproduction based on their fitness values. Reproduction is performed by the use of the crossover operators so that the generated child inherits some information from their parents. Mutation operators are incorporated to increase diversity among generations. In addition, to avoid the traditional GA from being stuck in a local optimum, we designed two local search algorithms called combining local search (CLS) and decomposing local search into the GA to improve its exploration capability.

5.1. Solution Encoding

A solution to the PTCP is encoded as a vector of integer numbers defining the index of containers to which instruments are assigned. The cardinality of the vector is equal to the total number of instances of instruments stocked in an inventory (i.e., i = 1 I max k { A i k } , which is the number of columns in the PIIM). The instruments are also presented in the order of the columns of the PIIM. For example, given the PIIM of Table 2 and 13 possible containers, including trays and peel packs, the vector [4,12,4,3,4,2,6,2,4,6,2,3,3] means that the first, the third, the fifth, and the ninth instruments are placed in container 4. The second instrument is in container 12 by itself. The fourth, twelfth and thirteenth instruments are placed in container 3. The sixth, eighth and eleventh instruments are in container 2. Finally, the seventh and the tenth instruments are in container 6. Therefore, containers 4, 3, 2, and 6 can be labeled as trays of instruments since they contain more than one instrument, while container 12 is a peel pack since it holds a single instrument. Thus, the above vector indicates that the instruments can be grouped into four tray types and one peel pack. A sample of solution encoding is illustrated in Figure 3.

5.2. Fitness Function

The fitness function can be used to assess the quality of chromosomes, of which chromosomes with better fitness values are more likely to survive. For a given chromosome which is a solution to the PTCP and determines values of z i j t , x t , and e t , the PTAP is solved to determine the values of the y t k decision variables. It can be recalled that the PTAP is an integer linear programming model of the objective function (1) subject to the constraint set (3). Finally, the fitness value corresponding to each chromosome is defined as 1 f , where f is the value of the objective function (1).

5.3. Genetic Operators

5.3.1. Selection Operator

This operator selects individuals in the current population for reproduction. There exist various selection operators such as tournament selection, random selection, and roulette-wheel selection. In this study, the roulette-wheel selection operator was chosen, in which the chance of individuals being selected for reproduction was proportional to their fitness value.

5.3.2. Crossover Operators

Two cut-point crossover operators: this operator randomly selects two genes as the cut-points within the selected parents and then interchanges the sections before and after these points to produce two new solutions [37]. The mechanism of this operator is shown in Figure 4.
Uniform crossover operator: this operator fills each gene of the offspring by copying the corresponding gene from one or the other parent with a 50% probability [38]. To do so, a binary vector with the same length as the parents can be generated using a uniform distribution. Each gene’s value in the offspring is copied from the first parent if the value of the corresponding gene in the binary vector is “1”; otherwise, the gene’s value of the offspring is copied from the second parent (see Figure 5).

5.3.3. Mutation Operators

Swap Operator: As illustrated in Figure 6, this operator swaps two arbitrary genes in a parent.
Inversion operator: this operator randomly selects two genes and then reverses the sequence of the section between the selected genes (Figure 7).
Shifts operator: this operator selects one gene to transfer it from its current position to the first position and shifts the rest of the genes forward (Figure 8).
Shuffle operator: this operator randomly selects two genes in a solution and then randomly updates the values of the section between these two genes (Figure 9).
If any of the mutation and crossover operators result in an infeasible solution by violating the weight limit constraint, a repairing mechanism can be performed to bring the solution into the feasible region. In the repairing mechanism, as long as there are trays that exceed the weight limit, the lightest instruments from those trays can be moved to the lightest trays to ensure that the weight limit constraint is met.

5.4. Combining Local Search (CLS)

The CLS attempts to reduce the number of containers in a given solution. For this purpose, some containers should be selected, and then their instruments should be combined together to become a single container. The selection of containers occurs based on their contribution to the objective function. Given that a container t holds I instruments, its contribution to the objective function is equal to the cost component f 1 (i.e., expression 1a) for the instrument trays and is equal to the cost component f 2 (i.e., expression 1b) for the peel packs, where the container is assumed to be assigned to all procedures (i.e.,   y t k = 1 ,     t { 1 , , T } and k { 1 , , K } ). Therefore, the contribution of a tray t , denoted as O t , can be calculated as follows:
O t = C 1 [ k = 1 K y t k F k ( i = 1 I j = 1 N i z i j t ) ( 1 i = 1 I j = 1 N i ( 1 z i j t P i j k ) ) ] ,       t { 1 , , T } |   t   is   a   tray ,
O t = C 2 k = 1 K F k i = 1 I j = 1 N i z i j t P i j k ,       t { 1 , , T } |   t   is   a   peel   pack .
In the next step, three walk points are considered in the search scheme as follows. Walk 1: combining the two containers with the lowest O t ; Walk 2: combining the two containers with the highest O t ; Walk 3: combining the container with the lowest O t and the container with the highest O t . The repairing mechanism can also be applied when these explorations yield a tray that violates the maximum weight of a tray. Among the visited solutions, the best one in terms of the objective function’s value can be added to the current population. This selected solution cab also be considered as the starting point for the next walk. To further explore the solution space and escape from the local optima, the CLS is enabled for random walks in a feasible solution space with the probability α . These steps are successively performed until the number of containers is reduced by γ percent. Algorithm A2 of Appendix A presents the CLS pseudo-code.
For example, it can be assumed that a solution [4,4,10,3,2,1,6,2,4,1,7,4,3] is selected for application to the CLS. This solution represents seven container types, which are indexed as 4, 10, 3, 2, 1, 6 and 7. Given the probabilities of instrument usage in Table 3, C 1 = 3 , C 2 = 2 , and F   k { 1 , ,   6 } = 10 , the following table shows the instruments in each container and the contribution of each container to the objective function (Table 4).
Now it can be assumed that walk 1 is intended to be explored. Since containers 7 and 2 contribute the lowest costs to the objective function, they are combined to reduce the number of containers. In this case, all instruments in container 7 can be moved to container 2, which results in a new solution S 1 as [4,4,10,3,2,1,6,2,4,1,2,4,3]. This new solution no longer includes container 7.

5.5. Decomposing Local Search (DLS)

In contrast with CLS, the decomposing local search (DLS) explores the solution space by creating new containers. Similar to the steps described for performing CLS, the DLS also consists of two steps. In the first step, the contribution of containers to the objective function can be calculated according to expressions (12) and expressions (13). The containers with the lowest and the highest contributions are the candidates to be decomposed. If the candidate container is a peel pack, then it can be excluded from the decomposition process. In the next step, given a container t , the contribution of each instrument to the objective function is calculated according to the expression (14).
O i = C 1 k = 1 K F k P i j k ,       i { 1 , , I }
Finally, the instruments with the lowest and the highest O i are removed from their current tray, and each is placed in a new container by itself (i.e., a peel pack). For the DLS algorithm, the same random walk strategy as described for the CLS algorithm can be implemented. The pseudo-code for the DLS algorithm is presented in Algorithm A3 of Appendix A. Since both CLS and DLS algorithms take the contributions of the trays and instruments into account and have random components, they are expected to exploit prior knowledge while exploring the solution space.

6. Experimental Design

6.1. Benchmark Problems

The performance of the proposed algorithms can be evaluated using three datasets introduced by Ahmadi et al. [2], as well as a new VLD dataset introduced in this paper for validation purposes. Each dataset contains five instances with the same number of surgeons and procedures. The usage probabilities of the instruments are randomly generated from a uniform distribution [0, 1]. The specifications of the datasets are presented in Table 5.
The small size instances in the VLD dataset are five variations of the example depicted in Table 1, which were constructed upon two surgeons, three procedures, five unique instruments, and thirteen copies of instruments. The weight of each instrument for this VLD dataset was assumed to be one pound with a tray capacity of 5 pounds. Additionally, the frequency of each procedure can be assumed to be one time per year. The unit cost of reprocessing an instrument in a tray (i.e., C 1 ) and the unit cost of reprocessing an instrument in a peel pack (i.e., C 2 ) were retrieved from Mhlaba et al. [7] at $0.40 and $0.80, respectively. The costs associated with the administration and handling of a tray type (i.e., C 3 ) and a new peel pack (i.e., C 4 ) were estimated as follows. Handling a tray was assumed to take 200 s, which included the time needed for picking the tray to build the case cart by logistical staff, the handling of the tray in the OR by clinical staff, and the handling of the tray by SPD staff for re-sterilization. Likewise, the time needed to perform these activities was assumed to be 120 s for peel packs. Considering $18 and $45 per hour as the labor costs for logistical and clinical staff [5], respectively, the average labor cost for handling a container would be $31.50 per hour. With this consideration, C 3 and C 4 were estimated to be $1.75 and $1.05, respectively. The PTOP mathematical model was implemented in GAMS 28.2.0 and ran using the BONMIN branch-and-bound (B-BB) algorithm. All other algorithms were coded and executed in MATLAB R2020a on a personal Intel® Core™ i7 10th GEN CPU @1.10 GHz with 16 GB RAM. The number of generations to terminate the GA for the VLD dataset was fixed at 50 generations, and for all other datasets was fixed at 500. The reason for this difference in the number of generations was that the VLD dataset was relatively simple and required fewer computational resources for convergence. However, the other datasets were more complex and required more time to reach convergence. Therefore, a higher number of generations was necessary to ensure that the GA could adequately explore the solution space. Additionally, a time limit of one hour was considered to be a stopping criterion for the BONMIN solver.

6.2. Parameter Settings

In order to tune the parameters of the developed GA, CLS, and DLS algorithms, the Taguchi method was used. Three parameters of n P o p , μ c , and μ m were pertinent to the GA, and two parameters of α and γ belonged to both CLS and DLS algorithms. In the GA, n P o p is the population size, μ c is the crossover rate, and μ m is the mutation rate. For each parameter, four possible levels were considered, as reported in Table 6. Thus, an orthogonal array of L 16 ( 4 5 ) was designed to perform 16 experiments for the calibration of each parameter. Each experiment is run ten times over the VLD dataset independently.
The Taguchi method optimizes the value of each parameter by expressing the variability of the response factor as a signal-to-noise (S/N) ratio [39]. Since the goal of tuning the parameters is concerned with obtaining both the optimized computational time and solution quality, the response factor is defined as the sum of the normalized execution time (i.e., N t i m e ) and normalized solution cost (i.e., N c o s t ) . The Taguchi design for the considered parameters is presented in Table A2 of Appendix A. The Taguchi analysis is performed on MINITAB 18. According to the main effects plot for S/N ratios depicted in Figure A1 of Appendix A, the selected values of the parameters are n P o p = 70 , μ c = 0.6 , μ m = 0.8 , α = 0.6 , and γ = 0.8 .

6.3. Computational Results

To evaluate the performance of the proposed methodology for solving the PTOP, the results of the developed algorithms were compared with that of the B-BB algorithm. Four approaches for solving the PTOP were considered. In the first approach, the classical crossover and mutation operators of the GA described in Section 5 were utilized. The second approach, called H-GA, incorporated the p-median heuristic in the initialization process by generating half of the population using the p-median heuristic. For this purpose, the p -median heuristic was successively applied to all feasible numbers of containers, and then the best n P o p 2 solutions were selected to be entered into the initial population. In the third approach, denoted as GA-CD, the GA enabled us to use two local search algorithms of CLS and DLS. Finally, in the fourth approach, denoted as H-GA-CD, the CLS and the DLS algorithms were used along with the p-median heuristic.
Solving the PTOP mathematical model, which belongs to the nonlinear mixed integer programming classes, using the B-BB algorithm reached an “out of memory” status without finding a feasible solution to the instances of 1S-7P, 2S-7P, and 5S-7P datasets. Thus, the performance of the developed H-GA-CD algorithm against the B-BB algorithm was reported only over the VLD dataset. Given that these algorithms were stochastic in nature, ten independent runs were performed, and the results are reported in Table 7. This table reports the average, standard deviation, computational time in seconds, and the best cost of solutions obtained over ten replications for each instance.
Considering the average and standard deviation, the H-GA-CD outperforms the B-BB algorithms in all instances. The H-GA-CD algorithm is shown to be highly robust as it yielded very small solution variation over multiple runs. The H-GA-CD algorithm could generate better solutions in considerably less computational time compared to the B-BB algorithm. The average computational time over five instances of this validation dataset was 1618 s and 25 s for B-BB and H-GA-CD algorithms, respectively. Given the computational resources that the commercial B-BB algorithm needs for such small-size instances and its failure to find even a feasible solution to larger instances, incorporating the designed metaheuristic algorithm was justified. Next, four solution approaches, including GA, H-GA, GA-CD and H-GA-CD, were employed to solve larger instances of the PTOP. The overall results of the five instances in each dataset are presented in Table 8. The computational time is reported in seconds.
The results in Table 8 show that incorporating the proposed CLS and DLS algorithms into the GA (i.e., GA-CD approach) resulted in solutions with, on average, 11.2%, 12.3%, and 2.6% lower costs than the traditional GA in the 1S-7P, 2S-7P, and 5S-7P datasets, respectively. The H-GA-CD approach also showed a similar trend as it yielded 11.4%, 5.3%, and 3.7% better solution quality than the GA approach in the studied datasets. The H-GA approach, however, only showed better performances in the 1S-7P and 2S-7P datasets, with a 9.7% and 5.4% improvement, respectively, compared to the GA approach. In the 5S-7P dataset, which is a large-size dataset and the traditional GA algorithm performs 5.9% better than H-GA. This is likely due to its becoming trapped at the local optima generated by the p-median heuristic, from which classic genetic operators could not assist with escaping from these local optima. It is worth noting that while the computational time required for larger instances is lengthy, once the trays are configured, they would not need frequent updates. Therefore, this is a long-term planning problem where solution quality is prioritized over computational time.
While the CLS and DLS algorithms require longer computational times, the higher-quality solutions they provide justify the increased time. In order to compare the results of these algorithms from a statistical viewpoint, non-parametric Wilcoxon rank-sum tests were examined with a significance level of 0.05. The non-parametric Wilcoxon rank-sum test was selected because it is unlikely that the distribution of the results generated by the algorithms follows the normal distribution. The p-value of the statistical test, along with the superior algorithm in terms of the solution quality, is reported in Table 9.
In order to have an insight into how these four approaches performed in each generation, the average of convergence plots over ten replications; for instance, #1 in the 2S-7P dataset, is illustrated in Figure A2 of Appendix A. This figure proves that using the CLS and the DLS algorithms, along with GA, can significantly improve the quality of the solutions. Figure A2 also demonstrates that incorporating the p-median heuristic into the initial population and utilizing the CLS and the DLS methods can make the GA consistently perform better than the other three approaches. The detailed results of the best solution obtained from the H-GA-CD method for the PTCP are presented in Table A3 of Appendix A. It shows how many of each instrument are placed in 14 configured containers, of which 13 containers are instrument trays, and one container is a peel pack. As mentioned earlier, given a solution to the PTCP, the PTAP can easily be solved. The results of the PTAP, along with the likelihood of each container to be used for a surgeon and procedure, are presented in Table 10.
As can be seen in Table 10, container 2, which is a tray, includes 23 instruments that are required for all procedures that are performed by both surgeon 1 and surgeon 2. It is expected that instruments within this tray are used during all procedures. Therefore, this tray can be opened before the procedures. There are some trays, however, that are unlikely to be used during certain procedures. For example, given the lap cholecystectomy procedure, surgeon 1 would use instruments within tray 5 only 50% of the time. Every time this tray is opened, it would cost $4.40 for re-sterilization. Since this procedure is performed 70 times a year, not opening the tray before the procedure could extrapolate to a yearly cost savings of $154.

7. Managerial Insights

The number of containers is one of the crucial factors for surgeons, clinical, and operational staff [22]. The results reported in Table 10 show that some trays (e.g., container 12 and container 13) include only a few instruments. From a practical point of view, it might not be efficient to configure such trays due to limited physical space. In this case, the number of containers to be created can be restricted by adding a new constraint t = 1 T x t R , where R is the maximum allowed number of containers. In Table 11, we conducted a sensitivity analysis on R to find how changes in the number of containers affected the expected cost of reprocessing instruments and handling containers.
Table 11 shows how a minimum of three containers was required and a configuration of 14 containers, 13 of which were trays and one of which was a peel pack, resulted in the lowest total cost. The details of this configuration are presented in Table 11 and Table A3 in Appendix A. By configuring fewer tray types the surgeons were exposed to more unrequested instruments, which resulted in higher re-sterilization costs. In contrast, configuring fewer tray types means less effort involved in handling them, resulting in a lower administrative and handling cost. Such a trade-off is depicted in Figure 10. In Table 11, the last row reports on the cost savings resulting from increasing the containers by one unit. Configuring four containers instead of three and five containers instead of four reduced the total costs by $1866 and $790, respectively.
Further cost savings can be made by providing surgeons with information about the re-sterilization costs of instruments and the potential annual savings. This will could them to better decide whether to open a tray with a low likelihood of usage before the procedure or to only open it during the procedure. Within our model, surgeons and inventory managers can estimate the potential cost savings associated with not opening trays before procedures. Figure 11 correlates the probability of containers being used and the potential cost savings if these containers are not opened prior to the procedures. The figure shows that if instruments are optimally grouped into four containers and surgeons choose not to open trays and peel packs that have less than a 50% likelihood of being used before their procedures, this results in a 3% annual cost savings in the reprocessing of instruments. Such cost savings can be increased to 10% and 16% if the instruments are configured into 5 and 10 containers, respectively. An analysis of this kind can identify surgeons, procedures, and containers that could be targeted to be a part of a pilot implementation of this study.
To examine the impact of uncertainty in a patient’s surgical instrument needs on the cost estimation of a prescribed tray configuration during a procedure, a Monte Carlo simulation was conducted. In this simulation, instance #1 from the 2S-7P dataset was selected with four containers configured, and the instrument used for each procedure was modeled using experimental probability distributions. This simulation also accounted for the fact that if the first copy of an instrument was not used, the other copies of the same instrument would not be used either. Therefore, if P i j k = 0 , then P i j k = 0 for j > j . The results of the simulation are presented in Figure 12, which shows the distribution of the total realized cost over 5000 simulation runs. The data indicates that there is a 23% chance that the actual costs exceed the estimated cost of $20,649 generated by the PTOP model, as reported in Table 11.
.

8. Conclusions and Discussion

This study dealt with the problem of configuring surgical trays based on the probabilities of instrument usage. This problem is called the probabilistic tray optimization problem (PTOP). The PTOP is formulated as an integer non-linear programming model. Since this problem is in the class of NP-hard problems, it is decomposed into two sub-problems of the probabilistic tray configuration problem (PTCP) and the probabilistic tray assignment problem (PTAP). Given a solution to the PTCP, the PTAP becomes a trivial problem of an integer linear programming model, which can be easily solved. A heuristic and a metaheuristic were developed to solve the PTCP. The heuristic algorithm provides an approximate solution to the PTCP in a short computational time by formulating the PTCP as a well-known p -median model. However, this method may result in a low-quality solution. Therefore, a genetic algorithm (GA) was also designed to improve the quality of solutions obtained by the p -median heuristic. Moreover, in order to enable the GA to better explore the solution space, two local search algorithms—combining local search (CLS) and decomposing local search (DLS)—were invented and embedded in the GA. The CLS and DLS algorithms attempted to decrease and increase the number of containers, respectively, in a given solution. The applicability of the model is demonstrated through the historical data available in the literature.
Close collaboration with surgeons to improve their conventional tray configurations is key for a successful implementation. The PTOP model can facilitate this collaboration, as the potential cost savings derived from the reconfiguration of trays and not opening trays with a low probability of usage become readily apparent through executing the PTOP model. Given the required data, the PTOP model can be applied in any hospital to reveal how much they can benefit from the optimization of the instrument trays.
Our approach provides the minimum number of trays and instruments that need to be stocked for performing procedures that are not scheduled for the same day. However, since, in practice, the re-sterilization process usually takes one day to be completed, the trays undergoing re-sterilization would not be available for use until the next day [17,19,22,25,40]. Therefore, multiple copies of a tray are needed when multiple procedures requiring the same tray type are scheduled to be performed in a day. Such an approach provides a maximum number of trays and instruments that need to be stocked for performing procedures. However, the optimal number of instruments and trays highly depends on the OR schedule. For example, if two procedures are scheduled on the same day only a few times a year, keeping two copies of a tray would not be optimal. Instead, by modifying the schedule of those procedures not to take place concurrently, the stocked instruments could be reduced.
It should be noted that our model does not eliminate any instruments that may be used rarely, as we understand their importance in specific procedures. Our model only rearranges all instruments that surgeons request into different tray types by separating frequently used instruments from rarely used ones. This model ensures that the surgeon’s preferred instruments are readily available in the operating room. By optimally configuring the instrument trays, the time it takes for surgeons to find the instruments they need during surgery can also be reduced. In this research, the cost savings resulting from our model are associated only with avoiding unnecessary labor-intensive re-sterilization processes rather than potential savings on capital and holding costs. To estimate capital savings, a standardization project, such as the work of Koyle et al. [8], which resulted in a reduction in surgical instruments, needs to be conducted. It is important to note that removing a rarely used instrument tray from stock may compromise patient safety, as the tray may be needed in rare occurrences. Therefore, in such a standardization project, safety stock for these trays should be considered to prevent the risk of tray non-availability.
Since many hospitals currently do not record the use of instruments, one of the considerations toward implementing the findings of this model is that direct observations are required to collect accurate instrument usage data. These observations could be performed after a procedure in the sterile processing department or in the operating rooms. Some surgical tool companies offer RFID-tagged instruments, which facilitate data collection. Since collecting accurate usage data, particularly for instruments without RFID tags, can be difficult, an alternative would be to have surgeons label instruments as “always needed”, “sometimes needed”, or “rarely needed” and then estimate how likely each label is to be used. Additionally, a time and motion study could help with a more accurate estimation of the labor costs associated with the sterilization process and handling of trays. One may benefit from the works of John-Baptiste et al. [41], and Mhlaba et al. [7] in estimating these labor costs, but analysts should be conscious that the sterilization process in hospitals may be different from each other. Another limitation of this study is that it assumes the probabilities of instrument usage to remain constant over time. However, in reality, these probabilities may change due to variations in surgical techniques or patient demographics. As a result, it may be necessary to periodically update the PTOP model to account for such changes.
Future studies could explore how instrument usages are interdependent and develop a model that captures these relationships. Another interesting future research direction is to incorporate patient demographic variables such as gender, race, age and patient clinical variables such as Body Mass Index (BMI), past surgical history, and past medical history to predict the necessity of an instrument for a given case. Furthermore, studying other heuristics and metaheuristics and evaluating their performance against our solution algorithms would provide another interesting future study.

Author Contributions

Conceptualization, E.A., D.T.M. and S.H.; Methodology, E.A. and D.T.M.; Software, E.A.; Resources, S.H.; Writing—original draft, E.A.; Writing—review & editing, E.A. and D.T.M.; Visualization, E.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The datasets used in this research are available at https://doi.org/10.17632/6vr7fp4npy.1 (accessed on 13 January 2022).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Algorithm A1: p-median based heuristic.
Inputs: A set of vertexes V . A set of edge E .
The distance matrix corresponding to the d i j m n
1: Set P = i = 1 I max k { A i k }
2: While there exist feasible solutions Do
3:       Solve the p -median model
4:      Solve the PTAP
5:        P P 1
6: End While
Output: The solution corresponding to the P that resulted in
the minimum objective function value in the PTAP
Algorithm A2: Combining Local Search (CLS).
1: Extract the initial number of containers, I n i N u m C o n , for a given solution S
2: N u m C o n I n i N u m C o n
3: While N u m C o n > I n i N u m C o n   ( 1 γ )
4:     For  t = 1  To  N u m C o n  Do
5:       Calculate O t using Equations (12) and (13)
6:     End For
7:     For  w = 1  To 3 Do
8:     If R a n d ( ) α    Then
9:      Perform random walk on S and obtain local optimum S w
10:     Else If w = 1
11:       Walk 1: combine the two containers with the lowest O t and obtain local optimum S 1
12:     Else If w = 2
13:       Walk 2: combine the two containers with the highest O t and obtain local optimum S 2
14:     Else If w = 3
15:       Walk 3: combine the containers with the lowest and the highest O t and obtain local
          optimum S 3
16:     End If
17:     If t h e   w e i g h t   o f   t h e   c o m b i n e d   c o n t a i n e r s > L  Then
18:         Perform the repairing mechanism
19:           End If
20:       End For
21:      Add the best solution among S 1 , S 2 , and S 3 to the current population
22:       N u m C o n N u m C o n 1
23: End While
Algorithm A3: Decomposing Local Search (DLS).
1: Extract the initial number of containers I n i N u m C o n for a given solution S
2: N u m C o n = I n i N u m C o n
3: While N u m C o n < I n i N u m C o n   ( 1 + γ )
4:      For  t = 1  To  N u m C o n except peel packs Do
5:          Calculate O t using Equation (12)
6:      End For
7:      Select the container t with the lowest O t
8:      For all instruments in tray t
9:        Calculate O i using Equation (14)
10:      End For
11:      For  w = 1  To 2 Do
12:     If R a n d ( ) α    Then
13:            Random walk: Generate a new solution S w by randomly selecting an instrument and putting it in a peel pack
14:     Else If w = 1
15:      Walk 1: Generate a new solution S 1 by putting the instrument with the lowest O i
in a new container as a peel pack
16:     Else If w = 2
17:      Walk 2: Generate a new solution S 2 by putting the instrument with the highest O i
in a new container as a peel pack
18:     End If
19:      End For
20:     Select the container t with the highest O t
21:     For all instruments in tray t
22:        Calculate O i using Equation (14)
23:     End For
24:     For  w = 3  To 4 Do
25:      If R a n d ( ) α    Then
26:            Random walk: Generate a new solution S w by randomly selecting an instrument and putting it in a peel pack
27:     Else If w = 1
28:      Walk 1: Generate a new solution S 3 by putting the instrument with the lowest O i
in a new container as a peel pack
29:      Else If w = 2
30:      Walk 2: Generate a new solution S 4 by putting the instrument with the highest O i
in a new container as a peel pack
31:     End If
32:      End For
33:      Add the best solution among S 1 , S 2 , S 3 , and S 4 to the current population
34:      N u m C o n N u m C o n + 1
35: End While
Table A1. The values of d i j m n corresponding to the PTOP example presented in Table 1.
Table A1. The values of d i j m n corresponding to the PTOP example presented in Table 1.
Instrument
12345678910111213
Instrument13.006.299.557.676.007.666.676.247.936.506.3610.696.30
26.290.506.273.391.025.292.231.244.531.501.368.581.82
39.556.272.765.995.538.226.365.637.435.555.548.896.09
47.673.395.991.282.576.593.752.805.002.712.678.373.31
56.001.025.532.570.014.691.400.263.990.520.388.010.94
67.665.298.226.594.692.345.404.807.515.185.049.615.19
76.672.236.363.751.405.400.691.564.811.881.748.692.25
86.241.245.632.800.264.801.560.124.220.740.608.051.16
97.934.537.435.003.997.514.814.221.994.104.079.864.77
106.501.505.552.710.525.181.880.744.100.250.778.031.42
116.361.365.542.670.385.041.740.604.070.770.188.021.28
1210.698.588.898.378.019.618.698.059.868.038.024.008.25
136.301.826.093.310.945.192.251.164.771.421.288.250.46
Table A2. Taguchi analysis for tuning parameters of the GA, CLS, and DLS algorithms.
Table A2. Taguchi analysis for tuning parameters of the GA, CLS, and DLS algorithms.
ExperimentParameter’s LevelCost ($)Time (s)NcostNtimeResponse
n P o p μ c μ m α γ
11111140.2412.861.000.001.00
21222240.1921.600.800.321.12
31333340.0416.030.250.120.36
41444440.0118.190.150.190.35
52123440.1220.210.540.270.80
62214340.0425.840.250.470.73
72341240.0917.040.440.150.60
82432140.0916.240.440.120.56
93134240.0624.110.350.410.76
103243140.0028.370.090.570.65
113312439.9921.550.050.320.37
123421339.9725.000.000.440.44
134142339.9940.250.051.001.05
144231440.1327.200.580.521.10
154324140.0626.970.350.520.86
164413240.0527.700.290.540.84
Figure A1. Main effects plot for S/N ratios.
Figure A1. Main effects plot for S/N ratios.
Mathematics 11 02219 g0a1
Figure A2. Convergence plot of the four approaches to solve the PTOP (instance #1 in the 2S-7P dataset).
Figure A2. Convergence plot of the four approaches to solve the PTOP (instance #1 in the 2S-7P dataset).
Mathematics 11 02219 g0a2
Table A3. A solution to the PTCP.
Table A3. A solution to the PTCP.
InstrumentsContainers
1234567891011121314
1 1
21
31
4 1
51
6 1
71
8 1
9 1
10 1
111
121
131
141
15 11
16 11
171 1
18 1 1
19 1 1
2011
211
22 1 1
23 1 1
241 1
25 1 1
26 1 1
271 1
281 1
291 1
30 11
311 1
32 1 1
331 1
34 1 1
35 1 1
36 2
371 1
38 1 1
39 2
40 1 1
41 1 1
421 1
43 1 1
44 1 1
45 1 1
46 1 1
47 11
48 1 1
49 11
50 11
51 11
52 1 1
531 1
542
55 1 1
56 2
571 1
581 1
59 2
60 1 1
611 1
62 11
63 1 1
641 1
65 2
66 1
67 1 1
681 1
69 1 1
70 1 1
71 1 1
72 11
73 11
74 1 1
75 11
76 1 1

References

  1. CMS. National Health Expenditure Projections 2019–2028, Office. 2020, 1–17. Available online: https://www.cms.gov/Research-Statistics-Data-and-Systems/Statistics-Trends-and-Reports/NationalHealthExpendData/NationalHealthAccountsProjected (accessed on 13 January 2022).
  2. Ahmadi, E.; Masel, D.T.; Schwerha, D.; Hostetler, S. A bi-objective optimization approach for configuring surgical trays with ergonomic risk consideration. IISE Trans. Healthc. Syst. Eng. 2019, 9, 327–341. [Google Scholar] [CrossRef]
  3. Schoenfelder, J.; Kohl, S.; Glaser, M.; McRae, S.; Brunner, J.O.; Koperna, T. Simulation-based evaluation of operating room management policies. BMC Health Serv. Res. 2021, 21, 271. [Google Scholar] [CrossRef]
  4. Görgülü, B.; Sarhangian, V. Newsvendor Approach to Design of Surgical Preference Cards. Serv. Sci. 2022, 14, 213–230. [Google Scholar] [CrossRef]
  5. Ahmadi, E.; Masel, D.T.; Hostetler, S. A robust stochastic decision-making model for inventory allocation of surgical supplies to reduce logistics costs in hospitals: A case study. Oper. Oper. Res. Health Care 2019, 20, 33–44. [Google Scholar] [CrossRef]
  6. Stockert, E.W.; Langerman, A. Assessing the magnitude and costs of intraoperative inefficiencies attributable to surgical instrument trays. J. Am. Coll. Surg. 2014, 219, 646–655. [Google Scholar] [CrossRef]
  7. Mhlaba, J.M.; Stockert, E.W.; Coronel, M.; Langerman, A.J. Surgical instrumentation: The true cost of instrument trays and a potential strategy for optimization. J. Hosp. Adm. 2015, 4, 82. [Google Scholar] [CrossRef]
  8. Koyle, M.A.; AlQarni, N.; Odeh, R.; Butt, H.; Alkahtani, M.M.; Konstant, L.; Pendergast, L.; Koyle, L.C.; Baker, G.R. Reduction and standardization of surgical instruments in pediatric inguinal hernia repair. J. Pediatr. Urol. 2018, 14, 20–24. [Google Scholar] [CrossRef]
  9. Farrokhi, F.R.; Gunther, M.; Williams, B.; Blackmore, C.C. Application of lean methodology for improved quality and efficiency in operating room instrument availability. J. Healthc. Qual. 2015, 37, 277–286. [Google Scholar] [CrossRef] [PubMed]
  10. Harvey, L.; Slocum, P.; Heft, J.; Van Meter, M.; Lovett, B.; Adam, R. Gynecologic Surgery Instrument Trays: Leveraging Surgeon Knowledge to Improve Supply Chain Efficiency. J. Gynecol. Surg. 2017, 33, 180–183. [Google Scholar] [CrossRef]
  11. Nast, K.; Swords, K. Decreasing operating room costs via reduction of surgical instruments. J. Pediatr. Urol. 2019, 15, 153.e1–153.e6. [Google Scholar] [CrossRef]
  12. Choobineh, A.; Movahed, M.; Tabatabaie, S.H.; Kumashiro, M. Perceived Demands and Musculoskeletal Disorders in Operating Room Nurses of Shiraz City Hospitals. Ind. Health 2010, 48, 74–84. [Google Scholar] [CrossRef] [PubMed]
  13. Seavey, R. Packaging for Sterilization: It’s not just a box, a bag or gift wrapping! Manag. Infect. Control. 2008, 8, 82–97. [Google Scholar]
  14. Sheikhzadeh, A.; Gore, C.; Zuckerman, J.D.; Nordin, M. Perioperating nurses and technicians’ perceptions of ergonomic risk factors in the surgical environment. Appl. Ergon. 2009, 40, 833–839. [Google Scholar] [CrossRef] [PubMed]
  15. Choi, S.D. A Review of the Ergonomic Issues in the Laparoscopic Operating Room. J. Healthc. Eng. 2012, 3, 587–604. [Google Scholar] [CrossRef]
  16. Morris, L.F.; Arenas, M.A.R.; Cerny, J.; Berger, J.S.; Borror, C.M.; Ong, M.; Cayo, A.K.; Graham, P.H.; Grubbs, E.G.; Lee, J.E.; et al. Streamlining variability in hospital charges for standard thyroidectomy: Developing a strategy to decrease waste. Surgery 2014, 156, 1441–1449. [Google Scholar] [CrossRef] [PubMed]
  17. Dollevoet, T.; van Essen, J.T.; Glorie, K.M. Solution methods for the tray optimization problem. Eur. J. Oper. Res. 2018, 271, 1070–1084. [Google Scholar] [CrossRef]
  18. Ahmadi, E.; Masel, D.T.; Metcalf, A.Y.; Schuller, K. Inventory management of surgical supplies and instruments in hospitals: A literature review. Health Syst. 2018, 8, 134–151. [Google Scholar] [CrossRef]
  19. Van de Klundert, J.; Muls, P.; Schadd, M. Optimizing sterilization logistics in hospitals. Health Care Manag. Sci. 2008, 11, 23–33. [Google Scholar] [CrossRef]
  20. Florijn, E.P. Optimisation of the Distribution of Surgical Instruments over Trays: Cost Effectiveness and Quality Improvement of an Operating Theatre. Master’s Thesis, University of Twente, Enschede, The Netherlands, 2008. [Google Scholar]
  21. Reymondon, F.; Pellet, B.; Marcon, E. Optimization of hospital sterilization costs proposing new grouping choices of medical devices into packages. Int. J. Prod. Econ. 2008, 112, 326–335. [Google Scholar] [CrossRef]
  22. Dobson, G.; Seidmann, A.; Tilson, V.; Froix, A. Configuring surgical instrument trays to reduce costs. IIE Trans. Healthc. Syst. Eng. 2015, 5, 225–237. [Google Scholar] [CrossRef]
  23. Dos Santos, B.M.; Fogliatto, F.S.; Zani, C.M.; Peres, F.A.P. Approaches to the rationalization of surgical instrument trays: Scoping review and research agenda. BMC Health Serv. Res. 2021, 21, 163. [Google Scholar] [CrossRef]
  24. Harris, S.; Claudio, D. A Goal Programming Approach to the Tray Optimization Problem. In IIE Annual Conference. Proceedings; Institute of Industrial and Systems Engineers (IISE): Peachtree Corners, GA, USA, 2021; pp. 680–685. [Google Scholar]
  25. Deshpande, V.; Mundru, N.; Rath, S.; Knowles, M.; Rowe, D.; Wood, B.C. Data-Driven Surgical Tray Optimization to Improve Operating Room Efficiency. Oper. Res. 2023. [Google Scholar] [CrossRef]
  26. AAMI. Comprehensive Guide to Steam Sterilzation and Sterility Assurance in Health Care Facilities; AAMI: Melbourne, VIC, Australia, 2006. [Google Scholar]
  27. Karp, R.M. Reducibility among combinatorial problems. In Complexity of Computer Computations; Springer: Berlin/Heidelberg, Germany, 1972; pp. 85–103. [Google Scholar]
  28. Berman, O.; Wang, J. The network p-median problem with discrete probabilistic demand weights. Comput. Oper. Res. 2010, 37, 1455–1463. [Google Scholar] [CrossRef]
  29. Mladenović, N.; Brimberg, J.; Hansen, P.; Moreno-Pérez, J.A. The p-median problem: A survey of metaheuristic approaches. Eur. J. Oper. Res. 2007, 179, 927–939. [Google Scholar] [CrossRef]
  30. Kusiak, A. The generalized group technology concept. Int. J. Prod. Res. 1987, 25, 561–569. [Google Scholar] [CrossRef]
  31. Wang, J.; Roze, C. Formation of Machine Cells and Part Families in Cellular Manufacturing: A Linear Integer Programming Approach. In Proceedings of the 1994 IEEE International Conference on Industrial Technology—ICIT ’94, Guangzhou, China, 5–9 December 1994; pp. 350–354. [Google Scholar]
  32. Hassanzadeh, R.; Mahdavi, I.; Mahdavi-Amiri, N.; Tajdin, A. A genetic algorithm for solving fuzzy shortest path problems with mixed fuzzy arc lengths. Math. Comput. Model. 2013, 57, 84–99. [Google Scholar] [CrossRef]
  33. Ahmadi, E.; Süer, G.A.; Al-Ogaili, F. Solving Stochastic Shortest Distance Path Problem by Using Genetic Algorithms. Procedia Comput. Sci. 2018, 140, 79–86. [Google Scholar] [CrossRef]
  34. Ahmadi, E.; Masel, D.T.; Hostetler, S.; Maihami, R.; Ghalehkhondabi, I. A centralized stochastic inventory control model for perishable products considering age-dependent purchase price and lead time. TOP 2020, 28, 231–269. [Google Scholar] [CrossRef]
  35. Vahdani, B.; Soltani, M.; Yazdani, M.; Mousavi, S.M. A three level joint location-inventory problem with correlated demand, shortages and periodic review system: Robust meta-heuristics. Comput. Ind. Eng. 2017, 109, 113–129. [Google Scholar] [CrossRef]
  36. Ahmadi, E.; Mosadegh, H.; Maihami, R.; Ghalehkhondabi, I.; Sun, M.; Süer, G.A. Intelligent inventory management approaches for perishable pharmaceutical products in a healthcare supply chain. Comput. Oper. Res. 2022, 147, 105968. [Google Scholar] [CrossRef]
  37. Talbi, E.-G. Metaheuristics: From Design to Implementation; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  38. Asefi, H.; Jolai, F.; Rabiee, M.; Araghi, M.E.T. A hybrid NSGA-II and VNS for solving a bi-objective no-wait flexible flowshop scheduling problem. Int. J. Adv. Manuf. Technol. 2014, 75, 1017–1033. [Google Scholar] [CrossRef]
  39. Aghaaminiha, M.; Ghanadian, S.A.; Ahmadi, E.; Farnoud, A.M. A machine learning approach to estimation of phase diagrams for three- component lipid mixtures. Biochim. Biophys. Acta BBA-Biomembr. 2020, 1862, 183350. [Google Scholar] [CrossRef] [PubMed]
  40. Diamant, A.; Milner, J.; Quereshy, F.; Xu, B. Inventory management of reusable surgical supplies. Health Care Manag. Sci. 2017, 21, 439–459. [Google Scholar] [CrossRef] [PubMed]
  41. John-Baptiste, A.; Sowerby, L.; Chin, C.; Martin, J.; Rotenberg, B. Comparing surgical trays with redundant instruments with trays with reduced instruments: A cost analysis. CMAJ Open 2016, 4, E404–E408. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Illustration of used and unused instruments for the Plastic tray. The used instruments are listed on the towel, and the unused instruments are collected in the tray.
Figure 1. Illustration of used and unused instruments for the Plastic tray. The used instruments are listed on the towel, and the unused instruments are collected in the tray.
Mathematics 11 02219 g001
Figure 2. The solution of PTCP obtained from p-median heuristic.
Figure 2. The solution of PTCP obtained from p-median heuristic.
Mathematics 11 02219 g002
Figure 3. Solution encoding for a sample PTCP.
Figure 3. Solution encoding for a sample PTCP.
Mathematics 11 02219 g003
Figure 4. Two cut-point crossover operator.
Figure 4. Two cut-point crossover operator.
Mathematics 11 02219 g004
Figure 5. Uniform crossover operator.
Figure 5. Uniform crossover operator.
Mathematics 11 02219 g005
Figure 6. Swap mutation operator.
Figure 6. Swap mutation operator.
Mathematics 11 02219 g006
Figure 7. Inversion mutation operator.
Figure 7. Inversion mutation operator.
Mathematics 11 02219 g007
Figure 8. Shifts mutation operator.
Figure 8. Shifts mutation operator.
Mathematics 11 02219 g008
Figure 9. Shuffle mutation operator.
Figure 9. Shuffle mutation operator.
Mathematics 11 02219 g009
Figure 10. The impact of number of containers on the cost components.
Figure 10. The impact of number of containers on the cost components.
Mathematics 11 02219 g010
Figure 11. Percentage of cost savings in the re-sterilization costs (i.e., f 1 + f 2 ) due to not opening the trays before procedures.
Figure 11. Percentage of cost savings in the re-sterilization costs (i.e., f 1 + f 2 ) due to not opening the trays before procedures.
Mathematics 11 02219 g011
Figure 12. Simulated solution for instance #1 in the 2S-7P dataset with R = 4 .
Figure 12. Simulated solution for instance #1 in the 2S-7P dataset with R = 4 .
Mathematics 11 02219 g012
Table 1. An example of surgeons’ preference cards.
Table 1. An example of surgeons’ preference cards.
Surgeon-ProceduresInstrument
12345
102031
201301
323212
411102
522102
620210
Table 2. The binary representation of the surgeons’ preference cards.
Table 2. The binary representation of the surgeons’ preference cards.
Surgeon-ProceduresInstrument
12345
100110000111 10
200100 111000 10
311111 110100 11
410100 100000 11
511110 100000 11
611000 11010000
Table 3. The probability distribution of instruments’ usage.
Table 3. The probability distribution of instruments’ usage.
Proc.Instrument
12345
10.000.000.950.700.00 0.000.000.00 0.760.250.18 0.950.00
20.000.000.560.000.00 0.480.260.12 0.000.000.00 0.800.00
30.800.150.700.530.01 0.600.180.00 0.450.000.00 0.750.14
40.900.000.400.000.00 0.680.000.00 0.000.000.00 0.850.12
50.450.150.150.050.00 0.190.000.00 0.000.000.00 0.650.20
60.850.200.000.000.00 0.390.250.00 0.780.000.00 0.000.00
Table 4. An example of calculating contribution of each container to the objective function.
Table 4. An example of calculating contribution of each container to the objective function.
Container TypeIncluded Instruments in the Container *Container TypeContainer’s Contribution ( O t )
1 3 1 , 4 2 Tray155.4
2 2 3 , 3 3 Tray7.8
3 2 2 , 5 2 Tray99.4
4 1 1 , 1 2 ,   4 1 ,   5 1 Tray667.1
6 3 2 Peel pack13.8
7 4 3 Peel pack3.6
10 2 1 Peel pack55.2
* i j indicates the jth copy of instrument i
Table 5. Specifications of the datasets.
Table 5. Specifications of the datasets.
Dataset NameNumber of InstancesNumber of SurgeonsNumber of ProceduresNumber of Unique InstrumentsTotal Number of Instruments
VLD523513
1S-7P51776136
2S-7P52776136
5S-7P55776250
Table 6. The values of the parameters corresponding to each level.
Table 6. The values of the parameters corresponding to each level.
Level n P o p μ c μ m α γ
1300.400.500.200.20
2500.500.600.400.40
3700.600.700.600.60
4900.700.800.800.80
Table 7. Validation of the proposed solution algorithms.
Table 7. Validation of the proposed solution algorithms.
DatasetInstance #B-BBH-GA-CD
AverageS.D.BestTimeAverageS.D.BestTime
VLD146.74.344.72064.840.00.239.924.3
258.46.846.6915.1037.40.137.323.0
343.92.442.82095.931.10.430.926.7
440.48.934.01773.635.00.534.724.8
543.45.237.41240.032.90.032.928.4
Table 8. Computational results of the solution approaches.
Table 8. Computational results of the solution approaches.
DatasetInstance #GAH-GAGA-CDH-GA-CD
Ave.S.D.BestTimeAve.S.D.BestTimeAve.S.D.BestTimeAve.S.D.BestTime
1S-7P110,0574669560189787393387162093874419585201772875718985382147
298893649341168589732238765223687811398573176387207885821621
3971735791311712881324585092058877617284891721883523286291682
4834133078551723751428472292051744615272601994731615171682134
5837025579661711784948675002443742512272212206744810273101981
2S-7P119,95144919,505507218,98412518,841525618,04625217,769787118,86411418,7867649
220,19034519,861509519,42025619,105494517,30936616,875896619,28527219,0488218
319,92138119,380491818,35722118,093486517,50230917,225794018,43914318,2237956
420,47539119,869505719,87860019,348481617,89333617,506803119,50214719,3327824
520,19837619,620492718,61714518,485459217,61924517,443775319,35114219,1887302
5S-7P184,102106882,81079,02490,673416586,75681,11182,950416178,972165,62083,409426977,557166,985
285,551124384,53678,90190,216179688,39884,09382,792232680,335172,94082,432215379,260165,187
385,042189582,72878,86989,550174487,57082,35982,269369878,203169,84480,998280076,270167,017
485,270184883,24977,98591,259409386,61880,32182,731286878,724171,44479,654159678,274165,616
585,411105284,43277,89388,867343485,06481,97483,674260779,609168,05883,248252780,420163,162
Table 9. Pairwise statistical results of the Wilcoxon rank-sum test for the average cost.
Table 9. Pairwise statistical results of the Wilcoxon rank-sum test for the average cost.
DatasetGAH-GAGA-CDH-GA-CD
GA----
H-GA0.585
(no significant difference)
---
GA-CD0.013
(GA-CD performs better)
0.021
(GA-CD performs better)
--
H-GA-CD0.018
(H-GA-CD performs better)
0.252
(no significant difference)
0.358
(no significant difference)
-
Table 10. A solution to the PTAP and the probability of containers to be used during the procedure.
Table 10. A solution to the PTAP and the probability of containers to be used during the procedure.
SurgeonProcedureContainers
1234567891011121314
Surgeon 1Mediport insertion100%100%53%19%99%9%77%25% 51% 82%
Excision of small lesion100%100% 76%
Lap appendectomy100%100%100%32%95%92%99%2% 72%
Lap cholecystectomy100%100%67%99%50%
Lap ventral hernia repair100%100%99%49%67%60%97%75%12%39% 11%
Open hernia repair100%100%100% 99%100%100%
Bowel resection100%100%100%86%100%96%100%48%35% 56%20%55%4%
Surgeon 2Mediport insertion100%100%41%26%69%6%98% 39%
Excision of small lesion 100%100%
Lap appendectomy97%100%100%16% 11%84%
Lap cholecystectomy 100% 100% 2%
Lap ventral hernia repair100%100%59%33%33%95%75% 55%
Open hernia repair100%100%99% 100% 68%99%
Bowel resection100%100%100%37%98%100%98% 21%74%98%9%
Number of instruments2723131111998765331
Weight (lb.)14.1711.568.065.676.63.165.433.623.72.521.41.60.890.56
Resterilization cost per procedure
if the container is opened
$10.80 $9.20 $5.20 $4.40 $4.40 $3.60 $3.60 $3.20 $2.80 $2.40 $2.00 $1.20 $1.20 $0.80
Table 11. Results of the sensitivity analysis on the number of containers for instance #1 in the 2S-7P dataset.
Table 11. Results of the sensitivity analysis on the number of containers for instance #1 in the 2S-7P dataset.
Value of R 3456789101112131415
Number of trays 34567889911121313
Number of peel packs0000001121111
f 1 $20,383$18,393$17,009$15,997$15,326$14,836$14,594$14,027$13,842$13,712$13,292$13,120$13,120
f 2 $0$0$0$0$0$0$1$16$60$1$3$1$1
f 3 $2132$2256$2849$3509$3948$3976$3982$4006$4067$4179$4552$4641$4641
f 4 $0$0$0$0$0$0$18$118$79$17$7$7$7
Total Cost$22,514$20,649$19,858$19,506$19,274$18,812$18,595$18,167$18,048$17,909$17,853$17,769$17,769
Cost Savings/ $1866$790$353$232$462$218$428$118$139$56$84$0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ahmadi, E.; Masel, D.T.; Hostetler, S. A Data-Driven Decision-Making Model for Configuring Surgical Trays Based on the Likelihood of Instrument Usages. Mathematics 2023, 11, 2219. https://doi.org/10.3390/math11092219

AMA Style

Ahmadi E, Masel DT, Hostetler S. A Data-Driven Decision-Making Model for Configuring Surgical Trays Based on the Likelihood of Instrument Usages. Mathematics. 2023; 11(9):2219. https://doi.org/10.3390/math11092219

Chicago/Turabian Style

Ahmadi, Ehsan, Dale T. Masel, and Seth Hostetler. 2023. "A Data-Driven Decision-Making Model for Configuring Surgical Trays Based on the Likelihood of Instrument Usages" Mathematics 11, no. 9: 2219. https://doi.org/10.3390/math11092219

APA Style

Ahmadi, E., Masel, D. T., & Hostetler, S. (2023). A Data-Driven Decision-Making Model for Configuring Surgical Trays Based on the Likelihood of Instrument Usages. Mathematics, 11(9), 2219. https://doi.org/10.3390/math11092219

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop