Healthcare Management Engineering Insights
Healthcare Management Engineering Insights
Healthcare
Management
Engineering
In Action
Applying Fundamental Management
Principles for Operational Decision
Making in Healthcare
Second Edition
Business Guides on the Go
“Business Guides on the Go” presents cutting-edge insights from practice
on particular topics within the fields of business, management, and
finance. Written by practitioners and experts in a concise and accessible
form the series provides professionals with a general understanding and a
first practical approach to latest developments in business strategy, leader-
ship, operations, HR management, innovation and technology manage-
ment, marketing or digitalization. Students of business administration or
management will also benefit from these practical guides for their future
occupation/careers.
These Guides suit the needs of today’s fast reader.
Alexander Kolker
Healthcare
Management
Engineering
In Action
Applying Fundamental Management
Principles for Operational Decision
Making in Healthcare
Alexander Kolker
MILWAUKEE, WI, USA
© The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer Nature Switzerland
AG 2012, 2024
This work is subject to copyright. All rights are solely and exclusively licensed by the Publisher, whether
the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of
illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and trans-
mission or information storage and retrieval, electronic adaptation, computer software, or by similar or
dissimilar methodology now known or hereafter developed.
The use of general descriptive names, registered names, trademarks, service marks, etc. in this publication
does not imply, even in the absence of a specific statement, that such names are exempt from the relevant
protective laws and regulations and therefore free for general use.
The publisher, the authors, and the editors are safe to assume that the advice and information in this book
are believed to be true and accurate at the date of publication. Neither the publisher nor the authors or
the editors give a warranty, expressed or implied, with respect to the material contained herein or for any
errors or omissions that may have been made. The publisher remains neutral with regard to jurisdictional
claims in published maps and institutional affiliations.
This Springer imprint is published by the registered company Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Paper in this product is recyclable.
Distinct Book Features
v
Preface to the Updated and Expanded 2nd
Edition
vii
viii Preface to the Updated and Expanded 2nd Edition
• Capacity: How many beds are required for a department or unit? How
many procedure rooms, operating rooms, or pieces of equipment are
needed for different services?
• Staffing: How many nurses, physicians, and other providers are needed
for a particular shift in a unit (department) in order to best achieve
operational and service performance objectives?
• Scheduling: What are the optimized staff schedules that help not only
in delivering safe and efficient care for patients but also take into
account staff preferences and convenience?
• Patient flow: What patient wait time is acceptable (if any at all) at the
service stations to achieve the system throughput goals?
• Resource allocation: Is it more efficient to use specialized resources or
pooled (interchangeable) resources (operating/procedure rooms, beds,
equipment, and staff)?
x Preface to the Updated and Expanded 2nd Edition
• Does it make economic sense to keep some patient service lines (or
drop them at all)?
• Forecasting: How to forecast the future patient volumes (demand) or
patient unit census forshort- and medium-term planning purposes?
• Optimized geographic location of facilities and facilities layout.
• Design of the facility-optimized workflow.
• Defining and comparing the staff or unit productivity with multiple
inputs and outputs.
• Optimizing a supply chain distribution network and inventory
management.
• Performing advanced multivariate statistical data analysis for market-
ing and budget planning.
This list can easily be extended to include any other area of operational
management that requires quantitative analysis to justify
decision-making.
The ultimate goal and the holy grail of management engineering meth-
odology is to provide aid and guidance to efficiently manage hospital
operations, i.e., reducing the costs of using resources for delivery of care
while keeping high safety and outcomes standards for patients.
It can be said that the entire hospital is a patient; management engi-
neering methodology plays the role of a medical field with different spe-
cialties developed to address different conditions and operational
problems, and the management engineer serves as a doctor who diagno-
ses the operational disease and develops a treatment plan for an ailing
hospital and its operations.
No concept or methodology can truly be convincing without multiple
concrete and practically relevant examples of its application. Kopach-
Konrad et al. (2007) state, “…we believe that widespread success will
only come when a critical mass of healthcare organizations recognize its
[healthcare engineering] value through concrete examples. Only then
will these organizations promote changes needed for its adoption.” The
author of this book shares this belief.
In this book, traditional managerial decision-making and manage-
ment engineering methodology are applied side-by-side to the same
problems to illustrate the deficiencies of the former and the power of the
Preface to the Updated and Expanded 2nd Edition xi
the right staff in the right shift, cross-trained staffing, and the optimal
staffing of an Eye Institute.
Chapter 4. Staffing Problems with Random Patient Demand: Solutions
Using Analytic Optimization Techniques.
This chapter includes examples of optimization problems using ana-
lytic techniques rather than numeric optimization. The problem of staff-
ing with minimal cost of under- and overstaffing with random patient
demand is presented in Sect. 4.1 using the so-called “newsvendor” frame-
work. A somewhat related problem of PACU (Post Anesthesia Care Unit)
optimal staffing is presented in Sect. 4.2 for highly variable day-to-day
and hour-to-hour patient census using multi-criteria optimization.
Section 4.3 describes an approach for minimizing the number of tests per
specimen at relatively low disease prevalence by pooling several speci-
mens into one batch.
Chapter 5. Resource Allocation and Service Line Problems: Solutions
Using Linear Optimization.
Linear optimization is a powerful methodology for finding the mini-
mal cost staffing and resource allocation solution when the objective
function and constraints can be presented as linear combinations of the
decision variables. There are plenty of problems amenable to this meth-
odology in healthcare settings. Excel solver is used for setting up and
running all examples presented in this chapter. The role of the sensitivity
report produced by the Excel solver is highlighted for analyzing the opti-
mal solution. The use of the sensitivity report is explained and demon-
strated in detail in multiple examples. Five problems are presented
including deciding to keep or drop a service line, minimal cost mix of
various skilled staff, optimizing the diagnostic lab tests offering, and the
minimal cost supplier’s distribution network with various numbers of
manufacturing facilities, distribution centers, and receiving hospitals.
Chapter 6. Staffing and Scheduling Problems with Sliding Days Off:
Solutions Using Integer Linear Optimization.
Integer linear optimization is a subset of general linear optimization.
Five problems are presented including staffing and scheduling for full-
and part-time employees with sliding mandatory days off, personal
employees’ preference, three-shift operations, and physician residents’
Preface to the Updated and Expanded 2nd Edition xiii
restricted work hours. The Excel solver is used for setting up and solving
these problems.
Chapter 7. Comparative Performance of Units and Organizations with
Multiple Inputs and Outputs: Introduction to Data Envelopment
Analysis.
The preferred methodology for evaluating the relative performance of
units and organizations with multiple inputs and outputs is data envelop-
ment analysis (DEA). The setting up of a DEA problem as a linear opti-
mization problem is explained in detail. The use of Excel solver is
demonstrated for setting up and solving several DEA problems. The use
of the shadow prices produced by the solver sensitivity report is also dem-
onstrated for the estimation of the improvement targets for lower-
performing units.
Chapter 8. Statistical Data Analytics for Decision-Making.
This chapter presents examples of the interconnection between man-
agement engineering and advanced statistical data analytics. Principal
component decomposition of the multi-variable data set to identify the
main contributing independent variables as well as cluster analysis of ZIP
codes in terms of patient contribution margin are presented. Also, a
method of forecasting the patient census time series using a recursive
linear digital filtering technique is included. An application of Bayesian
inference is presented for single and multiple diagnostic tests to demon-
strate that not all commonly used methods are effective for containing
the spread of COVID-19. In the last section, a method is presented for
calculating the expected number of patients discharged from ED using
the concept of conditional probability.
Chapter 9. Allocating Cost Savings (Gains) Between Cooperating
Providers: Introduction to the Concept of the Shapley Value.
The problem of cost allocation between cooperating parties arises in
the accounting of practically every organization. There is an advantage of
using the Shapley value for cost (savings) allocation between cooperating
providers. A simplified case of the Shapley value is also included based on
the airport landing fee framework.
Chapter 10. Summary of Fundamental Management Principles.
Concluding Remarks.
xiv Preface to the Updated and Expanded 2nd Edition
References
Berwick, D. (2011). Observations on initiating systems change in health-
care: Challenges to overcome. In Engineering a learning healthcare sys-
tem: A look at the future (Workshop summary, p. 58). Institute of
Medicine (IOM)/The National Academies Press. [Link]
[Link]?record_id=12213
Butler, T. (1995). Management science/operations research projects in
health care: The administrator’s perspective. Health Care Management
Review, 20(1), 19–25.
Buttell Crane, A. (2007, April). Management engineers (p. 50).
Hospitals & Health Networks (H&HN), Management/Governance.
[Link]
jsp?dcrpath=HHNMAG/Article/data/04APR2007/0704HHN_
FEA_Management&domain=HHNMAG
Carter, M. (2002). Health care management. Diagnosis: Mismanagement
of resources. Operation Research/Management Science (OR/MS)
Today, 2, 26–32.
Compton, W. D., & Reid, P. (2008, Spring). Engineering and health care
delivery system (editorial). The Bridge, 38(1), 3–5. National Academy
of Engineering. [Link]
Fabri, P. (2008). Can health care engineering fix health care? American
Medical Association Journal of Ethics, Virtual Mentor, 10(5), 317–319.
Preface to the Updated and Expanded 2nd Edition xix
Fowler, J., Benneyan, J., Carayon, P., Denton, B., Keskinocak, P., &
Runger, G. (2011). An introduction to a new journal for healthcare
systems engineering. IIE Transactions on Healthcare Systems
Engineering, 1, 1–5.
IOM (Institute of Medicine). (2011). Engineering a learning healthcare
system: A look at the future (Workshop summary). The National
Academies Press. [Link]
Kopach-Konrad, R., Lawley, M., Criswell, M., Hasan, I., Chakraborty,
S., Pekny, J., & Doebbeling, B. (2007). Applying systems engineering
principles in improving health care delivery. Journal of General Internal
Medicine, 22(3), 431–437.
Lawrence, D. (2010). Chapter 1: Healthcare: How did we get here and
where are we going? In Y. Yuehwern (Ed.), Handbook of healthcare
delivery systems (pp. 1.1–1.15). CRC Press.
Mayo Clinic. (2011, January 25). Mayo Clinic launches new center to focus
on how healthcare is delivered. Mayo Clinic News. [Link].
org/news2011-rst/[Link]
Reid, P., Compton, W., Grossman, J., & Fanjiang, G. (Eds.). (2005).
Building a better delivery system: A new engineering/healthcare partner-
ship. National Academy of Engineering and Institute of Medicine/The
National Academy Press.
Story, P. (2009, January). Are we thinking systems yet? American Society
for Quality (ASQ). [Link]
thinking-systems-[Link]
Valdez, R. S., & Brennan, P. F. (2009, May 10). Industrial and system
engineering and health care: Critical areas of research (Background
report). University of Wisconsin-Madison/Agency for Healthcare
Research and Quality (AHRQ) Publication, #09-0094-EF. http://
[Link]/engineeringhealthfinalreport
Vissers, J. M. H. (1998). Health care management modeling: A process
perspective. Health Care Management Science, 1, 77–85.
1 Traditional
Management, Management Engineering,
and Selection of a Portfolio of Projects 1
1.1 Traditional Management and Management
Engineering 1
1.2 Selection of a Portfolio of Projects in Healthcare
Organizations 7
References 12
2 Capacity
and Patient Flow Problems: Solutions Using
Queuing Analytics and Discrete Event Simulation 13
2.1 Queuing Analytic Modeling (QA): Its Use and
Limitations 14
2.2 Discrete Event Simulation Methodology: What Is a
Discrete Event Simulation Model and How Does a
Simple Model Work? 18
2.3 The Number of Exam Rooms and Staffing
in the Radiology Department 24
2.4 Outpatient Clinic: Centralized or Separate Locations? 32
2.5 Outpatient Clinic: Non-steady-State Operations 36
2.6 Outpatient Clinic: Limited Queue Size with Leaving
“Inpatient” Patient 38
xxiii
xxiv Contents
3 Staffing
and Scheduling Problems with Random Patient
Demand: Solutions Using Discrete Event Simulation 99
3.1 Scheduling Order for Appointments with Different
Duration Variability: Is 30 min + 60 min Always Equal
to 60 min + 30 min? 100
3.2 Centralized Discharge Versus Individual Units
Discharges103
3.3 Staffing of the Hospital Receiving Center 107
3.4 Staffing of the Unit with Cross-Trained Staff 110
3.5 Outpatient Clinic Costs and Staffing: Is the Right
Staff Used at the Right Time? 115
3.6 Staffing and Performance Analysis of the Eye Institute 121
References127
Contents xxv
4 Staffing
Problems with Random Patient Demand:
Solutions Using Analytic Optimization Techniques129
4.1 Staffing with Minimal Cost of Under- and Overstaffing
with Random Patient Demand: The “Newsvendor”
Framework130
4.2 PACU Hour-to-Hour Staffing Using Multicriteria
Optimization146
4.3 Optimized Pooled Screening Testing: Get Done
More by Doing Less 156
References159
5 Resource
Allocation and Service Line Problems: Solutions
Using Linear Optimization161
5.1 Patient Service Volumes for Three Service Lines: Keep
or Drop a Service Line? 165
5.1.1 The Solver Sensitivity Report: What Is It and
How to Read It? 168
5.2 Staffing Various Skill Mix Nursing Personnel
in the Outpatient Clinic 172
5.3 Diagnostic Test Lab: Which Tests and in What
Volumes Should Be Offered? 179
5.4 Medical Equipment Supplier’s Distribution Network.
Is the Total Shipping Cost Minimal? 182
References194
6 Staffing
and Scheduling Problems with Sliding Days Off:
Solutions Using Integer Linear Optimization195
6.1 Clinical Unit Staffing and Scheduling with Five
Consecutive Workdays and 2 Mandatory Sliding
Days Off 196
6.2 Weekly Staffing and Scheduling with Employee’s
Preferences and Seniority 200
6.3 Including Part-Time Employees: Optimization of
Full-Time and Part-Time Staffing and Scheduling
with Sliding Consecutive Days Off 206
xxvi Contents
7 Comparative
Performance of Units and Organizations
with Multiple Inputs and Outputs: Introduction to Data
Envelopment Analysis229
7.1 Eight Units with Two Inputs and Two Outputs 235
7.2 Twelve Medical Labs with Two Inputs and Three
Outputs240
7.3 A Network of Ten Hospitals with Five Inputs
and Three Outputs 244
References250
8 Statistical
Data Analytics for Decision-Making251
8.1 What Population Demographic Factors Are the
Biggest Contributors to Hospital Contribution
Margin? Introduction to Principal Components
Decomposition252
8.2 Cluster Analysis: Which Zip Codes Form Distinct
Contribution Margin Groups? 260
8.3 Forecasting of Patient Census Using Recursive Linear
Digital Filtering Technique 266
8.4 Your COVID Test Is Positive. Are You Infected?
Bayesian Inference for Single and Multiple Diagnostic
Tests277
8.5 The Expected Number of Patients Discharged from
the Emergency Department: Is the Average Length
of Stay a Good Predictor? 293
References298
Contents xxvii
9 Allocating
Cost Savings (Gains) Between Cooperating
Providers: Introduction to the Concept of the Shapley
Value299
9.1 Is Allocation of Savings Between Cooperating
Providers Fair? The Use of the Shapley Value Concept 300
9.1.1 Marginal Contribution 302
9.1.2 The Shapley Value Definition 303
9.1.3 Three-Member Group 304
9.1.4 Four-Member Group 308
9.2 Simplified Case of the Shapley Value: The Airport
Landing Fee Framework 312
References318
10 Summary
of Fundamental Management Principles.
Concluding Remarks321
References328
I ndex329
1
Traditional Management, Management
Engineering, and Selection
of a Portfolio of Projects
black box. The purpose of these advanced methods in their current state
is to infer some insights from a big empirical collection of various data
(big data) obtained “naturally” usually without any preliminary planning,
regardless of whether or not the collected data contain information
needed to address the specific management/business problem. Data col-
lected this way typically contain little useful signal and a lot of noise and
spurious correlations. Noise is a random pattern that is unstable and dis-
appears in the long run, but it is mistaken for a stable signal. The noise
can also be produced by overfitting the data model. However, it is diffi-
cult in general to tell the signal from the noise. It requires a lot of subject
matter knowledge. Data must be placed in a business context. Without
it, there is no sure way to differentiate the useful signal from the decep-
tive noise.
Current DS and ML help to explore merely some patterns and associa-
tions/correlations in the big datasets and attempt to use these patterns
and correlations to conclude on the current state of affairs, make a man-
agement decision, or expect some future patterns. However, all collected
data inevitably represent the past, i.e., what was already happened.
Although the future trend could include some statistically stable trend
from the past that was extrapolated into a short-term time horizon using
various techniques (see, e.g., Chap. 8, Sect. 8.3), in general, the future
trend is not defined solely by its past and a long-term extrapolation is
rarely if ever sustainable regardless of the amount of the past data. When
new data are collected, the data model should usually be retrained
(re-mapped).
Thus, to be sustainable, a managerial decision should include an
understanding of some hidden internal interconnections of various fac-
tors and components within a process or a system rather than merely
relying on a pattern in past data. This is where management engineering
comes into play.
In contrast to the DS and ML approach, management engineering is a
methodology for making managerial decisions that are based on the out-
comes of a validated mathematical model of the process/system with vari-
ous inputs. This model includes the structure and interconnections of the
process/system components and variables, not only input/output data.
Such a model is sometimes called the “digital twin” of the process/system.
4 A. Kolker
Hospitals tend to have more key stakeholders than other types of organi-
zations in which control rests with a relatively small number of individu-
als. Any hospital’s key stakeholder could potentially prevent or delay the
implementation of changes going forward. For a change to be imple-
mented (no matter how well it is justified), all key stakeholders must
perceive themselves to be better off with the change than without it.
Successful implementation of a well-justified management engineer-
ing solution requires that (i) a convincing need for change is effectively
communicated within the organization, (ii) the reservations and concerns
of all key stakeholders are recognized and addressed, and (iii) an imple-
mentation plan is properly executed and is sustainable.
These steps include making a convincing case to all key stakeholders
with veto power that the proposed change will offer them some net ben-
efits over implementation costs.
That’s why change within a department (say, the Emergency
Department or ICU alone) is easier than a systemic change within the
entire organization because organizational change involves more diverse
key stakeholders.
In summary, even a convincing justification of a management engi-
neering solution is not enough to make it practically implementable. If
the implementation of a solution (or change) makes at least one key
stakeholder better off and makes no other key stakeholder worse off, then
the solution is called Pareto efficient. As a consequence, a Pareto-efficient
solution has the highest chance of being implemented (Hopp & Lovejoy,
2013). Not taking this fundamental principle into account is, perhaps,
one of the main reasons for the slow implementation of management
engineering methodology for decision-making in healthcare settings
(along with the educational and technical reasons outlined above in the
report by Reid et al. (2005)).
An example of successful implementation of change based on manage-
ment engineering is presented, for example, by Ryckman et al. (2009).
The authors of this work write: “A system for smoothing flow, based on
an advanced predictive model (DES) for need, occupancy, and length of
stay, coupled with an active daily strategy for demand/capacity matching
of resources, allowed much better early planning, predictions, and capac-
ity management, thereby ensuring that all patients are placed in suitable
1 Traditional Management, Management Engineering… 7
Supply
chain Mass Purchasing Reducing supplier bullwhip effect.
Pharmaceuticals Optimized distribution network.
Reduce ordering Bullwhip effect Reorder time and amount.
Clinical resource optimization
One way to organize setting priorities and selecting areas of focus for
cost reduction portfolio of projects is to use a driver diagram to identify
and detail the core strategies (primary drivers or the primary type of
problems) and associated processes (secondary drivers or more specific
type of problems) that are most likely to lead to the cost reduction goal.
Figure 1.1 depicts an example of the driver diagram that follows Martin
et al. (2009).
The primary drivers include (i) clinical quality outcomes, (ii) staffing
and scheduling, (iii) capacity and patient flow, (iv) supply chain, and (v)
mismatched services. They are mostly applicable to almost any healthcare
organization. The primary drivers directly affect the desired business out-
come—in this case, generating “dark green” or better to say “hard” dol-
lars by reducing the annual operating budget typically by 1–3%. “Hard”
dollars are well-documented cost reductions in budgeted items. In con-
trast, the claimed savings that cannot be directly connected to the finan-
cial bottom line are referred to as “light green” or “soft” dollars.
1 Traditional Management, Management Engineering… 9
By breaking down the total wages per admission into its separate com-
ponents, different potential sources of cost reduction could be identified
(see, e.g., Chap. 4, Sect. 4.1 in this book):
great frustration for both the finance department and the improvement
team. The finance becomes skeptical of claims made by improvement
teams, while clinicians and management engineers feel that their efforts
go unrecognized and unrewarded.
Martin et al. (2009) give an example of improvement in pneumonia
care processes resulting in a switch from IV to oral antibiotics for 370
patients, for supply-cost “savings” of $319,000. However, it was difficult
to identify these savings because the year-on-year increase in drug costs
meant that the savings were lost in annual cost increases. Without this
improvement, total drug spending might indeed have been $319,000
higher than it actually was in the subsequent year. To move from “soft”
dollars to “hard” dollars, the amount saved from improvement needs to
be identified and isolated within the budget. This allows decision-makers
to reallocate funds or remove this cost from the budget. Distinguishing
between “soft” dollars and “hard” dollars—and understanding the pro-
cess by which the former is converted to the latter—is crucial in an orga-
nization that wants to enhance value. Doing so requires, according to
Martin et al. (2009), (i) a clear method for making that conversion and
for tracking costs across the organization; (ii) a clear sense of potential
savings before the project begins; and (iii) an expectation of how savings
will ultimately be applied and allocated (an approach for such allocation
is presented, e.g., in this book in Chap. 9, Sect. 9.1 “Is Distributing of
Savings between Cooperating Providers Fair? The Use of the Shapley
Value Concept”).
Martin et al. (2009) provided examples of two alternatives (along with
Excel spreadsheet templates) for tracking the cost reduction efforts. Both
use adjusted versions of standard financial figures to demonstrate if
“hard” dollars have been achieved.
In summary, the principles of selection of the problem types (primary
drivers) and project portfolio (secondary drivers) outlined in this section
were the basis for including them in this book (with removed proprietary
information and some distracting details that are not important for
this book).
12 A. Kolker
References
Chambers, C., Dada, M., & Williams, K. (2022). Improving processes for health
care delivery (p. 347). Springer.
Green, L. (2006). Queuing analysis in healthcare. In R. Hall (Ed.), Patient flow:
Reducing delay in healthcare delivery (pp. 281–307). Springer.
Hopp, W., & Lovejoy, W. (2013). Hospital operations: Principles of high efficiency
health care (p. 623). FT Press.
Joustra, P., de Witt, J., Van Dijk, N., & Bakker, P. (2011). How to juggle priori-
ties? In An interactive tool to provide quantitative support for strategic patient-
mix decisions: An ophthalmology case. Health Care Management Science.
[Link]
Kolker, A. (2014). Chapter 131: Management science for healthcare applica-
tions. In J. Wang (Ed.), Encyclopedia of business analytics and optimization
(pp. 1446–1456). IGI-Global. [Link]
management-science-for-healthcare-applications/107339
Martin, L. A., Neumann, C. W., Mountford, J., Bisognano, M., & Nolan,
T. W. (2009). Increasing efficiency and enhancing value in health care: Ways to
achieve savings in operating costs per year (IHI Innovation Series White Paper).
Institute for Healthcare Improvement. [Link]
Reid, P., Compton, W., Grossman, J., & Fanjiang, G. (Eds.). (2005). Building a
better delivery system: A new engineering/healthcare partnership. National
Academy of Engineering and Institute of Medicine. The National
Academy Press.
Ryckman, F., Yelton, P., Anneken, A., Kissling, P., Schoettker, P., & Kotagal,
U. (2009). Redesigning intensive care unit flow using variability manage-
ment to improve access and safety. The Joint Commission Journal on Quality
and Patient Safety, 35(11), 535–543.
Smiseth, O. A. (2023). Managing a hospital (p. 149). Springer. Business Guides
on the Go.
2
Capacity and Patient Flow Problems:
Solutions Using Queuing Analytics
and Discrete Event Simulation
Capacity and patient flow are one of the most practically important and
widespread types of problems in healthcare operations management. It is
widely acknowledged that the most powerful and versatile methodology
for analyzing this kind of problems is discrete event simulation (DES). At
the same time, queuing analytic modeling (QA) is often recommended as
a means of analyzing hospital capacity, patient flow, and staffing issues
(Litvak, 2007; McManus et al., 2004; Haraden et al., 2003). However,
such a recommendation underestimates some serious practical limita-
tions of QA for hospital applications but overestimates the difficulties of
using DES. In this chapter, these two methodologies (QA and DES) are
applied side-by-side to the same problems (along with a traditional man-
agerial approach) to demonstrate the pros and cons of each. The discrete
event simulation software package used in this book is ProcessModel v6.0
(the trial version can be downloaded from Process Model, Inc. website
[Link]. DES models presented in this chapter are avail-
able as supplemental electronic material SEM once the trial version is
downloaded and installed).
However, practically any other high-level DES package can also be
used such as ProModel, Arena, Simul8, AnyLogic, Simio, FlexSim, etc.
All of them provide a user-friendly graphical interface that makes the
efforts of building realistic simulation models no more demanding than
the efforts to make simplifications, adjustments, and calibrations to
develop rather complex but limited analytic queuing models. Dias et al.
(2016) and Forbus and Berleant (2022) provided a review and a com-
parative study of dozens of commercially available and open-source dis-
crete event simulation packages.
Before proceeding with practical examples of application, a brief over-
view of the above methodologies is provided in Sects. 2.1 and 2.2.
for example, limited queue size, customers leaving the system after wait-
ing a specified amount of time, multiple queues with different average
service times and different providers’ types, different service priorities,
etc. However, the use of these cumbersome formulas even built in Excel
spreadsheet functions or coded in R programming is rather limited
because they cannot capture the complexity of most healthcare systems of
practical interest.
Assumptions that allow deriving most queuing formulas are not always
valid for many healthcare processes. For instance, several patients some-
times arrive in the Emergency Department at the same time (several
people were injured in the same auto accident), and/or the probability of
new patient arrivals could depend on the previous arrivals when the ED
is close to its capacity, or the average arrival rate varies during a day, etc.
These possibilities alone make the arrival process a non-ordinary, nonsta-
tionary with after-effect, i.e., a non-Poisson process for which queuing
formulas are not valid. Therefore, it is important to properly apply statis-
tical goodness-of-fit tests to verify that the null hypothesis that actual
arrival data follow a Poisson distribution cannot be rejected at some level
of significance.
An example of a conclusion from the goodness-of-fit statistical test
that is not convincing enough can be found, for instance, in Harrison
et al. (2005). The authors tried to justify the use of a Poisson process by
using a chi-square goodness-of-fit test. They obtained the test p-value in
the range from 0.136 to 0.802 for different days of the week. Because
p-value was greater than the 0.05 level of significance, the authors did not
reject the null hypothesis of Poisson distribution (accepted the null
hypothesis). On the other hand, the fundamental property of a Poisson
distribution is that its mean value is equal to its variance (squared stan-
dard deviation). However, the authors’ own data indicated that the mean
value was not even close to the variance for at least 4 days of the week.
Thus, the use of a Poisson distribution was not actually convincingly jus-
tified for patient arrivals. Apparently, the chi-square test p-values were
not large enough to accept the null hypothesis with high enough confi-
dence (alternatively, the power of the statistical test was likely too low).
Despite its rather limited applicability to many actual patient arrival
patterns, a Poisson process is widely used in operation research as a
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 17
(ii) Discrete event simulation (DES) which operates mostly with indi-
vidual patients and their attributes. It is most appropriate for analyz-
ing business operations of separate hospitals and clinics.
(iii) Agent-based simulation (ABS) which operates mostly with the
actions and interactions of autonomous entities. It includes emerg-
ing rules of behavior that did not exist in the original model design,
and it is most appropriate for analyzing the effect of emerging indi-
vidual behavior on the system response as a whole. All three meth-
odologies can be merged if it is warranted by the problem to be
solved. However, the most powerful and versatile simulation meth-
odology used for healthcare business operations is discrete event
simulation (DES).
Table 2.1 A few random numbers sampled from a uniform inter-arrival time and
exponential service time
Inter-arrival time (min) Service time (min)
2.6 1.4
2.2 8.8
1.4 9.1
2.4 1.8
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 21
model, in which discrete events (changes) in the system are tracked when
they occur over time. In this particular example, we were tracking events
at discrete points in time t = 2.6, 4.0, 4.8, 6.2, 8.6, 13.6, 22.7.
Once the simulation is completed for any length of time, another set
of random numbers from the same distributions is generated, and the
procedure (called replication) is repeated. Usually, multiple replications
are needed to properly capture the system’s variability. In the end, the
system’s output statistics are calculated, e.g., the average patient and
server waiting time, its standard deviation, the average number of patients
in the queue, the confidence intervals, and so on.
In this example, only two patients out of four waited in the queue.
Patient 3 waited 13.6–6.2 = 7.4 min and patient 4 waited
22.7–8.6 = 14.1 min, so the simple average waiting time for all four
patients is (0 + 0 + 7.4 + 14.1)/4 = 5.4 min. Notice, however, that the first
two patients did not wait at all, while patient 4 waited 2.6 times longer
than the average. This illustrates that the simple average could be rather
misleading as a performance metric for highly variable processes without
some additional information about the spread of data around the average.
Similarly, the simple arithmetic average of the number of waiting
patients (average queue length) is 0.5. However, a more informative met-
ric of the queue length is the time-weighted average that takes into
account the length of time each patient was in the queue. In this case, it
is (1 × 7.4 + 1 × 14.1)/22.7 = 0.95. Sometimes the time-weighted average
is a better system performance metric than the simple average.
22 A. Kolker
• Flow chart of the process, i.e., a diagram that depicts the logical flow
of a process from its inception to its completion.
• Entities, i.e., items to be processed, e.g., patients, documents, cus-
tomers, etc.
• Activities, i.e., tasks performed on entities, e.g., medical procedures,
exams, document approval, customer check-in, etc.
• Resources, i.e., agents used to perform activities and move entities,
e.g., service personnel, equipment, nurses, and physicians.
• Entity routings that define directions and logical conditions flow
for entities.
1. Start with a process flow chart similar to drawing a house floor plan by
an architect. This step includes the purpose of the model, making a list
of the entities, activities, main routes, and logic of moving entities
through the activities.
2. Collect the required information to fill the model, such as an arrival
pattern, service time statistical distributions, capacity of each activity,
resources availability and work shifts, etc.
3. Collect the standard model’s building blocks and connect them
according to the flow chart similarly to connecting the house building
blocks, such as bricks, panels, drywalls, etc.
4. Similar to a house inspection, test run the model to identify some
errors in logic and connections usually using an animation feature.
How does one decide whether DES methodology is the right approach
for solving a particular problem? This is a methodology of choice for
analyzing the dynamic behavior of complex systems and processes with
random inputs and interactions of the system’s elements and feedback
loops. It should especially be preferred if there is a big decision to make
with a high potential for failure or reward. It provides a framework for
experimentation with the system and testing various management sce-
narios to reveal unintended consequences of management solutions.
Problem Description
The radiology department of the large community hospital had experi-
enced capacity and staffing issues that resulted in long patient wait times,
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 25
28
26
Average wait time, min
24
22 onset of steady-state regime ->
20
18
16
14 warm-up period
12
10
0 100 200 300 400 500 600 700 800
Simulation time, hours
Fig. 2.1 Illustration of the warm-up period estimate and the onset of the steady-
state regime with small fluctuations. Simulation with 300 replications. The dotted
lines are 95% confidence intervals, upper end, and low end, respectively, for each
data point. The warm-up period is ~240 h. The steady-state average wait time
within the confidence interval width is ~25.4 min
After spending 2 full days, a few ideas were generated at the end of this
workout, such as increasing the number of exam rooms, decreasing the
exam time, limits on the number of exam requests to match the depart-
ment’s capability, and making sure X-ray technicians are available on
time without delay. Nonetheless, the radiology director was not satisfied
28 A. Kolker
with the workout results because the generated ideas were on the surface
and self-evident but some basic capacity questions remained unanswered.
After all, no brainstorming could generate specific answers on capacity
and staffing without using some hard data and calculations that were
beyond the scope and capability of the workout participants.
-----Original Message-----
From: John Dow*
Sent: Thursday, July 20, 2006 12:16 PM
To: Alex Kolker
Subject: RE: Radiology prediction model
Alex,
I composed this email yesterday with this initial data. Here goes, I would
use 10 minutes as the mean exam time and 20 exams per hour as the
maximum # of exams with the 4 rooms we have.
What’s the average queuing for 4 rooms as well as expected queuing if we
add one more room?
Thanks, John Dow*, Director, Radiology Administration
-----Original Message-----
From: Alex Kolker
Sent: Thursday, July 20, 2006 3:59 PM
To: John Dow*
Subject: Radiology prediction model
John,
These are results but notice that I’ve applied the simplest queuing M/M/s
type model. Although some assumptions of this model are not accu-
rate, the model provides a valid relationship between resource utiliza-
tion, patient wait time and length of queue, and the number of exam
rooms. A more accurate model will follow based on discreet event
simulation of the patient flow in the radiology department.
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 29
Moreover, make a note that (i) 10 min is the average exam time that is
exponentially distributed. Actual exam time could much deviate from the
average; (ii) 20 exams per hour is the maximal peak demand that typically
happens in the afternoon hours. Nonetheless, solely for the simplicity of
the queuing model, I’ve assumed that this value is the average daily
Poisson random demand.
Queueing results with the above assumptions were obtained using the
file “Queuing spreadsheet_M_M_s_model” (SEM 2.1).
Four rooms:
Five rooms:
As you see, adding one more room makes a big difference and conve-
nience for patients significantly reducing their wait time and the queue
length. At the same time, adding one more room results in a decrease in
room utilization.
This is a general result: adding a resource results in a decrease in patient
wait time (good for patients) but at the expense of a lower resource utili-
zation for the service provider (bad for the provider). A balance could be
established between patients’ convenience and the service provider’s
resource use.
30 A. Kolker
-----Original Message-----
From: John Dow*
Sent: Thursday, July 20, 2006 4:10 PM
To: Alex Kolker
Subject: RE: Radiology prediction model
Alex,
Could you tell me what number of exams per hour with 4 rooms equates
to a queue of 1 or less?
Thanks, John*
From: Alex Kolker
Sent: Thursday, July 20, 2006 5:00 PM
To: John Dow*
Subject: Radiology prediction model
John, trial a few values for the number of exams gives the mean 16.7
exam requests per hour (for 4 rooms) to get the queue size ≤ ~ 1, and
the mean waiting time ~ 3.5 min.
Naturally, if the queue size drops (from 3 to 1) and the mean exam length
stays the same (10 min) the number of exam requests should also drop
(from 20 per hour down to 16.7 per hour).
In summary, the whole situation is another illustration of the supply and
demand problem.
Demand for service is the number of exam requests (per hour), i.e. the
patient’s arrival rate. Supply is the number of exams that can be per-
formed per hour using available resources (rooms and techs) given the
mean exam time.
The queueing depends on the balance of these two characteristics (actu-
ally their ratio): If demand significantly exceeds supply then the unlim-
ited queue grows with the unlimited increase of the waiting time. If
supply is adequate, i.e. the ratio demand/(supply × resources) < 1 then
a steady state queue is possible, the queue does not grow with time,
and waiting time is limited.
From: John Dow*
Sent: Thursday, July 20, 2006 6:10 PM
To: Alex Kolker
Subject: Re:
Radiology prediction model
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 31
Alex,
This is great—the spreadsheet as well.
Thanks much!
John Dow*, Director, Radiology Administration
(*Real name was changed)
1
s
Problem Description
A busy community clinic provides flu shots during a flu season peak on a
walk-in basis (no appointment necessary). The clinic is staffed with four
nurses. The average patient arrival rate is about 28 patients/h. Giving a
shot (including the time for necessary paperwork) takes on average
about 8 min.
During a typical clinic operation day, it was observed quite a long
waiting line and some patients complained about long waiting times. The
clinic’s operational performance should be improved.
rate λ = 28 patients/h, the average flu shot time t = 8 min, and the num-
ber of nurses (servers) N = 4.
Substituting λ = 28 patients/h, and t = 8 min = 0.1333 h in the Excel
queuing spreadsheet (available using the link to the book’s SEM in Sect.
2.3), we get the average number of patients in the queue Lq = 11.9 and
the average waiting time of about 25.5 min. The clinic’s average utiliza-
tion is 93%.
The new proposed system of two separately located clinics that was
supposed to be more convenient for patients and perform better consists
of two separate systems with N = 2. The arrival rate for each separate
clinic was going to be λ = 28/2 = 14 patients/h. Using the same Excel
spreadsheet, we get the average number of patients in each clinic queue
Lq = 12.6. The average waiting time in the queue will be about 54 min!
Thus, in the proposed “improved” clinics, the number of patients in
each clinic queue will remain about the same, while the average waiting
time will be about twice (!) that of the original operations.
It should be concluded that the proposed improvement change that
might look reasonable on the surface does not stand the scrutiny of the
simple quantitative analysis.
Although the ratio of the number of patients per server remains the
same in both cases (28:4 = 7 vs. 14:2 = 7), separating one random patient
flow on two equally divided separate random flows with the proportion-
ally divided number of servers does not result in efficiency improvement;
it makes things worse. This is because separate servers cannot help each
other if some of them become overworked for some time due to a surge
in patient arrival because of random variability of the patient flow.
The use of a discrete event simulation model (DES) to analyze the
same process without resorting to queuing analytics is given in the next
section.
elements are simply dragged down from the symbols’ pallet and then
connected in the right order.
The next step is to fill in the process information: patients arrive peri-
odically, one patient at a time. The inter-arrival time is exponentially
distributed with the average inter-arrival time of 60 min/28 = 2.14 min,
designated E(2.14), as indicated on the data arrival panel in Fig. 2.2. This
arrival pattern corresponds to the Poisson arrival rate of 28 patients/h (E
stands for exponential distribution). In the Flu_Clinic data panel, the
capacity input was 4 (four patients served concurrently by four nurses),
and the service time was an exponential random variable with an average
value of 8 min, i.e., E(8). This completes the model setup (the model is
available as the book’s supplemental electronic material (SEM 2.2)).
The model was simulated with 300 replications to capture the vari-
ability of patient arrivals and service time and get reasonably narrow con-
fidence intervals of the outputs. To identify the onset of the steady-state
output, the model was simulated for various time lengths, as it was dis-
cussed at the end of Sect. 2.2 in “A comment on the comparison of the
queuing analytic and DES models outputs” and presented in Fig. 2.1.
Fig. 2.2 Layout of the simulation model of the flu clinic. Information on the
panel indicates patient arrival type (Periodic) that repeats on average every
E(2.14) min
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 35
As was demonstrated earlier, the warm-up period in this case was about
240 simulation hours. The average steady-state number of patients in the
queue is ~12.2 (95% CI is 11.3–13.2), and the average steady-state wait-
ing time (non-value added time) is ~25.4 min (95% CI is 24.8–26.1 min).
This is the average of all fluctuating output values beyond the warm-up
of 240 h. If only ten replications were run with, say 600 h, to save the
overall simulation duration time, then the average number in the queue
and the wait time would be naturally somewhat different from the above
values, ~12.5 (95% CI 9.9–15.0) and ~26.6 min (95% CI,
21.3–32.1 min), respectively. Notice that the confidence intervals become
much wider with fewer replications. Less accurate simulation results
(wider confidence intervals and larger fluctuations) are always a price to
be paid for shorter simulation duration time.
Thus, with the simple DES model, we got practically the same steady-
state results as with the QA model using no analytic formulas at all.
Similarly, for the separately located clinics with the arrival rate of
14 patients/h and 2 nurses (capacity), the same DES model (with 600 h
of the simulation run and 300 replications) results in the average steady-
state 12.3 patients in the queue (95% CI 11.8–12.8) and the average
waiting time ~52.5 min (95% CI 50.4–54.6 min), i.e., again practically
the same results as given above by QA formulas. Of course, these steady-
state values will somewhat fluctuate depending on the simulation run
length beyond the warm-up period of 240 h and the number of
replications.
Thus, both QA and DES models illustrate a fundamental management
engineering principle: combined resources with random patient flow and
with unlimited queue size and no leaving patients are more efficient than
separate resources with the same total workload. If specialized (dedicated)
resources are needed due to patient privacy, infection control, non-
movable equipment, or other special factors, then some additional capac-
ity should be planned and budgeted to cover the loss of resources’
efficiency. Specialized resources (staff, operating or procedure rooms,
beds, etc.) typically cost more than mutually interchangeable (pooled) or
shared resources (see also Sect. 2.12).
36 A. Kolker
Problem Description
The manager of the same clinic decided to verify that the clinic would
operate smoothly enough with a new team of less experienced nurses
(servers) who would work only a little slower than the previous one. The
average time to give a shot will be about 8.6 min (instead of the average
8 min for more experienced staff, as in the previous section).
275
250
Patient wait time, min
225
200
175
150
125
100
75
50
25
0
0 100 200 300 400 500 600 700 800
Simulation time, hours
Fig. 2.3 DES simulated average patient wait time. Upper curve is the non-steady-
state process (average patient service time is 8.6 min). Bottom curve is the steady-
state process (average patient service time is 8.0 min)
38 A. Kolker
obtained from the queuing formulas. (A similar plot for the average num-
ber of patients in the queue for the steady-state and non-steady-state
clinic operations is not shown here to save space.)
This example also illustrates an important principle of “unintended
consequences.” An intuition that is not supported by the objective quan-
titative analysis says that a small change in the system input (service time
from the average of 8 min to 8.6 min) would result in a small change in
the output (small increase in the number of waiting patients and their
waiting time). For some systems, this is indeed true. Systems in which the
output is always directly proportional to the input are called linear
systems.
However, there are quite a few systems in which this simple reasoning
breaks down: a small change in the value of the system’s input parameter(s)
results in a dramatic change in the system’s output (behavior), e.g., from
a steady-state regime to a non-steady-state regime with unlimited growth
(sometimes, an opposite effect happens: a large change in the system
input results in a very small and weak system’s output showing low sensi-
tivity to the input). Such systems are called nonlinear. They often exhibit
complex behavior that is difficult to foresee without some type of numeric
analysis even though such systems can consist of only a few elements
(Schuster, 1998).
Problem Description
Unlimited queue size is not always a good model of real systems. In many
cases, patients wait in a waiting lounge that has usually a limited number
of chairs (space). QA models designated M/M/s/K are available that
include a limited queue size, K (Green, 2006; Lawrence & Pasternak,
2002; Hall, 1990). However analytic formulas become very cumber-
some. If the QA model should also include some patients who leave the
system after waiting some time in the queue before the service starts
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 39
Fig. 2.4 Layout of the model with patients reneging after random wait time in
the queue with the triangle distribution T(10,15, 25) min indicated on the input
data panel. The queue size limited capacity is 15
6.5
a
Average patient wait time, min
5.5
4.5
4
0 10 20 30 40 50 60 70 80 90 100 110 120
Simulation time, hours
3.4
b
Average number of patients in
3.2
3
2.8
queue
2.6
2.4
2.2
0 10 20 30 40 50 60 70 80 90 100 110 120
Simulation time, hours
Fig. 2.5 Average patient wait time (a) and average number of patients in the
queue (b) for a clinic with limited queue size (15) and “inpatient” patients leaving
the queue after waiting from 10 to 25 min with the most likely 15 min
Problem Description
In Sect. 2.5, the queuing system with unlimited queue size (average
arrival rate 28 patients/h and average shot time 8.6 min) was the non-
steady-state one, for which QA model formulas could not be used.
42 A. Kolker
However, the clinic’s manager realized that the average patient arrival
rate varies significantly during the day and that 28 patients/h was actually
a peak arrival rate, from noon to 3 p.m. (similar to the input peak demand
of 20 exams/h in Sect. 2.3). In the morning hours from 8 am to 10 am,
the arrival rate was lower, 15 patients/h. From 10 am to noon, it was
20 patients/h, and in the afternoon from 3 pm to 6 pm, it was about
23 patients/h.
How does the time-varying arrival rate affect the clinic’s operational
performance?
Traditional Approach
Because the arrival rate before noon is lower than the peak value in the
midday and then patient arrival slows down again, this should somewhat
compensate for the peak value. Therefore, the average arrival rate should
be lower than the peak one. It can be calculated as (15 + 20 + 28 + 23)/
4 = 21.5 patients/h. Thus, the clinic’s operational performance should be
close to the steady-state and acceptable.
Problem Description
There is a ten-bed ICU unit. The average daily patient arrival rate is two
patients per day (but the actual daily arrival rate varies from one to three
patients depending on the day of the week). Patient length of stay (LOS)
46 A. Kolker
is in the range from 1 to 3 days, with 2.5 days being the most likely. It is
observed that the daily average number of ICU patients (occupied beds)
is 5. The manager believes that the average daily utilization of 5 beds out
of the total available 10 beds, i.e., 50%, is too low. The manager wants to
improve the use of bed efficiency of the unit.
What is the probability that more than six beds will be needed with the
average daily five patients, P(#beds > 6)?
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 47
n 6
5n
P # beds 6 1 exp 5 24%
n 0 n!
Thus, more than 6 beds will be needed for about a quarter of the time
which is close to the above point estimate of ~28%. Therefore, the tradi-
tional management approach would create a regular bed shortage.
To apply QA formulas to estimate an average wait time and the num-
ber of waiting patients, a usual QA assumption should be made that the
length of stay in the ICU is exponentially distributed and, of course, that
patient arrival is a Poisson process.
Using the M/M/s model with the average inter-arrival rate
2/24 = 0.08333 patients/h and the average service time (LOS) of
2.5 days × 24 = 60 h (SEM 2.1), it is easy to calculate that for the six-bed
ICU, the average number of patients in the queue is 2.9, the average
waiting time is 35 h (about 1.5 days), and the average ICU utilization is
about 83%.
Notice that the weekly average arrival rate remains the same as above:
(4 + 4 + 1 + 3 + 1 + 1 + 0)/7 = 2 patients/day.
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 49
The simulation model layout is the same as above except that the input
daily arrival pattern is used instead of the periodic exponential random
arrival pattern. The steady-state simulation (warm-up 3600 h, 300 repli-
cations) resulted in a weekly average number of patients in the queue of
~0.29 and an average wait time of ~3.4 h. About 19.7% of patients would
wait more than 6 h. This is much better than the above simulation sce-
narios, but the wait time still might not be good enough to meet the
safety and quality standards for many patients waiting for admission into
ICU. Other scenarios could be tested.
This example illustrates a fundamental management engineering prin-
ciple: because of the variability of the number of patient arrivals and
length of stay, some degree of reserved capacity (sometimes up to 50%) is
needed to avoid regular operational problems.
Another possible way of improving the unit’s operations is the daily
load leveling (smoothing) of the number of elective procedures scheduled
for patients with required postsurgical ICU admission, as illustrated in
Sect. 2.11 and also in (Kolker, 2009).
In summary, capacity and staffing decisions/planning based only on
averages without taking into account particular arrival and service time
distributions usually result in significant miscalculations. This is called
the flaw of averages (see also Costa et al., 2003; Marshall et al., 2005; de
Bruin et al., 2007; Savage, 2009).
An illustration of the above principle in practice was also provided by
Green (2004). The American College of Obstetrics and Gynecology rec-
ommended that the occupancy level of obstetrics units should not exceed
75%. However, many hospitals had obstetrics units operating below this
level. According to the 1997 Institutional Cost Report, 117 out 148 of
New York state hospitals had average utilization (occupancy) below 75%.
Therefore, some hospitals have eliminated beds to reduce “excess” capac-
ity and costs. However, in a few years, these hospitals had to consider
restoring the beds because of multiple complaints of excessive admis-
sion delays.
50 A. Kolker
Problem Description
Frequently mixed patient arrival patterns exist, i.e., some patients are
scheduled to arrive at a specific time, while other patients arrive unex-
pectedly at random points in time. For instance, some clinics accept
patients who make an appointment but also accept urgent random walk-
in patients. Operating room suites schedule elective surgeries when sud-
denly a trauma patient arrives and emergency surgery is required. Such
mixed patient arrival patterns with different degrees of variability require
special treatment.
Suppose that there is one procedure room and there are eight proce-
dures scheduled for a day, at 6 am, 8 am, 10 am, 12 pm, 2 pm, 4 pm,
6 pm, and 8 pm. On this day six random emergency patients also arrived
with the average inter-arrival time of 4 h. The total average number of
patients for one day is 14. The procedure time is exponentially distrib-
uted with an average of 1 h.
What is the operational performance of this simple unit in terms of the
number of patients in the queue and their wait time?
Fig. 2.6 Mixed patient arrival pattern: random and scheduled. The panel demon-
strates a scheduled arrival pattern from 6 am to 8 pm every 2 h
52 A. Kolker
Table 2.6 Summary of three scenarios with various degrees of randomness for
mixed patient arrivals patterns: simultaneous random and scheduled arrivals
Average
number of
patients in Average
the queue wait time
Scenario (Lq) Wq (h) Note
Queuing model M/M/s 0.82 1.4 All arrivals are fully
exponentially
random all day
around
DES. Two random arrivals: 0.69–0.87 0.87–1.08 Intermediate
One with average every (95% CI) (95% CI) randomness: One
4 h all day around; arrival is limited to
another one with average the range 6 am to
every 2 h from 6 am to 8 pm
8 pm
DES. Two arrivals: One with 0.49–0.64 0.7–0.9 Low randomness; one
average every 4 h all day (95% CI) (95% CI) arrival is scheduled
around; another one is non-random
scheduled non-random
exactly every 2 h from
6 am to 8 pm
Problem Description
Hospitals or hospital departments often try to benchmark their opera-
tional performance against some other hospitals that are considered
industry best practices. For instance, if one hospital has a low patient
admission wait time, then another hospital with a higher admission wait
often launches an improvement project to meet that low admission target
that another hospital already achieved. That hospital studies admission
and other processes that help the best achiever to get its impressive per-
formance and then tries to implement what it learned hoping to get simi-
lar results.
Suppose that there are two hospitals: one is a small rural hospital with
a total of 25 beds, and another one is a large community hospital with a
total of 250 beds. The large hospital claims that it admits practically all
patients with no waiting time. The small rural hospital acknowledges that
its admission time is typically at least about 2 h and that a significant
percentage of patients wait more than 2 h. The small hospital wants to
learn and implement the best admission practices from the large counter-
part to get a similar low admission time.
It was observed that the small hospital typically scheduled one patient
a day and had about four emergency admissions. The large community
hospital typically admitted 10 scheduled patients a day and had about 40
emergency admissions.
It turns out that the ratio of the total number of admitted patients to
the number of hospital beds is exactly the same: 5 over 25 for the small
hospital and 50 over 250 for the large hospital. The management of the
small hospital reasoned that because those ratios are the same, the patient
length of stay is also the same, and the implemented admission process
steps are similar for both hospitals, their operational performance in
terms of admission wait time should also have been the same, i.e., very
low or none at all.
Yet, despite their best efforts, the waiting time for small rural hospitals
was still high enough, about 2 h. The management of the small hospital
management was puzzled.
practice, this is not the case. Empty hospitals would take a lot of time for
simulation warm-up to get a steady-state output for a valid comparison.
To eliminate or reduce a warm-up period, it makes sense to pre-fill the
hospitals with some initial number of patients. Let’s pre-fill a large hospi-
tal with 200 patients, and a small one with 20 patients, so the ratios pre-
fill to full capacity will stay the same: 200/250 and 20/25, accordingly.
This can be done by making one-time periodic arrivals of 200 or 20
patients (with 0 repeats), respectively. Two simulation models were run
with 100 replications for a full 45 weeks, i.e., 7560 h (SEM 2.8 and 2.9).
The steady-state simulation results presented in Table 2.7 are
instructive.
The large 250-bed hospital has indeed much lower (~10 times) admis-
sion wait time and a much lower % of patients waiting longer than 2 h
than the small hospital, while both have practically the same utilization.
This is an illustration of an important fundamental management engi-
neering principle: the performance of scale. Size does matter. Large hos-
pitals (units) always have better operational performance characteristics
(lower wait time and number of patients in the queue, higher utilization)
than small units with the same input relative to their size. This is a reflec-
tion of a nonlinear scale effect that is typical for healthcare organizations:
random arrival fluctuations are easier to spread out or smooth out over
larger available capacity. Failure to take this effect into account usually
invalidates simple linear benchmarking and proportional adjustments.
If the small hospital wants to reduce its admission waiting time, but
the patient arrival rate is beyond its control and the number of beds is
fixed, then the only option is reducing patient length of stay. For instance,
DES modeling results presented in the last line in Table 2.7 indicate that
the admission wait time as well as the percent of patients waiting more
than 2 h for the small 25-bed hospital will be low and in line with the
large 250-bed hospital if the maximal LOS does not exceed 5 days
(120 h). However, in this case, the average capacity utilization is also
significantly reduced, which is usually an undesired effect. Lower utiliza-
tion is the price paid for lower admission wait time similar to the take-
away in Sect. 2.3. Attempts to increase utilization (the desired effect) will
inevitably increase of the admission wait time (undesired effect).
In summary, process improvement efforts based on linear benchmark-
ing and simple linear proportional adjustments of input values could be
misguided and short-lived if the scale effect (organization size) is not
taken into account.
Problem Description
In most hospitals, random (emergency) surgeries compete for the same
operating room (OR) resources as scheduled (elective) surgeries. While
the variable number of daily emergency surgeries is beyond hospital con-
trol (this is a natural variability), there is a significant variation in the
number of daily scheduled elective surgical cases that could be actively
managed using a hospital scheduling system (Litvak & Long, 2000;
Kolker, 2009).
It is possible to manage the scheduling of elective cases in such a way
that smoothens the overall patient flow variability. A daily load leveling of
elective cases would reduce the chances of excessive peak demand for the
system’s capacity (operating rooms and ICU) and, consequently, would
reduce patient waiting time.
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 57
In 2011, the Leapfrog Group Hospital Survey included in the new sec-
tion Patient Experience of Care the use of operational management
methodology (management engineering) and, specifically, smoothing
elective patient scheduling (Leapfrog, 2011).
What is the quantitative effect of the daily load leveling of elective
surgeries on delay to start the case in the presence of the competing
demand from random emergency surgeries for OR resources?
Table 2.8 Elective, emergency, and daily load-leveled admissions for the
4-week period
Number of Number of Number of daily-
Day of elective emergency leveled (smoothed)
Week week admissions admissions elective admissions
1 Monday 9 16 6
1 Tuesday 11 14 7
1 Wednesday 8 14 7
1 Thursday 5 20 7
1 Friday 5 15 7
2 Monday 10 18 7
2 Tuesday 13 20 7
2 Wednesday 11 9 7
2 Thursday 8 11 7
2 Friday 3 20 7
3 Monday 5 17 7
3 Tuesday 9 11 7
3 Wednesday 8 15 7
3 Thursday 6 15 7
3 Friday 6 20 7
4 Monday 7 15 7
4 Tuesday 4 13 7
4 Wednesday 3 12 7
4 Thursday 4 11 7
4 Friday 3 20 6
Total 138 306 138
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 59
A smoothed schedule (last column in Table 2.8) has the same number
of cases over the 4-week period (138) as the original un-smoothed sched-
ule, i.e., not a single case was dropped but they are rather re-arranged
over the period to smooth daily peaks and valleys.
It was assumed that three interchangeable operating rooms (ORs) are
available in this case. The emergency surgery duration is in the range of
1.5–2.5 h, with this most likely time of 2.1 h. Elective surgery duration
is in the range from 1.5 to 3 h, with the most likely time of 2.4 h. Both
are represented by corresponding triangle distributions similar to those
used in the previous sections.
Emergency patients were assumed to arrive uniformly randomly 24 h
a day. Elective surgeries were scheduled daily from 7 am to 3 pm.
The simulation model layout is similar to the one presented in the
previous section (SEM 2.8).
Simulation for 672 h (4 weeks) and 100 replications for the original
un-smoothed elective schedule along with competing emergency cases
resulted in the average patient waiting time of 0.74 h for emergency cases
(95% CI is 0.72–0.76 h) and 1.23 h for elective cases (95% CI is
1.20–1.26 h).
The same simulation modeling with a smoothed (load-leveled) elective
schedule along with the same competing emergency cases resulted in the
average patient waiting time of 0.59 h for emergency (95% CI is
0.57–0.61 h) and 0.86 h for elective cases (95% CI is 0.83–0.88 h) (SEM
2.10 and 2.11).
Thus, in this particular example, elective daily load leveling results in
about a 25% reduction in waiting time for emergency surgeries and
about 39% reduction in waiting time for elective cases.
Elective schedule smoothing (daily load leveling) is indeed a very pow-
erful approach to reducing patient waiting time and improving efficiency.
This is a fundamental management engineering principle.
A simple simulation model also allows testing the effect of the different
smoothing schemes. For instance, if nearly the same daily number of
elective cases is not possible due to some practical limitations, it is possi-
ble to test another less perfect smoothing scheme to make sure that the
result is still worth the efforts of its implementation (or maybe not). No
60 A. Kolker
Problem Description
The issue of specialized vs. pooled (interchangeable) operating rooms
(ORs) caused a controversy in the literature on healthcare improvement.
Specifically, if surgical cases include both scheduled (elective) and
unscheduled emergency surgeries starting at random times, is it more
efficient to reserve specialized operating rooms (ORs) dedicated sepa-
rately for scheduled and emergency surgeries, or is it better to perform
both types of surgeries in any available OR using pooled or interchange-
able OR arrangement?
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 61
consider an OR suite with two operating rooms, OR1 and OR2. Both
emergency (random) and scheduled surgeries are included.
Let’s first consider the situation when the majority of cases are sched-
uled elective surgeries. Six elective surgeries are scheduled 5 days a week:
Monday to Friday at 7 am, 9 am, 11 am, 1 pm, 3 pm, and 5 pm.
Emergency surgeries are assumed to start independently randomly
24 h a day with the average patient inter-arrival time of 6 h (a Poisson
arrival rate is 0.166 patients/h).
Similar to the values used in the previous section, the most likely
scheduled surgery duration is assumed to be 2.4 h, and emergency sur-
gery duration is assumed to be 2.1 h (Wullink et al., 2007). The assumed
variability is from 1.5 to 3 h for scheduled surgeries and from 1.5 to 2.5 h
for emergency surgeries.
Using these arrival and service time data, let us consider two scenarios.
Scenario 1: there are two ORs, one is specialized and dedicated only for
scheduled surgeries (OR1) and stays open from 7 am to 7 pm; another
one is specialized and dedicated only for emergency surgeries (OR2)
and is open 24 h a day, as shown on Fig. 2.7. If the dedicated OR is
Fig. 2.7 Simulation model layout of the dedicated operating room for scheduled
elective OR1 and emergency surgery OR2
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 63
not available, then the new patient waits in the queue area until the
corresponding dedicated OR becomes available (SEM 2.12).
Scenario 2: There are also two ORs. However, they are pooled, i.e., fully
interchangeable: both emergency and scheduled patients go to any
available OR, as indicated in Fig. 2.8. If both ORs are not available,
then patients wait in the queue area until one of the ORs becomes
available. Emergency patients have higher priority and move first to
the available OR (SEM 2.13).
Table 2.9 Specialized ORs vs. pooled ORs. Most surgeries are scheduled electives
Dedicated OR Pooled OR
Characteristics Elective Emergency Elective Emergency
95% CI of the number of 28–29 20–21 30 20–21
performed surgeries
95% CI of the average wait time (h) 2.5–2.9 0.4–0.49 0.14–0.16 0.18–0.21
95% CI of the average number of 0.7–0.8 0.08–0.1 0.03–0.04 0.03–0.04
patients in queue
Average weekly OR utilization (%) 99 35 46
Table 2.10 Specialized ORs vs. pooled ORs. Most surgeries are unscheduled
emergency
Dedicated OR Pooled OR
Characteristics Elective Emergency Elective Emergency
95% CI of the number of 12 53–54 12 59–60
performed surgeries
95% CI of the average wait time (h) 0.17–0.19 6.6–7.6 0.8–0.9 0.62–0.71
95% CI of the average number of 0.02–0.03 3.6–4.3 0.11–0.12 0.32–0.38
patients in queue
Average weekly OR utilization (%) 55 92 62
Simulation results for 120 h (5 days) using 300 replications are given
in Table 2.10.
Notice that in this case, the number of performed emergency surgeries
is higher for pooled ORs. Although the average waiting time for sched-
uled elective surgeries increased up to about 0.8–0.9 h for pooled ORs,
the average wait time for emergencies dropped dramatically by more than
an order of magnitude, from about 7 h down to about 0.6 h (compare
this dramatic drop to Wullink et al.’s (2007) result). A similar picture is
observed for the average number of patients in the queue.
Overall, these results support the conclusions of Wullink et al. (2007)
that performing emergency surgeries in the interchangeable OR scenario
is more efficient than in the reserved dedicated emergency OR. These
authors provided a detailed instructive discussion on why the dedicated
OR scenario performs worse, especially for emergency surgeries, while
intuitively it seems that it should perform better, as Haraden et al. (2003)
assumed.
Wullink et al. (2007) pointed out that besides reserving OR capacity
for emergency surgery arrivals, ORs need to reserve capacity to cope with
the variability of surgery duration. In the pooled OR scenario, the reser-
vation might be shared to increase the flexibility for dealing with unex-
pected long case duration and emergency surgery, whereas the dedicated
scenario does not offer the opportunity to use the overflow principle.
66 A. Kolker
Problem Description
A specialized procedure operating rooms (SPR) unit is being planned to
unload the volume of outpatient day surgeries from the main surgical
department’s general operating rooms. At the end of the current year,
1813 special surgical procedures of type 1 and type 2 were performed.
The projected additional procedure volume for the next year is 179 pro-
cedures of type 1 (gastroenterology) and 13 procedures of type 2
(pulmonary).
Patients get prepared for procedures in the preparation bed area; then
they move to the available SPR for the procedures and move back into
the same bed area (not necessarily into the same bed) for post-procedure
recovery. After initial post-anesthesia recovery time (phase 1), inpatients
are moved to a regular nursing unit for full recovery. Outpatients stay in
the bed area for the full recovery time (phase 1 and phase 2).
The average time to perform a procedure is 0.82 h with the time dis-
tribution presented in Fig. 2.9 (panel a). The average post-procedure
recovery time is 3.3 h with time distribution presented in Fig. 2.9 (panel
b). SPR turnover time ranges from 10 to 20 min (cleaning the room, re-
stocking it with supplies, and making it ready for the next patient). Bed
turnover time (cleaning, changing linens, and making the bed ready for
the next patient) is also from 10 to 20 min. The SPR unit is supposed to
work annually 255 days from 7 am to 5 pm (10 h daily, no weekends).
Utilization was targeted at about 85%.
The following operational SPR performance criteria have been estab-
lished: (i) patient wait time to get into SPR for a procedure is less than
1 h for 98% patients and (ii) post-procedure patient wait time to get back
into the recovery bed from SPR is less than 5 min for 98% patients.
68 A. Kolker
Fig. 2.9 Special procedure room (SPR) procedure time (a) and recovery
bedtime (b)
Beds annual patient volume average bed time including turn over
/ total annual service hours utilization
1813 179 13 3.3 0.25 / 2550 0.85 3.3 ~ 3
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 69
SPRs annual patient volume average procedure time including turn over
/ total annual service hours utilization
1813 179 13 0.82 0.25 / 2550 0.85 0.99 ~ 1
Fig. 2.10 Simulation model layout of patient flow for special procedure operat-
ing rooms (SPR) and preparation and post-procedure recovery beds
Table 2.11 The effect of the number of preparation and post-procedure recovery
beds and special procedure rooms (SPR) on operational performance
characteristics
Percent
of
patients
Percent Average waiting
of post- longer
patients procedure than
Average waiting wait time 5 min to
patient longer to get back get back
wait time than 1 h to bed for for
to get to to get to recovery recovery
Number Number SPR (min) SPR (%) (min) (95% (%) (95% Performance
of beds of SPR (95% CI) (95% CI) CI) CI) criteria met?
3 1 35.4–37.2 33.9–34.7 22.2–23.4 23.6–24.6 No
4 1 42.0–43.2 27.4–28.1 3.6–3.7 9.7–10.2 No
5 2 3.0–3.1 0.8–0.95 1.8–2.4 5.9–6.3 No
6 2 3.0–3.6 0.9–1.1 0.6–0.7 1.7–1.95 Yes
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 71
will wait longer than 5 min to get back into bed from SPR for recovery,
respectively (vs. acceptable limit 2%).
Even five beds and two SPRs are not enough because the percentage of
patients waiting longer than 5 min to get back to post-procedure beds
exceeds accepted 2%. A simulation model indicates that the minimal 6
beds and 2 SPR are required to meet operational performance criteria. In
this case, only close to 2% of patients will wait longer than the accept-
able limits.
Although the average wait time is acceptable for all scenarios (much
less than 1 h to get to SPR and less than 5 min post-procedure to get back
to beds for all but the first scenario), the percentage of patients exceeding
1 h and 5 min limits of 2%, respectively, is high for all but the last one.
This is a consequence of wide time distributions around the averages.
Thus, the correct amount of resources in healthcare settings with high
input variability can only be calculated using process simulation method-
ology that takes into account statistical distributions of the inputs. The
use of only averages, as inputs, results in highly inaccurate calculations
and should be avoided. This is another illustration of the flaw of averages,
as it was already discussed earlier in Sect. 2.8.
Problem Description
A typical large community hospital includes four main units: an emer-
gency department (ED) with 25 beds; an intensive care unit (ICU) with
49 beds; a surgical department with 12 operating rooms (OR); and regu-
lar nursing units (NU) with a total capacity of 360 beds.
The total typical monthly patient volume was about 4478 including
ED patients transported by ambulance, ED walk-in patients, and surgi-
cal scheduled patients. There was also a significant percentage (~19%) of
patient readmission within 30 days of discharge.
72 A. Kolker
After a few months of a pilot project that tested most of these ideas, the
ED team reported a significant reduction in the average length of stay
(ALOS) and a reduction in ED diversion time. The management praised
a great ED improvement success.
However, both the OR and ICU started reporting increased patient
wait time to get in due to “no available OR or no ICU beds.” This, in
turn, resulted in an increase in the cancellation rate for scheduled surger-
ies and kept more postsurgical patients boarded in OR waiting for ICU
beds. Hence, more surgical cases were delayed.
The hospital management wanted to repeat an ED improvement suc-
cess by initiating the ICU or surgical department process improvement
projects. However, the management was not sure anymore that depart-
mental process improvement would be translated into the overall hospital
patient flow improvement. They were looking for a better analysis and a
sustainable hospital-wide solution.
It turns out that patient flow is a property of the entire hospital system
rather than the property of the separate departments/units. A detailed
analysis was required of the overall hospital system patient flow and the
interdependency of subsystems/units to establish the system’s weak link
and the right units for process improvement projects. Such an analysis
could only be performed using the DES methodology (Kolker, 2013).
Large systems are typically deconstructed into smaller subsystems
using natural breaks in the system: emergency, surgical, ICU, floor nurs-
ing units, etc. The subsystems can be modeled and analyzed separately.
However, they should be reconnected in a way that recaptures the most
74 A. Kolker
Fig. 2.11 Simulation model layout of patient flow for a typical entire hospi-
tal system
surgery, about 20% of admitted patients move into the ICU, and about
10% of patients are admitted from ED into the floor nursing units (NU).
When ED, OR, ICU, or NU is completely full (at full capacity), a
diversion status is declared. The units stay on diversion until at least one
bed in the unit is freed. Total unit diversion is defined here as the percent-
age of operational time when the unit is at full capacity and can no longer
accept new patients.
OR suite has 12 interchangeable operating rooms used both for ED
emergency and scheduled surgeries. There are four daily scheduled OR
cases at 6 am, 9 am, 12 pm, and 3 pm, Monday to Friday (there are no
scheduled surgeries on weekends). Scheduled cases form a separate OR
admissions flow, as indicated in Fig. 2.11.
Elective surgery duration depends on the surgical service type, such as
general surgery, orthopedics, neurosurgery, etc. For simplicity of this par-
ticular model, elective surgery duration was weighted by each service per-
centage, and the best statistical distribution fit was identified (inverse
Gaussian in this case). Emergency surgery duration was the best fit by
Pearson 6 statistical distribution.
About 40% of postsurgery patients were admitted from the OR into
the ICU (direct ICU admission), while 60% were admitted into the floor
NU. However, some patients (about 5%) are readmitted from the floor
NU back to the ICU (indirect ICU admission from the OR). The model
included 1–30-day patient readmission feedback loops. It was reported
(Jencks et al., 2009) that almost one-fifth (19.6%) of patients nationwide
who had been discharged from a hospital were re-hospitalized within
30 days, including ~13% back to ICU and ~6% back to floor NU.
ICU length of stay is assumed to range from 1 to 3 days, with the most
likely of 1.5 days, represented by a triangle distribution. Kolker (2009)
developed a detailed ICU simulation model and analysis.
Patient length of stay (LOS) in NU was assumed to range from 2 to
10 days, with the most likely of 5 days, also represented by a triangle
distribution. At the simulation start, the ED, ICU, and NU were pre-
filled with the midnight census of 15, 46, and 350 patients, respectively.
Pre-filling is used instead of a simulation warm-up period in this example
to approximate the steady-state regime. The DES model is SEM 2.15.
76 A. Kolker
(continued)
77
Table 2.12 (continued)
78
1 2 3 4 5 6
% patients 49.7–57.2 65.2–71.7 Much worse 54.5–60.6 Worse but not
waiting much
admission to
ICU from ED
A. Kolker
longer than 1 h
% patients 9.1–10.4 9.5–10.9 The same 10.5–11.7 About the
waiting same
admission to
floor_NU from
ED longer than
1h
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 79
Problem Description
The Institute of Medicine has advised state health departments to prepare
for overwhelming mass casualty disasters (IOM, 2009, 2012). While
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 81
patients who are either will not benefit from care or too sick to survive
regardless of the amount of used resources.
• A total capacity of 280 PICU beds was available in the event of CSC
activation, as determined by data included in the disaster management
plan of each hospital.
• Pattern of daily patient arrivals and patient volumes during the pan-
demic. To derive estimates of CSC patient volume, we used projec-
tions published by the DHHS (2005) describing a typical severe
pandemic, which estimates that approximately 6600 children per one
million population will require intensive care. Since this value is only
an approximation, the patient volumes were actually ranging from
5000 to 10,000+. Daily numbers of patients in PICUs were estimated
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 85
Fig. 2.14 The pattern of random daily arrival of children presenting with respira-
tory failure
life time, T1/2. The percentage of those who died while waiting for a
PICU admission was calculated as:
Given the input variables described above, the goal of the simulation
was to determine the specific thresholds T_POD and T_DV above which
patients would not be admitted that would yield the maximal survival
and bed occupancy over a range of patient volumes. The procedure for
maximizing the total survival is described below.
where w1, w2, and w3 are relative weights of the corresponding compo-
nents of the OF.
To place a higher priority on minimizing mortality over maximizing
bed occupancy, the weights w2 and w3 were assigned values twice the
weight of w1. The OF optimization is always a trade-off between the
desired and undesired trends of the contributing components. The three
components of the OF move its value in different directions: higher aver-
aged occupancy makes a higher OF value which is a desired trend, but
higher mortality reduces the OF value which is an undesired trend.
1. POD and DV values for each patient were generated as random vari-
ables from the statistical distributions derived from the VPS dataset.
These statistical distributions were produced using prediction equa-
tions with clinical and demographic variable inputs used from actual
patients in the VPS dataset. These best-fit statistical distributions
were built into the simulation model action logic.
2. The first trial set of the thresholds T_POD and T_DV was generated
by the model’s optimization module Sim_Runner.
3. POD and DV values were compared to the trial T_POD and T_DV
values from step 2. If POD < T_POD and DV < T_DV, the decision
was made to admit for treatment. All others were rejected from
treatment.
4. Patients triaged into the admitted group with no PICU bed immedi-
ately available were tested for exceeding the maximum wait time in
the queue. Those exceeding the corresponding limits were counted as
expired before receiving treatment (reneging route on the model lay-
out in Fig. 2.13).
5. The treatment group (those that passed triage and were admitted to
a PICU bed and existing patients occupying beds at the start of CSC
activation) was tested for mortality using best-fit equations to quan-
tify the segment of the treatment group that died.
6. Occupancy rates and the percentage of deaths (treated and untreated)
were calculated and incorporated into the simulation model to cal-
culate the objective function OF.
7. Steps 1–6 were replicated using a new set of random variables gener-
ated from POD and DV statistical distributions used in step 1 to
capture statistical variability and produce outputs with reasonably
narrow confidence intervals.
8. A new set of thresholds T_POD and T_DV was generated by the
optimization module Sim_Runner.
9. Steps 1–7 were repeated, comparing the new value of the objective
function (OF) with previous results. Threshold combinations that
yielded greater OF results were prioritized preferentially compared to
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 89
Simulation Results
The simulated optimal specific thresholds as a function of patient volume
are displayed in Figs. 2.16 and 2.17.
Because several of the model’s input parameters were fixed (PICU
surge bed capacity = 280, wait time limit = 1 h, T1/2 = 2 h), simulation
optimization results were reduced to a function of patient volume. (The
fixed parameters could, of course, be changed for another model applica-
tion.) The FCFS admission was simulated by assigning the triage thresh-
olds the values T_POD ≤ 100% and T_DV ≤ 18 days, which removed
admission selection and effectively made admitted all randomly arrived
patients provided a PICU bed was available.
Percent survival for the optimal admission triage and the correspond-
ing FCFS survival are presented in Tables 2.13 and 2.14.
It follows from Tables 2.13 and 2.14 that the optimal triage admission
resulted in improved survival of those admitted to the PICU and greater
survival among all patients compared to FCFS. In addition, the survival
advantage increased as patient volume increased.
The effect of the triage thresholds on PICU occupancy was also
examined.
PICU occupancy was averaging between 79% and 88%. Daily occu-
pancy values are displayed in Fig. 2.18 for 5000 and 10,000 patient
volumes.
With the initiation of the CSC algorithm, occupancy drops, indicat-
ing more selective criteria for PICU entry are deployed compared with
FCFS. As crisis standard care (CSC) progresses, the triage algorithm pro-
duces nearly 100% occupancy, which is in alignment with the FCFS
approach that fills each available bed with the next presenting patient.
In the final days of CSC, the expected decline in patients requiring
critical care occurs, which signals the end of CSC and the return to con-
ventional care. This pattern is distinct from that produced by the FCFS
strategy which assigns an ICU bed based on availability alone, thus main-
taining nearly 100% occupancy from the start. Despite this, the triage
scheme resulted in significantly improved overall population survival by
assigning beds to the most appropriate patients.
2
Table 2.13 Optimal triage thresholds and survival percent for various patient volumes (99% confidence interval in
parenthesis)
Optimal triage
Total patient Optimal triage days on Average 6-week Survival in treated Total survival (entire
volume T_POD % ventilator (T_DV) bed occupancy (admitted) subgroup patient pool)
5273 (5255–5290) 42.3% 4.90 78.6% (78.3–78.8%) 94.9% (94.8–95.0%) 78.5% (77.9–79.0%)
6275 (6255–6295) 40.8% 4.64 83.5% (83.3–83.7%) 95% (94.8–95.0%) 71.1% (70.6–71.6%)
7282 (7259–7305) 38.9% 4.12 85.2% (85.0–85.4%) 95.1% (95–95.2%) 65.8% (65.4–66.2%)
8271 (8242–8301) 35.5% 3.72 85.7% (85.5–85.9%) 95.35% 62% (61.6–62.4%)
(95.3–95.4%)
10,293 27.3% 2.73 88.0% (87.8–88.2%) 95.8% (95.7–95.9%) 54.7% (54.3–55.0%)
(10,265–10,322)
Capacity and Patient Flow Problems: Solutions Using Queuing…
91
92
A. Kolker
Table 2.14 Survival percent for first come-first served (FCFS) admission (99% confidence interval in parenthesis)
Total patient Optimal triage Optimal triage days Average 6-week Survival in treated Total survival (entire
volume T_POD % on ventilator (T_DV) bed occupancy (admitted) subgroup patient pool)
5273 (5255–5290) 42.3% 4.90 83.1% 91.1% (91.0–91.2%) 73.1% (70.9–73.7%)
(82.8–83.3%)
6275 (6255–6295) 40.8% 4.64 87.4% 91.1% (90.9–91.2%) 64.9% (64.4–65.4%)
(87.2–87.6%)
7282 (7259–7305) 38.9% 4.12 90.7% 91.1% (91.0–91.2%) 58.2% (57.8–58.7%)
(90.6–90.9%)
8271 (8242–8301) 35.5% 3.72 93.3% 91.2% (91.1–91.3%) 53% (52.6–53.4%)
(93.1–93.5%)
10,293 27.3% 2.73 96.5% 91.1% (91.0–91.2%) 44.1% (43.9–44.3%)
(10,265–10,322) (96.4–96.6%)
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 93
Model Validation
To make sure that mortality with no thresholds (FCFS) is always greater
than that for admitted patients with the optimal thresholds for the given
patient volume, the following steps were performed:
2. For each patient in this random sample, the actual POD % and DV
were calculated.
3. The optimal thresholds were applied and the patients who did not
meet both thresholds (they were not admitted) were excluded.
4. Percent mortality of admitted patients from this sample was calculated.
Takeaway
This simulation modeling approach presents a rational framework to
guide resource allocation in the event of an overwhelming pandemic
developed with actual pediatric critical care data. In this example, the
discrete event simulation model included also a nonlinear optimization
procedure. The non-linear optimization problem arises if the objective
function is nonlinearly related to the decision variables such as T_POD
and T_DV in this case. In contrast to the problems that can be treated
using linear optimization methodology (presented below in Chaps. 5 and
6), a nonlinear optimization based on the evolutionary algorithm is
essentially a heuristic methodology. It cannot guarantee finding the global
optimum rather than a local one. To increase the chance of finding the
global optimum, several searches should be performed using several dif-
ferent starting points within the decision variable constraints. The cost of
this strategy is, of course, an increase in computational time. Nonetheless,
combining DES with nonlinear optimization is a powerful methodology
worth considering for some complex optimization problems.
2 Capacity and Patient Flow Problems: Solutions Using Queuing… 95
References
AHA. (2011). Trend watch. Examining the drivers of readmissions and reducing
unnecessary readmissions for better patient care. American Hospital Association
(AHA). September.
Cho, S. J., Jeong, J., Han, S., Yeom, S., Park, S. W., Kim, H. H., et al. (2011).
Decreased emergency department length of stay by application of a comput-
erized consultation management system. Academic Emergency Medicine,
18(4), 398–402. [Link]
Christian, M. D., Toltzis, P., & Kanter, R. K. (2011). Task force for Pediatric
emergency mass critical care: Treatment and triage recommendations for
pediatric emergency mass critical care. Pediatric Critical Care Medicine, 12,
S109–S119.
Costa, A., Ridley, S., Shahani, A., Harper, P., De Senna, V., & Nielsen,
M. (2003). Mathematical modeling and simulation for planning critical care
capacity. Anesthesia, 58, 320–327.
De Bruin, A., van Rossum, A., Visser, M., & Koole, G. (2007). Modeling the
emergency cardiac in-patient flow: An application of queuing theory. Health
Care Management Science, 10, 125–137.
DHHS. (2005). Department of health and human services. Health and Human
Services Pandemic Influenza Plan.
Dias, L., Vieira, A., Pereira, G., & Oliveira, J. (2016). Discrete event simulation
software ranking – A top list of the worldwide most popular and used tools.
In Winter simulation conference (WSC). IEEE Xplore. [Link]
WSC.2016.7822165
Forbus, J. J., & Berleant, D. (2022). Discrete-event simulation in healthcare
settings: A review. Modelling, 3, 417–433. [Link]
modelling3040027
Gall, C., Wetzel, R., Kolker, A., Kanter, R., & Toltzis, P. (2016). Pediatric triage
in a severe pandemic: Maximizing survival by establishing triage thresholds.
Critical Care Medicine, 44(9), 1762–1768. [Link]
nal/Abstract/2016/09000/Pediatric_Triage_in_a_Severe_Pandemic_
[Link]
Gallivan, S., Utley, M., Treasure, T., & Valencia, O. (2002). Booked inpatient
admissions and hospital capacity: Mathematical modeling study. British
Medical Journal, 324, 280–282.
Goldratt, E., & Cox, J. (2004). The goal (3rd ed., p. 384). North River Press.
96 A. Kolker
Toltzis, P., Soto-Campos, G., Kuhn, E., & Wetzel, R. A. (2014). A pediatric
scheme to guide resource allocation in a mass casualty. Society of Critical Care
Medicine, 41, A148.
Wang, J. (2012). Reducing length of stay in emergency department: A simula-
tion study at a community hospital. IEEE Transactions on Systems, Man, and
Cybernetics. Part A: Systems and Humans, 42(6), C1–C1309.
Weber, D. O. (2006). Queue fever: Part 1 and Part 2. May 10. Hospitals &
Health Networks, Health Forum. [Link]
j s p / a r t i c l e d i s p l a y. j s p ? d c r p a t h = H H N M A G / Pu b s Ne w s A r t i c l e /
data/2006May/060509HHN_Online_Weber&domain=HHNMAG
Wullink, G., Van Houdenhoven, M., Hans, E., van Oostrum, J., van der Lans,
M., & Kazemier, G. (2007). Closing emergency operating rooms improves
efficiency. Journal of Medical Systems, 31, 543–546.
3
Staffing and Scheduling Problems
with Random Patient Demand:
Solutions Using Discrete Event
Simulation
Problem Description
An outpatient clinic has to schedule two groups of patient appointments:
one group is new patient appointments and another group is follow-up
patient appointments. On average, a new patient appointment takes
about 60 min but could be in the range from 45 to 90 min. The follow-
up appointment takes an average of 30 min and its duration is usually
from 25 to 35 min.
On a typical day, five new and six follow-up appointments are sched-
uled. Using the average appointment time of 60 and 30 min accordingly,
it is estimated that the clinic’s total work time will be
5 × 60 + 6 × 30 = 480 min, or 8 h.
The following managerial problem arises: does appointment order
affect patient wait time and clinic daily total work time?
In other words, is it better to schedule new appointments first and
then the follow-up appointments, or the other way around? Or does the
appointment order make a difference at all?
Indeed, it looks self-evident that the total clinic time for the average
60-min appointment followed by a 30-min appointment would be the
same as that for a 30-min appointment followed by a 60-min appointment.
Table 3.1 Simulation of scheduling rules for one clinic operational day
95% CI of the 95% CI of the number 95% CI of the
average patient wait of patients with NO total clinic time
Scheduling rule time (min) wait (h)
Smallest 6.7–8.5 2.9–3.6 9.6–10.1
variability
first
Random order 10.3–13.2 2.5–3.1 10.3–10.9
Largest 18.9–23.8 1.5–2.02 11.8–12.9
variability
first
102 A. Kolker
Table 3.2 Summary of 25 days of clinic operations using IHI Monte Carlo
simulation
Standard
Standard Average deviation of
Average deviation of number of number of
Scheduling wait time wait time patients with patients with NO
rule (min) (min) NO wait wait
Smallest 6.3 4.4 3.2 2
variability
first
Random 12.2 12.1 2.5 1.9
order
Largest 19.2 17.4 1.4 1.6
variability
first
3 Staffing and Scheduling Problems with Random Patient… 103
Problem Description
The inability to discharge patients in a timely fashion is a typical hospital
problem.
Let’s consider four nursing units. Discharges occur daily Monday to
Friday, usually in the afternoon between 2 pm and 9 pm. The number of
discharges from each unit is a random quantity with the average four
daily discharges, i.e., the number of discharges on a particular day could
be less or more than four. The most likely time to complete discharge was
about 30 min but could be in the range of 20 min to as high as 60 min.
So, on average about 2 h of nursing time was spent on the discharge pro-
cess for each of 4 units. This nursing time was taken away from direct
patient care. The nursing management wanted nurses to spend more time
on direct patient care rather than making the discharge paperwork. What
could be done about that?
104 A. Kolker
(i) The unit nurse would be free from a paper-intensive discharge pro-
cess to maximize their time for direct patient care.
(ii) Bed management could rely on one contact person who would
notify cleaning services to clean beds after discharge.
(iii) Case management would get a resource to help coordinate and plan
their activities.
(iv) Patient flow and throughput would improve because of a more
timely discharge process.
Scenario 1 Current process: there are four independent units. Each unit
has its own nurse who handles discharges as ordered. Four random dis-
charges take place daily in the afternoon from 2 pm to 9 pm. A model
layout is depicted in Fig. 3.1 (SEM 3.2). Scenario 2 Proposed discharge
process: one dedicated nurse who performs discharges for all units as
ordered. The nurse’s shift duration is 8 h, from 1 pm to 9:30 pm with a
30-min lunch break, as indicated on the panel layout in Fig. 3.2 (SEM
3.3). Simulation results are summarized in Table 3.3.
This is a rather unexpected result: a newly hired dedicated discharge
nurse would make significantly fewer discharges than nurses in the cur-
rent process creating an additional discharge backlog and increasing
patient wait time.
Why does this seemingly good management solution actually create a
worse problem compared to the existing discharge process?
3 Staffing and Scheduling Problems with Random Patient… 105
Fig. 3.1 Simulation model layout of current discharge arrangement. Each nurse
RN_A, RN_B, RN_C, and RN_D discharges patients individually from each sepa-
rate unit
It turns out that in the existing process, four nurses perform discharges
independently. Each nurse has her own path to discharge. If there is a
delay with a particular patient in a particular unit, none of the other
nurses in other units would be impacted.
With the new centralized discharge process, all discharges from all
units form a so-called series of dependent events (Goldratt & Cox, 2004;
Motwani et al., 1996). In a series of dependent events, a random delay
with the discharge of a particular patient would inevitably impact the rest
of the patients. In other words, hiring a dedicated resource would create
a system bottleneck. By definition, the system’s bottleneck is a resource
whose capacity is less than or equal to the demand placed on it (Goldratt
& Cox, 2004).
Based only on averages, the dedicated nurse would be able to perform
80 discharges for a week (4 discharges × 4 units × 5 days). If each dis-
charge takes on average 30 min, then 40 work hours for a week would be
needed to complete them. This would make the capacity of the dedicated
106 A. Kolker
Fig. 3.2 Simulation model layout of the proposed discharge arrangement. One
centralized nurse (central discharge) discharges patients from all units. The cen-
tral discharge nurse shift is from 1 pm to 8:30 pm with a 30-min lunch time (red
band) indicated on the shift editor panel
discharge nurse equal to the demand placed on her/him, i.e., would cre-
ate a bottleneck according to the above definition. In reality, the situation
would be even worse because of the demand variability and the nurse
availability that is less than 100% because of additional breaks, huddles,
meetings, etc.
On the “good” days, when the demand for discharges is less than the
nurse’s capacity to perform them (fewer than four discharges per unit that
take less than 30 min per discharge), the extra capacity cannot be “stored”
to serve the next day’s demand. Such an extra capacity on a particular
3 Staffing and Scheduling Problems with Random Patient… 107
“good” day would be lost. On the other hand, on the “bad” days, when
demand for discharges exceeds the nurse’s capacity, the unserved demand
would not be lost; it will have to be fulfilled the next day forming a back-
log (unless there is overtime). Thus, the effect of “bad” days is accumu-
lated, while the effect of “good” days is not. In other words, unsynchronized
variability of demand and capacity causes queues to grow.
This illustrates the following fundamental management engineering
principles: (i) in a series of dependent events, only the bottleneck defines
the throughput of the entire system regardless of the throughput and
capacity of non-bottlenecks and (ii) an unfulfilled service request backlog
(appointments, discharges, document processing, etc.) could remain or
grow even if the average high variability demand is less than the service
capacity (Savin, 2006).
There are possible solutions that can elevate the bottleneck’s constraint
of the centralized discharge nurse, e.g., scheduling easy discharges with
lower variability first (see Sect. 3.1) or the use of the shared staff (as ana-
lyzed in Sect. 3.4).
Problem Description
The hospital receiving center receives annually 127,139 packages. The
time to process one package (screen, scan, store, etc.) could range from 7
to 15 min with the typical time of 10 min. The department works
Monday to Friday from 8 am to 4:30 pm, with a 30-min lunch time and
two 15-min breaks for each staff member during a typical day. (There is a
total of 255 annual workdays.) It is required to develop a staffing plan for
the receiving department with no overtime.
Using the typical time of 10 min to handle a package and the available
daily work time of 7.5 h × 60 = 450 min, each staff member can handle
450/10 = 45 packages per day. Hence, 498.6/45 = 11.1 staff FTE would be
needed to handle the projected work volume without delay and overtime.
Fig. 3.3 Simulation model layout for staffing the hospital receiving center
Problem Description
A hospital case management department performs three types of transac-
tions: reservation, urgent admissions, and pre-registration. It is expected
3 Staffing and Scheduling Problems with Random Patient… 111
that the total annual transaction volume will be 53,855. About 14% of
all transactions are reservation, 32% is urgent admissions, and 52% is
pre-registration. Case management specialists for urgent admissions and
pre-registration are cross-trained and could substitute for each other if
needed. Reservation specialists work independently and they are not
involved in performing other types of transactions.
The department works Monday to Friday from 8 am to 4:30 pm, with
a 30-min lunch time and two 15-min breaks for each staff member dur-
ing a typical day. There are a total of 255 annual workdays. It was esti-
mated that total nonproductive time (unscheduled breaks, distractions,
local building trips, personal time, meetings, etc.) for reservation special-
ists was about 13%, for urgent admission specialists was 15%, and for
pre-registration specialists was about 22%.
Transaction time data were collected over some representative period
and summarized in Table 3.4.
It was needed to develop the minimal staffing (FTE) requirement for
each transaction type that allowed performing the annual transaction vol-
ume without overtime (one FTE-full time equivalent is assumed to
be 8 h).
30
A
P e r ce n t
15
0
0 2 4 6 8 10 12 14 16 18
Reservation time, min
20
B
P e r ce nt
10
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38
Urgent Admission time, min
40
C
P e r ce n t
20
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28
Pre-Registration time, min
Fig. 3.4 Transaction time distribution for reservation staff (a), urgent admission
staff (b), and pre-registration staff (c)
Fig. 3.5 Simulation model layout for cross-trained (shared) staffing of the case
management department. The dotted lines indicate shared (alternate) workload
assignments for urgent admission and pre-registration staff
Table 3.5 Simulation output for the case management department and
required FTE
Simulated annual transaction Required
Transaction type volume, 99% CI FTE
Reservation 7460–7575 0.5
Urgent admission 17,194–17,308 1.5
Pre-registration 29,020–29,153 2.5
Total 53,854–53,855 4.5
Target total annual 53,855
transaction volume
FTE. The reason for this discrepancy is already mentioned in long trans-
action times (tails from distribution time) that are not reflected in the
mean transaction time. Once again, this is one more illustration of the
flaw of averages discussed in detail in the previous sections.
In addition, as it was discussed in Sect. 3.3, it is easy to use the simula-
tion model to find out how much reduction in transaction time
3 Staffing and Scheduling Problems with Random Patient… 115
variability would be needed if the total staffing is strictly budgeted to, say,
3 FTE or 4 FTE. Simulation modeling (discrete event simulation mod-
els) is indeed indispensable to answering such type of questions (or other
questions mentioned at the end of the previous section).
Problem Description
An outpatient flu clinic is open during the flu season to provide flu vac-
cine shots on a walk-in basis. The clinic stays open from 8 am to 6 pm.
Giving a shot (including filling out the paperwork) takes on average
about 8 min but could be in the range from 6 to 10 min. Medical provid-
ers (staff) have unpaid 30-min lunch time and two additional paid
15-min short breaks (not overlapping with each other). The total sched-
uled daily time off for each medical provider is 1 h. The government pays
the clinic $30 for each flu shot. The staffing pay rate is $28.6/h.
The average weekly patient volumes are collected in Table 3.6.
Thus, patient arrival rate is highly variable not only during a typical
day but also during days of the week. Because the Table values are aver-
ages, actual arrivals were modeled as discrete random integer values from
the Poisson distribution with the Table averages.
For instance, the Excel Poisson random numbers generator provided
for the average value 17 the whole values as high as 26 and as low as 8
with various whole values in between. Or for the average value 45 the
actual whole values could be as high as 70 and as low as 30, and so on.
Table 3.6 The average weekly patient volumes for the outpatient clinic
Monday Tuesday Wednesday Thursday Friday Total
8 am to 10 am 17 11 20 20 12 80
10 am to 2 pm 44 36 46 62 52 240
2 pm to 4 pm 18 17 20 16 19 90
4 pm to 6 pm 36 32 37 45 40 190
Total 115 96 123 143 123 600
116 A. Kolker
Thus, the actual total weekly volume could be either higher or lower
than the nominal 600. The clinic’s management should decide how many
medical providers are needed to staff the clinic on a typical day and what
will the projected net revenue be weekly. The net revenue is defined here as
the difference between the gross revenue generated by served patients and
the operational staffing and supplies costs required to provide this service.
The model design and layout are a combination of the models described
in Sects. 2.5, 2.6, and 2.7 with patients leaving (reneging) after waiting
more than 20 min. Of course, any other numbers and input information
can be used (SEM 3.6). The simulation model layout is depicted in
Fig. 3.6.
A baseline (Scenario 1) simulation input data included one provider
working from 8 am to 5 pm (1 full-time equivalent (FTE)) and another
one starting later from 9 am until the end of the day at 6 pm (1 FTE).
On top of scheduled 30-min time off and two 15-min breaks, providers
are assumed to spend additionally about 5% of their time (20–25 min)
for personal needs, i.e., their practical availability is about 95%.
Simulation output for the 95% confidence interval (CI) for the weekly
number of served patients was only 474–477. This is much lower than
the expected average value of 600. On top of that, 142–145 patients
(about 23%) left weekly without a shot because of waiting longer than
20 min. The 95% CI for the weekly net revenue was from $11,830 to
118 A. Kolker
$11,899. This is much less than the expected average value of $15,569.
Such a significant decrease in the net revenue resulted from the inevitable
process variability as well as many patients left without the shot. Thus, in
contrast to the above traditional management estimates, staffing of only
two providers was not enough to meet the clinic’s performance targets.
The next step might be to find out how much it would help to increase
the net revenue and reduce the number of leaving patients if both provid-
ers work, for example, extended hours, from 8 am to 6 pm (Scenario 2).
From a traditional management standpoint, this should not be needed
because 16 h of working time on average should be enough to meet the
average patient demand for service time. Nonetheless, the simulation of
this scenario indicated that the 95% CI for the number of served patients
increased to 520–523, and 96–99 (~15–16%) patients left without a
shot. The 95% CI for net revenue has become markedly higher than that
for regular hours, $12,734–12,819. An additional work hour for each
provider (overtime) is paid at a 50% higher rate. This resulted in addi-
tional operating expenses (increased labor costs). However, these costs
would be amply offset by the clinic’s higher net revenue because more
revenue-generating patients would be served.
Despite the improved clinic’s performance compared to the scenario
with the regular hours, it is still sometimes desired to avoid overtime.
Therefore, another option could be, for example, hiring a part-time pro-
vider. Suppose, for example, that an additional 0.6 FTE with a total of
30 min paid time off is placed in the morning shift from 8 am to 1 pm
(Scenario 3). How much will it help?
Simulation of this scenario indicated that the number of served patients
will be lower than in Scenario 2 with the regular staff overtime: 95% CI
is only 507–509, and 95% CI net revenue is $12,117–12,170.
Does this mean that more than 0.6 FTE is needed? Not necessarily.
What if the same 0.6 FTE is placed in the second shift, from 1 pm to
6 pm? (Scenario 4). This scenario resulted in 95% CI for the number of
served patients, i.e., 558–561, only 57–60 (9–10%) patients left without
a shot and 95% CI for net revenue increases to $13,596–13,684. This
performance is much better than that for the baseline and Scenario 2.
Thus, simulation modeling demonstrates that a third provider (0.6
FTE) placed in the right shift does help to serve more patients and
3 Staffing and Scheduling Problems with Random Patient… 119
increases the net revenue despite the higher costs of keeping an additional
provider. Nonetheless, this is still lower than the expected nominal aver-
age number of served patients (600) and the estimated net revenue of
$15,569.
Of course, many other scenarios of staffing shifts and the clinic’s opera-
tion modes are possible to analyze using this simulation model. The clin-
ic’s manager might be interested, for example, in knowing what staffing
is needed to serve all 600 nominal weekly patients. Simulation modeling
indicates that practically all weekly patient volume will be served using
3.6 FTE (three full-time FTE and a part-time 0.6 FTE in the afternoon).
This scenario gives the 95% CI for the number of patients 608–610, and
9–11 patients left (~2%). The 95% CI for net revenue has become
$13,877–13,936. This is better than all of the above scenarios so far. Even
more patients are served than the nominal 600.
A graphical summary of all simulated scenarios is presented in Fig. 3.7.
It follows from Fig. 3.7 that the highest weekly net revenue (~$13,900)
is generated by 3.6 FTE (1 provider 8 am to 5 pm, 2 providers 9 am to
6 pm, and a part-time provider 0.6 FTE in the afternoon 1 pm to 6 pm).
Patient volume served was 608–610 with 9–11 patients left because of
3.6 FTE
2.6 FTE:
13900 0.6 FTE 1 pm - 6 pm
13700
weekly net revenue, $
13500
13300
2.25 FTE:
13100 3 FTE 4 FTE
2 hours
12900
12700
12500
12300
12100 2 FTE 2.6 FTE:
11900 0.6 FTE 8 am - 1 pm
11700
470 480 490 500 510 520 530 540 550 560 570 580 590 600 610 620
Fig. 3.7 Weekly net revenue vs. the number of served patients and correspond-
ing FTE for 7 simulated scenarios
120 A. Kolker
waiting more than 20 min. The cost of an additional 1.6 FTE (vs. base-
line 2 FTE) was offset by the additional revenue generated by serving
more patients.
The next best scenario in terms of net revenue (~$13,640) is 2.6 FTE
with 0.6 FTE in the afternoon. However, in this scenario, only about 560
patients were served. On the other hand, if there is a shortage of providers
or the staffing budget is strictly limited, then 2.6 FTE could be a via-
ble option.
One more scenario was simulated with 4 FTE (1 provider works 8 am
to 5 pm and 3 providers work 9 am to 6 pm). In this scenario, the num-
ber of served patients was ~604–606. It is not statistically significantly
different from that for 3.6 FTE. However, the net revenue (~$12,860)
was much lower than that for 3.6 FTE because in this case the cost of an
additional 0.4 FTE was not offset by the additional revenue from serving
about the same patient volume. A similar picture is observed comparing
3 FTE and 2.6 FTE: the cost of an additional 0.4 FTE was not offset by
the revenue from serving a few more patients, the value that is statistically
indistinguishable from that for 2.6 FTE.
Even the best net revenue (~$13,900) was still much lower than the
traditional management estimate of $15,569. No scenario demonstrated
that the latter value could be reached even serving all patient volumes.
The net revenue overestimation was a surprise for the clinic’s manage-
ment. The reason for such a miscalculation was a wrong (too low) staffing
level estimate ignoring the patient volume and service time variabilities
and, hence, miscalculated (too low) labor cost. The use of a lower staffing
level (lower labor costs) resulted in serving fewer patients and, hence,
resulted in lower generated revenue. On the other hand, to serve more
patients and, hence, generate more revenue, a higher staffing level is
needed. But this, in turn, increases the labor cost.
This analysis illustrates an important general trade-off between the
number of served patients, the cost of resources to serve them, and the
net revenue. The net revenue increases if the revenue growth from serving
more patients offsets the growth of staffing and supply costs (resources)
to serve them. However, at some point the growth of the resource costs
exceeds the growth of the revenue generated by serving only a few more
or the same number of patients; hence, the net revenue goes down. The
3 Staffing and Scheduling Problems with Random Patient… 121
right staffing that balances the number of served patients, the cost of
resources to serve them, and the net revenue could be found only using
simulation modeling.
Problem Description
This section presents an example of consulting work that the author of
this book had performed some time ago. This example is instructive
because it illustrates a typical request for analysis written using fuzzy lan-
guage by a healthcare manager who is an analytic layman. This is a rather
typical example of a problem in a healthcare setting that requires, first,
some interpretation and making reasonable assumptions. Below is the
text of the original request:
This problem description should have been first converted into a more
concise and specific form that would make it amenable to quantitative
analysis.
Below is a more concise problem description summary:
Thus, 120 patients cannot be seen for 4 h of clinic time even with five
doctors, unless significant clinic overtime is used. However, clinic over-
time should be avoided. Apparently, 120 patients could be seen without
overtime for 1 week (5 days) rather than daily. In this case, the average
number of daily patients would be about 120/5 ~ 24. This is consistent
with the above estimate.
Fig. 3.8 A baseline (an ideal case) simulation model layout that includes two
preliminary exam doctors and three final exam doctors
3 Staffing and Scheduling Problems with Random Patient… 125
• At the end of the simulation cycle, the 95% confidence interval of the
percentage of regular patients discharged within 2 h was ~43.7–50%.
• The average wait time for regular patients was about 85–92 min.
• Doctors were very busy: preliminary doctors’ non-stop utilization was
about 96–97%, while exam doctors’ utilization was about 88–89%.
• The number of the patients that were seen (both regular and compli-
cated) was in the range of only 16–17.
This is in line with the above estimate of the number of daily seen
patients that ignored the time variability. The target of 80% or more
regular patients discharged within 2 h was not even closely met.
Next, several possible staffing scenarios should be analyzed to find out
the staffing level that would meet the operation performance targets. The
first approach would be the use of trial and error to find out the required
staffing level. Another approach would be the use of the automated opti-
mization module Sim_Runner (similar to the one used in Sect. 2.15) to
maximize the patient volume objective function given the constraint that
more than 80% of regular patients should be discharged within 2 h. In
this case, the trial and error approach was simpler and it would save total
simulation time.
The first trial and error scenario was the use of three nurses instead of
two preliminary doctors and three exam doctors. The model is exactly the
same as the baseline except for making preliminary exam capacity equal
to three to accommodate the concurrent work of three nurses. The out-
put of this model was:
• At the end of the simulation cycle, the 95% confidence interval of the
percentage of regular patients discharged within 2 h was ~39–47%.
• The average wait time for regular patients was about 89–97 min;
126 A. Kolker
Thus, adding one more nurse for a short (4–6 min) preliminary exam
does not make a statistically significant difference compared to the base
model. Apparently, a bottleneck is doctors performing long final exams.
Now, let us test a scenario with one additional full-time (FT) and three
part-time (PT) doctors working a shorter shift from 9:30 am to 1 pm.
The staffing set is now 3 nurses + 4 final exam doctors (FT) + 3 PT doc-
tors (9:30 am to 1 pm). This scenario output was:
• At the end of the simulation cycle, the 95% confidence interval of the
percentage of regular patients discharged within 2 h increased almost
by a factor of 2 to ~76–80%, i.e., almost on the desired target of 80%.
• The average wait time for regular patients dropped down to about
49–52 min.
• The number of the patients that were seen (both regular and compli-
cated) was in the range of only 40–41.
Thus, adding full and part-time doctors does make a difference com-
pared to the base model. Apparently, the doctor’s bottleneck was elevated
in this scenario.
It is possible to simulate many other scenarios using this type of model
to find out the most satisfactory staffing level and shifts for meeting the
operation performance targets. Among such scenarios could be, for
example, (i) the effect of additional part-time exam doctors scheduled for
different shifts with various availability; (ii) restriction of the total patient
volume; (iii) changing the percentage of regular and complicated patients;
(iv) various patient arrival patterns; (v) extension of clinic hours or allow-
ing overtime; or (vi) combinations of the above.
3 Staffing and Scheduling Problems with Random Patient… 127
References
Ben-Arieh, D., & Wu, C.-H. (2011). Reducing patient waiting time at an
ambulatory surgical center. Chapter 12. In A. Kolker & P. Story (Eds.),
Management engineering for effective healthcare delivery: Principles and applica-
tions (pp. 1–512). IGI-Global. [Link]
Cayirli, T., Veral, E., & Rosen, H. (2006). Designing appointment scheduling
systems for ambulatory care services. Health Care Management
Science, 9, 47–58.
Costa, A., Ridley, S., Shahani, A., Harper, P., De Senna, V., & Nielsen,
M. (2003). Mathematical modeling and simulation for planning critical care
capacity. Anesthesia, 58, 320–327.
Goldratt, E., & Cox, J. (2004). The goal (3rd ed., p. 384). North River Press.
IHI. (2011). Institute for Healthcare Improvement (IHI). [Link]
knowledge/Pages/Tools/[Link]
Klassen, K. J., & Rohleder, T. R. (1996). Scheduling outpatient appointment in
a dynamic environment. Journal of Operations Management, 14, 83–101.
Motwani, J., Klein, D., & Harowitz, R. (1996). The theory of constraints in
services: Part 2-examples from health care. Managing Service Quality,
6(2), 30–34.
Savage, S. (2009). The flaw of averages (p. 392). Wiley.
Savin, S. (2006). Managing patient appointments in primary care. Chapter 5. In
R. W. Hall (Ed.), Patient flow: Reducing delay in healthcare delivery
(pp. 123–150). Springer.
Teow, K. L. (2009). Practical operations research applications for healthcare
managers. Annals Academy of Medicine Singapore, 38(6), 564–566.
4
Staffing Problems with Random Patient
Demand: Solutions Using Analytic
Optimization Techniques
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2024 129
A. Kolker, Healthcare Management Engineering In Action, Business Guides on the Go,
[Link]
130 A. Kolker
Problem Description
The current demand for more efficient use of healthcare resources, cost
reduction, and patient safety improvement stimulates the development
of novel operations management solutions. Workforce management has
emerged as the one with the biggest impact on cost efficiency and quality
of care in the USA. Indeed, labor cost typically absorbs about 54% of a
total hospital’s operating revenue (Herman, 2013). Therefore, an accu-
rate assessment of the required staffing that matches the highly variable
patient demand is an integral part of the hospital’s budgeting and plan-
ning process. However, random fluctuations in patient demand present
staffing challenges to many hospitals. Given the dynamic nature of the
healthcare supply and demand, the inevitable variability in patient
demand creates two types of problems: (i) over-staffing, which results in
some idle time and/or pay under contractual obligation for no patient
care activity which hurts operation margins, and (ii) under-staffing,
which requires overtime and/or premium pay for call-in from the nursing
extra staffing pool or external staffing agency that was not budgeted ahead
that also hurts operating margins and causes lower quality of care. The
latter problem affects also patients’ and staff’s satisfaction.
This section provides an overview and examples of the application of
the data analytics methodology called the “newsvendor” framework. This
methodology helps to determine the optimal staffing solutions for the
specified time periods for hospital units with randomly fluctuating daily
patient census.
4 Staffing Problems with Random Patient Demand: Solutions… 131
F S Cu / Cu Co (4.1)
Fig. 4.1 Annual daily midnight census and the corresponding staffing
but are kept fixed here for simplicity; other nursing units can
have different acuity percentages and PNRs). Each census value
was multiplied by the weighted sum of inverted PNR values to con-
vert it to the corresponding staffing, e.g., staffing = cen-
sus × (0.28/1 + 0.67/2 + 0.05/3). Results are presented in Fig. 4.1
(upper panel). Note that dividing each census value by the average
PNR = (0.28 × 1 + 0.67 × 2 + 0.05 × 3) = 1.77 would not be correct
because dividing by the average is a nonlinear operation that produces
a biased staffing value. (Dividing census by PNR would be correct if
PNR is an exact number, not the average value).
2. The use of Patient Classification System (PCS): In general, PCS is a
system that quantifies categories of care in order to estimate the
required nursing hours for direct patient care (Malloch & Meisel,
2013). PCS-generated hours of care for each patient are summed up,
and the sum is then divided by the shift length hours. This gives the
required staffing for the particular period. An example of the variable
staffing generated by the PCS is presented, for comparison, in Fig. 4.1
136 A. Kolker
for the same daily midnight census used for the above PNR example
(middle panel).
3. The use of the staffing grid. A staffing grid reflects the perceived nurs-
ing workload by the nursing management. It is usually budget target-
based and serves as a baseline guide. The skills and competency of the
staff are also a consideration. An example of a simplified staffing grid
for a medical-surgical unit is presented in Table 4.1.
The Hourly Costs of Under- and Overstaffing per Nurse and the
Cost Ratio
If too few nurses are planned and staffed vs. the actual nursing demand
(understaffing case), then the additional nurses can be called in from the
internal float pool at no extra cost (if a trained nurse is available), or from
the off-duty nursing pool, or a staffing agency at a premium above the
base pay rate. To illustrate, suppose that a trained nurse is available in the
internal float pool 50% of the time. Hence, a nurse from a staffing agency
should be called on 50% of the time at a premium pay 150% above the
base pay. Thus, the additional cost per called-on nurse over the base pay
rate P ($/h) is Cu = (1–0.5) × 1.5 × P = 0.75 × P ($/h). This cost can be
somewhat underestimated because the float nurses are usually less effi-
cient than the staff nurses.
If too many nurses are staffed (overstaffing), then the extra nurses can
be called off either into the internal float nursing pool or offered to take
paid/unpaid vacation or put on call. If the nurse is put on call, then she/
he does not usually get paid (in the majority of US hospitals). This makes
an impression that calling off a nurse does not cost anything (except in
the case of a contractual obligation to pay for putting on call, if any).
Hence, the unit always has a financial incentive to staff to the maximum.
While this might be true for the unit’s base payroll, there is a fixed cost
for the called-off nurse. KPMG reported (KPMG, 2011) that the average
nursing hourly cost in the USA is about P = $45/h. This includes the base
payroll cost of $35/h, and the total of insurance, recruiting, training, and
other costs per nurse of about $10–11/h, on average. (Currently, the base
payroll and the fixed nursing costs are higher. However, their ratio
remains about the same as it was reported by the KPMG.) Therefore,
even if the called-off nurse is not get paid through the unit’s payroll, she/
he still costs the hospital about 22–24% of the base pay, i.e., 0.24 × P
($/h), but brings no patient care value. Using this fixed cost of overstaff-
ing, the cost ratio in this example is r = 0.75/(0.75 + 0.24) = 0.76. Notice
that the base pay rate is canceled out in this ratio.
In most cases, the understaffing costs more than the overstaffing.
Therefore, the ratio r is usually greater than 0.5. In some cases the under-
and overstaffing costs are the same, making r = 0.5.
Note that an approach was suggested in which the cost parameters Co
and Cu were imputed from the observed actual staffing rather than
treated as the model inputs (Olivares et al., 2008). However, this approach
assumes that the actual staffing was already close to the optimal.
Otherwise, the imputed cost parameters would be significantly biased.
This ECDF curve can be directly used for graphically finding the optimal
staffing solution.
However, if it is desired to build a numerical procedure for staffing
software to find a more accurate optimal solution using an Excel spread-
sheet or any coding language, such as R, Python, or .NET, then the
numerical conversion of staffing into ECDF is needed. The following
generic algorithm can be used.
If a random set of n data points S1, S2, …, Sn contains m data points
that are less than or equal to some value Sk, then ECDF Fk(S) is defined
as the fraction m/n. The ECDF is calculated using the rule:
0, if Sk Smin
Fk S {m / n, if Smin Sk Smax , k 1, 2,
1, if Sk Smax
(i) Make the FTE bins from minimal to maximal staffing values with
the increment of 0.25 FTE. (From a practical standpoint, it is usu-
ally enough to get the staffing solution accuracy ±0.125 FTE. Hence,
a 0.25 FTE staffing increment is usually sufficient.)
4 Staffing Problems with Random Patient Demand: Solutions… 139
Fig. 4.2 ECDF for variable staffing. The dotted lines illustrate the graphical solu-
tion of Eq. (4.1)
(ii) For each incremented FTE value S, test if the corresponding ECDF
value F(S) < r. If this condition is true, then go to the next value of
F(S), and so on.
(iii) If F(S) ≥ r then stop. The optimal staffing will approximately be the
mean value
Example 4.1 Annual and monthly optimal nursing staffing (St. John
hospital, medical-surgical unit). The daily variable required nursing
staffing for the medical-surgical unit is presented in Fig. 4.3a, as well as
the corresponding ECDF Fig. 4.3b.
Using the unit’s available financial data, the extra cost of understaffing
was $30/h/nurse, while the overstaffing per nurse was assumed to be the
fixed rate of $11/h. This makes r = 0.73. The staffing bins and ECDF
values are presented in Table 4.2 (it is abbreviated to save space).
It is seen that r = 0.73 lies in the interval between ECDF values 0.715
and 0.759. Thus, the annual optimal nursing staffing level is the mean of
the corresponding FTE bins, i.e., (6.5 + 6.75)/2 = 6.625 FTE as shown
in the third column in Table 4.2. The graphical solution is indicated by
the dotted lines in Fig. 4.3b.
Fig. 4.3 The daily nursing staffing demand for the medical-surgical unit (a) and
the corresponding ECDF (b). The dotted line on (b) is the graphical optimal solu-
tion ~6.6 FTE for r = 0.73
4 Staffing Problems with Random Patient Demand: Solutions… 141
Table 4.2 FTE bins, ECDF, and the optimal annual nursing staffing level
Staffing FTE bins ECDF, F(S) Optimal staffing, FTE
6.25 0.674
6.50 0.715
6.75 0.759 6.625
7.00 0.775
7.25 0.811
Fig. 4.4 The optimal monthly staffing and the annual optimal staffing level (flat
solid line)
The total costs of under- and overstaffing for the monthly optimal
staffing FTE were compared to that for the average staffing FTE in
Table 4.3.
Example 4.2 The optimal staffing for caregivers’ skill mix. Previous
examples demonstrated the optimal staffing calculation only for one cat-
egory of caregivers—the nurses. However, most nursing units include
other staff skills mix, such as nursing assistants (NA), clinical nurse spe-
cialists (CNS), and nurse practitioners (NP), among others.
The typical cost of under- and overstaffing for RNs in this unit was
similar to the one used in the previous example. However, the costs for
NA and CNS were not readily available. It was assumed that the costs of
4 Staffing Problems with Random Patient Demand: Solutions… 143
Table 4.3 Monthly optimal and average FTE and the total costs of under- and
overstaffing
Total costs of under- and overstaffing
Optimal Monthly For optimal For average Cost saving
Month FTE average FTE FTE ($) FTE ($) (%)
Jan 6.6 5.1 7,766 9,697 20
Feb 8.1 6.5 5,200 7,486 31
Mar 6.6 5.7 4,988 7,094 30
Apr 4.9 4.6 5,093 5,893 14
May 6.9 5.3 7,870 9,887 20
Jun 6.9 6.2 7,165 8,167 12
Jul 6.4 5.0 5,104 7,232 29
Aug 7.1 5.5 8,924 10,670 16
Sep 5.9 4.7 5,545 7,477 26
Oct 4.4 3.4 6,511 7,195 10
Nov 5.4 4.0 5,087 6,749 25
Dec 9.9 8.0 13,625 15,351 11
Table 4.4 Standard hours of care per patient for various levels of patient acuity
Acuity level Standard hours of care
Level 03 2.8
Level 04 3.4
Level 05 4.0
Level 06 4.6
Level 07 5.5
Level 08 6.6
Admissions 0.5
Discharges 0.5
under- and overstaffing for NA were about the same ($13–15/h) result-
ing in the NA cost ratio r = 0.5. The understaffing costs for CNS were
assumed twice the cost of overstaffing resulting in r = 0.67. At the same
time, the calculated optimal staffing values for these caregiver categories
turned out not too sensitive to the r values. Due to the particular staffing
demand and, hence, the shape of ECDF, the NA optimal staffing is not
changed with r values in the range of r = 0.43–0.56. CNS optimal staff-
ing is not changed with r values in an even wider range of r = 0.56–0.75.
Thus, some variations of the poorly known costs were not too critical for
the calculated optimal staffing values presented in Table 4.6, which also
includes the total costs of under- and overstaffing for three categories of
Table 4.5 The number of patients at various levels of acuity, admissions, and discharges (ADT) and the daily required staffing FTE for
three categories of caregivers—NA, RN, and CNS
NA CNS
share of RN share share
total of total of total
Admissions Total hours of hours of hours NA RN CNS
& discharges hours care (28 care of care staffing, staffing, staffing,
Date Level03 Level04 Level 05 Level 06 Level07 (ADT) of care %) (57%) (15%) FTE FTE FTE
1/1/12 1 16 3 93.9 26.8 53.7 13.4 3.3 6.7 1.7
1/2/12 1 17 4 3 107.0 30.5 61.2 15.3 3.8 7.6 1.9
1/3/12 3 11 7 6 107.3 30.6 61.3 15.3 3.8 7.7 1.9
1/4/12 1 21 5 7 134.9 38.5 77.1 19.3 4.8 9.6 2.4
1/5/12 2 22 3 11 133.9 38.2 76.5 19.1 4.8 9.6 2.4
1/6/12 23 5 4 137.0 39.1 78.3 19.6 4.9 9.8 2.4
1/7/12 23 4 7 134.5 38.3 76.8 19.2 4.8 9.6 2.4
1/8/12 1 18 4 7 115.0 32.8 65.7 16.4 4.1 8.2 2.1
1/9/12 18 4 3 107.6 30.7 61.5 15.4 3.8 7.7 1.9
1/10/12 14 9 2 116.0 33.1 66.3 16.6 4.1 8.3 2.1
1/11/12 23 6 7 145.5 41.5 83.2 20.8 5.2 10.4 2.6
1/12/12 20 3 1 109.2 31.1 62.4 15.6 3.9 7.8 2.0
1/13/12 1 21 4 7 129.3 36.9 73.9 18.5 4.6 9.2 2.3
1/14/12 1 20 4 2 119.8 34.1 68.4 17.1 4.3 8.6 2.1
1/15/12 17 3 94.4 26.9 54.0 13.5 3.4 6.7 1.7
1/16/12 5 11 1 5 81.1 23.1 46.3 11.6 2.9 5.8 1.4
1/17/12 1 25 1 5 129.0 36.8 73.7 18.4 4.6 9.2 2.3
1/18/12 1 26 2 8 142.1 40.5 81.2 20.3 5.1 10.2 2.5
1/19/12 1 28 3 9 157.2 44.8 89.8 22.5 5.6 11.2 2.8
1/20/12 1 14 5 7 102.3 29.2 58.4 14.6 3.6 7.3 1.8
1/21/12 14 7 4 106.9 30.5 61.1 15.3 3.8 7.6 1.9
4 Staffing Problems with Random Patient Demand: Solutions… 145
Table 4.6 Optimal and average staffing FTE for three caregiver categories and
the total costs of under- and overstaffing
Total cost of under- and
overstaffing
Optimal Average for average
Caregiver category FTE FTE for optimal FTE FTE
Nursing assistant (NA) 4.4 4.3 $20,705 $25,681
Registered nurse (RN) 9.9 8.6 $49,701 $75,245
Certified nursing 2.4 2.2 $8,414 $15,146
specialist (CNS)
caregivers both for the optimal FTE and for the suboptimal average
FTE. It is seen that the costs of the former are much lower than that of
the latter, similar to the results demonstrated in Example 4.1.
In summary, these examples demonstrated that the optimal staffing
levels vary with the time periods for which they are calculated (annual,
monthly, or any other). They also demonstrated that the mathematically
proved minimal total costs of under and overstaffing for the optimal staff-
ing values are indeed less than that for any suboptimal staffing levels such
as the average. It is recommended that the newsvendor framework is
widely applied to justify a staffing plan whenever it is being developed
and budgeted for.
At the same time, some issues should be further resolved within the
“newsvendor” framework. One is determining the optimal staffing for
the future time periods with some projected gradual non-random trend
in demand. This trend should be incorporated into the current ECDF to
account for some shifts in the future staffing demand distribution.
Another problem is related to the true costs of occurrences of under-
and overstaffing. The called-in float or agency nurses are usually new to
the unit; thus, they are less productive than the regular staff nurses.
Similarly, calling off nurses by putting them on call or floating into the
pool after scheduling up for work usually negatively impacts their morale
and productivity. These non-monetary factors could result in underesti-
mating the true costs of under- and overstaffing. It seems that this prob-
lem has not been adequately addressed yet. Moreover, the classic
“newsvendor” framework assumes that the costs of under- and overstaff-
ing are constant within the time period. However, these costs could be
146 A. Kolker
Problem Description
In the large PACU (post-anesthesia care unit), patients’ census and nurs-
ing hourly and daily workload are highly variable. Census often changes
fast from the minimal value to the peak value within an hour or two. The
required adequate number of nurses to care for the volume of patients
entering the PACU from the operating rooms (OR) is not known ahead
of time. The anesthesiologist and the OR nurses are required to take care
of the patient until a PACU nurse becomes available. This, in turn, results
in (i) delays in OR because anesthesiologist and OR nurse are not
4 Staffing Problems with Random Patient Demand: Solutions… 147
Fig. 4.5 Layout of the simulation model for calculating PACU census and its
change when an admission or discharge takes place
An example of the census values outputs for 1-h time slots is given in
Table 4.7. This table includes only three 1-h time slots to save space,
and it was much shortened for illustration purposes. To filter in the
census values from the entire time series column “Census” into the
corresponding time slot columns, say, 7:30–8:30, the Excel function
was used if(and(‘time of day’>=7.5, ‘time of day’<8.5), ‘Census’, “”). For
the next time slot, the values 7.5 and 8.5 should be changed according
to the boundaries of that time slot indicated in this table as 8.5–9.5,
9.5–10,5, and so on up to 22.5–23.5 and 23.5–24.0. To filter in any
other time slots, the corresponding boundaries can easily be changed,
e.g., for 30 min slots as 7.5–8, 8–8.5, and so on.
Step 4. The number of nurses N that maximized the OF(N) for each 1-h
time slot was taken as the optimal staffing along with the correspond-
ing percent of patient coverage and nursing usage. (Staffing was also
calculated for each 30-min time slot from 7:30 to 22:30 using the same
method. Results are not shown in this section because of limited space.)
The record of PACU annual patient volume included 10,011 admissions
and 10,011 discharges from Monday week 1 to Monday week 53.
Input data file into the DES model included week number, day of
week, and time of day for each admission and discharge.
4 Staffing Problems with Random Patient Demand: Solutions… 151
Table 4.7 An example of the output census values for 1-h time slots
Time of day (h) 7:30–8:30 8:30–9:30 9:30–10:30
Census (in decimals) (7.5–8.5) (8.5–9.5) (9.5–10.5)
1 7.90 1
2 8.25 2
1 8.50 1
2 8.70 2
3 8.87 3
4 8.88 4
5 9.05 5
4 9.05 4
3 9.12 3
4 9.22 4
5 9.33 5
4 9.47 4
3 9.67 3
2 9.73 2
1 9.83 1
0 10.07 0
1 10.13 1
2 10.47 2
3 10.48 3
A short sample of this record (out of total records for 53 weeks) is pre-
sented in Table 4.8. The entire file is presented in SEM 4.2.
152 A. Kolker
The DES model layout for calculating census each time moment it
changes is presented in Fig. 4.5. This is an example of the DES model
that is actually used only for counting the number of admissions and
discharges at specific moments. The simulation length was 8760 h
(52 weeks × 168 h + 24 h for one more Monday on week 53). There were
no random inputs here, so only 1 replication was enough to run. The
DES model is presented in SEM 4.3.
The model used two inputs: admissions and discharges (week, day of
the week, and time of the day) with a total of 20,022 entries. The “attach”
routing made the admitted patient wait for the corresponding time from
the discharge input and then left the PACU system at this time. The cen-
sus change counting occurred at this moment. The output contained a
total of 20,022 census time series data points in the range from 0 to 12.
An example of the simulated PACU census is presented for the first
5 days of week 1 (to save space) in Fig. 4.6. Each group of peaks repre-
sents the times of day when the census changed, from day 1 to day 5
(Monday to Friday), up to 120 h (5 days) with no patients in PACU at
night time and on Saturday and Sunday.
9.00
Week 1; Day 1 Week1; Day 2 Week 1; Day 3 Week 1; Day 4 Week 1; Day 5
8.00
7.00
C
6.00
e
n 5.00
s 4.00
u 3.00
s
2.00
1.00
0.00
0.00 12.00 24.00 36.00 48.00 60.00 72.00 84.00 96.00 108.00 120.00
Time, hours
Fig. 4.6 Simulated 1-week census: day 1 to day 5 (Monday to Friday), 7:30 am to
11:30 pm
4 Staffing Problems with Random Patient Demand: Solutions… 153
0.610
0.605
Weighted OF(N)
0.600
0.595
0.590
0.585
0.580
5 6 7 8 9 10 11 12
Number of nurses, N
Fig. 4.7 An example of the linear weighted sum objective function for various
numbers of nurses for the time slot 16:30–17:30. The maximum corresponds to the
optimal staffing N = 7 for this time slot
12
11
10
9
8
7
Census
6
5 Optimal staffing
4 Current staffing
3
2
1
0
Fig. 4.8 The optimal daily nursing staffing and the currently used staffing for
each 1-h time slot from 7:30 to 24:00. Note: the national standard requires that a
minimum of two nurses be present at all times when a patient (even only one) is
in the PACU. Thus, all optimal staffing values of one must be increased to two
Table 4.9 PACU Staffing plan for seasonal variability and different days of
the week
Monday January to May June to August September to December
7:30–9:30 2 1a 2
9:30–11:30 6 8 5
11:30–13:30 6 8 7
13:30–15:30 7 7 5
15:30–17:30 6 7 7
17:30–19:30 3 6 5
19:30–21:30 2 2 2
Tuesday to
Thursday
7:30–9:30 5 5 5
9:30–11:30 7 8 8
11:30–13:30 8 8 8
13:30–15:30 11 9 9
15:30–17:30 6 8 7
17:30–19:30 3 4 4
19:30–21:30 1a 3 1a
Friday
7:30–9:30 4 6 7
9:30–11:30 10 8 7
11:30–13:30 7 9 9
13:30–15:30 6 9 8
15:30–17:30 7 6 9
17:30–19:30 4 8 5
19:30–21:30 3 2 2
a
Note: To comply with the ASPAN guidelines, staff value of one should be
treated as two
This methodology and the PACU staffing plans were implemented for
planning surgical services at the Children’s Hospital.
Problem Description
The US Center for Disease Control and Prevention (CDC) has revised its
recommendations for screening for human immunodeficiency virus
(HIV) and now recommends HIV screening for all patients aged 13–64 in
all healthcare settings, including hospital emergency departments, urgent
care clinics, inpatient services, sexually transmitted disease clinics, tuber-
culosis clinics, and primary care offices (Armstrong & Taege, 2007;
Bozzette, 2005).
A large testing laboratory is staffed and equipped with a testing capac-
ity of 60 HIV specimens per day. Due to new CDC recommendations,
the specimen daily volume has increased to about 100 per day. This
results in a testing backlog and frequent staff overtime.
The management is eager to increase the testing capacity to reduce the
backlog and staff overtime.
1 1 P n n 1 1 P n
N 1 1 P 1 / n
n
P 1 n 1/ n
158 A. Kolker
The maximum of the right-hand side for this inequality for integer
numbers is about 0.306 for n = 3. Therefore, a reduction of the number
of tests per specimen is theoretically possible only if the probability of a
positive specimen (disease prevalence), P, is less than about 30%. For
each P in this range, N has its minimal value. If 12.4% < P < 30.6%, then
the optimal batch size which minimizes N is exactly n = 3. If P < 11.1%,
then the optimal batch size which minimizes N can be approximated by
the formula (rounded to the nearest integer number)
1
n 0.5
P
References
Armstrong, W., & Taege, A. (2007). HIV screening for all: The new standard of
care. Cleveland Clinic Journal of Medicine, 74(4), 297–301.
ASPAN. (2004). Standards of perianesthesia nursing practice. ASPAN.
Bozzette, S. (2005). Routine screening for HIV infection-timely and cost effec-
tive. New England Journal of Medicine, 352, 620–621.
CDC Report. (2008). HIV prevalence estimates-USA, 2006. MMWR Weekly,
57(39), 1073–1076.
Choi, T. M. (Ed.). (2012). Handbook of newsvendor problems: Models, extensions
and applications. Springer.
Dexter, F. (2007). Why calculating PACU staffing is so hard and why/how oper-
ations research specialists can help. Journal of Perianesthesia Nursing,
22(5), 357–359.
Herman, B. (2013, December 10). Ten statistics on hospital labor costs as a
percentage of operating revenue. Becker’s Hospital Review. Retrieved from
[Link]
on-hospital-labor-costs-as-a-percentage-of-operating-[Link]
Hopp, W., & Lovejoy, W. (2013). Hospital operations: Principles of high efficiency
health care (p. 623). FT Press.
Kolker, A. (2017). The optimal workforce staffing solutions with random patient
demand. Chapter 322. In Encyclopedia of information science and technology
(4th ed., pp. 3711–3724). IGI-Global.
KPMG. (2011). US hospital nursing labor costs study. Retrieved from http://
[Link]/pdfs/KPMG_2011_Nursing_Labor_CostStudy.pdf
160 A. Kolker
Malloch, K., & Meisel, M. (2013). Patient classification system: State of the sci-
ence. Nurse Leader, 11(6), 35–37. [Link]
Olivares, M., Terwiesch, C., & Cassorla, L. (2008). Structural estimation of the
newsvendor model: An application to reserving operating room time.
Management Science, 54(1), 41–45.
Saraniti, B. (2006). Optimal pooled testing. Health Care Management Science,
9, 143–149.
5
Resource Allocation and Service Line
Problems: Solutions Using Linear
Optimization
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2024 161
A. Kolker, Healthcare Management Engineering In Action, Business Guides on the Go,
[Link]
162 A. Kolker
Problem Description
A hospital is in the process of evaluating the financial viability of three
patient service lines. These service lines require the use of five types of
hospital resources. The average values of these resources (both per patient
and the annual total limits) are presented in Table 5.1. The average net
revenue per patient for each service line is also included in this table.
Zero values mean that this type of resource is not required at all for the
corresponding service line.
The hospital management would like to maximize the total net reve-
nue from these services by deciding which services should be offered and
in what volumes (a somewhat similar problem was considered by
Ozcan, 2009).
Table 5.1 Resources required for service lines and total annual limits
Service line Service line Service line Total annual
Resource type 1 MDC-2 2 MDC-19 3 MDC-21 limit
Length of stay 3 days/ 6 days/ 4 days/ 19,000
patient patient patient patient-days
Nursing time 3 h/patient 5 h/patient 4.5 h/ 16,000 h
patient
Interventional 0.5 h/ 1 h/patient 0 4000 h
radiology patient
Laboratory 1 h/patient 2 h/patient 3.5 h/ 6000 h
procedures patient
Operating rooms 2 h/patient 0 4 h/patient 1100 h
Net revenue, $ $560/ $790/ $1100/ Maximize
patient patient patient
166 A. Kolker
net revenue per patient, $1100, followed by volumes for service lines 2
(MDC-19) and 1 (MDC-2). However, it is practically very difficult to
make sure that the total annual resource limits will be met without some
kind of quantitative analysis.
Binding constraint: A constraint that forms the optimal point of the fea-
sible solution, i.e., calculated constraint is equal to the input value.
Slack: A difference of the calculated constraint using the optimal values
of the decision variables and the right-hand side value of the ≤ type
constraint.
Surplus: A difference of the calculated constraint using the optimal values
of the decision variables and the right-hand side value of ≥ type
constraint.
Shadow Prices: How much a one-unit increase in the right-hand side of
a constraint would increase the value of the objective function (OF).
Reduced Cost: How much the coefficient for a variable in the objective
function could be varied holding the same optimal solution.
The sensitivity report for the above problem is presented in Table 5.2.
The report includes two types of data for (i) variable cells and (ii)
constraints.
Table 5.2 Sensitivity report for three service lines’ problem
Variable cells
Cell Name Final Reduced Objective Allowable Allowable
value cost coefficient increase decrease
$B$21 Solution (Xi) 550 0 560 1E+30 165
MDC-2
$C$21 Solution (Xi) 2725 0 790 330 790
MDC-19
$D$21 Solution (Xi) 0 −612.5 1100 612.5 1E+30
MDC-21
Constraints
Cell Name Final Shadow Constraint Allowable Allowable
value price R.H. side increase decrease
$E$16 LOS RHS 18,000 0 19,000 1E+30 1000
$E$17 Nursing hours 15,275 0 16,000 1E+30 725
RHS
$E$18 Radiology 3000 0 4000 1E+30 1000
procedures
RHS
$E$19 Laboratory 6000 395 6000 290 5450
procedures
RHS
$E$20 Operating 1100 82.5 1100 2900 1100
room RHS
170 A. Kolker
with the $1100 coefficient because the OF is not affected by the term in
which the patient charge is multiplied by zero volume, i.e., by X3 = 0.
Suppose now that we vary the coefficient of X2 in the objective
function.
How will the optimal solution change? The current solution X2 is
2725. The allowable increase is $330 while the allowable decrease is $790.
The optimal solution remains unchanged if the coefficient X2 varies
within the range from $1120 to 0. If the net revenue for patient line 2 is
actually outside this range, say, $1121, then the optimal solution becomes
X1 = 0, X2 = 3000, X3 = 0 with the OF = $3,363,000. Thus, dropping
two service lines X1 and X3, and keeping only line 2 with the patient
charge of $1121 would become financially preferred (ignoring that such
a solution would limit patient access to care). It is easy to verify this result
directly using the model in Fig. 5.1.
Problem Description
An outpatient clinic of the Department of Pediatric Internal Medicine
was staffed with five RNs, three LPNs, and five CNAs. The clinic’s man-
agement was interested in justifying the current staffing level and reduc-
ing, if it is possible, the total labor cost.
CNA (Certified Nurse Assistant) hourly pay rate was $11.5. CNAs
perform jobs that require minimal training, such as getting vital signs,
blood pressure, glucometer/finger sticks, cleaning and restocking the
exam rooms between patients, etc.
LPN (Licensed Practitioner Nurse) hourly pay rate was $14. In addi-
tion to the tasks of a CNA, an LPN performs preparing and giving injec-
tions, monitoring catheters, applying dressings, and other more
complicated tasks.
RN (Registered Nurse) hourly pay rate was $24. In addition to CNA
and LPN tasks, RNs perform more complex tasks, such as preparing and
giving IVs for patients, assessing patients’ current health status, assisting
physicians in performing procedures, etc.
Over 8 months, each nursing task was timed as it took place on the
clinic floor using a digital clock. A tally was then obtained of the fre-
quency of each task and its average duration along with the standard
5 Resource Allocation and Service Line Problems… 173
errors of the average. Using adjustments to account for time for miscel-
laneous tasks, it was determined that the clinic weekly required on aver-
age 134 h of CNA time, 88 h of LPN time, and 138 h of RN time. The
total mean weekly nursing time was 360 h with the standard error of the
mean (SEM) of 17.3 h (= standard deviation/sqrt(N)).
1. Decision variables: the number of FTEs for CNA, LPN, and RN,
respectively.
Let X1 be the number of FTE for CNA, X2 be the number of FTE
for LPN, and X3 be the number of FTE for RN.
2. Objective function is total labor cost using hourly pay rate, $/h.
The objective function OF was the total cost that should be mini-
mized, i.e.
OF 11.5 X1 14 X2 24 X3 MIN
3. Constraints:
(i) An LPN can perform a CNA job but a CNA is not allowed and
cannot perform an LPN’s job.
(ii) Based on the collected data, the minimal total nurse time is 360 h.
(iii) The minimal RN time is 138 h. No CNA or LPN can perform an
RN’s job.
(iv) The maximal RN time is 360 h which potentially covers all the
nursing work.
(v) The maximal CNA time is 134 h. A CNA cannot perform any
other nursing job (neither LPN’s nor RN’s job).
(vi) The maximal LPN time is 222 h (88 + 134) which is the time that
covers his/her own weekly time (88 h) and the CNA time (134 h).
174 A. Kolker
The model setup using Excel solver is presented in Fig. 5.2. (The Excel
file is available in supplemental electronic material SEM 5.2.)
Parameters of the OF and constraints are in cells B15–D15 and B16–
D16, respectively. Solution X1, X2, and X3 are in cells B21–D21. The
objective function is calculated in cell G15 using the Excel function
sumproduct(B15:D15, B$21:D$21) which calculates the sum of consecu-
tive products of the pay rate parameters and solutions X1, X2, and X3.
The left-hand sides (LHS) of constraints are in cells in column E16:E20
taken as solutions X1, X2, and X3 from solution cells B$21:D$21. The
right-hand sides of constraint values are in column G16:G20. The objec-
tive function in $G$15 was set to be minimized by changing variable cells
$B$21:$D$21 subject to the set of the above constraints.
Fig. 5.2 Excel solver model setup. Staffing the various skill mix nursing personnel
in the outpatient clinic
5 Resource Allocation and Service Line Problems… 175
Table 5.3 Base model optimal solution with “ignore integer constraints”
box checked
X1 X2 X3
CNA LPN RN LHS Z
Objective function (Cj) 11.5 14.0 24.0 152.1 OF—total weekly
cost
Weekly nurse hours 134.0 88.0 138.0 9.00 ≥ 9.00 Total # of FTE
CNA 3.35 ≤ 3.35 RN min # of FTE
LPN 2.20 ≤ 5.55 CNA FTE
RN 3.45 ≤ 9.00 LPN FTE
RN 3.5 ≥ 3.45 RN max # of FTE
Solution (Xi)-number 3.35 2.20 3.45
of FTE
The base model solution was obtained with checked boxes and “Make
unconstrained variables non-negative” and “Select a solving method-
Simplex LP” in Solver options. Thus, it produced fractional FTE as a
solution. It is presented in Table 5.3.
If the solution is rounded off, then X1 = 3, X2 = 2, and X3 = 3.
If the additional constraint was added for searching the solution that
should be only integer (B$21:D$21 = int), then no rounding off was
needed. The integer solution is presented in Table 5.4.
Thus, the rounded-off FTE solution for RN is lower, i.e., 3 vs. 4.
However, more constraints become binding, so the nursing resources will
be tighter for the rounded-off solution. More information on the optimal
solution and possible cost reduction comes from a sensitivity report.
However, the sensitivity report is possible only for non-integer solutions
(it will be considered later in this section). The sensitivity report cannot
be generated for the integer solution.
Table 5.4 Base model integer optimal solution with “ignore integer constraints”
box unchecked
X1 X2 X3
CNA LPN RN LHS Z
Objective function (Cj) 11.5 14.0 24.0 158.5 OF—total weekly
cost
Weekly nurse hours 134.0 88.0 138.0 9.00 ≥ 9.00 Total # of FTE
CNA 3.00 ≤ 3.35 RN min # of FTE
LPN 2.00 ≤ 5.55 CNA FTE
RN 4.00 ≤ 9.00 LPN FTE
RN 4.0 ≥ 3.45 RN max # of FTE
Solution (Xi)-number 3.00 2.00 4.00
of FTE
turns out that the FTE optimal solutions for the mean nurse time (360 h)
and lower limit nurse time (326 h) constraints are the same: X1 = 3,
X2 = 2, X3 = 4. However, the X2 solution for the upper limit nurse time
constraint (394 h) is different: X1 = 3, X2 = 3, X3 = 4. Therefore, to be
on the safe side as a hedge against a sudden spike in nurse time, it is better
to accept the optimal solution for the upper limit of the total nurse time.
A summary of the current and optimal staffing is presented in Table 5.5.
Thus, potential cost saving with the optimal staffing was
(219.5 − 172.5)/219.5 ~ 21%.
To summarize, the optimal staffing solution suggests that CNA could
be reduced by 2 FTE (from 5 FTE) because their time demand could be
handled by 3 FTE. Indeed, the constraint for CNA is not binding—there
is a slack in it. Thus, cutting the CNA payroll from 5 FTE down to 3
FTE would save ~$46,000 annually. The optimal solution suggests that
the number of LPN should not be changed. This recommendation was
also supported by the clinic observation. The LPNs on the floor are
5 Resource Allocation and Service Line Problems… 177
always on the move with little time to relax unless there were patient
no-shows.
The optimal solution suggests also reducing the current payroll by 1
RN FTE. Although RNs are the most trained nurses, they are not too
efficient during increased patient demand unless they assume in addition
the role of LPN to cover some lower skill tasks. But this does not fit their
level of expertise and is not cost-efficient given their pay rate. Instead of
RN, hiring an additional LPN could prove more cost-efficient staffing
solution. This scenario can also be considered using the LO model frame-
work. As it often happens, low-level skill (CNA) and high-level skill (RN)
staff are somewhat less versatile than the medium-level skill staff, such as
LPN that could cover most (but not all) of the weekly tasks.
Table 5.6 The sensitivity report for the non-integer optimal solution
Variable cells
Cell Name Final Reduced Objective Allowable Allowable
value cost coefficient increase decrease
$B$21 Solution 3.35 0 11.5 2.5 1E+30
(Xi)-number
of FTE CNA
×1
$C$21 Solution 2.2 0 14 10 2.5
(Xi)-number
of FTE LPN
×2
$D$21 Solution 3.45 0 24 1E+30 10
(Xi)-number
of FTE RN ×3
Constraints
Cell Name Final Shadow Constraint Allowable Allowable
value price R.H. side increase decrease
$E$16 Weekly nurse 9 14 9 3.35 2.2
hours
$E$17 CNA 3.35 −2.5 3.35 2.2 3.35
$E$18 LPN 2.2 0 5.55 1E+30 3.35
$E$19 RN 3.45 0 9 1E+30 5.55
$E$20 RN 3.45 10 3.45 2.2 3.35
The next section of the sensitivity report includes the shadow prices.
For the total nursing hours, it is 14 and the allowable increase is 3.35
while the allowable decrease is 2.2. If we elevate this constraint by 3.35
making total nursing hours 9 + 3.35 = 12.35 h, then we expect an increase
of the OF by 14 × 3.35 = 46.9 making it 152.1 + 46.9 = $199. This could
be directly verified using the model. The optimal solution becomes
(rounded off): X1 = 3, X2 = 6, X3 = 3, i.e., the number of LPN is signifi-
cantly increased, from 2 to 6. This is consistent with an earlier conclusion
of the LPN role’s importance.
The next shadow price value is 10 for RNs with an allowable increase
of 2.2 and an allowable decrease of 3.35. If we decrease RN by, say, 1 FTE
(out of allowable 3.35), then the OF is decreased thereby by 10 making
it $142.1 instead of $152.1. The optimal rounded-off solution will
become X1 = 3, X2 = 3, and X3 = 2. Thus, the number of CNAs and
LPNs remains unchanged when the RN is reduced by one FTE.
5 Resource Allocation and Service Line Problems… 179
Problem Description
A lab is planning to offer some new diagnostic tests. Management assessed
the required staff and analysis time, as well as the materials budget
(monthly) presented in Table 5.7.
The management wants to decide which tests and in what volumes should
be offered by maximizing the total lab’s monthly net revenue.
Fig. 5.3 Excel solver model setup. The diagnostic test lab
Table 5.8 The optimal number of each test that maximizes the total monthly
net revenue
Total
Test type → 1 2 3 4 5 RHS Z
Objective function Cj $8.00 $10.00 $8.00 $7.00 $10.00 $1821
Staff time per test, min 15.0 15.0 15.0 20.0 25.0 3400 ≥ 3400
Testing time per 20.0 40.0 40.0 60.0 45.0 6000 ≥ 6000
specimen, min
Material cost per test, $ $12.0 $15.0 $16.0 $14.0 $14.0 $2700 ≥ 2700
Solution (Xi) 120 55 0 0 31
procedures per month
Table 5.9 The sensitivity report for diagnostic test lab problem
Variable cells
Cell Name Final Reduced Objective Allowable Allowable
value cost coefficient increase decrease
$B$14 Solution (Xi) 119.7 0 8 0.18 1.04
procedures/
month
$C$14 Solution (Xi) 55.2 0 10 0.46 0.33
procedures/
month
$D$14 Solution (Xi) 0 −2.58 8 2.58 1E+30
procedures/
month
$E$14 Solution (Xi) 0 −2.93 7 2.93 1E+30
procedures/
month
$F$14 Solution (Xi) 30.9 0 10 4.5 0.66
procedures/
month
Constraints
Cell Name Final Shadow Constraint Allowable Allowable
value price R.H. side increase decrease
$I$11 ≥Z 3400 0.0516 3400 494.23 400
$I$12 ≥Z 6000 0.0129 6000 1687.8 1427.77
$I$13 ≥Z 2700 0.5806 2700 240 734.28
The allowable increase of the reduced cost is $2.58. Hence, the revenue
should be at least $8 + 2.58, i.e., about $10.6 or greater. Indeed, substi-
tuting $10.6 instead of $8 per test 3 makes its volume 51 rather than
zero. Similarly, test 4 volume can be made nonzero if its net revenue per
test becomes at least 7 + 2.93 = 9.93 or, say, 9.94. This makes its volume
34 rather than zero.
Another option to keep tests 3 or 4 at non-zero volume (if they are
clinically much needed) is sacrificing some total lab’s net revenue. Indeed,
the reduced cost for test 3 is negative, −$2.58. If the minimal desired test
3 volume is needed at the level of, say, 30 per month, then the financial
loss for keeping this volume will be −$2.58 × 30 = $77.4. This could be
verified by tightening the non-negativity constraint X3 ≥ 0 to X3 ≥ 30
making the additional constraint in the cell $d$16 ≥ 30. The solution
with this additional constraint is X1 = 116, X2 = 23, X3 = 30, X4 = 0,
182 A. Kolker
X5 = 35, and the OF value is $1743 which is less than the original
OF = $1821 by the expected value of $77.4. Notice also an additional
drop in volumes for tests 1 and 2, from 120 and 55 to 116 and 23,
respectively.
Next, the management identified some way of getting an additional
$200 of the materials budget. How much more the net revenue would be
expected? This question can be answered using shadow price data in the
section “constraints” of the sensitivity report.
Indeed, the shadow price for materials is 0.581 as shown in the last line
of Table 5.9. Hence, the net revenue will be increased by
0.581 × $200 = $116.2, i.e., becomes $1937.
Another question from the management was if an additional 6 h of
monthly time becomes available, which resource should this time be
applied to staff time or testing time? Staff time has a much higher shadow
price; therefore, additional 6 h should be applied to this resource making
it 3400 + 360 = 3760 min with the additional net revenue
0.0516 × 360 = $18.58.
Problem Description
This section presents another important application of LO methodology,
namely, the optimization of a supplier’s transportation cost over the dis-
tribution network.
A common supply chain distribution network structure is a three-level
system that consists of manufacturing facilities (factories), distribution
centers (DC), and consumers (clients). A transportation company deliv-
ers products from several factory locations, accumulates them at several
distribution centers (DC) located at various distances from the factories,
and delivers the products to several customers located at various distances
from the distribution centers. The company is interested in a methodol-
ogy that helps minimize the total transportation cost over a distribution
5 Resource Allocation and Service Line Problems… 183
Fig. 5.4 A three-level shipping distribution network: two factories, F1 and F2;
three distribution centers, DC1, DC2, and DC; four hospitals, H1, H2, H3, and H4.
Arrowhead lines are allowable distribution shipping routes from F(i) → DC(j) and
from DC(j) → H(k). Excluded routes are F1 → DC3; DC1 → H3 and H4; DC2 → H4;
and DC3 → H1
network while meeting customers’ orders and available factories and dis-
tribution centers’ capacity.
As an example, Fig. 5.4 illustrates the network representing the distri-
bution system that consists of two manufacturing facilities (factories), F1
and F2; three DCs, DC1, DC2 and DC3; and four hospitals, H1, H2,
H3, and H4.
It follows from this figure that, e.g., factory F1 ships to DC1 and DC2
but there is no shipment from F1 to DC3. Similarly, DC1 ships to hos-
pitals H1 and H2 but provides no shipment to hospitals H3 and H4
and so on.
The distances between the corresponding factories, DCs, and hospitals
as well as shipping rates are given in Tables 5.10 and 5.11, respectively.
Table 5.10 includes the column and the row for DC4 and H5. These data
are not used for the network analyzed in this section but will be used for
an expanded network in the following sub-section.
184 A. Kolker
Table 5.10 Distance (in miles) between factories and DCs and DCs and hospitals in
the distribution network
DC1 DC2 DC3 DC4
Factory 1 500 300 250 155
Factory 2 600 150 450 300
Hospital 1 340 550 240 1560
Hospital 2 120 190 234 1660
Hospital 3 400 430 340 890
Hospital 4 235 560 1900 1540
Hospital 5 678 987 287 1670
Table 5.12 Shipping costs, $/unit, based on inputs from Tables 5.10 and 5.11.
Symbols in parenthesis are the corresponding designations of the shipping costs
routes “from → to”
DC1 DC2 DC3 DC4
Factory 1 $1500 $900 (CF1DC2) $750 (CF1DC3) $698 (CF1DC4)
(CF1DC1)
Factory 2 $3000 $675 (CF2DC2) $1350 $900 (CF2DC4)
(CF2DC1) (CF2DC3)
Hospital 1 Hospital 2 Hospital 3 Hospital 4 Hospital 5
DC1 $1020 $540 (CDC1H2) $1200 $705 (CDC1H4) $3390
(CDC1H1) (CDC1H3) (CDC1H5)
DC2 $2750 $855 (CDC2H2) $1290 $2800 $4935
(CDC2H1) (CDC2H3) (CDC2H4) (CDC2H5)
DC3 $720 (CDC3H1) $702 (CDC3H2) $1020 $9500 $861 (CDC3H5)
(CDC3H3) (CDC3H4)
DC4 $7800 $8300 $4450 $7700 $8350
(CDC4H1) (CDC4H2) (CDC$H3) (CDC4H4) (CDC4H5)
The transportation costs from factories to DCs and from DCs to hos-
pitals based on these input data are presented in Table 5.12. This table
also includes the column and the row for DC4 and H5 which are not
used for the network in this section but will be used for an expanded
network in the following sub-section.
5 Resource Allocation and Service Line Problems… 185
The first constraint means that the total shipment out of each factory
i = 1, 2 cannot exceed their corresponding manufacturing capacity,
respectively. Constraint 2 means that the total shipment out of both fac-
tories to each DC(j) cannot exceed DC(j) capacity. Constraint 3 means
that the total shipment out of DC(j) to customers H(k) over the corre-
sponding allowable routes cannot exceed the corresponding total ship-
ment to DC(j) from factories F(i). Constraint 4 means that shipments
out of the corresponding DC(j) or their combinations must be exactly
equal to the orders placed by the customers H(k).
Constraints 1–3 are of the type “shipment ‘out’ ≤ shipment ‘in’” or
“shipment out ≤ capacity.” Sometimes it is recommended to use, instead,
constraints of the type “shipment ‘out’ = shipment ‘in’” or “shipment
‘out’ = capacity” (Schrage, 1997). It means that all units manufactured in
the factory must be shipped out to DC or, similarly, all units shipped to
DC must leave DC and be shipped out to the customers. Nonetheless, it
seems that these types of “ = ” constraints are too restrictive. Indeed, some
units could be manufactured in the factory but not necessarily shipped in
the current truckload or some units shipped to DC could stay in DC
waiting for an additional order. Therefore, the models presented in this
section include the constraints of the type “out ≤ in” rather than “out = in”.
The Excel solver setup is presented in Fig. 5.5. (The model file is SEM
5.4. Table 2F 4DC 5H).
Cells C33:E34 were reserved for decision variables XF(i)DC(j). Cells
C38:F40 were reserved for decision variables XDC(j)H(k).
The OF was calculated in cell B46 as the sum of two functions:
sumproduct(C22:E23, C33:E34) to account for shipping costs from fac-
tories F(i) to DC(j) and sumproduct(C26:F28, C38:F40) to account for
shipping costs from DC(j) to H(k).
In search of the solution, to force the units to be shipped only over the
allowable routes and keep zeroes for routes that must be excluded (accord-
ing to the network structure in Fig. 5.4), the changing variable cells
ranges were made C33:D33, C34:E34, C38:D38, C39:E39, and
D40:F40 ($ sign is omitted here for brevity).
Alternatively, it was possible to use only two changing variable cell
ranges C33:E34 and C38:F40 but to add the additional constraints such
as E33 = 0, E38:F38 = 0, F39 = 0, C40 = 0 that force the decision
5 Resource Allocation and Service Line Problems… 187
Fig. 5.5 Solver setup and solution for minimal transportation cost over the net-
work of two factories, three distribution centers (DC), and four hospitals with
restricted selected routes according to Fig. 5.4
variables for the excluded routes to be zeros and, hence, not counted in
the calculation the OF that should reflect the network structure.
The optimal number of the shipped units (solution) is presented in
Table 5.13.
The total minimal shipping cost (OF) was $78,445. All constraints
and customer orders were satisfied with this solution according to the
network structure.
As indicated in Table 5.13, the total number of units shipped out of
factory 1 (8 units) was less than its full manufacturing capacity (10 units).
All units manufactured in factory 2 were at capacity (9 units) and were
shipped out to DCs.
The solution exhibits a useful feature that is applied to the solution of
any network problem: if the capacities and requirement coefficients are
integer values, then the solution will also be integer without the addi-
tional integer-forcing variables constraints. This feature allows generating
a sensitivity report that otherwise could not be generated for integer
solutions.
188 A. Kolker
Table 5.13 The optimal number of units shipped over the distribution network
DC1 DC2 DC3 Total units Capacity of
shipped F(i)
from F(i)
to DC(j)
Factory 1 3 5 0 8 10
Factory 2 0 4 5 9 9
Hospital 1 Hospital 2 Hospital 3 Hospital 4 Total
shipment
from DC(j)
DC1 3 0 0 0 3
DC2 0 5 4 0 9
DC3 0 0 0 5 5
Table 5.14 The sensitivity report for the optimal solution of the distribu-
tion network
Final Reduced Objective Allowable Allowable
Cell Name value cost coefficient increase decrease
Variable cells
$C$33 Factory 1 DC1 3 0 1500 1130 285
$D$33 Factory 1 DC2 5 0 900 285 225
$C$34 Factory 2 DC1 0 1725 3000 1E+30 1725
$D$34 Factory 2 DC2 4 0 675 225 1E+30
$E$34 Factory 2 DC3 5 0 1350 1E+30 405
$C$38 DC1 3 0 1020 1130 1E+30
hospital 1
$D$38 DC1 0 285 540 1E+30 285
hospital 2
$C$39 DC2 0 1130 2750 1E+30 1130
hospital 1
$D$39 DC2 5 0 855 285 1E+30
hospital 2
$E$39 DC2 4 0 1290 405 1E+30
hospital 3
$D$40 DC3 0 522 702 1E+30 522
hospital 2
$E$40 DC3 0 405 1020 1E+30 405
hospital 3
$F$40 DC3 5 0 9500 1E+30 1E+30
hospital 4
Constraints
$G$33 Factory 1 8 0 10 1E+30 2
Total
shipped to
DC
$G$34 Factory 2 9 −225 9 5 2
Total
shipped to
DC
to DC2 while one more unit out of F2 to DC2 compared to the solution
shown in Table 5.13. Thus, increasing the capacity of F2 resulted in a
somewhat lower total shipping cost with some rearrangement of the
number of units shipped over the allowable routes. Note that the shadow
price for factory 1 was zero because its capacity constraint was not bind-
ing, and so the capacity change does not affect the OF.
190 A. Kolker
Fig. 5.6 Solver setup and the minimal transportation cost solution over the
expanded distribution network (two factories, four DC, and five hospitals) with
no restrictions on possible routes
reserved for the solution. This was in contrast to the previous network
structure in which some ranges of cells reserved for solution were fixed to
zero to exclude counting the cost of shipping over the corresponding
excluded routes.
The total cost of shipping OF (cell B$46) was calculated similarly to
that in the previous sub-section as the sum of two sumproduct functions
that included the cells to account for DC4 and H5. Constraints were also
similar to those in the previous sub-section including the expanded
ranges of cells to account for DC4 and H5.
The optimal solution was the total minimal shipping cost of $39,536
and the number of the shipped units is shown in Table 5.16.
It is seen from Table 5.16 that although the factories’ manufacturing
capacities remained the same as in the previous sub-section, the number
of units shipped out of F1 using any route is now 10, while only 8 units
were shipped from F2, making total 18. At the same time, the total mini-
mal shipping cost with the use of all possible combinations of routes was
now almost twice as low compared to the previous case with the restricted
network.
192 A. Kolker
Table 5.16 Shipping solution over the expanded distribution network with no
restricted routes
DC1 DC2 DC3 DC4 Total Capacity
shipped of factory
from (i)
factories
(i) to all
DC
Factory 1 5 0 5 0 10 10
Factory 2 0 5 1 2 8 9
Total 5 5 6 2
shipment
from all
factories
to DC(j)
DC(j) 5 5 6 6
capacity
Hospital 1 Hospital 2 Hospital 3 Hospital 4 Hospital 5 Total
shipment
from
DC(j) to
all
hospitals
DC1 0 0 0 5 0 5
DC2 0 5 0 0 0 5
DC3 3 0 2 0 1 6
DC4 0 0 2 0 0 2
Total 3 5 4 5 1 18
shipment
from all
DC(j) to
the H(k)
Hospital 3 5 4 5 1 18
order
Table 5.17 The sensitivity report for the optimal shipping network with unre-
stricted routes
Variable Name Final Reduced Objective Allowable Allowable
cell value cost coefficient increase decrease
$C$33 Factory 1 5 0 1500 900 1E+30
DC1
$D$33 Factory 1 0 825 900 1E+30 825
DC2
$E$33 Factory 1 5 0 750 397.5 900
DC3
$F$33 Factory 1 0 397.5 697.5 1E+30 397.5
DC4
$C$34 Factory 2 0 900 3000 1E+30 900
DC1
$D$34 Factory 2 5 0 675 825 1E+30
DC2
$E$34 Factory 2 1 0 1350 900 397.5
DC3
$F$34 Factory 2 2 0 900 397.5 900
DC4
Constraints Name Final Shadow Constraint Allowable Allowable
cell value price R.H. side increase decrease
$G$33 Total 10 −600 10 1 1
shipped
from
factory 1
to all DC
$G$34 Total 8 0 9 1E+30 1
shipped
from
factory 2
to all DC
the total shipping cost down to $38,936. The number of units shipped by
factory 1 is 11 (equal to capacity), but the number of units shipped by
factory 2 dropped down to 7 making total 18 units shipped.
On the other hand, decreasing the capacity by 1 (allowable decrease)
making it 9 units resulted in the increase of the shipping cost by $600
making it $40,136. In this case all 9 units from both factories were
shipped making the same total 18 units. Thus, changing the capacity
within allowable increase and decrease results in rearrangement of the
number of shipped units keeping total number the same.
194 A. Kolker
References
Kall, P., & Mayer, J. (2005). Stochastic linear programming. Models, theory and
computation (p. 397). Springer Science & Business Media.
Ozcan, Y. (2009). Quantitative methods in health care management (2nd ed.,
p. 438). Jossey-Bass. A Wiley Imprint.
Prekopa, A. (1995). Stochastic programming. Kluwer Academic Publisher.
Schrage, L. (1997). Optimization modeling with LINDO (5th ed.). Brooks/Cole
Publishing.
Sen, S., & Higle, J. (1999). An introductory tutorial on stochastic linear pro-
gramming models. Interfaces, 29(2), 33–61.
6
Staffing and Scheduling Problems
with Sliding Days Off: Solutions Using
Integer Linear Optimization
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2024 195
A. Kolker, Healthcare Management Engineering In Action, Business Guides on the Go,
[Link]
196 A. Kolker
Problem Description
Scheduling staff is a typical but difficult and time-consuming task in
operations management. A typical full-time unit nurse (such as an ICU
nurse) usually works 5 days a week with 2 consecutive mandatory rotat-
ing days off and rotating shifts. In this section, the optimal daily shift
staffing schedule is presented for a typical ICU unit.
Minimal staffing coverage is given in Table 6.1.
The problem was determining the minimal total nursing pool size and
the number of nurses scheduled to work for each day given the minimal
daily coverage.
Table 6.1 Minimal daily staffing coverage for the ICU unit
Day of
week Monday Tuesday Wednesday Thursday Friday Saturday Sunday
Minimal 4 5 3 5 4 3 6
staffing
coverage
6 Staffing and Scheduling Problems with Sliding Days… 197
OF X1 X2 X3 X 4 X5 X6 X7 min
Fig. 6.1 Integer linear optimization model setup. Weekly scheduling 5 workdays
with 2 mandatory sliding days off
In addition, X1, X2, X3, X4, X5, X6, and X7 must be positive and integer.
Step 4. Parameters: There are no parameters in this problem.
The Excel problem setup is presented in Fig. 6.1. (The model file is SEM
6.1; Tab Min staffing pool)
The key to setting up this type of problem is making an index table
that indicates the days the staff can be scheduled to work according to
his/her schedule ID (but not necessarily actually scheduled) and the days
when the staff cannot be scheduled to work. There are a total of seven
schedule IDs in this example numbered 1–7, as indicated in the left-hand
panel on Fig. 6.1.
The index table is presented as the middle panel titled “index:
1-working; 0-not working” on Fig. 6.1. Allowable work days are desig-
nated by index 1, otherwise by 0. For instance, for schedule ID 1 with
days off Saturday and Sunday, the index table shows 1 for Monday to
Friday and 0 for Saturday and Sunday. Similarly, for schedule ID 2 with
Sunday and Monday days off, the index table shows 1 for Tuesday to
Saturday and 0 for Monday and Sunday and so on for all other
schedule IDs.
6 Staffing and Scheduling Problems with Sliding Days… 199
Table 6.2 Breakdown of the number of scheduled staff per day per schedule ID
Solution
Schedule Days minimal # of
ID off staff Mon Tue Wed Thu Fri Sat Sun
1 Sat, 0
Sun
2 Sun, 1 1 1 1 1 1
Mon
3 Mon, 1 1 1 1 1 1
Tue
4 Tue, 0
Wed
5 Wed, 2 2 2 2 2 2
Thu
6 Thu, Fri 0
7 Fri, Sat 3 3 3 3 3 3
Total 7 Total 5 6 5 5 4 4 6
staff daily
→ →
Problem Description
This problem presents an extension of the previous problem presented in
Sect. 6.1.
In many hospital units, it is required to determine (in addition to the
previous constraints) which staff to assign to which schedule if their per-
sonal preferences and seniority factors are taken into account.
6 Staffing and Scheduling Problems with Sliding Days… 201
Table 6.5 The total score as a weighted sum of normalized seniority (weight 0.65)
and normalized preferences (weight 0.35)
Total score: weighted normalized seniority + weighted normalized preference
ID_1 score ID_2 score ID_3 score ID_4 score ID_5 score ID_6 score ID_7 score
1.00 0.70 0.75 0.85 0.90 0.95 0.80
0.64 0.59 0.54 0.69 0.84 0.79 0.74
0.43 0.48 0.68 0.38 0.53 0.58 0.63
0.55 0.30 0.40 0.25 0.50 0.35 0.45
0.72 0.77 0.92 0.67 0.62 0.82 0.87
0.63 0.48 0.43 0.38 0.33 0.58 0.53
0.61 0.76 0.56 0.46 0.51 0.66 0.71
The decision variables were binary variables with the value 1 if the
schedule ID was assigned to the employee and, otherwise, 0 if the par-
ticular schedule ID was not assigned.
Hence, for each employee, the schedule decision score was calculated
as the sumproduct function of the product of the weighted normalized
score (from Table 6.5) and the corresponding value of the binary variable
1 or 0. If the binary variable value is 1, then that schedule ID is desig-
nated as “assigned”; otherwise, it is 0 and is not included in the total
decision score.
The objective function (OF) was maximizing the total decision scores
for all employees as a team. This maximal total decision score was searched
by varying the binary variables for each employee with a trial of each
schedule ID (1 or 0, i.e., “assigned” or “not assigned”) until the maximal
total score value was found.
6 Staffing and Scheduling Problems with Sliding Days… 203
Constraints
Fig. 6.2 Integer linear optimization model setup. Weekly staffing and schedul-
ing with employee’s preferences and seniority
204 A. Kolker
This solution gave all employees 2 consecutive days off every week and
assigned them to either their first, second, or third choice of schedule ID
weighted by their seniority.
Notice that even this rather simple LO problem includes a total of 56
decision variables: the number of staff for 7 possible schedules and 49
binary variables to assign 7 employees to 7 possible schedules. There are
also 63 constraints: 7 constraints for the sum of binary variables equal to
1; 7 constraints for making the actual number of employees equal to the
required number from the optimal solution in Sect. 6.1; and 49 con-
straints for making sure that the variables (0 or 1) to assign 7 employees
to 7 possible schedules are indeed binary.
Sometimes, employees might not be happy with the assigned sched-
uled ID given by the optimal solution. Indeed, Alex might not like
assigning him his third choice ID schedule despite his highest seniority.
Counterintuitive schedule assignments like this happen because the goal
(objective function) was maximizing the total score for the entire team
rather than maximizing personal employee seniority and preference score.
Given the constraints, satisfying personal scores for each team member
was not possible. Therefore, a team player is supposed to accept the out-
come that he/she does not like when it happens.
Nonetheless, one may wish to relax or tighten some constraints to get
another more intuitively satisfying solution. However, changing con-
straints without some idea of the solution’s sensitivity to such a change
could result in an unfeasible solution. (Unfortunately, there is no help
here from the sensitivity report discussed in detail in Chap. 5. As it was
already mentioned, integer LO does not produce such reports.)
To illustrate, let us tighten the daily constraints used in the problem
setup in Sect. 6.1 making them all, say, 7 as indicated in Table 6.6 in the
last row.
6 Staffing and Scheduling Problems with Sliding Days… 205
Fig. 6.3 Weekly staffing and scheduling with employee’s preferences and senior-
ity. Unfeasible solution with elevated daily constraints
Using these tightened constraints and the same model as in Fig. 6.1,
the minimal total staffing pool size becomes 10 rather than 7 with the
schedules ID breakdown: 1 employee with ID_1; 2 with ID_2; 1 with
ID_3; 2 with ID_4; 1 with ID_5; 1 with ID_6; and 2 with ID_7 (total
staffing pool 10).
However, the extension of this model that takes into account seniority
and schedule preference described in this section results in an unfeasible
solution with these new constraints, as shown in Fig. 6.3. Three con-
straints cannot be satisfied (actual number of employees vs. required opti-
mal number for ID_4, ID_5, and ID_6) meaning that there is no schedule
ID assignments scheme with the total staffing pool 10 and the current
206 A. Kolker
staff schedule preference choice and seniority. So, either the weights of
seniority and schedule preference should be changed or some employees’
schedule preferences should be modified and the model runs again to see
if all the problem’s constraints could be satisfied in the new scenario.
Problem Description
The same nursing unit discussed in Sect. 6.1 was exploring the weekly
work schedule for full-time nurses along with some part-time nurses. The
total minimal number of nurses required per day was the same as was
indicated in Table 6.1 in Sect. 6.1. It was assumed that the same staffing
would be required week after week for at least a few months. Full-time
nurses work 5 days with 2 consecutive days off. The management wanted
to explore a minimal staffing cost option for combined full-time and
part-time staff if the part-time nurses would be working only 2 consecu-
tive days a week with the 5 other days of the week off.
The daily full-time nursing pay per shift was $280 (=$35/h × 8 h shift),
while the part-time was paid daily the amount of $320 (=$40/h × 8 h).
Fig. 6.4 Excel solver model setup for full- and part-time nursing staffing
PT1, PT2, …, PT7 with only 2 consecutive work days and, hence, 5 days
off corresponding to two values of index 1 and 5 values of zero index in
each row.
Constraints for FT and PT were calculated separately using the sum-
product function and then added up to make the total number of employ-
ees, as indicated in line 22 on Fig. 6.4. Objective function was the total
weekly employees’ pay to be minimized. The pay for FT and PT was
calculated separately for each day and then added up.
The model outcome is presented in Fig. 6.4. The minimal weekly pay
is $8920 (weekly $7000 full-time and $1920 PT pay). The FT pool size
is 5 and the PT pool size is 3 making the total pool size 8. The breakdown
of the number of scheduled staff per day per ID is presented in the right-
hand panel in Fig. 6.4. It follows that two nurses with ID 5 (days off Wed
and Thu) should be scheduled on Monday, Tuesday and Friday, Saturday,
and Sunday. In addition, three nurses with FT7 (Fri, Sat days off) should
be scheduled each day except their days off Friday and Saturday.
On top of that, one part-time nurse (with schedule ID PT1) should be
scheduled on Saturday and Sunday, and two part-time nurses with sched-
ule ID PT6 should be scheduled on Thursday and Friday.
208 A. Kolker
Table 6.7 Full-time and part-time schedule with FT daily pay of $320 and PT daily
pay of $280
Schedule Days
ID off Solution Mon Tue Wed Thu Fri Sat Sun
FT7 Fri, Sat 1 1 1 1 1 1
PT1 Mon– 3 3 3
Fri
PT2 Tue– 2 2 2
Sat
PT3 Wed– 1 1 1
Sun
PT4 Thu– 3 3 3
Mon
PT5 Fri–Tue 0
PT6 Sat– 4 4 4
Wed
PT7 Sun– 0
Thu
Total FT → 1 Total 4 5 4 5 4 3 6
pool size FT + PT
per day →
Total PT → 13
pool size
Next, let us explore the effect of the daily pay if it increased for FT
nurses to $320 ($40/h) while the daily pay for part-time decreased to
$280 ($35/h). In this case, the situation dramatically changed (the full
Excel panel is not shown here to save space): most of the full-time staff
was essentially priced out of the schedule while the minimal total weekly
pay became lower than in the above case, $8880. Only 1 FT nurse with
schedule ID FT7 is left to be scheduled each day except Friday and
Saturday, while the rest of the staff becomes part-timers with the PT pool
size 13 to cover missing FT staff. These schedules are presented in
Table 6.7.
Full-time staff costs more, so the management may prefer part-timers.
Because part-time employees’ pay is lower and they don’t receive benefits,
employers may view them as more valuable than FT. This case can be
compared to the case considered in Chap. 5 in Sect. 5.2 when CNAs are
essentially priced out if their pay rate exceeds $14/h.
6 Staffing and Scheduling Problems with Sliding Days… 209
Problem Description
A typical clinical unit has some minimal staffing requirements based on
the average shift patient census and assumed nurse-to-patient ratio. This
ratio is based on assessed patient acuity level or external regulations.
It is assumed here that the pay rate is $35/h (base wages and overhead)
with a 50% pay rate increase for Saturday and Sunday shifts. Usually,
three 8-h shifts per day should be covered.
For a 7-day week and three shifts per day, there are a total of 21 differ-
ent schedules possible. These schedules are presented in Table 6.8 along
210
Table 6.8 Binary index variable Is,ds for s = 21 schedules and days and shifts: 1, on shift; 0, off shift
Days Mon_ Mon_ Mon_ Tue_ Tue_ Tue_ Wed_ Wed_ Wed_ Thu_ Thu_ Thu_
Schedule, s off shift_1 shift_2 shift_3 shift_1 shift_2 shift_3 shift_1 shift_2 shift_3 shift_1 shift_2 shift_3
1 Sat, 1 0 0 1 0 0 1 0 0 1 0 0
Sun
2 Sat, 0 1 0 0 1 0 0 1 0 0 1 0
A. Kolker
Sun
3 Sat, 0 0 1 0 0 1 0 0 1 0 0 1
Sun
4 Sun, 0 0 0 1 0 0 1 0 0 1 0 0
Mon
5 Sun, 0 0 0 0 1 0 0 1 0 0 1 0
Mon
6 Sun, 0 0 0 0 0 1 0 0 1 0 0 1
Mon
7 Mon, 0 0 0 0 0 0 1 0 0 1 0 0
Tue
8 Mon, 0 0 0 0 0 0 0 1 0 0 1 0
Tue
9 Mon, 0 0 0 0 0 0 0 0 1 0 0 1
Tue
10 Tue, 1 0 0 0 0 0 0 0 0 1 0 0
Wed
11 Tue, 0 1 0 0 0 0 0 0 0 0 1 0
Wed
12 Tue, 0 0 1 0 0 0 0 0 0 0 0 1
Wed
13 Wed, 1 0 0 1 0 0 0 0 0 0 0 0
Thu
14 Wed, 0 1 0 0 1 0 0 0 0 0 0 0
Thu
15 Wed, 0 0 1 0 0 1 0 0 0 0 0 0
6
Thu
16 Thu, 1 0 0 1 0 0 1 0 0 0 0 0
Fri
17 Thu, 0 1 0 0 1 0 0 1 0 0 0 0
Fir
18 Thu, 0 0 1 0 0 1 0 0 1 0 0 0
Fri
19 Fri, 1 0 0 1 0 0 1 0 0 1 0 0
Sat
20 Fri, 0 1 0 0 1 0 0 1 0 0 1 0
Sat
21 Fri, 0 0 1 0 0 1 0 0 1 0 0 1
Sat
Average census per 10 24 20 21 18 15 8 15 24 18 20 21
shift
Minimal staff 5 12 10 10 9 7 4 7 12 9 10 10
demand per day
and shift Nds, Min
Demand
with the average shift census and the minimal staff demands for each shift
assuming nurse to patient ratio of 1:2 for all shifts (only Monday to
Thursday are shown in this table. Friday, Saturday, and Sunday are struc-
tured similarly but not shown here to save space).
The management is supposed to develop a staffing schedule to meet
the minimum coverage for each day and shift with 5 workdays and 2
consecutive days off for each staff member in such a way that the total
weekly staffing cost is minimized.
21
N ds X s I s ,ds for each ds 1,, 21
s 1
21 21 21
C N ds Pds Lds Pds Lds X s I s ,ds min
ds 1 ds 1 s 1
21
N ds X s I s ,ds N ds ,min Demand , for each ds 1,, 21
s 1
Fig. 6.5 Excel solver model setup to minimize the total weekly labor cost for
three daily shifts staffing coverage with 2 consecutive days off
The setup structure is similar to that described in Sects. 6.1, 6.2, and
6.3 except that the index table now includes three shifts per each day,
such as Mon_shift_1; Mon_shift_2; Mon_shift_3; Tue_shift_1; Tue_
shift_2; Tue_shift_3 and so on. The objective function to minimize is in
cell L24 and calculated constraints are in cells E26:Y26. Constraints are
calculated using the sumproduct function of solution values times corre-
sponding 1 or 0 values for each day and shift from the Index table.
The optimal solution is presented in Table 6.9 for Monday to Thursday
shifts (Friday to Sunday are not shown to save space). Thus, to imple-
ment the schedule with three daily shifts and rotating consecutive 2 days
21
off, the total required nursing pool size should be N tot X s 36 . The
minimal total weekly cost is $77,800. s 1
Table 6.9 provides information for the number of nurses that should
be assigned on the particular day for the particular shift. For instance, for
shift 1 on Monday, two nurses should be assigned with schedule s = 1 (Sat
and Sun off), one nurse should be assigned with schedule s = 13 (Wed
and Thu off), and two nurses should be assigned with schedule s = 19 (Fri
and Sat off). Thus, a total of five nurses should be assigned for shift 1 on
Monday. Similarly, say, for shift 3 on Wednesday, eight nurses should be
Table 6.9 Optimal solution to scheduling problem (Friday to Sunday are not shown to save space)
Schedule Optimized
s = 1, Days solution, Mon_ Mon_ Mon_ Tue_ Tue_ Tue_ Wed_ Wed_ Wed_ Thu_ Thu_ Thu_
…, 21 off Xs: shift_1 shift_2 shift_3 shift_1 shift_2 shift_3 shift_1 shift_2 shift_3 shift_1 shift_2 shift_3
1 Sat, 2 2 2 2 2
Sun
2 Sat, 7 7 7 7 7
Sun
3 Sat, 8 8 8 8 8
Sun
4 Sun, 5 5 5 5
Mon
5 Sun, 1 1 1 1
Mon
6 Sun, 1 1 1 1
Mon
7 Mon, 1 1 1
Tue
8 Mon, 0 0 0
Tue
9 Mon, 1 1 1
Tue
10 Tue, 0 0 0
Wed
11 Tue, 4 4 4
Wed
12 Tue, 0 0 0
Wed
13 Wed, 1 1 1
Thu
(continued)
Table 6.9 (continued)
Schedule Optimized
s = 1, Days solution, Mon_ Mon_ Mon_ Tue_ Tue_ Tue_ Wed_ Wed_ Wed_ Thu_ Thu_ Thu_
…, 21 off Xs: shift_1 shift_2 shift_3 shift_1 shift_2 shift_3 shift_1 shift_2 shift_3 shift_1 shift_2 shift_3
14 Wed, 1 1 1
Thu
15 Wed, 0 0 0
Thu
16 Thu, 0 0 0 0
Fri
17 Thu, 0 0 0 0
Fri
18 Thu, 2 2 2 2
Fri
19 Fri, 2 2 2 2 2
Sat
20 Fri, 0 0 0 0 0
Sat
21 Fri, 0 0 0 0 0
Sat
Actual staff demand per 5 12 10 10 9 11 10 8 12 10 12 10
shift, Nds
Minimal staff demand per 5 12 10 10 9 7 4 7 12 9 10 10
shift, Nds, Min Demand
6 Staffing and Scheduling Problems with Sliding Days… 217
assigned with schedule s = 3 (Sat and Sun off), one nurse with schedule
s = 6 (Sun and Mon off), one nurse with schedule s = 9 (Mon and Tue
off), and two nurses with schedule s = 18 (Thu and Fri off). The total staff-
ing for shift 3 on Wednesday is 12.
Many other scenarios are possible to analyze using LO methodology.
For instance, nurse-to-patient ratio and pay rate could be different from
shift to shift and day to day, or part-time staff can be added to the
schedule.
At the same time, there is an important caveat for LO model result
interpretation similar to that mentioned in Chap. 5 at the end of Sect.
5.1. Census values are treated here as fixed parameters although these
numbers are just the averages over the variable census described by some
statistical distributions. Therefore, this problem is, strictly speaking, a sto-
chastic linear optimization problem. However, describing and using sto-
chastic linear optimization is beyond the scope of this book. As was
mentioned in the previous section, stochastic linear optimization coverage
was provided, for example, by Prekopa (1995) and Kall and Mayer (2005).
Daytime Scheduling
LO approach to scheduling in this case is somewhat similar to the one
used in the Sects. 6.1, 6.2, 6.3, and 6.4. However, this scheduling prob-
lem differs because the IOM recommendations require actually the devel-
opment of a monthly (4 weeks) schedule rather than a weekly one to
implement at least one monthly 2-day off period (golden weekend) along
with at least 1 weekly day off for all other weeks in this month. These
requirements result in 196 possible schedules! Indeed, it is possible to
make seven schedules with any 2 consecutive days off for any 1 of 4 weeks
for the month along with the other weeks with a single day off, say, on
Monday, i.e., 7 × 4 = 28 schedules for 1 month. This pattern can be
repeated for the weeks with other single days off, i.e., Tuesday, Wednesday,
and so on. Thus, there are a total of 28 × 7 = 196 schedules possible.
Decision variables for this problem are the number of residents, Xs,
assigned to each of s = 196 possible schedules. The objective function is
the total number of residents with all schedules Ntot that should be mini-
mized, i.e.
196
N tot X s min
s 1
Constraints are the minimal total staff demand for each of the 7 days
of each of the 4 weeks (dw = 1, …, 28). Monday week 1 (Mon1) is
indexed as dw = 1, Tuesday week 1 (Tue1) is indexed as dw = 2, and so on
to Sunday week 4 (Sun4) indexed as dw = 28. Thus, Ndw, min = 4 for each
weekday (Monday to Friday for each of 4 weeks), and Ndw, min = 3 for each
of four weekends (Saturday and Sunday for each week). Thus, the opti-
mal solution Xs (s = 1, …, 196) must satisfy the inequalities
196
N dw X s I s ,dw N dw,min , foreach dw 1,, 28
s 1
The values of the binary index variable Is,dw = 0 or Is,dw = 1 are indicated
in Table 6.10. Only the first s = 7 × 6 = 42 schedules for the first 2 weeks
are shown in this table to save space. Mon234 designates Monday off on
Table 6.10 Binary index variable for daytime shift schedule development: 1, on shift; 0, off shift (first 42 schedules for
220
Nighttime Scheduling
The overall LO approach for the development of an optimal night shift
schedule is similar to the one described above. However, there is a com-
plication due to the requirement of the cycle of 2 consecutive days off
after three or four consecutive nights. Total cycle length (on-off) is less
than 7 days/week. Therefore, there are no fixed days off for different
weeks (as it was in the previous daytime schedule development). Assuming
here 4 consecutive work nights, a pattern with 2 consecutive days off will
be sliding, and the total cycle (on-off) is repeated every 6 full weeks. The
binary index variables (on-off) for s = 6 schedules are presented in
Table 6.12.
6 Staffing and Scheduling Problems with Sliding Days… 223
There are only six different schedules in this case that start on week 1.
For instance, schedule 1 starts with 4 consecutive work nights on week 1
(Mon1–Thu1) with 2 consecutive days off (Fri1 and Sat1). Then there is
another string of 4 work nights (Sun1 to Wed2) with another 2 days off
(Thu2 and Fri2), and so on. Schedule 2 starts on week 1 with Mon1 day
224 A. Kolker
Table 6.12 Binary index variables for night time schedule development: 1, on
shift; 0, off shift
Week Schedule, s → 1 2 3 4 5 6
Two sliding
consecutive days-off Fri1, Sat1, Mon1, Tue1, Wed1, Thu1,
↓ starting on week 1 → Sat1 Sun1 Tue1 Wed1 Thu1 Fri1
1 Mon1 1 0 0 1 1 1
1 Tue1 1 1 0 0 1 1
1 Wed1 1 1 1 0 0 1
1 Thu1 1 1 1 1 0 0
1 Fri1 0 1 1 1 1 0
1 Sat1 0 0 1 1 1 1
1 Sun1 1 0 0 1 1 1
2 Mon2 1 1 0 0 1 1
2 Tue2 1 1 1 0 0 1
2 Wed2 1 1 1 1 0 0
2 Thu2 0 1 1 1 1 0
2 Fri2 0 0 1 1 1 1
2 Sat2 1 0 0 1 1 1
2 Sun2 1 1 0 0 1 1
3 Mon3 1 1 1 0 0 1
3 Tue3 1 1 1 1 0 0
3 Wed3 0 1 1 1 1 0
3 Thu3 0 0 1 1 1 1
3 Fri3 1 0 0 1 1 1
3 Sat3 1 1 0 0 1 1
3 Sun3 1 1 1 0 0 1
4 Mon4 1 1 1 1 0 0
4 Tue4 0 1 1 1 1 0
4 Wed4 0 0 1 1 1 1
4 Thu4 1 0 0 1 1 1
4 Fri4 1 1 0 0 1 1
4 Sat4 1 1 1 0 0 1
4 Sun4 1 1 1 1 0 0
5 Mon5 0 1 1 1 1 0
5 Tue5 0 0 1 1 1 1
5 Wed5 1 0 0 1 1 1
5 Thu5 1 1 0 0 1 1
5 Fri5 1 1 1 0 0 1
5 Sat5 1 1 1 1 0 0
5 Sun5 0 1 1 1 1 0
6 Mon6 0 0 1 1 1 1
(continued)
6 Staffing and Scheduling Problems with Sliding Days… 225
Week Schedule, s → 1 2 3 4 5 6
Two sliding
consecutive days-off Fri1, Sat1, Mon1, Tue1, Wed1, Thu1,
↓ starting on week 1 → Sat1 Sun1 Tue1 Wed1 Thu1 Fri1
6 Tue6 1 0 0 1 1 1
6 Wed6 1 1 0 0 1 1
6 Thu6 1 1 1 0 0 1
6 Fri6 1 1 1 1 0 0
6 Sat6 0 1 1 1 1 0
6 Sun6 0 0 1 1 1 1
off, which is consecutive to Sun6 day off because the entire pattern is
repeated every 6 weeks, and so on.
The optimal solution that minimizes the total number of residents
assigned to all schedules subject to a set of constraints (a minimum of
three residents on Mondays to Fridays and a minimum of two residents
on Saturdays and Sundays) was obtained using Microsoft Excel solver.
(The model file is SEM 6.6.)
The optimal solution is given in Table 6.13.
It follows from this table that, for example, three residents should be
scheduled on Monday, week 1 (Mon1): one resident starting on Monday
with Friday and Saturday days off this week (Fri1 and Sat1) and two resi-
dents starting on Monday (Mon1) with Wednesday and Thursday days-
off this week (Wed1 and Thu1). Other schedules are interpreted similarly.
Thus, LO is indeed a very powerful and resourceful methodology for
developing rather complex optimized time coverage schedules.
Table 6.13 Optimal solution to nighttime resident scheduling problem
226
Week Schedule, s → 1 2 3 4 5 6
Two sliding
consecutive
days-off starting Wed1, Thu1, Total
↓ on week 1 → Fri1, Sat1 Sat1, Sun1 Mon1, Tue1 Tue1, Wed1 Thu1 Fri1 residents Constraints
A. Kolker
1 Mon1 1 2 3 3
1 Tue1 1 1 2 4 3
1 Wed1 1 1 1 3 3
1 Thu1 1 1 1 3 3
1 Fri1 1 1 2 4 3
1 Sat1 1 2 3 2
1 Sun1 1 2 3 2
2 Mon2 1 1 2 4 3
2 Tue2 1 1 1 3 3
2 Wed2 1 1 1 3 3
2 Thu2 1 1 2 4 3
2 Fri2 1 2 3 3
2 Sat2 1 2 3 2
2 Sun2 1 1 2 4 2
3 Mon3 1 1 1 3 3
3 Tue3 1 1 1 3 3
3 Wed3 1 1 2 4 3
3 Thu3 1 2 3 3
3 Fri3 1 2 3 3
3 Sat3 1 1 2 4 2
3 Sun3 1 1 1 3 2
4 Mon4 1 1 1 3 3
4 Tue4 1 1 2 4 3
4 Wed4 1 2 3 3
4 Thu4 1 2 3 3
4 Fri4 1 1 2 4 3
4 Sat4 1 1 1 3 2
4 Sun4 1 1 1 3 2
5 Mon5 1 1 2 4 3
5 Tue5 1 2 3 3
5 Wed5 1 2 3 3
5 Thu5 1 1 2 4 3
5 Fri5 1 1 1 3 3
5 Sat5 1 1 1 3 2
5 Sun5 1 1 2 4 2
6 Mon6 1 2 3 3
6 Tue6 1 2 3 3
6 Wed6 1 1 2 4 3
6 Thu6 1 1 1 3 3
6 Fri6 1 1 1 3 3
6 Sat6 1 1 2 4 2
6 Sun6 1 2 3 2
6 Staffing and Scheduling Problems with Sliding Days…
227
228 A. Kolker
References
Blum, A., Shea, S., Czeisler, C., Landrigan, C., & Leape, L. (2011). Implementing
the 2009 IOM recommendations on resident physician work hours, supervi-
sion and safety. White paper. Dovepress. Nature and Science of Sleep, 3, 47–85.
[Link]
Kall, P., & Mayer, J. (2005). Stochastic linear programming. Models, theory and
computation (p. 397). Springer Science & Business Media.
McLaughlin, D., & Olson, J. (2012). Healthcare operations management
(AUPHA) (p. 441). Health Administration Press.
Prekopa, A. (1995). Stochastic programming. Kluwer Academic Publisher.
Ulmer, C., Wolman, D., & Johns, M. (Eds.). (2009). Resident duty hours:
Enhancing sleep, supervision, and safety. Institute of Medicine. National
Academic Press. [Link]
Enhancing-Sleep-Supervision-and-[Link]
7
Comparative Performance of Units and
Organizations with Multiple Inputs
and Outputs: Introduction to Data
Envelopment Analysis
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2024 229
A. Kolker, Healthcare Management Engineering In Action, Business Guides on the Go,
[Link]
230 A. Kolker
Traditional Approach
A commonly used method is calculating ratios. Typically, we take some
output measure and divide it by some input measure. We view units as
taking inputs and converting them into outputs. In this example, there is
a single input measure, i.e., the number of staff members, and a single
output measure, i.e., the number of transactions. Hence, we have:
Unit 1 has the highest ratio of personal transactions per staff member,
whereas unit 4 has the lowest ratio of personal transactions per staff mem-
ber. As unit 1 has the highest ratio of 6.9, we can compare all other units
to it and calculate their relative efficiency with respect to unit 1.
7 Comparative Performance of Units with Multiple Inputs… 231
To do this we divide the ratio for any unit by 6.9. This gives
However, typically there are more than one input and one output that
characterize a system’s performance. Hence, there are several different
ratios used in the traditional approach. One problem with comparison
using ratios is that different ratios can give different performance pic-
tures, and it is difficult to combine the entire set of ratios into a single
numeric score that gives a single picture of relative performance. The
problem of having different ratios giving different pictures would be
especially true for multiple numbers of units and many input/output
measures.
DEA Approach
Data envelopment analysis (DEA) is a method for evaluating the effi-
ciency of a set of peer entities called decision-making units (DMUs). The
definition of a DMU is generic and flexible. It could be any manufactur-
ing or service unit that consumes some input resources to produce some
output products/services. DEA allows multiple inputs and outputs to be
combined into a single efficiency score that characterizes the DMU’s
performance.
The term “envelopment” in DEA comes from the enveloping property
of the most efficient DMUs that form the frontier against which all other
DMUs are compared. The main idea of the method is that a set of inputs
and outputs can be identified that are common across a set of units to be
compared in terms of their efficiency. These inputs and outputs are used
to create a scoring function for each unit relative to other units in the set.
DEA’s empirical orientation and the absence of a need for the a priori
assumptions that accompany other approaches (such as standard forms of
statistical regression analysis) have resulted in its wide use in healthcare,
finance, and manufacturing (Cooper et al., 2011; Zhu, 2002).
232 A. Kolker
where output weights are u1,…,uM and input weights are ν1,…,νN.
The inputs and output (I and O values) are directly observed and serve
as inputs in the DEA model setup. They should not necessarily share a
common metric as long as they are on the same basis for all DMUs. The
fact that inputs and outputs may use different measure units (or be
dimensionless) is one of the great advantages of DEA methodology.
The values of weights (u, and ν values) need to be found as an optimal
solution that maximizes the score for each DMU. If inputs and outputs
bear some dimensions, then the solution weights u and ν have the cor-
responding inverse dimensions, so that each product uMOM and νNIN
remains dimensionless.
Thus, unit k may get a scoring function that makes it as good as pos-
sible (maximized) subject to no other unit getting a score >1. If unit k
gets a score of 1, it means that there is no other unit strictly dominating
that unit k.
234 A. Kolker
1 I1k 2 I 2 k N I Nk 1
In other words, the problem becomes the following: given the fixed
weighted sum of inputs, maximize the weighted sum of outputs with
respect to weights. This is a so-called output-oriented objective. It is pos-
sible to formulate an input-oriented objective with a fixed nominator
(weighted sum of outputs) and minimized weighted sum of inputs.
Details of these two approaches can be found in Cooper et al. (2011). In
this chapter, only the first formulation will be used.
7 Comparative Performance of Units with Multiple Inputs… 235
Problem Description
There are four performance metrics collected for eight hospital units
(DMUs). The first two metrics are treated as inputs: cost per patient (1)
and % of Medicaid patients for each DMU (2). The last two metrics are
treated as outputs: surgical quality score (1) and the length of medically
necessary patient stay-LOS (2) (Table 7.2).
Which units (DMU) can be considered efficient, and which ones are
less efficient?
236 A. Kolker
Which inputs or outputs and by how much does a less efficient unit
need to change to become as efficient as the best units?
DEA Model
Excel solver setup is presented in Fig. 7.1.
In the Excel setup for this problem, there are four decision variables:
two output weights and two input weights. The cells B21:E21 were
reserved for these decision variables. The objective function to be maxi-
mized is in cell C24. Using the INDEX function, it is taken from column
cells F13:F20 which are weighted outputs for each corresponding unit i
from number 1 to 8 in this case. The OF in cells F13:F20 for each unit
was calculated as the sumproduct function of the two corresponding out-
puts for each unit and the decision variables, e.g., =sumproduct(d13:e13,
d$21:e$21) for unit 1, and so on. The formula in cell C24
is = INDEX(F13:F20, C23,1). This function looks at the list of cells
F13:F20 and takes the value in the row indicated in the above cell C23.
If the cell in C23 contains the number 1, this formula returns the value
from the first row in the range F13:F20 or F13. Similarly, if we look at
units 2, 3, 4, etc., this function will return the second, third, fourth entry,
etc. from this range and so on.
The same function is used to make sure that we have the correct con-
straint for the unit under consideration. Thus, in cell C25 we
use = INDEX(H13:H20, C23,1). Whatever unit number is listed in cell
7 Comparative Performance of Units with Multiple Inputs… 237
Fig. 7.1 DEA Excel solver model setup for eight units with two inputs and
two outputs
C23, cell C25 represents the appropriate weighted input for this unit
from the range H13:H20 that must be constrained to be equal to 1.
Additional constraints are the sum of weighted outputs ≤ sum of
weighted inputs for all k = 8 units. These constraints are in cells F13:F20
and H13:H20, respectively, and the relationship indicator “≤” is placed in
column G for clarity. One more set of constraints was added to avoid
possible solution weight values equal to zero. Therefore, the non-
negativity solution constraints were tightened to B21:E21 ≥ ɛ, where ɛ is
a small positive number, such as 0.00001 (the model file is SEM 7.1; Tab
Units 1–8 Index funct).
Calculated efficiency scores for all units are presented in Table 7.3.
Thus, the best-performing unit is #6 with an efficiency score of 1 and
the least-performing unit is #7 with a score of 0.703. What can be done
to improve the efficiency score of unit 7?
One approach is improving the surgical quality score. But by how
much? For instance, let us make the surgical quality score 34 (instead of
24.7). This would make the efficiency score for unit 7 equal to 0.968, i.e.,
return this unit to an almost perfect efficiency score. Is it practically
238 A. Kolker
Problem Description
A company manages 12 medical labs in the region. The company is in the
process of evaluating how efficiently each lab was operated. The data col-
lected are shown below for each of the 12 DMUs in Table 7.5.
The outputs include each lab’s net profit (in $100,000 s), the average
customer satisfaction rating, and the average patients’ health out-
comes score.
The inputs include total labor hours (in 100,000 s) and total operating
costs (in $M).
It is needed to compare the efficiency of these labs (units). What fac-
tors can be leveraged to improve the efficiency of the least-performing unit?
DEA Model
The DEA model setup is presented in Fig. 7.2 (the model file is SEM 7.2;
Tab 12 units).
Similar to the previous model setup, we indicate which DMU is being
evaluated in cell C17.
Cell C18 contains the function that returns the weighed output for
this DMU from the list of weighted outputs in column G, i.e.,
INDEX($G$2:$G$13, C17,1). Maximizing this value for each DMU
from cell C17 is the objective function (OF).
Cell C19 contains the function INDEX($I$2:$I$13, C17,1). This
function looks at the list of cells I2:I13 and takes the value in the row
indicated in C17 and column 1. Whatever DMU number is indicated in
cell C17, cell C19 represents the appropriate weighted input that must be
constrained to equal 1.
Additional constraints are the weighted outputs for all lab units ≤
weighted inputs, e.g.,
for unit 1 sumproduct(D2:F2, D$14:F$14) ≤ sumproduct(B2:C2,
B$14:C$14);
Fig. 7.2 DEA Excel solver model setup for 12 lab units with 2 inputs and
3 outputs
242 A. Kolker
Table 7.7 The sensitivity report and shadow prices lab units for the least-
performing lab 8
Final Shadow
Cell Name value price
$C$20 Input DMU # constraint: input 2 total 1 0.7719
operating costs
$G$3 Lab 1 weighted outputs 0.7206 0
$G$4 Lab 2 weighted outputs 0.8610 0.561
$G$5 Lab 3 weighted outputs 0.6071 0
$G$6 Lab 4 weighted outputs 0.6431 0
$G$7 Lab 5 weighted outputs 0.4243 0
$G$8 Lab 6 weighted outputs 0.6010 0
$G$9 Lab 7 weighted outputs 0.3655 0.299
$G$10 Lab 8 weighted outputs 0.7719 0
$G$11 Lab 9 weighted outputs 0.7229 0
$G$12 Lab 10 weighted outputs 0.8297 0
$G$13 Lab 11 weighted outputs 0.7073 0
$G$14 Lab 12 weighted outputs 0.9945 0.180
244 A. Kolker
Let us test the target value input 1 “labor hours” at 5.47 instead of the
current 7.09. This gives a much better score of 0.842 for lab unit 8. Next,
let’s test input 2 total “operating cost” 6.7 M instead of 8.69 M. The score
is now 0.995, almost perfect. The current outputs for lab unit 8 and the
target outputs are rather close to each other. They can be tested but not
much improvement is expected. We could therefore combine the two
improvement target inputs: 5.47 h and $6.7 M and get the perfect
score of 1.
Of course, any combination of intermediate values could be used if the
improvement target values are too challenging to achieve in practice. The
efficiency score will not be perfect 1 but some improvement will be
achieved.
Problem Description
A network of ten regional hospitals was considered for comparative hos-
pitals’ efficiency.
A regional disparity was becoming an increasingly important growth
constraint. Policymakers needed quantitative knowledge to design effec-
tive and targeted policies. Summary of input and output data for ten
network hospitals is presented in Table 7.8.
When questioning hospital efficiency, it is important to take into
account the quality of the services and patients’ satisfaction. In this case,
five inputs were chosen as hospitals’ characteristics: (i) total number of
7
Table 7.8 Summary of input and output data for ten network hospitals
Input 5: # of Output 1: Output 2:
Input 2 total portable average bed average Output 3:
Input 1 # # of medical Input 3 Input 4 # ultrasound occupancy, nursing outcome
Input/output of beds staff # of CT of MRI devices days time, days quality score
Hosp 1 500 1100 5 4 1149 244 8.0 89
Hosp 2 395 996 6 5 1100 240 6.0 87
Hosp 3 350 993 5 4 800 240 7.2 87
Hosp 4 403 990 4 3 1277 246 8.6 96
Hosp 5 446 990 5 5 1100 200 6.5 75
Hosp 6 444 970 6 3 998 230 10.0 90
Hosp 7 439 987 6 5 1150 250 8.9 84
Hosp 8 390 900 5 3 900 235 9.9 80
Hosp 9 495 1090 5 4 1050 239 7.5 79
Hosp 10 455 980.0 6 4 995.0 236.0 6.8 88
Comparative Performance of Units with Multiple Inputs…
245
246 A. Kolker
beds, (ii) total number of medical staff, (iii) the number of magnetic reso-
nance imaging (MRI) scanners, (iv) the number of computer tomogra-
phy (CT) scanners, and (v) the number of portable ultrasound devices.
Three outputs were chosen, such as (i) the average bed-day occupancy,
(ii) the average annual nursing days, and (iii) patients’ outcomes qual-
ity score.
This choice of inputs and outputs was somewhat justified. For inputs,
the number of beds is one of the indicators that reflect the relative size of
the hospital. It could serve as a proxy for marginal profit to the hospital.
The number of beds is one of the commonly used indicators for compar-
ing a hospital’s size.
The number of the medical staff represents the total number of medi-
cal employees without subdivision into subgroups, such as doctors,
nurses, physician assistants, technicians, etc. This value often serves as a
basic proxy indicator of the patient outcome quality because it is assumed
that the patient outcome quality is directly related to the number of med-
ical staff. Of course, the total number of medical staff directly affects the
total payroll costs. The total medical staff of hospitals as the aggregate
variable—human resources—is frequently used as a factor that affects
total hospital efficiency.
The number of magnetic resonance (MR) and computer tomography
(CT) scanners often serves as a factor of how widely the hospital uses the
advanced high-tech medical technology that is related positively to
patient quality outcomes. These pieces of equipment are used for making
improved diagnoses, and they are the most expensive pieces of medical
technology (priced at $$M). Therefore, the number of this machinery is
counted separately as a possible efficiency and reputation factor.
The number of portable ultrasound devices is also frequently consid-
ered as an indicator of diagnostic advancement and the technology factor
that affects the hospital’s efficiency and reputation.
For outputs, the bed average annual occupancy and nursing time are
used by many experts as indicators that reflect the use of resources avail-
able to the hospital. If these values are too low, it could be a warning sign
for inefficient use of resources and hospital capacities, which could lead
to a reduction in the number of active staffed beds or nurses with unaf-
fected patient outcomes. Finding the optimal size of the bedding and
7 Comparative Performance of Units with Multiple Inputs… 247
DEA Model
DEA model setup and its description is similar to that presented in
Sect. 7.2. Therefore, it is not included here to save space (the model file
is SEM 7.3; Tab 10 Hosp).
However, the use of shadow prices to calculate the improvement tar-
gets is important and is demonstrated here one more time along with
some caveats to keep in mind. Calculated efficiency scores for all ten
hospitals are summarized in Table 7.9.
The least-performing unit is hospital 5 with a score of 0.809. The best-
performing units are hospitals 3, 4, 6, and 8. Therefore, these units can
be used to set up the target for improvement for hospital 5. The
References
Cooper, W., Seiford, L., & Zhu, J. (2011). Data envelopment analysis: History,
models, and interpretations. In Handbook on data envelopment analysis (2nd
ed.). Springer.
Jung, S., Son, J., Kim, C., & Chung, K. (2023). Efficiency measurement using
data envelopment analysis (DEA) in public healthcare: Research trends from
2017 to 2022. PRO, 11, 811. [Link]
Zhu, J. (2002). Quantitative models for performance evaluation and benchmark-
ing: Data envelopment analysis with spreadsheets and DEA excel solver. Kluwer
Academic Publishers.
8
Statistical Data Analytics
for Decision-Making
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2024 251
A. Kolker, Healthcare Management Engineering In Action, Business Guides on the Go,
[Link]
252 A. Kolker
Problem Description
A hospital plans a major market share expansion to improve its long-term
financial viability. The management would like to know what population
demographic variables specific to the local core service area zip codes are
the biggest contributors to the financial contribution margin (CM $).
The contribution margin is defined as the difference between payments
collected from a patient and the patient’s variable cost to the hospital.
A set of population demographic data was collected for 51 local area
zip codes, and corresponding median contribution margins were calcu-
lated for each zip code (CM $).
The following groups of demographic variables were collected for each
zip code as the percentage of the total zip code population. Actual data
are not shown here because of the file size and limited space (the data file
is available in SEM 8.1):
Thus, a total of 32 data variables were collected for each of the 51 zip
codes and included in the demographic database.
PC1 11 X1 21 X 2 p1 X p
PCr 1r X1 2 r X 2 pr X p
Notice also that the sum of all eigenvalues in Table 8.1 (31.997) is
practically equal to the number of the original data variables, 32 (within
the rounding-off error), as it is required by the correct PC decomposition
technique.
The next step is to perform a regression analysis to relate CM $
(response function) to mutually uncorrelated PC variables and then to
determine the original variables (factors) that contribute the most
to CM $.
Regression analysis with PC has substantial advantages over regular
regression analysis with the original data, and it is much more reliable
(Jobson, 1992; Glantz & Slinker, 2001).
Because PCs are mutually uncorrelated, the variation of the dependent
variable (CM $) is accounted for by each component independently of
other components, and their contributions are directly defined by the
coefficients of the regression equation. The presence of any PC does not
affect the regression coefficients of the other PC.
Usually, if the purpose of PCD analysis is just the reduction of the
number of independent variables and the search for the underlying data
pattern, the PCs with eigenvalues much less than one (minor PCs) can be
dropped because they retain little or no information. However, in mul-
tiple regression with PC aimed at determining the contributions of PCs
to the response variable, it is not a good idea to drop even the minor PCs
because doing so can introduce an uncertain bias in the regression coef-
ficients (Jobson, 1992; Glantz & Slinker, 2001). In this particular case,
the first PCs PC1–PC9 accounted for 95.7% of the original data.
8 Statistical Data Analytics for Decision-Making 257
At the same time, only one predictor, PC7, was statistically significant
at the 5% significance level with P-value = 0.003 ≪ 0.05, as shown in
Table 8.2 (highlighted in bold).
Age category: “18–24 years old” (0.566) followed by “60–64 years old”
(0.316) and “65–74 years old” (0.197)
8 Statistical Data Analytics for Decision-Making 259
Table 8.3 Coefficients for the statistically significant principal component PC7
Variable PC7 coefficient Variable PC7 coefficient
Age 18–24 0.566 Some HS 0.041
Age 25–34 −0.056 Income <$15 K 0.046
Age 35–44 −0.47 Income $15–$25 K 0.016
Age 45–54 −0.039 Income $25–$35 K 0.018
Age 55–59 0.217 Income $35–$50 K 0.058
Age 60–64 0.316 Income $50–$75 K 0.015
Age 65–74 0.197 Income $75–$100 K −0.038
Age 75+ −0.119 Income $100–$150 K −0.056
AD degree 0.188 Income $150–$250 K −0.038
BD degree −0.139 Income $250–$500 K −0.01
PhD 0.072 Income $500 K+ −0.029
HS 0.165 Occupation: Healthcare −0.229
Less than HS −0.172 Labor 0.038
MD −0.029 Professional/administrative 0.047
Pro degree −0.107 Public service 0.174
Some college 0.146 Service industry 0.046
Problem Description
A regional community hospital collected contribution margin (CM) data
for all patients from local area zip codes. Median CM (in descending
order), are shown in Fig. 8.1 (actual zip codes are replaced by dummy zip
codes; for figure clarity, only a limited number of zip codes is presented).
The hospital management wants to get an answer to the following
question using this data: is it possible to combine some zip codes into a
few relatively homogeneous groups (clusters) with similar CM to analyze
which factors make one group (cluster) different from others? Or each of
a few dozen zip codes should be analyzed separately to find out why its
CM is different from others?
10000
8000
CM, $
6000
4000
2000
0
3 9 4 3 0 7 8 4 5 7 9 5 1 1 6 2 7 4 5 2 7 0 2 3 5 1 2 6 8
40 02 02 22 22 02 21 22 21 22 21 04 05 21 21 21 21 21 00 02 03 21 07 21 22 22 22 22 22
13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13
ZIP codes
Fig. 8.1 Median contribution margins (CM, $) for some zip codes
8 Statistical Data Analytics for Decision-Making 261
Cluster Analysis
Cluster analysis is used to reduce the number of observations of the vari-
ables (the number of rows of this dataset). Cluster analysis is widely used
to identify the groups of objects (clusters) in such a way that the objects
within one cluster have much higher similarity compared to objects
within another cluster; thus, the total number of observational values of
the variable is replaced by a few clusters of observational values. Clustering
262 A. Kolker
The dendrogram shows how the clusters are formed either by joining two
individual observations or pairing an individual observation with an
existing cluster. It can be seen at what similarity levels the clusters are
formed and the composition of the clusters of the final partition. On the
other hand, if the dataset is too large, the dendrogram visualization capa-
bility is poor and it provides little help.
The partitional procedure uses nonhierarchical clustering of observa-
tions. K-means clustering, the most widely used partitional method,
works best when sufficient information is available to make good starting
cluster designation or expected number of clusters, k. K-means procedure
first randomly selects k centroids (objects) and then decomposes objects
into k disjoint groups by iteratively relocating objects based on the simi-
larity between the centroids and the objects. In K-means, a cluster cen-
troid is the mean value of objects in the cluster. The method is generally
more accurate than hierarchical clustering. However, it provides no visu-
alization aid. On the other hand, hierarchical procedure is much slower
in computational time than partitional procedure. The latter can practi-
cally handle much larger datasets than the former.
Cluster analysis of the data presented in Fig. 8.1 was performed with
the software package Minitab 17 using hierarchical and K-mean methods
with the centroid linkage and Euclidian distance. An example of a five-
cluster dendrogram and corresponding cluster partitioning (k = 5) are
presented in Fig. 8.2.
It follows from both graphs (top and bottom) that according to hierar-
chical clustering, there is one zip code 13228, which forms a 1-object
cluster with sharply distinct minimal CM, $5588 (cluster 1). Five zip
codes (13226, 13222, 13221, 13225, and 13213) form cluster 2 with
relatively close CM, from about $6700 to $7100. The next nine zip codes
(from 13072 to 13216) form cluster 3 with the CM range from $7500 to
$8600. Seven zip codes (from 13211 to 13224) form cluster 4 with the
CM range from $8900 to $9500, and the last six zip codes from 13218
to 13403 form cluster 5 ($10,200–$11,500).
Comparative results of the K-mean clustering procedure for k = 5 clus-
ters using the same original data are presented in Fig. 8.3.
In contrast to the hierarchical procedure, the K-mean procedure selects
zip code 13403 as a one-object cluster (cluster 1), which has the highest
264 A. Kolker
58.19
Similarity
72.12
cluster 1
86.06
cluster 5 cluster 4 cluster 2
cluster 3
100.00
3 9 4 3 0 7 8 4 5 7 9 5 1 1 6 2 7 4 5 2 7 0 2 3 5 1 2 6 8
40 0 2 02 22 2 2 02 21 22 2 1 22 21 0 4 05 21 2 1 21 21 21 0 0 02 03 2 1 07 21 22 2 2 22 22 2 2
13 13 13 13 13 13 1 3 13 13 1 3 13 13 13 1 3 13 13 1 3 13 13 13 13 13 13 1 3 13 13 1 3 13 13
ZIP codes
10000
7500
CM, $
5000
2500
0
3 9 4 3 0 7 8 4 5 7 9 5 1 1 6 2 7 4 5 2 7 0 2 3 5 1 2 6 8
40 02 02 2 2 22 02 21 2 2 21 22 21 0 4 05 21 21 2 1 21 21 00 02 0 3 21 07 21 2 2 22 22 22 2 2
13 1 3 13 13 13 1 3 13 13 13 1 3 13 13 13 1 3 13 13 13 13 1 3 13 13 13 1 3 13 13 13 1 3 13 13
ZIP codes
Fig. 8.2 Hierarchical clustering: 5-cluster dendrogram (top) and cluster partition-
ing (bottom) for CM data from Fig. 8.1
CM, $11,518. The next two zip codes (13029 and 13024) form a two-
member cluster 2 with CM from $11,312 to $11,376. A four-member
cluster 3 is formed by four zip codes from 13223 to 13027, with CM
from $10,873 to $10,567, respectively, and so on, as indicated in Fig. 8.3.
It appears that K-mean partitioning in this example is more intuitively
appealing and makes more sense than hierarchical partitioning.
Thus, depending on the required granularity of data, the variable CM
with a total of 29 observations (zip codes) can be reduced to the variable
with 5 (or another small number) observational groups (clusters) with
values that are similar to each other within the cluster with a given simi-
larity metric. Thus, if there is some known information available about
an object (zip code), one can take full advantage of this information to
understand the other objects (other zip codes) in the cluster to which this
object belongs. Analysis of a few relatively small groups of similar objects
(zip codes) is usually easier than that of a large number of original un-
clustered observations.
8 Statistical Data Analytics for Decision-Making 265
12000 cluster 1
cluster 2
cluster 3
10000
cluster 4
8000
cluster 5
CM, $
6000
4000
2000
0
3 9 4 3 0 7 8 4 5 7 9 5 1 1 6 2 7 4 5 2 7 0 2 3 5 1 2 6 8
40 02 02 22 22 02 21 22 21 22 21 04 05 21 21 21 21 21 00 02 03 21 07 21 22 22 22 22 22
13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13 13
ZIP codes
Fig. 8.3 K-means clustering: cluster partitioning for CM data from Fig. 8.1
Problem Description
Patient census forecasting is paramount for cost-efficient workforce labor
management and an accurate assessment of the nursing staffing level (see
Sects. 4.1 and 4.2). Various time series forecasting techniques are avail-
able in the literature such as moving averages of various orders, exponen-
tial and double exponential smoothing with trend, Holt-Winters with
seasonality and trends, autoregressive integrated moving average
(ARIMA), and many others. Their applications and comparative pros
and cons are reviewed elsewhere in the literature (Ozcan, 2009) and in a
more advanced treatment, e.g., by Shumway and Stoffer (2017).
The California nurse-to-patient staffing ratio assembly bill AB 394 in
2004 required the State Department of Health Services to establish mini-
mal nurse-to-patient ratios for acute care and specialty hospitals in this
state. To comply with the law, Jacobs and Salinas Medical Centers col-
lected census data for every day for every hour (a total of a year’s worth of
data) in a number of their nursing units. They requested to make a fore-
cast of patient census for, at least, the next day’s shifts for rapid adjusting
of the corresponding staffing levels. To make census forecasting practi-
cally useful for that purpose, the forecasting error should have been
within ±1 or 2 patients. This is a challenging task.
0.3, and 0.5, respectively, with higher weight applied to the most recent
data points. In the single exponential smoothing, the oldest piece of data
is eliminated once a new piece is added. The forecast is calculated by add-
ing the previous forecast and the difference between the actual previous
value and the previous forecast multiplied by a weighting or smoothing
factor. The values 0.3–0.6 for this factor were tested. The accuracy of
these forecasts varied from case to case and, in general, it was not enough.
There are also available specialized forecasting software packages.
Among them Forecast Pro ([Link]), SAS Business
Analytics ([Link]), and International Institute of Forecasters
([Link] These packages required commercial licenses
and were not tested in the medical centers specifically for patient census
forecasting capability.
points could be considered too “old” if they are weakly correlated or not
correlated at all to the newer points. These “old” data points should not
be used for forecasting. Otherwise, the forecast will most likely be skewed.
The data points that are strongly correlated to the newer ones can be used
for making the forecast. Note that a strong linear correlation implies a
strong interdependency of the data points. The opposite, however, is not
true: data points could be interdependent but not linearly correlated.
Thus, too many past census data points are not necessarily needed for
making an accurate forecast and could even be detrimental to it.
The maximum number of “steps back to the past” at which the older
data points are still strongly enough linearly correlated to the newer data
points can be estimated using an autocorrelation function (ACF) of a
time series. An autocorrelation function (ACF) of a time series is a mea-
sure of the linear interdependency between the data points separated by
k time units (time lag). ACF is represented as a plot of the linear correla-
tion coefficients of the time series data points with themselves separated
by k time units as a function of k.
A possible measure of a correlation cutoff value k is the first zero cross-
ing of the time-lag axis at which time series data points become totally
uncorrelated. However, these metrics will include quite a few data points
that are weakly correlated to the new data points, such as correlation
coefficients that are in the range of 0–0.5. Sometimes, the correlation
cutoff lag is defined in the statistical literature as the smallest value k that
makes ACF(k) < K, where K = e−1 = 0.37 or K = 0.5. At the same time, it
is usually accepted that a strong enough linear correlation coefficient
starts with ~0.6. Therefore, it was assumed here that the number of the
past data points that should be used for the purpose of forecasting (a
cutoff value of the time-lag) corresponded to the ACF value of about 0.6
of the census time series.
As an example, 1 year’s worth of census data collected every 4 h each
day of the year is presented in Fig. 8.4 (Jacobs Medical Center (JAMC),
Med-Surgical unit).
The autocorrelation function (ACF) of this census time series is pre-
sented in Fig. 8.5.
8 Statistical Data Analytics for Decision-Making 269
26.0
25.0
24.0
23.0
22.0
21.0
20.0
19.0
18.0
17.0
census
16.0
15.0
14.0
13.0
12.0
11.0
10.0
9.0
8.0
7.0
6.0
5.0
time, hours
Fig. 8.4 One year’s worth of census data (JAMC, med-surg unit)
Fig. 8.5 The autocorrelation function of the census time series presented in
Fig. 8.4. Each census data point was collected every 4 h; hence, the cutoff of 96 h
corresponds to 24 data points
calculation using the remaining 2159 data points. In this formula, the
census values were stored in cells C$2:C$2159. Column A contained the
index variable from 1 to 2158 to count the number of census values. The
formula should have stopped about 10 data points before the end of the
training part of the time series.
ACF(k) can, of course, be calculated in R or RStudio software using
the ACF () function from the “tseries” library or using the Python library
of functions.
It follows from this plot that the ACF goes down pretty fast with some
small periodic fluctuations and remains low in the range of 0.2–0.3 from
the time lag of 96 h. However, only 24 census data points remain rather
strongly correlated to each other (>~0.6) (census was collected here every
4 h, so 96 h correspond to 24 data points), as indicated in Fig. 8.5. Thus,
only the last 24 data points were used for forecasting (out of the total
available 2158 data points).
Another example is 1 year’s worth of census data from Salinas Medical
Center (unit 1.607) collected every hour presented in Fig. 8.6.
12.0
11.0
10.0
9.0
8.0
7.0
census
6.0
5.0
4.0
3.0
2.0
1.0
0.0
time, hours
Fig. 8.6 One year’s worth of census data (Salinas Medical Center, unit 1.607)
8 Statistical Data Analytics for Decision-Making 271
Fig. 8.7 The autocorrelation function of the census time series presented in
Fig. 8.6
272 A. Kolker
41
40
39
38
37
36
35
34
33
32
31
30
Census
29
28
27
26
25
24
23
22
21
20
19
18
17
16
Time, hours
Fig. 8.8 One year’s worth of census data (Salinas Medical Center, Medical-
Surgical unit)
Fig. 8.9 The autocorrelation function of the census time series presented in
Fig. 8.8
Many more nursing units from Salinas, JAMC, and Baptist Medical
Center were similarly analyzed.
8 Statistical Data Analytics for Decision-Making 273
m
yn d j yn j xn , (8.1)
j 1
n n
right set of coefficients dj, one can predict the future data points of a time
series from a record of its past.
Because the predicted future data points, yn, are generated using previ-
ously calculated data-points, yn−j, such a procedure is called the recursive
procedure; its predicting behavior pattern quickly becomes much more
complex than a straight line, or a high-order polynomial. It is especially
successful at predicting time series that are rather smooth and oscillatory
though not necessarily periodic (Press et al., 1988).
However, in order to achieve its full usefulness, recursive procedure
must be stable. Recursive procedures feed on their own output; therefore,
they are not always stable, i.e., some particular “bad” sets of coefficients,
dj, can generate an exponentially growing output.
274 A. Kolker
n
z n d j z n j 0 (8.2)
j 1
are inside the unit circle, i.e., satisfy the condition |z| ≤ 1. Press et al.
(1988) provide a detailed description of the computational procedure for
using recursive prediction of the future data points of a time series using
the above equations.
Parameter m which is the number of linear recursive coefficients is the
only model’s parameter that can be estimated by training data points by
comparing forecast with the actual data. This parameter should be chosen
to be rather small by experimenting with the data. Usually, m should be
in the range 5–40. Larger values require massaging all complex roots of
the above polynomial equation that can be quite sensitive to round-off
errors and should be avoided (Press et al., 1988). Notice that the maximal
number of coefficients, m, must be one less than the number of used data
points, n, i.e., m ≤ n − 1. Moreover, the sum of coefficients dj must be
equal to 1 (or practically as close to 1 as possible) to avoid forecasting
biasing (Press et al., 1988).
Forecasting Results
Census forecasting was performed for the next day’s time horizon: fore-
casting 6 data points for every 4 h for JAMC (forecast 6 data points for
the next 24 h), and 24 data points for every 1 h for all units for Salinas
medical center.
An example of the forecasting screenshot is presented in Fig. 8.10.
The most recent census data points equal to the next day’s number of
forecast horizon data points were not included in the calculation ACF(k)
and the number of linear recursive coefficients, m, in Eq. 8.1. These data
points were used only for the estimation of the forecast error.
The mean forecasting error %, the mean absolute forecasting error (vs.
actual census), the number of used past data points determined from the
8 Statistical Data Analytics for Decision-Making 275
Fig. 8.10 JAMC medical center. Forecast of 6 census data points (future) every 4 h
for the next 24 h (screenshot)
corresponding ACF of the census time series, and the number of coeffi-
cients, m, are presented in Table 8.4 for JAMC (Jacobs Medical Center)
and for the various nursing units of Salinas Medical Center. Mean fore-
casting error % was defined as the mean absolute difference between the
actual and forecasted census divided by the actual census. The mean abso-
lute forecasting error was defined as the mean absolute difference between
the actual and predicted census. The absolute difference between the
actual and forecasted census did not exceed 2 for the majority of units, so
the forecasting for the next day’s time horizon was reasonably accurate.
Forecasting accuracy for a longer time horizon (2 days or longer) sig-
nificantly deteriorated. This was expected. When the new census data
276 A. Kolker
points become available for the day, the short-term forecast for the fol-
lowing day must be reevaluated using the same procedure. Naturally, no
forecast can be very accurate all the time; forecasting errors are inevitable.
In general, it is impossible to accurately predict the future pattern (espe-
cially for a long-time horizon) based only on past data.
Nonetheless, if the underlying factors make the patient census time
series stable enough and it contains no rare “black swan” events, such as
an emergency patient admission due to a major industrial accident or a
mass terrorist attack, then a recursive linear digital filtering could serve as
a reliable short-term forecasting methodology.
8 Statistical Data Analytics for Decision-Making 277
Problem Description
As authorities try to manage the COVID-19 crisis and contain the spread
of the virus, they face the challenge of deciding who is really infected and
who is not to send the infected people to quarantine while allowing the
non-infected to go. From the authorities’ point of view, there could be
four main scenarios: (i) if one is infected and diagnosed with coronavirus,
this individual should be quarantined for the public benefit, (ii) if one is
not infected but diagnosed with coronavirus, the individual may be
wrongly quarantined causing personal inconvenience. The public will
suffer no major harm but the authorities will have to spend some
resources; (iii) if one is infected and not diagnosed with coronavirus, then
the individual will be wrongly released likely contributing to spreading
the virus. This puts the public in danger of an expanded outbreak; (iv) if
one is not infected and not diagnosed with coronavirus, this individual
will rightly be released. The errors that the authorities can make are in
scenarios (ii) and (iii). Scenario (ii) is a minor inconvenience (if not done
too often), but scenario (iii) is the major issue that can cascade into a
larger outbreak. If an outbreak occurs, the authorities will have to do
contact tracing for possibly hundreds of people given the contagiousness
of coronavirus. This will be costly, so the primary interest is in minimiz-
ing the probability of the third scenario.
Thus, the authorities need to make an accurate assessment, so that they
can avoid releasing the infected individuals and quarantining the non-
infected ones. To achieve this goal, the authorities should first make an
initial assessment of all suspected infections whether they are patients at
the medical clinic or travelers from places with active outbreaks.
There are clues which hint at a COVID-19 infection, such as:
Not all of these clues are immediately available. Authorities at the air-
port can only screen people for temperature and identify travel history if
needed. A doctor will know the symptoms disclosed by the patient, but
it is the doctor who is supposed to take the initiative to find and link the
clues together. And asymptomatic cases may not have any clues.
So, suppose a coronavirus diagnostic test was ordered to decide on the
individual infection status. The test result was reported as positive. Based
on this single test outcome, could it be concluded that the individual was
infected and counted in statistics as such, and should this individual be
quarantined? If the test was negative, could the individual be declared
infection-free and let go? Should population mass testing be mandated?
What decision on further action should be made by the authorities based
on the population mass testing results?
However, the urge for mass overboard testing does not take into
account the issues with a test interpretation. Why the authorities cannot
rely directly on mass testing outcomes with a clear-cut answer: infected or
not infected? Interpreting the result of the test for COVID-19 is typically
associated with two types of mistakes. The first one is that a positive test
result means a patient has been infected, while a negative result means
that a patient does not. But this is an incorrect interpretation because no
test is 100% accurate.
Usually, the companies that manufacture the test kits specify the accu-
racy of the test. If a company markets its test as, say, “90% accurate,”
another common mistake is that a positive result means a 90% chance of
being infected and a negative result means a 90% chance of not being
infected. This interpretation is also not correct, but it is actually surpris-
ingly common. How to avoid these misinterpretations?
By definition, the sensitivity (Sn) is the true positive rate, i.e., the frac-
tion of yielding a correct positive outcome when the condition being
tested is actually present. The specificity (Sp) of the test is the true nega-
tive rate, i.e., the fraction of yielding a correct negative result when the
condition being tested is absent.
Bayes’ formula includes three parameters: (i) initial assessment or the
disease prevalence (prior probability) in the general population or in the
target subpopulation, Pr, (ii) the test’s specificity, Sp, and (iii) the test
sensitivity, Sn. The test’s positive predictive value (PPV) or the posterior
probability that the condition is present (infected) given a positive test
outcome is
Sn Pr
PPV infected | positive test (8.3)
Sn Pr 1 Sp 1 Pr
The test negative predictive value (NPV) or the posterior probability that
the condition is absent (not infected) given a negative test outcome is
Sp 1 Pr
NPV not infected | negative test (8.4)
Sp 1 Pr 1 Sn Pr
There are at least two types of tests for mass testing: real-time reverse
transcriptase or RT-PCR and a rapid antigen test (Ag-RDTs). There were
quite a few studies published aimed at determining the specificity and
sensitivity of these tests. Unfortunately, there is no accepted “gold stan-
dard” to compare with. Nonetheless, according to one of the recent sys-
tematic reviews of the RT-PCR test accuracy (Watson et al., 2020), the
current estimate of its lower end for the sensitivity (Sn) is about 70% and
for specificity (Sp) 95%.
Soni et al. (2023) demonstrated that serial rapid antigen testing Ag-
RDT taken twice 48 h apart resulted in an aggregated sensitivity of
93.4% for symptomatic participants; two-time serial testing among
asymptomatic participants was lower at 62.7% but it improved to 79.0%
with testing three times at 48-h intervals. There are pros and cons for
8 Statistical Data Analytics for Decision-Making 281
each type of test, but RT-PCR (as a single test) is deemed in general more
accurate and reliable than a rapid antigen test.
One could face two issues in practical applications of formulas
(8.3–8.4):
ln PPV PPV 1 Pr
(8.5)
Sn Pr, Sp PPV Sn Pr, Sp Sn
Sn Pr 1 Sp 1 Pr
and
282 A. Kolker
ln PPV PPV 1 Pr
Sp (8.6)
Pr, Sn PPV Sp Pr, Sn Sn Pr 1 Sp 1 Pr
and
ln NPV NPV 1 1 Sn Sp
(8.8)
Pr Sn, Sp NPV Pr Sn, Sp 1 Pr
Sp 1 Pr 1 Sn Pr
B C D E F G H I J K L
Single test
Single test 2 tests PPV, 2 tests NPV, Relative Relative Relative Relative Relative Relative
Infection NPV,
PPV, both tests both tests PPV PPV NPV NPV PPV NPV
Prevalence negative-
positive- positive- negative-not sensitivity sensitivity sensitivity sensitivity sensitivity sensitivity
% not
infected infected infected to Sn, % to Sp, % to Sn, % to Sp, % to Pr, % to Pr, %
infected
1% 12.4% 99.7% 66.4% 99.9% 1.3 17.5 0.01 0.003 88.50 -0.32
5% 42.4% 98.4% 91.2% 99.5% 0.8 11.5 0.05 0.017 12.12 -0.34
15% 71.2% 94.7% 97.2% 98.3% 0.4 5.8 0.18 0.056 2.26 -0.41
25% 82.4% 90.5% 98.5% 96.8% 0.25 3.5 0.32 0.10 0.94 -0.51
50% 93.3% 76.0% 99.5% 90.9% 0.10 1.3 0.80 0.25 0.27 -0.96
75% 97.7% 51.4% 99.8% 77.0% 0.03 0.5 1.62 0.51 0.12 -2.59
85% 98.8% 35.8% 99.9% 63.9% 0.02 0.25 2.14 0.68 0.10 -5.03
90% 99.2% 26.0% 99.9% 52.7% 0.01 0.16 2.47 0.78 0.09 -8.22
Fig. 8.11 PPV and NPV values and their sensitivities to Sn, Sp, and Pr values for
single and twice repeated tests (Sn = 70%, Sp = 95%)
8 Statistical Data Analytics for Decision-Making 283
2
Sn Pr
PPV infected | two positive tests (8.9)
2 2
Sn Pr 1 Sp 1 Pr
For two repeated negative tests, the combined specificity and sensitivity
(in contrast to PPV) should be calculated as Sp(comb) = Sp2 and
Sn(comb) = 1 − (1 − Sn)2. (This is analogous to parallel testing). Hence,
2
Sp 1 Pr
NPV not infected | two negative tests (8.10)
2 2
Sp 1 Pr 1 Sn Pr
It was assumed here that the Sn and Sp of the repeated tests were the
same. It should be noted though that Sn and Sp could vary over time
because of underlying changes in the disease progress (and, consequently,
a change in viral load). Indeed, each test would likely be conducted at
different times (say, RT-PCR followed by a rapid Ag-test or chest CT if
warranted). This would make it difficult to get combined PPV or NPV
without determining first Sn and Sp specifically for the timing of each test
which is not practical in most cases. Therefore, the assumption of the
same values of Sn and Sp could be justified (or not) by taking into account
the PPV and NPV sensitivities to Sn and Sp using formulas (8.5)
and (8.6).
Formulas 8.3–8.10 were built in the Excel spreadsheet (the spread-
sheet file is available in SEM 8.2).
Using the point values Sn = 70% and Sp = 95% for RT-PCR, some
PPV and NPV values for single and twice repeated tests are presented in
Table in Fig. 8.11. The sensitivities of PPV and NPV to Sn and Sp values
284 A. Kolker
Sp comb 1 i Sp i and Sn comb i Sn i .
Hence,
Pr i Sn i
PPV infected all positive tests (8.11)
Pr i Sn i 1 Pr i 1 Sp i
For NPV,
Sp comb i Sp i and
Sn comb 1 i 1 Sn i . Hence,
8 Statistical Data Analytics for Decision-Making 285
1 Pr i Sp i
NPV not infected | all negative tests (8.12)
1 Pr i Sp i Pr i 1 Sn i
Fig. 8.12 The equivalence of combined Sp(comb) and Sn(comb) and updated PPV
and NPV as the new priors
negative) tests out of the total M tests given the probability of true posi-
tive (or true negative) outcome of a single test.
As a summary of the overall content of this section, it seems that mass
population testing to confirm or rule out the infection during an active
outbreak is not a feasible solution. There are simply not enough resources
to test everyone especially twice or more. This may explain why the
authorities are trying to implement the over-the-board large-scale con-
tainment measures: they assume that everyone in the general population
is, or will be, infected. While the assumption of high infection prevalence
could be justified for more susceptible subpopulations mentioned above,
it is a questionable practice when applied to the general population caus-
ing a lot of inconvenience for too many people and a lot of economic
distress for society, as was seen in the last few years since the
COVID-19 onset.
Despite the caveats discussed above, the Bayesian inference remains
the most powerful methodology used in interpretation and making sense
of the diagnostic testing results.
Table 8.6 Relation between the test outcome and its measure (OM matrix)
Independent of its association with the disease, the presence of the test
measure is also associated with positive or negative test outcome (O) as
shown in the OM matrix in Table 8.6. The term “test outcome” means
the result of applying the test procedure to capture the presence of the
test measure. It addresses the question of how frequently the test outcome
turns positive or negative in the presence or absence of the test measure
regardless of the presence or absence of the disease itself, e.g., the tester’s
ability to interpret and report the number of RT-PCR cycles as infected
or not. Thus, the small tpr(sn), fpr, fnr, and tnr(sp) in the OM matrix are
numerically different from and independent of the corresponding coun-
terparts in the MD matrix in Table 8.5.
As indicated earlier, the characteristics of a test are usually described
only in terms of relationships between the test outcome and the presence
or absence of the disease (OD matrix) without paying attention to the
hidden relationships between the test outcome and its measure.
Sonnenberg (2004) suggested the following treatment to establish such a
relationship.
To re-assemble the original OD matrix of test outcome vs. disease, the
OM matrix needs to be multiplied by the MD matrix, i.e.,
OD = OM × MD (in this particular order). This matrix multiplication
yields the desired OD matrix shown in Table 8.7.
290 A. Kolker
Table 8.7 Relationship between test outcome and disease (OD matrix)
sn fpr Sn FPR
OD OM MD
fnr sp FNR Sp
sn Sn fpr FNR sn FPR fpr Sp
fnr Sn sp FNR fnr FPR sp Sp
Pr Sn sn 1 sp 1 Sn
PPV
Pr Sn sn 1 sp 1 Sn (8.13)
1
Pr sn 1 Sp 1 Sp Sp
and
8 Statistical Data Analytics for Decision-Making 291
Table 8.8 PPV and NPV values for the single test with an imperfect test examiner
Prevalence, PPV (infected-positive test), NPV (not infected-negative test),
% % %
1 6.6 99.6
5 26.8 98.0
15 55.2 93.5
25 69.9 88.5
50 87.5 71.9
75 95.4 46.0
85 97.5 31.1
90 98.4 22.1
Test characteristics: Sn = 0.7 and Sp = 0.95. A test examiner performance
characteristics: sn = 0.9 and sp = 0.95
1 Pr 1 sn 1 Sp sp Sp
NPV
1 Pr 1 sn 1 Sp sp Sp (8.14)
Pr 1 sn Sn sp 1 Sn
These formulas were built in the Excel spreadsheet. Once again, with
perfect capture of the test measure sp = sn = 1 and fpr = fnr = 0, these
equations are reduced to the previous Bayes equations for PPV and NPV
(8.3) and (8.4).
It is assumed here that most of the test examiners are qualified and
trained in capturing the test measure, so their own sensitivity, sn, and
specificity, sp, are high enough. Assuming sn = 0.9 and sp = 0.95 along
with the earlier used test values Sn = 0.7 and test Sp = 0.95, the PPV and
NPV for various prevalences are presented in Table 8.8.
Comparing this table with results in Table in Fig. 8.11 for a single test
demonstrates that even a high enough (but not perfect) ability of the test
examiner to capture the test measure results in a significant change in the
diagnostic test outcome. Indeed, the PPV with the imperfect test
examiner drops down from ~12.4% to 6.6% for prevalence of 1% and
from ~42.4% to 26.8% for a prevalence of 5%. For the high-end preva-
lence values, PPV with the imperfect test examiner also becomes smaller
but not too much compared to values in Table on Fig. 8.11 because PPV
moves closer to 100% anyway. A similar drop in NPV values is also
292 A. Kolker
observed for the imperfect test examiner but it looks relatively less pro-
nounced than that for PPV.
In the extreme case of dropping the test examiner’s performance charac-
teristics to sn = sp = 50%, the positive predictive values PPV drops down
to the level of prevalence, i.e., PPV = Pr while negative predictive values
NPV become 1 − Pr irrespective of the test’s own high level of Sp and Sn.
According to Sonnenberg (2004), in this extreme case “…the test does not
help in diagnostics at all…and a random choice from the tested popula-
tion would serve as well as performing the test.” This demonstrates that it
is the test examiner who ultimately influences the test diagnostic outcome.
Sensitivities of PPV and NPV to the test examiner’s performance char-
acteristics could be evaluated similarly to those performed earlier.
To summarize, the characteristics of a diagnostic test are determined,
first, by the association between a diagnostic parameter and a disease and,
second, by the technical ability to measure the diagnostic parameter. In
general, the test examiner’s competence, technical performance, and asso-
ciation of the test measure with the disease are described in terms of
comprehensive test characteristics such as sensitivity and specificity.
However, the relationship between the test outcome and diagnosis can be
broken down into two separate underlying contributing terms: (i) the
relationship between the test measure and diagnosis and (ii) the relation-
ship between the test outcome and the test measure. It is the second term
that takes into account and quantifies the influence of the test examiner
(be it a human being or an artificial intelligence system, AI) in making
the diagnostic test conclusion. This highlights the ultimate role of staff
training and competence or the software quality and reliability in per-
forming mass testing when the diagnostic conclusion involves the test
examiner’s judgment.
The following fundamental management principles were illustrated in
this section:
(i) the lower the prevalence of a disease or infection, the more likely a
positive test will be the false positive. That is to say, the PPV value
will be too low to convincingly conclude that the disease or infection
is present. In this case, a repeated test could be warranted. Thus, mass
screening of general populations with low disease (infection) preva-
8 Statistical Data Analytics for Decision-Making 293
lence aimed at its early detection (irrespective of the risk status of the
individual) is not an effective strategy.
(ii) Mass screening of the selected subpopulations with moderate to high
disease (infection) prevalence (expectation) could be effective in
detecting the disease with a single positive test but a negative test
warrants to be repeated.
Problem Description
It was demonstrated in Sect. 2.14 that over-improvement in the upstream
unit, such as ED, can negatively impact the performance of downstream
units, such as operation rooms and ICU. This happens because the capac-
ity of the downstream units is often not enough to handle a patient vol-
ume increase coming from the ED, especially if that increase was not
anticipated with sufficient lead time. To improve the preparedness of the
downstream inpatient units to handle patient volume coming from ED,
it would be helpful to predict the number of patients that are expected to
be discharged from ED in the next period of time and admitted as
inpatients.
This information would give a sufficient lead time to the downstream
units’ staff to expedite cleaning and making the needed beds available or
to consider discharging or transferring some “old” patients to make the
appropriate room for the new incoming patients.
T t
p f T dT
0
ESI
T t T
T
fESI T dT fESI T dT 1 fESI T dT qESI
0 0 0
8 Statistical Data Analytics for Decision-Making 295
Hence,
T t T
fESI T dT fESI T dT
FESI T t FESI T
qESI T ,t 0 0
T
1 FESI T
1 fESI T dT
0
where FESI(T) is the cumulative length of stay distribution for each ESI,
respectively.
Littig and Isken (2007) used a similar formula (sometimes called the
hazard ratio) to compute the conditional probability that the patient will
leave in the next time period (24 h) given the current LOS and cumula-
tive LOS distributions (applied to the entire hospital occupancy predic-
tion model).
Examples of LOS distribution densities for four ESI along with their
cumulative distributions for a case study hospital are given in Figs. 8.13
and 8.14.
ESI=1 ESI=2
30 30
Percent
20 20
Percent
10 10
0 0
1 2 3 4 5 6 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Length of stay (LOS), hours Length of stay (LOS), hours
ESI=3 ESI=4
30
20
15 20
Percent
Percent
10
10
0 0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 2 3 4 5 6 7 8 9 10 11 12
Length of stay ( LO S) , hour s Length of stay (LOS), hours
ESI=1 ESI=2
100 100
Cumulative Percent
75 75
Cumulative Percent
50 50
25 25
0 0
1 2 3 4 5 6 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Length of stay (LOS), hours Length of stay (LOS), hours
ESI=3 ESI=4
100 100
75 75
Cumulative Percent
Cumulative Percent
50 50
25 25
0 0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 2 3 4 5 6 7 8 9 10 11 12
Length of stay (LOS), hours Length of stay (LOS), hours
INPUT OUTPUT
Expected
number of
Number of Lead Probability of patients to be
patients in Current Discharge discharge in the next discharged in
ESI ED LOS=T hrs Time, t hrs time period, t hrs the next t hrs
ESI=1 5 4 2 0.792 4
6 5 2 0.480 3
Total ESI=1
expected to be
11 7
Total ESI=1 discharged in the
in ED next 2 hours
ESI=2 6 8 2 0.654 4
7 9.5 2 0.607 4
Total ESI=2
13 expected to be 8
Total ESI=2 discharged in the
in ED next 2 hours
ESI=3 5 10 2 0.719 4
8 6 2 0.595 5
Total ESI=3
13 expected to be 9
Total ESI=3 discharged in the
in ED next 2 hours
ESI=4 3 8 2 0.873 3
4 7 2 0.521 2
Total ESI=4
7 expected to be 5
Total ESI=4 discharged in the
in ED next 2 hours
Fig. 8.15 An example of the calculator input and output panels for the antici-
pated number of discharged ED patients
References
Fuhs, P., Martin, J., & Hancock, W. (1979). The use of length of stay distribu-
tions to predict hospital discharges. Medical Care, XV11(4), 355–368.
Glantz, S., & Slinker, B. (2001). Applied regression & analysis of variance (2nd
ed.). McGraw-Hill, Inc.
Jobson, J. D. (1992). Applied multivariate data analysis (V. 2. Categorical and
multivariate methods). Springer.
Littig, S., & Isken, M. (2007). Short term hospital occupancy prediction. Health
Care Management Science, 10, 47–66.
Ozcan, Y. (2009). Quantitative methods in health care management (2nd ed.,
p. 438). Jossey-Bass. A Wiley Imprint.
Peek, N. (2010). Chapter 24: Data mining. In Y. Yih (Ed.), Handbook of health-
care delivery systems. CRC Press.
Press, W., Flannery, B., Teukolsky, S., Vetterling, W., (1988). Numerical recipes
in C: The art of scientific computing (2nd ed., pp. 735). Cambridge University
Press, Cambridge.
Shumway, R., & Stoffer, D. (2017). Time series analysis and its applications with
R examples (Springer texts in statistics) (4th ed.). Springer.
Soni, A., Herbert, C., & +44 authors. (2023). Performance of rapid antigen
tests to detect symptomatic and asymptomatic SARS-CoV-2 infection.
Annals of Internal Medicine. American College of Physicians. [Link]
org/10.7326/M23-0385
Sonnenberg, A. (2004). We only see what we already know—A modified Bayes’
formula to explain inherent limitations of diagnostic tests. Medical Hypothesis,
63, 759–763. Elsevier. [Link]
Tibshirani, R., Walther, G., & Hastie, T. (2001). Estimating the number of
clusters in a dataset via the gap statistic. Journal of Royal Statistical Society, B,
63, 411–423.
Watson, J., Whiting, P., & Brush, J. (2020). Interpreting a Covid-19 test result.
British Medical Journal, 369. [Link]
WHO. (2020). [Link]
who-head-our-key-message-is-test-test-test
Yoo, I., Alafaireet, P., Marinov, M., Pena-Hernandez, K., Gopidi, R., Chang,
J.-F., & Hua, L. (2011). Data mining in healthcare and biomedicine: A sur-
vey of the literature. The Journal of Medical Systems., Published Online.
[Link]
9
Allocating Cost Savings (Gains)
Between Cooperating Providers:
Introduction to the Concept
of the Shapley Value
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2024 299
A. Kolker, Healthcare Management Engineering In Action, Business Guides on the Go,
[Link]
300 A. Kolker
Equitable allocation: gives everyone the same satisfaction level, i.e., the
proportion each player receives by their own valuation is the same for
all of them. This is a difficult aim as players might not be truthful if
asked about their valuation.
Proportional allocation: guarantees that each player gets his share. For
instance, if three people divide up an asset, then each gets at least a
third by their own valuation.
Envy-free allocation: everyone prefers his own share to the others. No one
is jealous of anyone else. No one would trade his share with any-
one else’s.
An efficient or Pareto optimal allocation: ensures that no other allocation
would make someone better off without making someone else worse
off. The term efficiency comes from the economic idea of an effi-
cient market.
Merit-based allocation: the more one brings to the coalition, the more one
gets out of the division of the accumulated gains.
The focus of this chapter is the Shapley value because it exhibits some
attractive properties (including monotonicity) and it is easier to compute
than the nucleolus. Therefore, it is more widely used for practical applica-
tions (Young, 1994).
The Shapley value aims at the “fair” allocation of collective costs or
gained savings between the collaborating participants based on the rela-
tive (marginal) contribution of each participant to their cooperative
activities. Basically, the more one brings to the coalition, the more one
gets out of the division of the accumulated gains. This is an example of
meritocracy fairness.
An example of using the Shapley value for cost allocation of the adver-
tising budget among three marketing channels and their combinations is
presented by Nagadevara (2019) although without showing the internal
calculation details.
It should also be noted that an approach was developed in which a
player’s payoff depends on the worth of the coalitions to which one
belongs and not necessarily on one’s marginal contributions, the so-called
membership values (Klenberg & Weiss, 2013). This means that one
would get some share regardless of whether or not one makes a marginal
contribution to the welfare of the coalition. It was demonstrated that the
set of non-marginal membership values includes those that embody
widely held notions of fairness, such as partial “benefit equalization,”
individual rationality, and “greater rewards follow from greater contribu-
tions,” where one’s contributions are not measured marginally.
produce $85 worth of the product. The one’s added value to the group is
defined as
$85 × 4 − $75 × 3 − $80 × 1 = $35, i.e. $85 − $80 = $5 for the newly
joined member and $10 for each of the three others. This is also called
the “marginal contribution” of the participant. It is a measure of what
this group member’s presence is worth, above the minimum that he
would require for his service. Notice that if one cannot produce any-
thing working alone ($0), then the added value becomes greater, e.g.,
$85 × 4 − $75 × 3 = $115, i.e., $85 for the newly joined member and
$10 for each of three others. Thus, joining the group provides more
benefits if one cannot produce anything alone.
s 1! n s ! V
Sh k s V s k ,
S n n!
Problem Description
For a long time, policymakers have been increasingly frustrated with fee-
for-service payment system in healthcare settings. Fee-for-service rewards
volumes and encourages silos and fragmentation of care. Several provi-
sions of the Affordable Care Act (ACA) healthcare legislation in the
United States seek to shift provider payments to value-based approaches
that encourage quality improvement and cost reduction. Payment bun-
dling promoted by the Center for Medicare Services (CMS) is one such
9 Allocating Cost Savings (Gains) Between Cooperating… 305
Table 9.1 Standalone costs and the Shapley cost allocation for group size n = 3
Physician
# Participant (member) → Hospital {1} group {2} SNF {3} Total
1 Median Medicare cost $30,000 $6000 $4000 $40,000
for standalone
providers per
episode-of-care
2 2-member coalitions → {12} {13} {23}
3 Cost reduction for 3.0% 2.0% 1.5%
2-member coalition
vs. standalone
member cost, %
4 Cost for 2-member $34,920 $33,320 $9850
coalitions
5 3-member coalition → {123}
6 Cost reduction for 3.5%
3-member coalition
vs. standalone
member cost, %
7 Cost for 3-member $38,600
coalitions
8 Shapley value for Shapley value Shapley value Sum of
hospital: provider {1} for Physician for SNF: Shapley
group: provider {3} values
provider {2}
9 $29,290 $5555 $3755 $38,600
10 Saving {1} Saving {2} Saving {3} Sum of
savings
11 $710 $445 $245 $1400
12 Saving allocation % {1} Saving Saving
allocation allocation
% {2} % {3}
13 2.4% 7.4% 6.1%
14 Coalition → {12} {13} {23}
15 Condition 1 for OK OK OK
non-empty core:
Shapley costs is less
than standalone
members’ costs
16 Coalition → {12} {13} {23}
(continued)
9 Allocating Cost Savings (Gains) Between Cooperating… 307
Physician
# Participant (member) → Hospital {1} group {2} SNF {3} Total
17 Condition 2 for OK OK OK
non-empty core:
two-member coalition
costs are less than the
combined members’
costs
All cores are non-empty
where C1, C2, and C3 are the costs for standalone providers, respectively;
C12, C13, and C23 are the costs for two-member coalitions, respectively;
and C123 is the cost for the three-member coalition.
For s = 1 the weight w1 = (1 − 1)!(3 − 1)!/3! = 1/3; for s = 2 the weight
w2 = (2 − 1)!(3 − 2)!/3! = 1/6; and for s = 3 the weight w3 = (3 − 1)! × (3
− 3)!/3! = 1/3.
These formulas were built in the Excel spreadsheet (the spreadsheet file
is available in SEM 9.1).
308 A. Kolker
Problem Description
CMS’s comprehensive care for joint replacement (CJR) model went into
effect on April 1, 2016 (Federal Register, 2015). Under this model, acute
9 Allocating Cost Savings (Gains) Between Cooperating… 309
Table 9.2 Standalone costs and the Shapley value costs allocation for group
size n = 4
line
# Hospital {1} IPR {2} SNF {3} HHA {4} Total
1 $16,023 $1568 $5034 $2849 $25,474
2 {12} {13} {14} {23} {24} {34}
3 3% 4% 1% 2% 3% 2%
4 $17,063 $20,215 $18,683 $6503 $4284 $7725
5 {123} {124} {134} {234}
6 3% 3% 4% 2%
7 $22,059 $19,827 $22,950 $9262
8 {1234}
9 5.0%
10 $24,200
11 Shapley Shapley Shapley Shapley
value for value for value for value for
provider 1, provider 2, provider 3, provider 4,
H{1} IPR{2} SNF{3} HHA{4}
12 $15,476 $1410 $4705 $2608 $24,200
13 Saving {1} Saving {2} Saving {3} Saving {4} Total
savings
14 $547 $158 $329 $241 $1274
15 Saving Saving Saving Saving
allocation allocation allocation allocation
% {1} % {2} % {3} % {4}
16 3.4% 10.0% 6.5% 8.5%
17 {12} {23} {34} {13} {14} {24}
18 OK OK OK OK OK OK
19 {123} {124} {134} {234}
20 OK OK OK OK
All cores are not empty
The Shapley values for all four providers, their savings, and the percent
of savings allocations are indicated in lines 11–16. The cores for all coali-
tions are not empty. All four providers have financial incentives to stay in
coalitions because their standalone costs would be higher. The Shapley
values allocate higher percentage savings to providers 2 (IPR) and 4
(HHA) (line 16) than simple cost-proportional allocations. This will
likely encourage their cooperation with the hospital/physicians and SNF
9 Allocating Cost Savings (Gains) Between Cooperating… 311
Problem Description
The Shapley value formula is difficult to use when the coalition size
increases above 4. The reason is that the number of permutations increases
exponentially with the number of participants, and the cost data are
needed for all k-member coalitions (k = 2, 3, 4, …, n). These multiple-
cost data are difficult to get, or they are not known. Therefore, some
special cases were developed in which simplified or approximate expres-
sions of the Shapley value are used (Reinhardt & Dada, 2005). In par-
ticular, there are widespread situations in which participants need to use
only a portion of the available total capacity. The total capacity is required
to serve only the biggest participant. The users of the portions of the total
capacity want to pay only their corresponding share of the biggest user’s
cost. One can think of toll roads, computer networks, or shared health-
care facilities.
One such practically important case is historically called an “airport
game” developed originally for setting landing fees for different types of
aircraft (small and big) that share the same runway. The runway length
(and the cost) is defined essentially by the biggest type of aircraft that
requires the full runway length, while smaller-type aircraft can use only a
portion of the total runway length. Littlechild and Owen (1973) first
demonstrated that an empirical cost allocation rule for landing fees is the
Shapley value in this particular case.
C1 < C2 < C3 < … < Cm, where m is the number of types of players.
9 Allocating Cost Savings (Gains) Between Cooperating… 313
Thus, the pattern is that the Shapley value for each player is equal to that
for the previous one plus the corresponding price differential divided by
the total number of players but the sum of all previous smaller ones.
If only one player is present in each cost category, i.e., if
n1 = n2 = n3 = …1 then
Sh1 C1 / N
Sh2 Sh1 C2 C1 / N 1
Sh3 Sh2 C3 C2 / N 2
Sh4 Sh3 C4 C3 / N 3
Shm Shm 1 Cm Cm 1 / N m 1
Four Hospitals
Four hospitals are considering offering imaging services at a new Center
for Diagnostic Imaging (CDI). Annual startup and operating costs, as
314 A. Kolker
Table 9.4 Annual startup and operating costs for new CDI
Hospital 1 Hospital 2 Hospital 3 Hospital 4
Annual patient volumes 1000 950 800 1800
Estimated annual startup costs $5 M $4 M $2 M $9 M
Step 1.
Step 2.
Step 3.
Table 9.5 Dental clinic comparison chart: capacity and costs for fixed clinics and
clinics using portable equipment
Size 3 chairs 6 chairs 9 chairs 12 chairs
Number of dentists 1 2 3 4
Total estimated startup costs $624,452 $1,016,217 $1,268,330 $1,549,538
316 A. Kolker
should be realized that there is no one best “fair” cost allocation method.
The choice of a method depends on the participants’ ideas about “fair-
ness,” their economic and political views, and other subjective factors.
Moreover, the ease of calculation and interpretation play a significant role
for the practical use.
The attractive features of the Shapley value framework make it valuable
for various applications, especially in healthcare settings where providers’
cooperation and team-based care are currently a long-term strategy for
improving quality of care and cost reduction.
At the same time, the Shapley value’s axiomatic basis has some short-
comings (Tijs & Driessen, 1986; Young, 1994). The Shapley value for-
mula is difficult to use when the group size n is greater than 4 because the
number of double, triple, etc. permutations, P, increases exponentially
with the group size as P = 2n − (n + 1). For instance, for a modest group
size n = 10, more than 1000 subgroups and their marginal costs should
be considered. This is impractical. Although some special approximations
were developed, such as an airport game framework, further develop-
ment of other Shapley value approximations is needed for large group
sizes. Specifically, an approximation is needed that does not necessarily
require evaluating every possible coalition using, for example, the con-
straint generation technique.
Another problem is the possibility that some players may be more
likely to act together than others. Even though coalition structures arise
naturally in cost allocation problems, their applications have not been
explored sufficiently. These are some fruitful areas for further develop-
ment of the Shapley value concept.
There is one more problem: to apply the Shapley value, we need to
know what everyone brings to any coalition. Yet, for instance, if people
collaborate, we cannot observe what happens without the collaboration.
We could ask participants what they could do on their own, but they’d
have incentives to lie: this is an unreliable approach.
For instance, two hospitals are building their new imaging centers.
Construction costs for standalone centers are 15 M and 8 M, respectively.
If the hospitals cooperate and build one shared center, the cost will be
19 M. Hence, the Shapley values for each hospital will be:
318 A. Kolker
Sh1 = 1/2 × (15 + (19 − 8)) = $13 M and Sh2 = 1/2 × (8 + (19 − 1
5)) = $6 M.
However, the second hospital could say: “We re-assessed, and actually
our stand-alone cost will be lower, only 6 M.” Hence, now the first hos-
pital Sh1 = 1/2 × (15 + (19-6)) = $14 M, i.e., greater than the original
contribution of $13 M. The second hospital Sh2 = 1/2 × (6 + (19-15))
= $5 M, i.e., less than the original $6 M. Clearly, there is an incentive for
the second hospital to lower its own standalone cost. One has to be mind-
ful of this potential issue.
References
Barton, T. L. (1992). A unique solution for the nucleolus in accounting alloca-
tions. Decision Sciences, 23(2), 365–375.
Dyrda, L. (2016, March 30). CJR bundles to pay $25k per episode: 8 statistics on
cost breakdown. Becker’s Spine Review. [Link]
orthopedic-spine-practices-improving-profits/item/30556-cjr-bundles-to-
pay-25k-per-episode-8-statistics-on-cost-[Link]
Federal Register (2015). Medicare Program; Comprehensive Care for Joint
Replacement Payment Model for Acute Care Hospitals Furnishing Lower
Extremity Joint Replacement Services. [Link]
2015-29438, [Link] gov/a/2015-29438
Klenberg, N., & Weiss, J. (2013). On membership and marginal values.
International Journal of Game Theory, 42(2), 357–373.
Kolker, A. (2017). Chapter 182: The concept of the shapley value and the cost
allocation between cooperating participants. In Encyclopedia of information
science and technology (4th ed., pp. 2095–2107). IGI-Global. [Link]
igi-g [Link]/chapter/the-c oncept-o f-t he-s hapley-value-a nd-t he-c ost-
allocation-between-cooperating-participants/183923
Leng, M., & Parlar, M. (2010). Analytic solution for the nucleolus of a three-
player cooperative game. Naval Research Logistics, 57, 667–672.
Littlechild, S., & Owen, G. (1973). A simple expression for the Shapley value in
a special case. Management Science, 20(3), 370–372.
9 Allocating Cost Savings (Gains) Between Cooperating… 319
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2024 321
A. Kolker, Healthcare Management Engineering In Action, Business Guides on the Go,
[Link]
322 A. Kolker
–– Corollary 3. Including too old data points to fit the past pattern
usually results in overfitting that is detrimental to the forecast accu-
racy (illustrated in Sect. 8.3).
Implementation Principle
• For low disease prevalence (less than 30%), pooled screening specimen
testing is more efficient than individual specimen testing (illustrated in
Sect. 4.3).
• The lower the prevalence (expectation) of a disease or infection, the
more likely a positive test would be the false positive and a negative
test would be the true negative (illustrated in Sect. 8.4).
–– Corollary 1. Mass screening of the general population with low dis-
ease (infection) prevalence (expectation) aimed at its early detection
(irrespective of the risk status of the individuals) is not an effective
strategy. A positive test should be suspicious and repeated.
–– Corollary 2. Mass screening of the selected subpopulations with
moderate to high disease (infection) prevalence (expectation) could
be an effective way of detecting the disease. A positive test would
most likely be a true positive. A negative test would likely be a false
negative and should be repeated.
• The value of a diagnostic test is determined by how much the test
affects the patient care decision.
–– Corollary 1. If a test is unlikely to change the clinical decision that
was already made, then the test should not be ordered.
326 A. Kolker
Concluding Remarks
The fascinating journey into the world of healthcare management engi-
neering has come to an end. A lot of topics were covered with various
depths, from beds and operating rooms’ capacity and patient flow to
optimal staffing and scheduling, from hospital admission policy to data
envelopment analysis, from the minimal cost distribution network to
Bayesian inference for diagnostic testing, from patient census forecasting
to principal component decomposition of the multivariable dataset, from
cluster analysis to the Shapley values. These topics were illustrated with
multiple examples taken or adapted from healthcare practice.
Where do we go from here? How will it all help us? According to the
seminal novel Fahrenheit 451 by famous novelist Ray Bradbury (1953,
page 84–85): “Only if the three necessary things could be given: (i) qual-
ity of information, (ii) leisure to digest it, and (iii) the right to carry out
actions based on what we learn from the interactions of the first two.”
I hope, as the author of this book, that I contributed to the first two of
the above. However, the ultimate success of the management engineering
mission depends on the third one.
In the spirit of the above quotation, let us summarize the main points
of this Business Guide on the Go Series book. The main goal of the book
was helping to bridge the gap in mutual understanding and communica-
tion between management engineering professionals and hospital and
clinic administrators empowered to make managerial decisions, i.e., the
right to carry out actions based on what they learned from this book. In
this regard, it was highlighted the role of the key hospital stakeholders
with veto power. It was also emphasized that a Pareto efficient manage-
ment solution has the highest chance of being implemented.
This book was intended primarily for hospital/clinic leadership to con-
vince them that management engineering methodology is worth a slog of
getting through, that it is an indispensable aid in discovering some hid-
den interconnections in the maze of hospital operations making it
10 Summary of Fundamental Management Principles… 327
References
Bradbury, R. (1953). Fahrenheit 451 (50th Anniv. ed.). The Random House
Publishing Group.
Hopp, W., & Lovejoy, W. (2013). Hospital operations: Principles of high efficiency
health care (p. 623). FT Press.
Index
© The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer 329
Nature Switzerland AG 2024
A. Kolker, Healthcare Management Engineering In Action, Business Guides on the Go,
[Link]
330 Index
D H
Data envelopment analysis (DEA) Healthcare system
multiple inputs and complex, 4, 16, 17
outputs, 230–250 Hospital
scoring function, 231, 233 large, 53–55, 323
target for improvement, 238, small, 53–56, 323
243, 247
Data science (DS), 2, 3, 20
Diagnostic lab, xii, 179–182 L
Discrete event simulation (DES) Linear optimization
non-steady state, 24, constraints, 146, 162–164,
36–39, 41, 45 168, 171, 173, 179, 182,
steady state, 15, 18, 24, 30, 31, 188, 191, 197, 199, 203,
34, 36, 38, 39, 42, 48, 49, 207, 212–214, 218, 219,
51, 55, 75 222, 225
warm-up, 24, 35, 37, 44, 45, 48, decision variables, 94, 146,
49, 55, 75 162, 164, 166, 169, 170,
Distribution network 173, 179, 185–187, 197,
restricted routes, 190, 192 202, 204, 212, 213, 218,
unrestricted routes, 193 219, 222
general for resource
allocation, 162
E index table, 198, 199, 206,
Emergency department (ED) 214
discharge from, 73, 294 integer for staffing and
Excel solver, 164, 167, 168, 174, scheduling, 196–225
177, 179, 180, 186, 188, 199, objective function, 146, 162–164,
207, 214, 222, 225, 233, 166–171, 173, 174, 176,
235–238, 241, 249, 250 179, 185, 197, 199, 202,
204, 207, 212–214, 218,
219, 234
F reduced cost, 169, 170, 177, 181,
Flaw of averages, 49, 71, 108, 188, 194
114, 163 sensitivity report, 166, 168–172,
Forecasting 175, 177, 178, 180–182, 187,
patient census, 326 189, 192, 204
Fundamental management shadow price, 169, 171, 172, 178,
principles, 80, 292, 321–327 182, 188, 189, 192
Index 331