Return to search

Parameter estimation in a cardiovascular computational model using numerical optimization : Patient simulation, searching for a digital twin

Developing models of the cardiovascular system that simulates the dynamic behavior of a virtual patient’s condition is fundamental in the medical domain for predictive outcome and hypothesis generation. These models are usually described through Ordinary Differential Equation (ODE). To obtain a patient-specific representative model, it is crucial to have an accurate and rapid estimate of the hemodynamic model parameters. Moreover, when adequate model parameters are found, the resulting time series of state variables can be clinically used for predicting the response to treatments and for non-invasive monitoring. In the Thesis, we address the parameter estimation or inverse modeling, by solving an optimization problem, which aims at minimizing the error between the model output and the target data. In our case, the target data are a set of user-defined state variables, descriptive of a hospitalized specific patient and obtained from time-averaged state variables. The Thesis proposes a comparison of both state-of-the-art and novel methods for the estimation of the underlying model parameters of a cardiovascular simulator Aplysia. All the proposed algorithms are selected and implemented considering the constraints deriving from the interaction with Aplysia. In particular, given the inaccessibility of the ODE, we selected gradient-free methods, which do not need to estimate numerically the derivatives. Furthermore, we aim at having a small number of iterations and objective function calls, since these importantly impact the speed of the estimation procedure, and thus the applicability of the knowledge gained through the parameters at the bedside. Moreover, the Thesis addresses the most common problems encountered in the inverse modeling, among which are the non-convexity of the objective function and the identifiability problem. To assist in resolving the latter issue an identifiability analysis is proposed, after which the unidentifiable parameters are excluded. The selected methods are validated using heart failure data, representative of different pathologies commonly encountered in Intensive Care Unit (ICU) patients. The results show that the gradient-free global algorithms Enhanced Scatter Search and Particle Swarm estimate the parameters accurately at the price of a high number of function evaluations and CPU time. As such, they are not suitable for bedside applications. Besides, the local algorithms are not suitable to find an accurate solution given their dependency on the initial guess. To solve this problem, we propose two methods: the hybrid, and the prior-knowledge algorithms. These methods, by including prior domain knowledge, can find a good solution, escaping the basin of attraction of local minima and producing clinically significant parameters in a few minutes. / Utveckling av modeller av det kardiovaskulära systemet som simulerar det dynamiska beteendet hos en virtuell patients är grundläggande inom det medicinska området för att kunna förutsäga resultat och generera hypoteser. Dessa modeller beskrivs vanligtvis genom Ordinary Differential Equation (ODE). För att erhålla en patientspecifik representativ modell är det viktigt att ha en exakt och snabb uppskattning av de hemodynamiska modellparametrarna. När adekvata modellparametrar har hittats kan de resulterande tidsserierna av tillståndsvariabler dessutom användas kliniskt för att förutsäga svaret på behandlingar och för icke-invasiv övervakning. I avhandlingen behandlar vi parameteruppskattning eller invers modellering genom att lösa ett optimeringsproblem som syftar till att minimera följande felet mellan modellens utdata och måldata. I vårt fall är måldata en uppsättning användardefinierade tillståndsvariabler som beskriver en specifik patient som är inlagd på sjukhus och som erhålls från tidsgenomsnittliga tillståndsvariabler. I avhandlingen föreslås en jämförelse av befintlinga och nya metoder. för uppskattning av de underliggande modellparametrarna i en kardiovaskulär simulator, Aplysia. Alla föreslagna algoritmer är valts och implementerade med hänsyn tagna till de begränsningar som finnis i simulatorn Aplysia. Med tanke på att ODE är otillgänglig har vi valt gradientfria metoder som inte behöver uppskatta derivatorna numeriskt. Dessutom strävar vi efter att ha få interationer och funktionsanrop eftersom dessa påverkar hastigheten på estimeringen och därmed den kliniska användbartheten vid patientbehandling. Avhandlingen behandlas dessutom de vanligaste problemen vid inversmodellering som icke-konvexitet och identifierbarhetsproblem. För att lösa det sistnämnda problemet föreslås en identifierbarhetsanalys varefter de icke-identifierbara parametrarna utesluts. De valda metoderna valideras med hjälp av data om hjärtsvikt som är representativa för olika patologier som ofta förekommer hos Intensive Care Unit (ICU)-patienter. Resultaten visar att de gradientfria globala algoritmerna Enhanced Scatter Search och Particle Swarm uppskattar parametrarna korrekt till priset av ett stort antal funktionsutvärderingar och processortid. De är därför inte lämpliga för tillämpningar vid sängkanten. Dessutom är de lokala algoritmerna inte lämpliga för att hitta en exakt lösning eftersom de är beroende av den ursprungliga gissningen. För att lösa detta problem föreslår vi två metoder: hybridalgoritmer och algoritmer med förhandsinformation. Genom att inkludera tidigare domänkunskap kan dessa metoder hitta en bra lösning som undviker de lokala minimernas attraktionsområde och producerar kliniskt betydelsefulla parametrar på några minuter.

Identiferoai:union.ndltd.org:UPSALLA1/oai:DiVA.org:kth-320464
Date January 2022
CreatorsTuccio, Giulia
PublisherKTH, Skolan för elektroteknik och datavetenskap (EECS)
Source SetsDiVA Archive at Upsalla University
LanguageEnglish
Detected LanguageSwedish
TypeStudent thesis, info:eu-repo/semantics/bachelorThesis, text
Formatapplication/pdf
Rightsinfo:eu-repo/semantics/openAccess
RelationTRITA-EECS-EX ; 2022:412

Page generated in 0.0031 seconds