A fast generalizable solution method for glucose control algorithms

Abstract

In critical care tight control of blood glucose levels has been shown to lead to better clinical outcomes. The need to develop new protocols for tight glucose control, as well as the opportunity to optimize a variety of other drug therapies, has led to resurgence in model-based medical decision support in this area. One still valid hindrance to developing new model-based protocols using so-called virtual patients, retrospective clinical data, and Monte Carlo methods is the large amount of computational time and resources needed. This paper develops fast analytical-based methods for insulin–glucose system model that are generalizable to other similar systems. Exploiting the structure and partial solutions in a subset of the model is the key in finding accurate fast solutions to the full model. This approach successfully reduced computing time by factors of 5600–144 000 depending on the numerical error management method, for large (50–164 patients) virtual trials and Monte Carlo analysis. It thus allows new model-based or model-derived protocols to be rapidly developed via extensive simulation. The new method is rigorously compared to existing standard numerical solutions and is found to be highly accurate to within 0.2%.