Simulation Software for Assessment of Nonlinear and Adaptive Multivariable Control Algorithms: Glucose - Insulin Dynamics in Type 1 Diabetes

Recent advances in the precision control strategy of artificial pancreas

The scientific diagnosis and treatment of patients with diabetes require frequent blood glucose testing and insulin delivery to normoglycemia. Therefore, an artificial pancreas with a continuous blood glucose (BG) monitoring function is an urgent research target in the medical industry. The problem of closed-loop algorithmic control of the BG with a time delay is a key and difficult issue that needs to be overcome in the development of an artificial pancreas. Firstly, the composition, structure, and control characteristics of the artificial pancreas are introduced. Subsequently, the research progress of artificial pancreas control algorithms is reviewed, and the characteristics, advantages, and disadvantages of proportional-integral-differential control, model predictive control, and artificial intelligence control are compared and analyzed to determine whether they are suitable for the practical application of the artificial pancreas. Additionally, key advancements in areas such as blood glucose data monitoring, adaptive models, wearable devices, and fully automated artificial pancreas systems are also reviewed. Finally, this review highlights that meal prediction, control safety, integration, streamlining the optimization of control algorithms, constant temperature preservation of insulin, and dual-hormone artificial pancreas are issues that require further attention in the future.

An ensemble machine learning approach for the detection of unannounced meals to enhance postprandial glucose control

BACKGROUND

Hybrid automated insulin delivery systems enhance postprandial glucose control in type 1 diabetes, however, meal announcements are burdensome. To overcome this, we propose a machine learning-based automated meal detection approach; METHODS:: A heterogeneous ensemble method combining an artificial neural network, random forest, and logistic regression was employed. Trained and tested on data from two in-silico cohorts comprising 20 and 47 patients. It accounted for various meal sizes (moderate to high) and glucose appearance rates (slow and rapid absorbing). To produce an optimal prediction model, three ensemble configurations were used: logical AND, majority voting, and logical OR. In addition to the in-silico data, the proposed meal detector was also trained and tested using the OhioT1DM dataset. Finally, the meal detector is combined with a bolus insulin compensation scheme; RESULTS:: The ensemble majority voting obtained the best meal detector results for both the in-silico and OhioT1DM cohorts with a sensitivity of 77%, 94%, 61%, precision of 96%, 89%, 72%, F1-score of 85%, 91%, 66%, and with false positives per day values of 0.05, 0.19, 0.17, respectively. Automatic meal detection with insulin compensation has been performed in open-loop insulin therapy using the AND ensemble, chosen for its lower false positive rate. Time-in-range has significantly increased 10.48% and 16.03%, time above range was reduced by 5.16% and 11.85%, with a minimal time below range increase of 0.35% and 2.69% for both in-silico cohorts, respectively, compared to the results without a meal detector; CONCLUSION:: To increase the overall accuracy and robustness of the predictions, this ensemble methodology aims to take advantage of each base model's strengths. All of the results point to the potential application of the proposed meal detector as a separate module for the detection of meals in automated insulin delivery systems to achieve improved glycemic control.

Artificial Intelligence and Machine Learning for Improving Glycemic Control in Diabetes: Best Practices, Pitfalls, and Opportunities

OBJECTIVE

Artificial intelligence and machine learning are transforming many fields including medicine. In diabetes, robust biosensing technologies and automated insulin delivery therapies have created a substantial opportunity to improve health. While the number of manuscripts addressing the topic of applying machine learning to diabetes has grown in recent years, there has been a lack of consistency in the methods, metrics, and data used to train and evaluate these algorithms. This manuscript provides consensus guidelines for machine learning practitioners in the field of diabetes, including best practice recommended approaches and warnings about pitfalls to avoid.

METHODS

Algorithmic approaches are reviewed and benefits of different algorithms are discussed including importance of clinical accuracy, explainability, interpretability, and personalization. We review the most common features used in machine learning applications in diabetes glucose control and provide an open-source library of functions for calculating features, as well as a framework for specifying data sets using data sheets. A review of current data sets available for training algorithms is provided as well as an online repository of data sources.

SIGNIFICANCE

These consensus guidelines are designed to improve performance and translatability of new machine learning algorithms developed in the field of diabetes for engineers and data scientists.

NIH Fifth Artificial Pancreas Workshop 2023: Meeting Report: The Fifth Artificial Pancreas Workshop: Enabling Fully Automation, Access, and Adoption

The Fifth Artificial Pancreas Workshop: Enabling Fully Automation, Access, and Adoption was held at the National Institutes of Health (NIH) Campus in Bethesda, Maryland on May 1 to 2, 2023. The organizing Committee included representatives of NIH, the US Food and Drug Administration (FDA), Diabetes Technology Society, Juvenile Diabetes Research Foundation (JDRF), and the Leona M. and Harry B. Helmsley Charitable Trust. In previous years, the NIH Division of Diabetes, Endocrinology, and Metabolic Diseases along with other diabetes organizations had organized periodic workshops, and it had been seven years since the NIH hosted the Fourth Artificial Pancreas in July 2016. Since then, significant improvements in insulin delivery have occurred. Several automated insulin delivery (AID) systems are now commercially available. The workshop featured sessions on: (1) Lessons Learned from Recent Advanced Clinical Trials and Real-World Data Analysis, (2) Interoperability, Data Management, Integration of Systems, and Cybersecurity, Challenges and Regulatory Considerations, (3) Adaptation of Systems Through the Lifespan and Special Populations: Are Specific Algorithms Needed, (4) Development of Adaptive Algorithms for Insulin Only and for Multihormonal Systems or Combination with Adjuvant Therapies and Drugs: Clinical Expected Outcomes and Public Health Impact, (5) Novel Artificial Intelligence Strategies to Develop Smarter, More Automated, Personalized Diabetes Management Systems, (6) Novel Sensing Strategies, Hormone Formulations and Delivery to Optimize Close-loop Systems, (7) Special Topic: Clinical and Real-world Viability of IP-IP Systems. "Fully automated closed-loop insulin delivery using the IP route," (8) Round-table Panel: Closed-loop performance: What to Expect and What are the Best Metrics to Assess it, and (9) Round-table Discussion: What is Needed for More Adaptable, Accessible, and Usable Future Generation of Systems? How to Promote Equitable Innovation? This article summarizes the discussions of the Workshop.

Detection of Meals and Physical Activity Events From Free-Living Data of People With Diabetes

BACKGROUND

Predicting carbohydrate intake and physical activity in people with diabetes is crucial for improving blood glucose concentration regulation. Patterns of individual behavior can be detected from historical free-living data to predict meal and exercise times. Data collected in free-living may have missing values and forgotten manual entries. While machine learning (ML) can capture meal and exercise times, missing values, noise, and errors in data can reduce the accuracy of ML algorithms.

METHODS

Two recurrent neural networks (RNNs) are developed with original and imputed data sets to assess detection accuracy of meal and exercise events. Continuous glucose monitoring (CGM) data, insulin infused from pump data, and manual meal and exercise entries from free-living data are used to predict meals, exercise, and their concurrent occurrence. They contain missing values of various lengths in time, noise, and outliers.

RESULTS

The accuracy of RNN models range from 89.9% to 95.7% for identifying the state of event (meal, exercise, both, or neither) for various users. "No meal or exercise" state is determined with 94.58% accuracy by using the best RNN (long short-term memory [LSTM] with 1D Convolution). Detection accuracy with this RNN is 98.05% for meals, 93.42% for exercise, and 55.56% for concurrent meal-exercise events.

CONCLUSIONS

The meal and exercise times detected by the RNN models can be used to warn people for entering meal and exercise information to hybrid closed-loop automated insulin delivery systems. Reliable accuracy for event detection necessitates powerful ML and large data sets. The use of additional sensors and algorithms for detecting these events and their characteristics provides a more accurate alternative.

Multivariable Automated Insulin Delivery System for Handling Planned and Spontaneous Physical Activities

BACKGROUND

Hybrid closed-loop control of glucose levels in people with type 1 diabetes mellitus (T1D) is limited by the requirements on users to manually announce physical activity (PA) and meals to the artificial pancreas system. Multivariable automated insulin delivery (mvAID) systems that can handle unannounced PAs and meals without any manual announcements by the user can improve glycemic control by modulating insulin dosing in response to the occurrence and intensity of spontaneous physical activities.

METHODS

An mvAID system is developed to supplement the glucose measurements with additional physiological signals from a wristband device, with the signals analyzed using artificial intelligence algorithms to automatically detect the occurrence of PA and estimate its intensity. This additional information gained from the physiological signals enables more proactive insulin dosing adjustments in response to both planned exercise and spontaneous unanticipated physical activities.

RESULTS

In silico studies of the mvAID illustrate the safety and efficacy of the system. The mvAID is translated to pilot clinical studies to assess its performance, and the clinical experiments demonstrate an increased time in range and reduced risk of hypoglycemia following unannounced PA and meals.

CONCLUSIONS

The mvAID systems can increase the safety and efficacy of insulin delivery in the presence of unannounced physical activities and meals, leading to improved lives and less burden on people with T1D.

ReplayBG: A Digital Twin-Based Methodology to Identify a Personalized Model From Type 1 Diabetes Data and Simulate Glucose Concentrations to Assess Alternative Therapies

OBJECTIVE

Design and assessment of new therapies for type 1 diabetes (T1D) management can be greatly facilitated by in silico simulations. The ReplayBG simulation methodology here proposed allows "replaying" the scenario behind data already collected by simulating the glucose concentration obtained in response to alternative insulin/carbohydrate therapies and evaluate their efficacy leveraging the concept of digital twin.

METHODS

ReplayBG is based on two steps. First, a personalized model of glucose-insulin dynamics is identified using insulin, carbohydrate, and continuous glucose monitoring (CGM) data. Then, this model is used to simulate the glucose concentration that would have been obtained by "replaying" the same portion of data using a different therapy. The validity of the methodology was evaluated on 100 virtual subjects using the UVa/Padova T1D Simulator (T1DS). In particular, the glucose concentration traces simulated by ReplayBG are compared with those provided by T1DS in five different scenarios of insulin and carbohydrate treatment modifications. Furthermore, we compared ReplayBG with a state-of-the-art methodology for the scope. Finally, two case studies using real data are also presented.

RESULTS

ReplayBG simulates with high accuracy the effect of the considered insulin and carbohydrate treatment alterations, performing significantly better than state-of-art method in almost all considered situations.

CONCLUSION

ReplayBG proved to be a reliable and robust tool to retrospectively explore the effect of new treatments for T1D on the glucose dynamics. It is freely available as open source software at https://github.com/gcappon/replay-bg.

SIGNIFICANCE

ReplayBG offers a new approach to preliminary evaluate new therapies for T1D management before clinical trials.

A practical approach based on learning-based model predictive control with minimal prior knowledge of patients for artificial pancreas

BACKGROUND AND OBJECTIVES

Complete identification of the glucose dynamics for a patient generally requires prior clinical procedures and several measurements for the patient. However, these steps may not be always feasible. To address this limitation, we propose a practical approach integrating learning-based model predictive control (MPC), adaptive basal and bolus injections, and suspension with minimal requirements of prior knowledge of the patient.

METHODS

The glucose dynamic system matrices were periodically updated using only input values, without any pretrained models. The optimal insulin dose was calculated based on a learning-based MPC algorithm. Meal detection and estimation modules were also introduced. The basal and bolus insulin injections were fine-tuned using the performance of glucose control from the previous day. To validate the proposed method, evaluations with 20 virtual patients from a type 1 diabetes metabolic simulator were employed.

RESULTS

Time-in-range (TIR) and time-below-range (TBR) were 90.8% (84.1% - 95.6%) and 0.3% (0% - 0.8%), as represented by the median, first (Q1), and third quartiles (Q3), respectively, when meal intakes were fully announced. When one out of three meal intake announcements was missing, TIR and TBR were 85.2% (75.0% - 88.9%) and 0.9% (0.4% - 1.1%), respectively.

CONCLUSIONS

The proposed approach obviates the need for prior tests from patients and shows effective regulation of blood glucose levels. From the perspective of practical implementation in clinical environments, to deal with minimal prior information of the patient, our study demonstrates how essential clinical knowledge and learning-based modules can be integrated into a control framework for an artificial pancreas.

InsNET: Accurate Basal and Bolus Insulin Dose Prediction for Closed Loop Diabetes Management

It has been demonstrated that closed-loop diabetes management results in better glycemic control and greater compliance than open-loop diabetes management. Deep learning models have been used to implement different components of artifical pancreas. In this work, a novel deep learning model InsNET has been proposed to estimate the basal and bolus insulin level and insulin bolus in patients with type I diabetes utilizing subcutaneous insulin infusion pumps for closed loop diabetes management system. The proposed InsNET is formed with a Wide-Deep combination of LSTM and GRU layers. Additionally, physical activity level has been included as an input in comparison to previous models where only past glucose levels (CGM), meal intake (CHO) and past insulin dosage were used as inputs. The proposed model was tested on In-silico data, and it achieved a Mean Absolute Error (MAE) of 0.002 and Root Mean Squared Error (RMSE) of 0.007 for UVA/Padova Dataset and MAE of 0.001 and RMSE OF 0.003 for mGIPsim Dataset.Clinical relevance- Insulin dose determination is an important as aspect of artificial pancreas. This work describes a deep learning model to determine accurate basal and bolus insulin dosage.

Adaptive Personalized Prior-Knowledge-Informed Model Predictive Control for Type 1 Diabetes

This work considers the problem of adaptive prior-informed model predictive control (MPC) formulations that explicitly incorporate prior knowledge in the model development and is robust to missing data in the output measurements. The proposed prediction model is based on a latent variables model to extract glycemic dynamics from highly-correlated data and incorporates prior knowledge of exponential stability to improve the prediction ability. Missing data structures are formulated to enable model predictions when output measurements are missing for short periods of time. Based on the latent variables model, the MPC strategy and adaptive rules are developed to automatically tune the aggressiveness of the MPC. The adaptive prior-knowledge-informed MPC is evaluated with computer simulations for the control of blood glucose concentrations in people with Type 1 diabetes (T1D) using simulated virtual patients. Due to the variability among people with T1D, the hyperparameters of the prior-knowledge-informed model are personalized to individual subjects. The percentage of time spent in the target range is 76.48% when there are no missing data and 76.52% when there are missing data episodes lasting up to 30 mins (6 samples). Incorporating the adaptive rules further improves the percentage of time in target range to 84.58% and 84.88% for cases with no missing data and missing data, respectively. The proposed adaptive prior-informed MPC formulation provides robust, effective, and safe regulation of glucose concentration in T1D despite disturbances and missing measurements.

A simulator with realistic and challenging scenarios for virtual T1D patients undergoing CSII and MDI therapy

In silico simulations have become essential for the development of diabetes treatments. However, currently available simulators are not challenging enough and often suffer from limitations in insulin and meal absorption variability, which is unable to realistically reflect the dynamics of people with type 1 diabetes (T1D). Additionally, T1D simulators are mainly designed for the testing of continuous subcutaneous insulin infusion (CSII) therapies. In this work, a simulator is presented that includes a generated virtual patient (VP) cohort and both fast- and long-acting Glargine-100 U/ml (Gla-100), Glargine-300 U/ml (Gla-300), and Degludec-100 U/ml (Deg-100) insulin models. Therefore, in addition to CSII therapies, multiple daily injections (MDI) therapies can also be tested. The Hovorka model and its published parameter probability distributions were used to generate cohorts of VPs that represent a T1D population. Valid patients are filtered through restrictions that guarantee that they are physiologically acceptable. To obtain more realistic scenarios, basal insulin profile patterns from the literature have been used to identify variability in insulin sensitivity. A library of mixed meals identified from real data has also been included. This work presents and validates a methodology for the creation of realistic VP cohorts that include physiological variability and a simulator that includes challenging and realistic scenarios for in silico testing. A cohort of 47 VPs has been generated and in silico simulations of both CSII and MDI therapies were performed in open-loop. The simulation outcome metrics were contrasted with literature results.

Generation of post-meal insulin correction boluses in type 1 diabetes simulation models for in-silico clinical trials: More realistic scenarios obtained using a decision tree approach

BACKGROUND AND OBJECTIVE

In type 1 diabetes (T1D) research, in-silico clinical trials (ISCTs) notably facilitate the design/testing of new therapies. Published simulation tools embed mathematical models of blood glucose (BG) and insulin dynamics, continuous glucose monitoring (CGM) sensors, and insulin treatments, but lack a realistic description of some aspects of patient lifestyle impacting on glucose control. Specifically, to effectively simulate insulin correction boluses, required to treat post-meal hyperglycemia (BG > 180 mg/dL), the timing of the bolus may be influenced by subjects' behavioral attitudes. In this work, we develop an easily interpretable model of the variability of correction bolus timing observed in real data, and embed it into a popular simulation tool for ISCTs.

METHODS

Using data collected in 196 adults with T1D monitored in free-living conditions, we trained a decision tree (DT) model to classify whether a correction bolus is injected in a future time window, based on predictors collected back in time, related to CGM data, previous insulin boluses and subject's characteristics. The performance was compared to that of a logistic regression classifier with LASSO regularization (LC), trained on the same dataset. After validation, the DT was embedded within a popular T1D simulation tool and an ISCT was performed to compare the simulated correction boluses against those observed in a subset of data not used for model training.

RESULTS

The DT provided better classification performance (accuracy: 0.792, sensitivity: 0.430, specificity: 0.878, precision: 0.455) than the LC and presented good interpretability. The most predictive features were related to CGM (and its temporal variations), time since the last insulin bolus, and time of the day. The correction boluses simulated by the DT, after implementation in the simulation tool, showed a good agreement with real-world data.

CONCLUSIONS

The DT developed in this work represents a simple set of rules to mimic the same timing of correction boluses observed on real data. The inclusion of the model in simulation tools allows investigators to perform ISCTs that more realistically represent the patient behavior in taking correction boluses and the post-prandial BG response. In the future, more complex models can be investigated.

Integrating Multiple Inputs Into an Artificial Pancreas System: Narrative Literature Review

Background

Type 1 diabetes (T1D) is a chronic autoimmune disease in which a deficiency in insulin production impairs the glucose homeostasis of the body. Continuous subcutaneous infusion of insulin is a commonly used treatment method. Artificial pancreas systems (APS) use continuous glucose level monitoring and continuous subcutaneous infusion of insulin in a closed-loop mode incorporating a controller (or control algorithm). However, the operation of APS is challenging because of complexities arising during meals, exercise, stress, sleep, illnesses, glucose sensing and insulin action delays, and the cognitive burden. To overcome these challenges, options to augment APS through integration of additional inputs, creating multi-input APS (MAPS), are being investigated.

Objective

The aim of this survey is to identify and analyze input data, control architectures, and validation methods of MAPS to better understand the complexities and current state of such systems. This is expected to be valuable in developing improved systems to enhance the quality of life of people with T1D.

Methods

A literature survey was conducted using the Scopus, PubMed, and IEEE Xplore databases for the period January 1, 2005, to February 10, 2020. On the basis of the search criteria, 1092 articles were initially shortlisted, of which 11 (1.01%) were selected for an in-depth narrative analysis. In addition, 6 clinical studies associated with the selected studies were also analyzed.

Results

Signals such as heart rate, accelerometer readings, energy expenditure, and galvanic skin response captured by wearable devices were the most frequently used additional inputs. The use of invasive (blood or other body fluid analytes) inputs such as lactate and adrenaline were also simulated. These inputs were incorporated to switch the mode of the controller through activity detection, directly incorporated for decision-making and for the development of intermediate modules for the controller. The validation of the MAPS was carried out through the use of simulators based on different physiological models and clinical trials.

Conclusions

The integration of additional physiological signals with continuous glucose level monitoring has the potential to optimize glucose control in people with T1D through addressing the identified limitations of APS. Most of the identified additional inputs are related to wearable devices. The rapid growth in wearable technologies can be seen as a key motivator regarding MAPS. However, it is important to further evaluate the practical complexities and psychosocial aspects associated with such systems in real life.

Incorporating Prior Information in Adaptive Model Predictive Control for Multivariable Artificial Pancreas Systems

BACKGROUND

Adaptive model predictive control (MPC) algorithms that recursively update the glucose prediction model are shown to be promising in the development of fully automated multivariable artificial pancreas systems. However, the recursively updated glycemic prediction models do not explicitly consider prior knowledge in the identification of the model parameters. Prior information of the glycemic effects of meals and physical activity can improve model accuracy and yield better glycemic control algorithms.

METHODS

A glucose prediction model based on regularized partial least squares (rPLS) method where the prior information is encoded as the regularization term is developed to provide accurate predictions of the future glucose concentrations. An adaptive MPC is developed that incorporates dynamic trajectories for the glucose setpoint and insulin dosing constraints based on the estimated plasma insulin concentration (PIC). The proposed adaptive MPC algorithm is robust to disturbances caused by unannounced meals and physical activities even in cases with missing glucose measurements. The effectiveness of the proposed adaptive MPC based on rPLS is investigated with in silico subjects of the multivariable glucose-insulin-physiological variables simulator (mGIPsim).

RESULTS

The efficacy of the proposed adaptive MPC strategy in regulating the blood glucose concentration (BGC) of people with T1DM is assessed using the average percent time in range (TIR) for glucose, defined as 70 to 180 mg/dL inclusive, and the average percent time in hypoglycemia (<70 and >54 mg/dL) and level 2 hypoglycemia (≤54 mg/dL). The TIR for a cohort of 20 virtual subjects of mGIPsim is 81.9% ± 7.4% (with no hypoglycemia or severe hypoglycemia) for the proposed MPC compared with 73.9% ± 7.6% (0.2% ± 0.1% in hypoglycemia and 0.1% ± 0.1% in level 2 hypoglycemia) for an MPC based on a recursive autoregressive exogenous (ARX) model.

CONCLUSIONS

The adaptive MPC algorithm that incorporates prior knowledge in the recursive updating of the glucose prediction model can contribute to the development of fully automated artificial pancreas systems that can mitigate meal and physical activity disturbances.

Prior Informed Regularization of Recursively Updated Latent-Variables-Based Models with Missing Observations

Many data-driven modeling techniques identify locally valid, linear representations of time-varying or nonlinear systems, and thus the model parameters must be adaptively updated as the operating conditions of the system vary, though the model identification typically does not consider prior knowledge. In this work, we propose a new regularized partial least squares (rPLS) algorithm that incorporates prior knowledge in the model identification and can handle missing data in the independent covariates. This latent variable (LV) based modeling technique consists of three steps. First, a LV-based model is developed on the historical time series data. In the second step, the missing observations in the new incomplete data sample are estimated. Finally, the future values of the outputs are predicted as a linear combination of estimated scores and loadings. The model is recursively updated as new data are obtained from the system. The performance of the proposed rPLS and rPLS with exogenous inputs (rPLSX) algorithms are evaluated by modeling variations in glucose concentration (GC) of people with Type 1 diabetes (T1D) in response to meals and physical activities for prediction windows up to one hour, or 12 sampling instances, into the future. The proposed rPLS family of GC prediction models are evaluated with both in-silico and clinical experiment data and compared with the performance of recursive time series and kernel-based models. The root mean squared error (RMSE) with simulated subjects in the multivariable T1D simulator where physical activity effects are incorporated in GC variations are 2.52 and 5.81 mg/dL for 30 and 60 mins ahead predictions (respectively) when information for all meals and physical activities are used, increasing to 2.70 and 6.54 mg/dL (respectively) when meals and activities occurred, but the information is with-held from the modeling algorithms. The RMSE is 10.45 and 14.48 mg/dL for clinical study with prediction horizons of 30 and 60 mins, respectively. The low RMSE values demonstrate the effectiveness of the proposed rPLS approach compared to the conventional recursive modeling algorithms.

Mathematical Models of Meal Amount and Timing Variability With Implementation in the Type-1 Diabetes Patient Decision Simulator

BACKGROUND

In type 1 diabetes (T1D) research, in-silico clinical trials (ISCTs) have proven effective in accelerating the development of new therapies. However, published simulators lack a realistic description of some aspects of patient lifestyle which can remarkably affect glucose control. In this paper, we develop a mathematical description of meal carbohydrates (CHO) amount and timing, with the aim to improve the meal generation module in the T1D Patient Decision Simulator (T1D-PDS) published in Vettoretti et al.

METHODS

Data of 32 T1D subjects under free-living conditions for 4874 days were used. Univariate probability density function (PDF) parametric models with different candidate shapes were fitted, individually, against sample distributions of: CHO amounts of breakfast (CHOB), lunch (CHOL), dinner (CHOD), and snack (CHOS); breakfast timing (TB); and time between breakfast-lunch (TBL) and between lunch-dinner (TLD). Furthermore, a support vector machine (SVM) classifier was developed to predict the occurrence of a snack in future fixed-length time windows. Once embedded inside the T1D-PDS, an ISCT was performed.

RESULTS

Resulting PDF models were: gamma (CHOB, CHOS), lognormal (CHOL, TB), loglogistic (CHOD), and generalized-extreme-values (TBL, TLD). The SVM showed a classification accuracy of 0.8 over the test set. The distributions of simulated meal data were not statistically different from the distributions of the real data used to develop the models (α = 0.05).

CONCLUSIONS

The models of meal amount and timing variability developed are suitable for describing real data. Their inclusion in modules that describe patient behavior in the T1D-PDS can permit investigators to perform more realistic, reliable, and insightful ISCTs.

Method to generate a large cohort in-silico for type 1 diabetes

BACKGROUND AND OBJECTIVE

In the last decade, several technological solutions have been proposed as artificial pancreas systems able to treat type 1 diabetes; most often they are built based on a control algorithm that needs to be validated before it is used with real patients. Control algorithms are usually tested with simulation tools that integrate mathematical models related mainly to the glucose-insulin dynamics, but other variables can be considered as well. In general, the simulators have a limited set of subjects. The main goal of this paper is to propose a new computational method to increase the number of virtual subjects, with physiological characteristics, included in the original mathematical models.

METHODS

A subject is defined by a set of parameters given by a mathematical model. From the available reduced number of subjects in the model, the covariance of each parameter of every subject is obtained to establish a mathematical relationship. Then, new sets of parameters are calculated using linear regression methods; this generates larger cohorts, which allows for testing insulin therapies in open-loop or closed-loop scenarios. The new method proposed here increases the number of subjects in a virtual cohort using two versions of Hovorka's mathematical model.

RESULTS

Two covariant cohorts are obtained with linear regression. Both cohorts are clustered to avoid overlapping in the glucose-insulin dynamics and are compared in terms of their qualitative and quantitative behaviours in the normoglycemic range. As a result, there have been generated two larger cohorts (256 subjects) than the original population, which contributes to improving the variability in in-silico tests. In addition, for analysing the characteristics of the covariant generation method, two random cohorts have been generated, where the parameters are obtained individually and independently from each other, exhibiting only distribution limitations so that these cohorts do not have physiological subjects.

CONCLUSIONS

The proposed methodology has enabled the generation of a large cohort of 256 subjects, with different characteristics that are plausible in the T1DM population, significantly increasing the number of available subjects in existing mathematical models. The proposed methodology does not limit the number of subjects that can be generated and thus, it can be used to increase the number of cohorts provided by other mathematical models in diabetes, or even other scientific problems.