Open Access
Estimation of patient flow in hospitals using up-to-date data: application to bed demand prediction during pandemic waves
(Public Library of Science, 2023) García de Vicuña Bilbao, Daniel; López-Cheda, Ana; Jácome, María Amalia; Mallor Giménez, Fermín; Institute of Smart Cities - ISC
Hospital bed demand forecast is a first-order concern for public health action to avoid healthcare systems to be overwhelmed. Predictions are usually performed by estimating patients flow, that is, lengths of stay and branching probabilities. In most approaches in the literature, estimations rely on not updated published information or historical data. This may lead to unreliable estimates and biased forecasts during new or non-stationary situations. In this paper, we introduce a flexible adaptive procedure using only near-real-time information. Such method requires handling censored information from patients still in hospital. This approach allows the efficient estimation of the distributions of lengths of stay and probabilities used to represent the patient pathways. This is very relevant at the first stages of a pandemic, when there is much uncertainty and too few patients have completely observed pathways. Furthermore, the performance of the proposed method is assessed in an extensive simulation study in which the patient flow in a hospital during a pandemic wave is modelled. We further discuss the advantages and limitations of the method, as well as potential extensions.
Open Access
Coping with stress in emergency department physicians through improved patient-flow management
(Elsevier, 2020) Cildoz Esquíroz, Marta; Ibarra, Amaia; Mallor Giménez, Fermín; Estatistika, Informatika eta Matematika; Institute of Smart Cities - ISC; Estadística, Informática y Matemáticas
This paper provides a method for the real-time monitoring of job stress in emergency department (ED) physicians. It is implemented in a Decision Support System (DSS) designed for patient-to-physician assignment after triage. Our concept of job stress includes not only the workload but also time pressure and uncertainty. A job stress function is estimated based on the consensus views of ED physicians obtained through a novel methodology involving stress factor analysis, questionnaire design, and the statistical analysis of expert opinions. The resulting stress score enables the assessment of job stress using workload data from the ED physicians’ whiteboard. These data can be used for the real-time measurement and monitoring of ED physician job stress in a stochastic and dynamic environment, which is the main novelty of this method as compared to previous workload and stress measurement proposals. A further advantage of this methodology is that it is general enough to be adapted to physician job stress monitoring in any ED. The use of the DSS for ED patient-flow management reduces job stress and spreads it more evenly among the whole team of physicians, while also improving other important ED performance measures such as arrival-to-provider time and the percentage of compliance with patient waiting time targets. A case study illustrates the application of the methodology for the construction of a stress-score, the monitoring of physician stress levels, and ED patient-flow management.
Open Access
A GRASP-based algorithm for solving the emergency room physician scheduling problem
(Elsevier, 2021) Cildoz Esquíroz, Marta; Mallor Giménez, Fermín; Mateo, Pedro; Institute of Smart Cities - ISC
This paper addresses a physician scheduling problem in an Emergency Room (ER) requiring a long-term work calendar to allocate work days and types of shift among all the doctors. The mathematical model is created without simplifications, using the real calendar, including holidays. This precludes the possibility of cyclic-type solutions, and involves numerous and varied constraints (demand, workload, ergonomics, fairness, etc.). An effective solution to this very difficult practical problem cannot be obtained, for large instances, with exact solution methods. We formulate a mathematical representation of a real-world ER physician scheduling problem featuring a hybrid algorithm combining continuous linear programming with a greedy randomized adaptive search procedure (GRASP). Linear programming is used to model a general physician-demand covering problem, where the solution is used to guide the construction phase of the GRASP, to obtain initial full schedules for subsequent improvement by iterative application of Variable Neighborhood Descent Search (VNDS) and Network Flow Optimization (NFO). A computational study shows the superiority of our approach over the Integer Linear Programming method in a set of instances of varying size and difficulty inspired by a real setting. The methodology is embedded in a software tool for generating one-year-ahead physician schedules for a local ER. These solutions, which are now in use, outperform the manually-created schedules used previously. © 2021 Elsevier B.V.
Open Access
The problem of the last bed: contextualization and a new simulation framework for analyzing physician decisions
(Elsevier, 2019) Azcárate Camio, Cristina; Esparza, Laida; Mallor Giménez, Fermín; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika
Faced with a full Intensive Care Unit (ICU), physicians need to decide between turning away a new patient in need of critical care and creating a vacancy by prematurely discharging a current occupant. This dilemma is widely discussed in the medical literature, where the influencing factors are identified, the patient discharge process described and the patient health consequences analyzed. Nevertheless, the existing mathematical models of ICU management practices overlook many of the factors considered by physicians in real-world triage decisions. This paper offers a review of the medical and mathematical literature on patient discharge decisions, and a proposal for a new simulation framework to enable more realistic mathematical modeling of the real-world patient discharge process. Our model includes a) the times at which discharge decisions are made and setup times for patient transfer from the ICU to a general ward and preparation of an ICU bed for an incoming patient, in order to capture the impossibility of an immediate switch of patients; b) advance notice of the number of patients due to arrive from elective surgery requiring intensive postoperative care and potentially triggering the need for early discharges to avoid surgery cancelations; and c) patient health status (to reflect the dependency of physicians’ discharge decisions on health indicators) by modeling length of stay with a phase-type distribution in which a medical meaning is assigned to each state. A simulation-based optimization method is also proposed as a means to obtain optimal discharge decisions as a function of the health status of current patients, the bed occupancy level and the number of planned arrivals from elective surgery over the following days. Optimal decisions should strike a balance between patient rejection and LoS reduction. This new simulation framework generates an optimal discharge policy, which closely resembles real decision-making under a cautious discharge policy, where the frequency of early discharge increases with the ICU occupancy level. This is a contrast with previous simulation models, which consider only the triage of the last bed, disregarding the pressures on physicians faced with high bed occupancy levels.
Open Access
Cálculo de la distribución del tiempo de vida de componentes mediante autopsia en sistemas binarios aditivos, serie-paralelo y paralelo-serie
(Universitat Politècnica de Catalunya, 1997) Mallor Giménez, Fermín; Azcárate Camio, Cristina; Pérez Prados, Antonio; Estadística e Investigación Operativa; Estatistika eta Ikerketa Operatiboa
En este artículo se estudia el problema de determinar la función de distribución del tiempo de vida de las componentes de un sistema binario, a partir del conocimiento de las leyes que rigen el funcionamiento del sistema y del conjunto de componentes que causa su fallo (obtenida mediante autopsia del sistema en el momento de su deterioro). Se presentan los resultados de Meilijson (1981) y Nowik (1990) que proponen un sistema de ecuaciones impíıcito para obtener estas distribuciones. Sin embargo, se observa que este sistema es de muy difícil resolución práctica, por lo que nosotros consideramos un método cuya utilización es más restringida pero más sencilla, y estudiamos su aplicación a sistemas binarios aditivos, serie-paralelo y paralelo-serie.
Open Access
Dynamic mean absolute error as new measure for assessing forecasting errors
(Elsevier, 2018-02-14) Frías Paredes, Laura; Mallor Giménez, Fermín; Gastón Romeo, Martín; León, Teresa; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute of Smart Cities - ISC
Accurate wind power forecast is essential for grid integration, system planning, and electricity trading in certain electricity markets. Therefore, analyzing prediction errors is a critical task that allows a comparison of prediction models and the selection of the most suitable model. In this work, the temporal error and absolute magnitude error are simultaneously considered to assess the forecast error. The trade-off between both types of errors is computed, analyzed, and interpreted. Moreover, a new index, the dynamic mean absolute error, DMAE, is defined to measure the prediction accuracy. This index accounts for both error components: temporal and absolute. Real cases of wind energy forecasting are used to illustrate the use of the new DMAE index and show the advantages of this new index over other error indices.
Open Access
Motor unit action potential duration, II: a new automatic measurement method based on the wavelet transform
(Lippincott, Williams & Wilkins, 2007) Rodríguez Carreño, Ignacio; Gila Useros, Luis; Malanda Trigueros, Armando; García Gurtubay, Ignacio; Mallor Giménez, Fermín; Gómez Elvira, Sagrario; Rodríguez Falces, Javier; Navallas Irujo, Javier; Ingeniería Eléctrica y Electrónica; Estadística e Investigación Operativa; Ingeniaritza Elektrikoa eta Elektronikoa; Estatistika eta Ikerketa Operatiboa
To present and evaluate a new algorithm, based on the wavelet transform, for the automatic measurement of motor unit action potential (MUAP) duration. A total of 240 MUAPs were studied. The waveform of each MUAP was wavelet-transformed, and the start and end points were estimated by regarding the maxima and minima points in a particular scale of the wavelet transform. The results of the new method were compared with the gold standard of duration marker positions obtained by manual measurement. The new method was also compared with a conventional algorithm, which we had found to be best in a previous comparative study. To evaluate the new method against manual measurements, the dispersion of automatic and manual duration markers were analyzed in a set of 19 repeatedly recorded MUAPs. The differences between the new algorithm’s marker positions and the gold standard of duration marker positions were smaller than those observed with the conventional method. The dispersion of the new algorithm’s marker positions was slightly less than that of the manual one. Our new method for automatic measurement of MUAP duration is more accurate than other available algorithms and more consistent than manual measurements.
Open Access
Assessing the impact of physicians' behavior variability on performance indicators in emergency departments: an agent-based model
(IEEE, 2025-01-20) Baigorri Iguzquiaguirre, Miguel; Cildoz Esquíroz, Marta; Mallor Giménez, Fermín; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute of Smart Cities - ISC
In emergency departments (EDs), traditional simulation models often overlook the variability in physician practice, assuming uniform service provision. Our study introduces a hybrid agent-based discrete-event simulation (AB-DES) model to capture this variability. Through simulation scenarios based on real ED data, we assess the impact of physician behavior on key performance indicators such as patient waiting times and physician stress levels. Results show significant variability in both individual physician performance and average metrics across scenarios. By integrating physician agent modeling, informed by literature from medical and workplace psychology, our approach offers a more nuanced representation of ED dynamics. This model serves as a foundation for future developments towards digital twins, facilitating real-time ED management. Our findings emphasize the importance of considering physician behavior for accurate performance assessment and optimization.
Open Access
Acuity-based rotational patient-to-physician assignment in an emergency department using electronic health records in triage
(SAGE, 2023) Cildoz Esquíroz, Marta; Ibarra Bolt, Amaya; Mallor Giménez, Fermín; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute of Smart Cities - ISC
Emergency department (ED) operational metrics generated by a new acuity-based rotational patient-to-physician assignment (ARPA) algorithm are compared with those obtained with a simple rotational patient assignment (SRPA) system aimed only at an equitable patient distribution. The new ARPA method theoretically guarantees that no two physicians’ assigned patient loads can differ by more than one, either partially (by acuity levels) or in total; whereas SRPA guarantees only the latter. The performance of the ARPA method was assessed in practice in the ED of the main public hospital (Hospital Compound of Navarra) in the region of Navarre in Spain. This ED attends over 140 000 patients every year. Data analysis was conducted on 9,063 ED patients in the SRPA cohort, and 8,892 ED patients in the ARPA cohort. The metrics of interest are related both to patient access to healthcare and physician workload distribution: patient length of stay; arrival-to-provider time; ratio of patients exceeding the APT target threshold; and range of assigned patients across physicians by priority levels. The transition from SRPA to ARPA is associated with improvements in all ED operational metrics. This research demonstrates that ARPA is a simple and useful strategy for redesigning front-end ED processes.
Open Access
Early detection of new pandemic waves: control chart and a new surveillance index
(Public Library of Science, 2024) Cildoz Esquíroz, Marta; Gastón Romeo, Martín; Frías Paredes, Laura; García de Vicuña Bilbao, Daniel; Azcárate Camio, Cristina; Mallor Giménez, Fermín; Institute of Smart Cities - ISC
The COVID-19 pandemic highlights the pressing need for constant surveillance, updating of the response plan in post-peak periods and readiness for the possibility of new waves of the pandemic. A short initial period of steady rise in the number of new cases is sometimes followed by one of exponential growth. Systematic public health surveillance of the pandemic should signal an alert in the event of change in epidemic activity within the community to inform public health policy makers of the need to control a potential outbreak. The goal of this study is to improve infectious disease surveillance by complementing standardized metrics with a new surveillance metric to overcome some of their difficulties in capturing the changing dynamics of the pandemic. At statistically-founded threshold values, the new measure will trigger alert signals giving early warning of the onset of a new pandemic wave. We define a new index, the weighted cumulative incidence index, based on the daily new-case count. We model the infection spread rate at two levels, inside and outside homes, which explains the overdispersion observed in the data. The seasonal component of real data, due to the public surveillance system, is incorporated into the statistical analysis. Probabilistic analysis enables the construction of a Control Chart for monitoring index variability and setting automatic alert thresholds for new pandemic waves. Both the new index and the control chart have been implemented with the aid of a computational tool developed in R, and used daily by the Navarre Government (Spain) for virus propagation surveillance during post-peak periods. Automated monitoring generates daily reports showing the areas whose control charts issue an alert. The new index reacts sooner to data trend changes preluding new pandemic waves, than the standard surveillance index based on the 14-day notification rate of reported COVID-19 cases per 100,000 population.