scholarly journals Mathematical Modeling and Robustness Analysis to Unravel COVID-19 Transmission Dynamics: The Italy Case

Biology ◽  
2020 ◽  
Vol 9 (11) ◽  
pp. 394
Author(s):  
Chiara Antonini ◽  
Sara Calandrini ◽  
Fabrizio Stracci ◽  
Claudio Dario ◽  
Fortunato Bianconi

This study started from the request of providing predictions on hospitalization and Intensive Care Unit (ICU) rates that are caused by COVID-19 for the Umbria region in Italy. To this purpose, we propose the application of a computational framework to a SEIR-type (Susceptible, Exposed, Infected, Removed) epidemiological model describing the different stages of COVID-19 infection. The model discriminates between asymptomatic and symptomatic cases and it takes into account possible intervention measures in order to reduce the probability of transmission. As case studies, we analyze not only the epidemic situation in Umbria but also in Italy, in order to capture the evolution of the pandemic at a national level. First of all, we estimate model parameters through a Bayesian calibration method, called Conditional Robust Calibration (CRC), while using the official COVID-19 data of the Italian Civil Protection. Subsequently, Conditional Robustness Analysis (CRA) on the calibrated model is carried out in order to quantify the influence of epidemiological and intervention parameters on the hospitalization rates. The proposed pipeline properly describes the COVID-19 spread during the lock-down phase. It also reveals the underestimation of new positive cases and the need of promptly isolating asymptomatic and presymptomatic cases. The results emphasize the importance of the lock-down timeliness and provide accurate predictions on the current evolution of the pandemic.

2017 ◽  
Vol 140 (1) ◽  
Author(s):  
Na Qiu ◽  
Chanyoung Park ◽  
Yunkai Gao ◽  
Jianguang Fang ◽  
Guangyong Sun ◽  
...  

In calibrating model parameters, it is important to include the model discrepancy term in order to capture missing physics in simulation, which can result from numerical, measurement, and modeling errors. Ignoring the discrepancy may lead to biased calibration parameters and predictions, even with an increasing number of observations. In this paper, a simple yet efficient calibration method is proposed based on sensitivity information when the simulation model has a model error and/or numerical error but only a small number of observations are available. The sensitivity-based calibration method captures the trend of observation data by matching the slope of simulation predictions and observations at different designs and then utilizing a constant value to compensate for the model discrepancy. The sensitivity-based calibration is compared with the conventional least squares calibration method and Bayesian calibration method in terms of parameter estimation and model prediction accuracies. A cantilever beam example, as well as a honeycomb tube crush example, is used to illustrate the calibration process of these three methods. It turned out that the sensitivity-based method has a similar performance with the Bayesian calibration method and performs much better than the conventional method in parameter estimation and prediction accuracy.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Christopher Martin ◽  
Stuart McDonald ◽  
Steve Bale ◽  
Michiel Luteijn ◽  
Rahul Sarkar

Abstract Background This paper describes a model for estimating COVID-19 related excess deaths that are a direct consequence of insufficient hospital ward bed and intensive care unit (ICU) capacity. Methods Compartmental models were used to estimate deaths under different combinations of ICU and ward care required and received in England up to late April 2021. Model parameters were sourced from publicly available government information and organisations collating COVID-19 data. A sub-model was used to estimate the mortality scalars that represent increased mortality due to insufficient ICU or general ward bed capacity. Three illustrative scenarios for admissions numbers, ‘Optimistic’, ‘Middling’ and ‘Pessimistic’, were modelled and compared with the subsequent observations to the 3rd February. Results The key output was the demand and capacity model described. There were no excess deaths from a lack of capacity in the ‘Optimistic’ scenario. Several of the ‘Middling’ scenario applications resulted in excess deaths—up to 597 deaths (0.6% increase) with a 20% reduction compared to best estimate ICU capacity. All the ‘Pessimistic’ scenario applications resulted in excess deaths, ranging from 49,178 (17.0% increase) for a 20% increase in ward bed availability, to 103,735 (35.8% increase) for a 20% shortfall in ward bed availability. These scenarios took no account of the emergence of the new, more transmissible, variant of concern (b.1.1.7). Conclusions Mortality is increased when hospital demand exceeds available capacity. No excess deaths from breaching capacity would be expected under the ‘Optimistic’ scenario. The ‘Middling’ scenario could result in some excess deaths—up to a 0.7% increase relative to the total number of deaths. The ‘Pessimistic’ scenario would have resulted in significant excess deaths. Our sensitivity analysis indicated a range between 49,178 (17% increase) and 103,735 (35.8% increase). Given the new variant, the pessimistic scenario appeared increasingly likely and could have resulted in a substantial increase in the number of COVID-19 deaths. In the event, it would appear that capacity was not breached at any stage at a national level with no excess deaths. it will remain unclear if minor local capacity breaches resulted in any small number of excess deaths.


Author(s):  
Christopher J. Arthurs ◽  
Nan Xiao ◽  
Philippe Moireau ◽  
Tobias Schaeffter ◽  
C. Alberto Figueroa

AbstractA major challenge in constructing three dimensional patient specific hemodynamic models is the calibration of model parameters to match patient data on flow, pressure, wall motion, etc. acquired in the clinic. Current workflows are manual and time-consuming. This work presents a flexible computational framework for model parameter estimation in cardiovascular flows that relies on the following fundamental contributions. (i) A Reduced-Order Unscented Kalman Filter (ROUKF) model for data assimilation for wall material and simple lumped parameter network (LPN) boundary condition model parameters. (ii) A constrained least squares augmentation (ROUKF-CLS) for more complex LPNs. (iii) A “Netlist” implementation, supporting easy filtering of parameters in such complex LPNs. The ROUKF algorithm is demonstrated using non-invasive patient-specific data on anatomy, flow and pressure from a healthy volunteer. The ROUKF-CLS algorithm is demonstrated using synthetic data on a coronary LPN. The methods described in this paper have been implemented as part of the CRIMSON hemodynamics software package.


2021 ◽  
Vol 13 (1) ◽  
pp. 143-158
Author(s):  
Tomasz Zwęgliński

Poles are highly aware of the fact that Polish civil protection assets are being deployed abroad in order to assist other disaster- and crisis-stricken countries around the world. Such a type of urgent assistance provided from one country to another in an emergency response situation is regulated and organised by the European Union. Poland, as a state participating in the EU international system, is very active in sharing its civil protection assets, such as in the Beirut explosion case in Lebanon (2020), and during forest fires in Sweden (2018). Using its civil protection resources to assist other countries poses a question on the potential influence of such activities on the homeland’s internal security. Solving the problem has to be preceded by answering the following question — How is the process of international civil protection deployment organised in Poland? Responding to this question required utilising such methods as a review and analysis of literature sources on the European dimension of civil protection, EU legal regulations, statistics, reports as well as Polish standard operating procedures and internal regulations on the national level and was the key method applied in the research. Furthermore, semi-formal interviews with Polish and EU experts were done. The findings prove that operational planning in the researched area is well organised, which significantly diminishes the identified risk for internal security.


2017 ◽  
Vol 14 (18) ◽  
pp. 4295-4314 ◽  
Author(s):  
Dan Lu ◽  
Daniel Ricciuto ◽  
Anthony Walker ◽  
Cosmin Safta ◽  
William Munger

Abstract. Calibration of terrestrial ecosystem models is important but challenging. Bayesian inference implemented by Markov chain Monte Carlo (MCMC) sampling provides a comprehensive framework to estimate model parameters and associated uncertainties using their posterior distributions. The effectiveness and efficiency of the method strongly depend on the MCMC algorithm used. In this work, a differential evolution adaptive Metropolis (DREAM) algorithm is used to estimate posterior distributions of 21 parameters for the data assimilation linked ecosystem carbon (DALEC) model using 14 years of daily net ecosystem exchange data collected at the Harvard Forest Environmental Measurement Site eddy-flux tower. The calibration of DREAM results in a better model fit and predictive performance compared to the popular adaptive Metropolis (AM) scheme. Moreover, DREAM indicates that two parameters controlling autumn phenology have multiple modes in their posterior distributions while AM only identifies one mode. The application suggests that DREAM is very suitable to calibrate complex terrestrial ecosystem models, where the uncertain parameter size is usually large and existence of local optima is always a concern. In addition, this effort justifies the assumptions of the error model used in Bayesian calibration according to the residual analysis. The result indicates that a heteroscedastic, correlated, Gaussian error model is appropriate for the problem, and the consequent constructed likelihood function can alleviate the underestimation of parameter uncertainty that is usually caused by using uncorrelated error models.


2021 ◽  
Author(s):  
Sheng Zhang ◽  
Joan Ponce ◽  
Zhen Zhang ◽  
Guang Lin ◽  
George Karniadakis

AbstractEpidemiological models can provide the dynamic evolution of a pandemic but they are based on many assumptions and parameters that have to be adjusted over the time when the pandemic lasts. However, often the available data are not sufficient to identify the model parameters and hence infer the unobserved dynamics. Here, we develop a general framework for building a trustworthy data-driven epidemiological model, consisting of a workflow that integrates data acquisition and event timeline, model development, identifiability analysis, sensitivity analysis, model calibration, model robustness analysis, and forecasting with uncertainties in different scenarios. In particular, we apply this framework to propose a modified susceptible–exposed–infectious–recovered (SEIR) model, including new compartments and model vaccination in order to forecast the transmission dynamics of COVID-19 in New York City (NYC). We find that we can uniquely estimate the model parameters and accurately predict the daily new infection cases, hospitalizations, and deaths, in agreement with the available data from NYC’s government’s website. In addition, we employ the calibrated data-driven model to study the effects of vaccination and timing of reopening indoor dining in NYC.


Author(s):  
B. Sandeep Reddy ◽  
Ashitava Ghosal

This paper deals with the issue of robustness in control of robots using the proportional plus derivative (PD) controller and the augmented PD controller. In the literature, a variety of PD and model-based controllers for multilink serial manipulator have been claimed to be asymptotically stable for trajectory tracking, in the sense of Lyapunov, as long as the controller gains are positive. In this paper, we first establish that for simple PD controllers, the criteria of positive controller gains are insufficient to establish asymptotic stability, and second that for the augmented PD controller the criteria of positive controller gains are valid only when there is no uncertainty in the model parameters. We show both these results for a simple planar two-degrees-of-freedom (2DOFs) robot with two rotary (R) joints, following a desired periodic trajectory, using the Floquet theory. We provide numerical simulation results which conclusively demonstrate the same.


Sensors ◽  
2020 ◽  
Vol 20 (20) ◽  
pp. 5934
Author(s):  
Xiao Li ◽  
Wei Li ◽  
Xin’an Yuan ◽  
Xiaokang Yin ◽  
Xin Ma

Lens distortion is closely related to the spatial position of depth of field (DoF), especially in close-range photography. The accurate characterization and precise calibration of DoF-dependent distortion are very important to improve the accuracy of close-range vision measurements. In this paper, to meet the need of short-distance and small-focal-length photography, a DoF-dependent and equal-partition based lens distortion modeling and calibration method is proposed. Firstly, considering the direction along the optical axis, a DoF-dependent yet focusing-state-independent distortion model is proposed. By this method, manual adjustment of the focus and zoom rings is avoided, thus eliminating human errors. Secondly, considering the direction perpendicular to the optical axis, to solve the problem of insufficient distortion representations caused by using only one set of coefficients, a 2D-to-3D equal-increment partitioning method for lens distortion is proposed. Accurate characterization of DoF-dependent distortion is thus realized by fusing the distortion partitioning method and the DoF distortion model. Lastly, a calibration control field is designed. After extracting line segments within a partition, the de-coupling calibration of distortion parameters and other camera model parameters is realized. Experiment results shows that the maximum/average projection and angular reconstruction errors of equal-increment partition based DoF distortion model are 0.11 pixels/0.05 pixels and 0.013°/0.011°, respectively. This demonstrates the validity of the lens distortion model and calibration method proposed in this paper.


2020 ◽  
Vol 139 (3) ◽  
pp. 54-66
Author(s):  
Joanna Grzela

In Iceland, since 2008, there has been a centralisation and coordination of forces and resources needed to protect the population and civil security. Duties in the fi eld of civil protection at the national level are delegated to the National Police Chief, who acts in accordance with the government’s policy of civil protection and security. The system’s reliability has been achieved thanks to the close cooperation of offi cials from many levels of government and self-government administration, which is overseen by the National Police Chief.


2018 ◽  
Vol 2018 ◽  
pp. 1-10 ◽  
Author(s):  
Nelson Rangel-Valdez ◽  
Eduardo Fernandez ◽  
Laura Cruz-Reyes ◽  
Claudia Gomez-Santillan ◽  
Gilberto Rivera ◽  
...  

One of the main concerns in Multicriteria Decision Aid (MCDA) is robustness analysis. Some of the most important approaches to model decision maker preferences are based on fuzzy outranking models whose parameters (e.g., weights and veto thresholds) must be elicited. The so-called preference-disaggregation analysis (PDA) has been successfully carried out by means of metaheuristics, but this kind of works lacks a robustness analysis. Based on the above, the present research studies the robustness of a PDA metaheuristic method to estimate model parameters of an outranking-based relational system of preferences. The method is considered robust if the solutions obtained in the presence of noise can maintain the same performance in predicting preference judgments in a new reference set. The research shows experimental evidence that the PDA method keeps the same performance in situations with up to 10% of noise level, making it robust.


Sign in / Sign up

Export Citation Format

Share Document