Stochastic Seasonal Models for Glucose Prediction in the Artificial Pancreas

A simple modeling framework for prediction in the human glucose-insulin system

Forecasting blood glucose (BG) levels with routinely collected data is useful for glycemic management. BG dynamics are nonlinear, complex, and nonstationary, which can be represented by nonlinear models. However, the sparsity of routinely collected data creates parameter identifiability issues when high-fidelity complex models are used, thereby resulting in inaccurate forecasts. One can use models with reduced physiological fidelity for robust and accurate parameter estimation and forecasting with sparse data. For this purpose, we approximate the nonlinear dynamics of BG regulation by a linear stochastic differential equation: we develop a linear stochastic model, which can be specialized to different settings: type 2 diabetes mellitus (T2DM) and intensive care unit (ICU), with different choices of appropriate model functions. The model includes deterministic terms quantifying glucose removal from the bloodstream through the glycemic regulation system and representing the effect of nutrition and externally delivered insulin. The stochastic term encapsulates the BG oscillations. The model output is in the form of an expected value accompanied by a band around this value. The model parameters are estimated patient-specifically, leading to personalized models. The forecasts consist of values for BG mean and variation, quantifying possible high and low BG levels. Such predictions have potential use for glycemic management as part of control systems. We present experimental results on parameter estimation and forecasting in T2DM and ICU settings. We compare the model's predictive capability with two different nonlinear models built for T2DM and ICU contexts to have a sense of the level of prediction achieved by this model.

Assessment of Seasonal Stochastic Local Models for Glucose Prediction without Meal Size Information under Free-Living Conditions

Accurate blood glucose (BG) forecasting is key in diabetes management, as it allows preventive actions to mitigate harmful hypoglycemic/hyperglycemic episodes. Considering the encouraging results obtained by seasonal stochastic models in proof-of-concept studies, this work assesses the methodology in two datasets (open-loop and closed-loop) recorded in free-living conditions. First, similar postprandial glycemic profiles are grouped together with fuzzy C-means clustering. Then, a seasonal stochastic model is identified for each cluster. Finally, real-time BG forecasting is performed by weighting each model’s prediction. The proposed methodology (named C-SARIMA) is compared to other linear and nonlinear black-box methods: autoregressive integrated moving average (ARIMA), its variant with input (ARIMAX), a feed-forward neural network (NN), and its modified version (NN-X) fed by BG, insulin, and carbohydrates (timing and dosing) information for several prediction horizons (PHs). In the open-loop dataset, C-SARIMA grants a median root-mean-squared error (RMSE) of 20.13 mg/dL (PH = 30) and 27.23 mg/dL (PH = 45), not significantly different from ARIMA and NN. Over a longer PH, C-SARIMA achieves an RMSE = 31.96 mg/dL (PH = 60) and RMSE = 33.91 mg/dL (PH = 75), significantly outperforming the ARIMA and NN, without significant differences from the ARIMAX for PH ≥ 45 and the NN-X for PH ≥ 60. Similar results hold on the closed-loop dataset: for PH = 30 and 45 min, the C-SARIMA achieves an RMSE = 21.63 mg/dL and RMSE = 29.67 mg/dL, not significantly different from the ARIMA and NN. On longer PH, the C-SARIMA outperforms the ARIMA for PH > 45 and the NN for PH > 60 without significant differences from the ARIMAX for PH ≥ 45. Although using less input information, the C-SARIMA achieves similar performance to other prediction methods such as the ARIMAX and NN-X and outperforming the CGM-only approaches on PH > 45min.

Automatic Artificial Pancreas Systems Using an Intelligent Multiple-Model PID Strategy

In this paper, an individualized intelligent multiple-model technique is proposed to design automatic artificial pancreas (AP) systems for the glycemic regulation of type 1 diabetic patients. At first, using the multiple-model concept, the insulin-glucose regulatory system is mathematically identified by constructing some local models. In this step, trade-offs between the number of local models and the complexity of the overall closed-loop system are made by defining and solving a bi-objective optimization problem. Then, optimal AP systems are designed by tuning a bank of proportionalintegralderivative (PID) controllers via the genetic algorithm (GA). A fuzzy gain scheduling strategy is employed to determine the participation percentages of the PID controllers in the control action. Finally, two safety mechanisms, called insulin on board (IOB) constraint and pump shut-off, are installed in the AP systems to enhance their performance. To assess the proposed AP systems, in silico experiments are performed on virtual patients of the UVA/Padova metabolic simulator. The obtained results reveal that the proposed intelligent multiple-model methodology leads to AP systems with limited hyperglycemia and no severe hypoglycemia.

A review of biosensor technology and algorithms for glucose monitoring

Diabetes mellitus (DM) has become a serious illness in the whole world. Until now, there is no effective cure for patients with DM. It is well known that the glucose level is one key factor to determine the progress of DM. It is also an important reference to carry out the accurate and timely treatment for patients with DM. In this article, the related biosensors technology that can be utilized to identify and predict glucose level are reviewed in detail, including the algorithms that can help to achieve numerical value of glucose level. Firstly, the biosensor technology based on the physiological fluids are illustrated, including blood, sweat, interstitial fluid, ocular fluid, and other available fluids. Secondly, the algorithms for achieving numerical value of glucose level are investigated, including the physiological model-based method and the machine learning-based method. Finally, the future development trend and challenges of glucose level monitoring are given and the conclusions are drawn.

Glucose Prediction under Variable-Length Time-Stamped Daily Events: A Seasonal Stochastic Local Modeling Framework

Accurate glucose prediction along a long-enough time horizon is a key component for technology to improve type 1 diabetes treatment. Subjects with diabetes might benefit from supervision and control systems that accurately predict risks and trigger corrective actions early enough with improved mitigation. However, large intra-patient variability poses big challenges to glucose prediction. In previous works by the authors, clustering and local modeling techniques with seasonal stochastic models proved to be efficient, allowing for good glucose prediction accuracy for long prediction horizons. Continuous glucose monitoring (CGM) data were partitioned into fixed-length postprandial time subseries and clustered with Fuzzy C-Means to collect similar behaviors, enforcing seasonality at each cluster after subseries concatenation. Then, seasonal stochastic models were identified for each cluster and local predictions were integrated into a global prediction. However, free-living conditions do not support the fixed-length partition of CGM data since daily events duration is variable. In this work, a new algorithm is provided to overcome this constraint, allowing better coping with patient's variability under variable-length time-stamped daily events in supervision and control applications. Besides predicted glucose, two real-time indices are additionally provided-a crispness index, indicating good representation of current glucose behavior by a single model, and a normality index, allowing for the detection of an abnormal glucose behavior (unusual according to registered historical data). The framework is tested in a proof-of-concept in silico study with ten patients over four month training data and two independent two month validation datasets, with and without abnormal behaviors, from the distribution version of the UVA/Padova simulator extended with diverse sources of intra-patient variability.

Classification of postprandial glycemic patterns in type 1 diabetes subjects under closed-loop control: an in silico study

In this contribution we explore some alternatives in order to obtain filtered and low dimension CGM data to provide well processed CGM data to AP systems. The presented approach explores the possible association of certain patient behaviors with certain glucose patterns. We compare the classical clustering algorithms (K-means, and fuzzy C-means), which has shown some limitations for CGM data processing, with a new clustering algorithm (K-means ellipsoid algorithm) more suited to CGM data. We test this new algorithm in a variety of complex scenarios including variabilty in the amount of ingested carbohydrates, absorption time and intrapatient parameters. The new algorithm overcomes the perceived problems and is able to discriminate between normoglycaemic, moderate and severe hyperglycaemic post-prandial behaviour, even with similar amounts of carbohydrates contained in a meal.

Seasonal Local Models for Glucose Prediction in Type 1 Diabetes

Linear empirical dynamic models have been widely used for blood glucose prediction and risks prevention in people with type 1 diabetes. More accurate blood glucose prediction models with longer prediction horizon (PH) are desirable to enable warnings to patients about imminent blood glucose changes with enough time to take corrective actions. In this study, a blood glucose prediction method is developed by integrating the predictions of a set of seasonal local models (each of them corresponding to different glucose profiles observed along historical data). In the modeling step, the number of sets and their corresponding glucose profiles characteristics are obtained by clustering techniques (Fuzzy C-Means). Then, Box-Jenkins methodology is used to identify a seasonal model for each set. Finally, blood glucose predictions of local models are integrated using different techniques. The proposed method is tested by using 18 60-h closed-loop experiments (including different exercise types and artificial pancreas strategies) and achieving mean absolute percentage error (MAPE) of 2.94%, 3.89%, 5.41%, 6.29% and 8.66% for 15-, 30-, 45-, 60-, and 90-min PHs, respectively.

Long-Term Glucose Forecasting Using a Physiological Model and Deconvolution of the Continuous Glucose Monitoring Signal

(1) Objective: Blood glucose forecasting in type 1 diabetes (T1D) management is a maturing field with numerous algorithms being published and a few of them having reached the commercialisation stage. However, accurate long-term glucose predictions (e.g., >60 min), which are usually needed in applications such as precision insulin dosing (e.g., an artificial pancreas), still remain a challenge. In this paper, we present a novel glucose forecasting algorithm that is well-suited for long-term prediction horizons. The proposed algorithm is currently being used as the core component of a modular safety system for an insulin dose recommender developed within the EU-funded PEPPER (Patient Empowerment through Predictive PERsonalised decision support) project. (2) Methods: The proposed blood glucose forecasting algorithm is based on a compartmental composite model of glucose-insulin dynamics, which uses a deconvolution technique applied to the continuous glucose monitoring (CGM) signal for state estimation. In addition to commonly employed inputs by glucose forecasting methods (i.e., CGM data, insulin, carbohydrates), the proposed algorithm allows the optional input of meal absorption information to enhance prediction accuracy. Clinical data corresponding to 10 adult subjects with T1D were used for evaluation purposes. In addition, in silico data obtained with a modified version of the UVa-Padova simulator was used to further evaluate the impact of accounting for meal absorption information on prediction accuracy. Finally, a comparison with two well-established glucose forecasting algorithms, the autoregressive exogenous (ARX) model and the latent variable-based statistical (LVX) model, was carried out. (3) Results: For prediction horizons beyond 60 min, the performance of the proposed physiological model-based (PM) algorithm is superior to that of the LVX and ARX algorithms. When comparing the performance of PM against the secondly ranked method (ARX) on a 120 min prediction horizon, the percentage improvement on prediction accuracy measured with the root mean square error, A-region of error grid analysis (EGA), and hypoglycaemia prediction calculated by the Matthews correlation coefficient, was 18.8 % , 17.9 % , and 80.9 % , respectively. Although showing a trend towards improvement, the addition of meal absorption information did not provide clinically significant improvements. (4) Conclusion: The proposed glucose forecasting algorithm is potentially well-suited for T1D management applications which require long-term glucose predictions.