WO2024138024A2 - Identification de modèle avec des données d'entrée incomplètes de diabète de type i - Google Patents

Identification de modèle avec des données d'entrée incomplètes de diabète de type i Download PDF

Info

Publication number
WO2024138024A2
WO2024138024A2 PCT/US2023/085495 US2023085495W WO2024138024A2 WO 2024138024 A2 WO2024138024 A2 WO 2024138024A2 US 2023085495 W US2023085495 W US 2023085495W WO 2024138024 A2 WO2024138024 A2 WO 2024138024A2
Authority
WO
WIPO (PCT)
Prior art keywords
glucose
insulin
parameters
data
patient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/US2023/085495
Other languages
English (en)
Other versions
WO2024138024A3 (fr
Inventor
Basak OZASLAN
Eleonora Maria AIELLO
III Francis J. Doyle
Eyal Dassau
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harvard University
Original Assignee
Harvard University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harvard University filed Critical Harvard University
Publication of WO2024138024A2 publication Critical patent/WO2024138024A2/fr
Publication of WO2024138024A3 publication Critical patent/WO2024138024A3/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/20—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for electronic clinical trials or questionnaires
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
    • G16H20/17—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients delivered via infusion or injection
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/67—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for remote operation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients

Definitions

  • the present disclosure relates generally to treatments for diabetes and more specifically to a method of identification of individualized model parameters via determining corrections for incomplete input data aiming for a correct prediction of glucose outcomes.
  • Type 1 Diabetes is associated with the lack of endogenous insulin production, leading to high blood glucose levels. Short of a cure, strategic planning of insulin injection and carbohydrate intake are essential to keep blood glucose levels within a safe range. Yet, glucose management is a challenging task in daily life as (i) glucose metabolism is, in general, difficult to predict, and (ii) the insulin effect happens on a timescale that makes timely glucose regulation non-trivial. Therefore, mathematical models that capture the relationship between the main inputs (i.e., insulin and meal) and glucose outcomes are highly sought to guide immediate and long-term treatment decisions for achieving satisfactory glucose management.
  • free-living data may be affected by many confounding factors and may also be incomplete in many cases.
  • a major challenge in fitting models to glucose metabolism in people with Type 1 diabetes is incomplete data as its collection partially relies on self-reporting and does not include all relevant events.
  • the time and amount of these meals can be wrongly reported or not reported at all (e.g., the person may simply choose not to or forget to report meals or snacks). Such inaccuracies may lead to inaccurate models that result in inappropriate treatment decisions.
  • a method for determining a modified dynamic system model for predicting an output from a dynamic system based on an input is disclosed.
  • Output data is collected from the dynamic system over a period of time.
  • Input data is collected from the dynamic system over a period of time.
  • Parameters of a dynamic system model are optimized via solving a nonlinear least squares problem based on the collected data while fixing a correction value.
  • a modified correction value is determined based on a mixed integer quadratic program executed on the collected data while fixing a constant value for the parameters of the dynamic system model.
  • a final correction value and final optimized parameters are provided to the dynamic system model from iterating the optimizing the parameters of the dynamic system model and determining the modified correction value.
  • a further implementation of the example method is where the dynamic system is a patient and the input data relates to a drug taken by the patient. Another implementation is where the dynamic system is one of a biological system, an industrial system, a population system, a market system, or a resource distribution system. Another implementation is where the example method further includes determining a maximum number of correction values based the number of times an absolute relative difference value of collected output data and an output of an initial fit from the least squares method exceeds a threshold value.
  • Another disclosed example is method for determining an individualized model for predicting glucose levels of a patient.
  • Glucose concentration output data is collected from the patient over a period of time.
  • Carbohydrate size and insulin infusion input data is collected from the patient over a period of time.
  • Parameters of a glucose dynamics model are optimized via solving a nonlinear least squares problem based on the collected data while fixing a constant time and a constant amount of a correction value.
  • the time and the amount of a modified correction value is determined based on executing a mixed integer quadratic program on the collected data while incorporating the optimized parameters of the glucose dynamics model.
  • a final correction value and final optimized parameters of the glucose dynamics model are provided based on iterating the optimizing parameters of the glucose dynamics model and determining the time and the amount of the modified correction value.
  • a further implementation of the example method is where the parameters of the glucose dynamics model include glucose at equilibrium, insulin at equilibrium, insulin sensitivity, meal transport rate, real absorption rate, insulin utilization rate, and insulin clearance rate.
  • the collecting glucose data includes periodically collecting sets of recorded glucose measurements from the patient, and wherein the collecting the input data includes periodically collecting data on insulin injections and carbohydrate intakes from the patient.
  • the glucose measurements are collected from a glucose monitor, the carbohydrate intake is entered from a patient input to a computing device, and the insulin injections are determined from an insulin pump.
  • optimizing parameters of the glucose dynamics model includes minimizing the sum of the squared errors from the nonlinear least squares problem while the number of corrections does not exceed a predetermined value.
  • the predetermined value is determined by the number of times an absolute relative difference value of a recorded glucose trace and an output of an initial fit from the least squares method exceeds a threshold value.
  • the predetermined number of corrections is limited by a predetermined time period.
  • the correction value compensates for incomplete input data including at least one of missing meal reports, miscalculated meals, delay in eating a meal, and insulin pump occlusion.
  • the example method includes storing the model with the final optimized parameters and final correction value in a storage device.
  • Another disclosed example is a system for determining an individualized glucose dynamics model for a patient.
  • the system includes a memory storing machine executable code.
  • a control system is coupled to the memory having one or more processors.
  • the control system is configured to execute the machine executable code to cause the control system to collect glucose data from the patient over a period of time.
  • the control system is configured to collect carbohydrate input data from the patient over a period of time.
  • the control system is configured to optimize parameters of a glucose dynamics model via solving a nonlinear least squares problem based on the collected data while fixing a constant time and a constant mount of a correction value.
  • the control system is configured to determine the time and the amount of a modified correction value based on executing a mixed integer quadratic program on the collected data while incorporating the optimized parameters of the glucose dynamics model.
  • the control system is configured to provide a final correction signal value and final optimized parameters of the glucose dynamics model based on iterating the optimizing parameters of the glucose dynamics model and determining the time and the amount of the modified correction value.
  • a further implementation of the example system is where the parameters of the glucose dynamics model include glucose at equilibrium, insulin at equilibrium, insulin sensitivity, meal transport rate, real absorption rate, insulin utilization rate, and insulin clearance rate.
  • the collection of glucose data includes periodically collecting sets of recorded glucose measurements from the patient, and where the collection of the input data includes periodically collecting data on insulin injections and carbohydrate intakes from the patient.
  • the example system includes a glucose monitor recording the glucose measurements.
  • the example system includes a computing device with a patient input to enter the carbohydrate intakes.
  • the example system includes an insulin pump configured to collect the insulin injections.
  • the optimization of parameters of the glucose dynamics model includes minimizing the sum of the squared errors from the nonlinear least squares problem while the number of corrections does not exceed a predetermined value.
  • the predetermined value is determined by the number of times an absolute relative difference value of a recorded glucose trace and an output of an initial fit from the least squares method exceeds a threshold value.
  • the predetermined number of corrections is limited by a predetermined time period.
  • the correction value compensates for incomplete input data including at least one of missing meal reports, miscalculated meals, delay in eating a meal, and insulin pump occlusion.
  • the example system includes a storage device storing the model with the final optimized parameters and final correction value.
  • Another disclosed example is a non-transitory machine readable medium having stored thereon instructions for performing a method comprising machine executable code which when executed by at least one machine, causes the machine to collect glucose data from the patient over a period of time.
  • the instructions cause the machine to collect carbohydrate input data from the patient over a period of time.
  • the instructions cause the machine to optimize parameters of a glucose dynamics model via solving a nonlinear least squares problem based on the collected data while fixing a constant time and a constant mount of a correction value.
  • the instructions cause the machine to determine the time and the amount of a modified correction value based on executing a mixed integer quadratic program on the collected data while incorporating the optimized parameters of the glucose dynamics model.
  • the instructions cause the machine to provide a final correction value and final optimized parameters of the glucose dynamics model based on iterating the optimizing parameters of the glucose dynamics model and determining the time and the amount of the modified correction value.
  • the system includes an insulin pump configured to deliver insulin into a patient.
  • the system includes a memory storing machine executable code.
  • the system includes a control system coupled to the memory.
  • the control system has one or more processors and is configured to execute the machine executable code to cause the control system to accept a carbohydrate size input and accept an insulin infusion input from the patient.
  • the control system is configured to modify the carbohydrate size input with a correction value of an individualized glucose dynamics model including optimized parameters.
  • the correction value and optimized parameters are determined by collecting glucose concentration output training data from the patient and collecting carbohydrate size and insulin infusion input training data from the patient over a period of time.
  • the determining includes iteratively optimizing parameters of a glucose dynamics model via solving a nonlinear least squares problem based on the collected training data while fixing a constant time and amount of a correction value; and determining the time and amount of a modified correction value based on executing a mixed integer quadratic program on the collected training data while fixing a constant value for the parameters of the glucose dynamics model.
  • the control system provides the insulin infusion input and the modified carbohydrate size input to the individualized glucose dynamics model to determine a prediction of glucose outcome.
  • the control system is configured to send a command to the insulin pump to deliver a determined insulin dosage based on the predicted glucose outcome.
  • a further implementation of the example system is where the parameters of the glucose dynamics model include glucose at equilibrium, insulin at equilibrium, insulin sensitivity, meal transport rate, real absorption rate, insulin utilization rate, and insulin clearance rate.
  • the example system includes a computing device in communication with the control system.
  • the computing device includes a patient input to enter the carbohydrate intake.
  • the insulin pump is configured to collect the insulin infusion input.
  • the optimization of parameters of the glucose dynamics model includes minimizing the sum of the squared errors from the nonlinear least squares problem while the number of corrections does not exceed a predetermined value.
  • the predetermined value is determined by determining the number of times an absolute relative difference value of a recorded glucose trace and an output of an initial fit from the least squares method exceeds a threshold value.
  • the predetermined number of corrections is limited by a predetermined time period.
  • the correction value compensates for incomplete input data including at least one of missing meal reports, miscalculated meals, delay in eating a meal, and insulin pump occlusion.
  • Another disclosed example is a method of providing treatment for a diabetic patient.
  • a carbohydrate size input from a patient is accepted.
  • An insulin infusion input from the patient is accepted.
  • the carbohydrate size input is modified with a correction value of an individualized glucose dynamics model including optimized parameters.
  • the correction value and optimized parameters are determined by: collecting glucose concentration output training data from the patient and collecting carbohydrate size and insulin infusion input training data from the patient over a period of time; and iteratively optimizing parameters of a glucose dynamics model via solving a nonlinear least squares problem based on the collected training data while fixing a constant time and amount of a correction value; and determining the time and amount of a modified correction value based on executing a mixed integer quadratic program on the collected training data while fixing a constant value for the parameters of the glucose dynamics model.
  • the insulin infusion input and the modified carbohydrate size input to the model are provided to determine a prediction of glucose outcome.
  • a command to an insulin pump is sent to deliver a determined insulin dosage based on the predicted glucose outcome.
  • a further implementation of the example method is where the parameters of the glucose dynamics model include glucose at equilibrium, insulin at equilibrium, insulin sensitivity, meal transport rate, real absorption rate, insulin utilization rate, and insulin clearance rate.
  • the carbohydrate intake input is received from a computing device including a patient input to enter the carbohydrate intake, and where the insulin infusion input is collected from an insulin pump.
  • the optimization of parameters of the glucose dynamics model includes minimizing the sum of the squared errors from the nonlinear least squares problem while the number of corrections does not exceed a predetermined value.
  • the predetermined value is determined by determining the number of times an absolute relative difference value of a recorded glucose trace and an output of an initial fit from the least squares method exceeds a threshold value.
  • Another implementation is where the predetermined number of corrections is limited by a predetermined time period. Another implementation is where the correction value compensates for incomplete input data including at least one of missing meal reports, miscalculated meals, delay in eating a meal, and insulin pump occlusion.
  • Another example is a non-transitory machine readable medium having stored thereon instructions for performing a method comprising machine executable code which when executed by at least one machine, causes the machine to accept a carbohydrate size input from a patient.
  • the instructions cause the machine to accept an insulin infusion input from the patient.
  • the instructions cause the machine to modify the carbohydrate size input with a correction value of an individualized glucose dynamics model including optimized parameters.
  • the correction value and optimized parameters are determined by: collecting glucose concentration output training data from the patient and collecting carbohydrate size and insulin infusion input training data from the patient over a period of time; and iteratively optimizing parameters of a glucose dynamics model via solving a nonlinear least squares problem based on the collected training data while fixing a constant time and amount of a correction value; and determining the time and amount of a modified correction value based on executing a mixed integer quadratic program on the collected training data while fixing a constant value for the parameters of the glucose dynamics model.
  • the instructions cause the machine to provide the insulin infusion input and the modified carbohydrate size input to the model to determine a prediction of glucose outcome.
  • the instructions cause the machine to send a command to an insulin pump to deliver a determined insulin dosage based on the predicted glucose outcome.
  • FIG. 1 is an insulin delivery system with a controller for insulin doses with an individualized model for dosage prediction
  • FIG. 2 is a set of graphs showing the need for correction for inaccurate information towards the identification of a dynamic glucose model
  • FIG. 3 is an example algorithm to determine an estimated number of corrections to be applied to the example individualized model
  • FIG. 4 is an example algorithm to provide optimized parameters and corrections for an example individualized model
  • FIG. 5A is a table of corresponding symbols for the plots in FIG. 5B;
  • FIG. 5B is a series of plots that show different scenarios of accurate and inaccurate input data
  • FIG. 6 is a table of root mean square error values across subjects for a comparison of models including a model generated by the example method applied to the different scenarios of accurate and inaccurate input data;
  • FIG. 7A is a series of boxplots showing the estimated number of corrections required for each of the different scenarios;
  • FIG. 7B is a series of boxplots showing how the number of estimated corrections change with different ARD thresholds
  • FIG. 8A shows graphs of model fits compared with the original output and corresponding input data in a training dataset along with the identified input corrections by the example method
  • FIG. 8B shows graphs of model fits compared with the original output and corresponding input data in a testing dataset
  • FIG. 9 is a set of graphs showing an example application on real-life data to obtain individualized glucose model for pregnant women with diabetes.
  • FIG. 10 is a flow diagram of the example routine to obtain corrections and optimized model parameters for an individualized model.
  • words of direction such as “top,” “bottom,” “left,” “right,” “above,” and “below” are intended to relate to the equivalent direction as depicted in a reference illustration; as understood contextually from the object(s) or element(s) being referenced, such as from a commonly used position for the object(s) or element(s); or as otherwise described herein.
  • This disclosure relates to a system and method that uses a sequential identification algorithm that alternates between solving a nonlinear least squares problem to optimize dynamic system model parameters and executing Mixed-Integer Quadratic Programming (MIQP) to optimize the time and amount of input corrections.
  • a dynamic system model is mathematical representation of the dynamics between the inputs and outputs of a dynamic system.
  • the input corrections account for missing or incorrect inputs in a dataset for identifying correct dynamic system model parameters.
  • the unreported or misreported parts in the data are reconciled by adding sparse corrections via the MIQP leading to an improved identification of the model parameters.
  • the present disclosure relates to a model parameter estimation method that is robust against incomplete or incorrect input data.
  • the example method uses available inputoutput data such as a drug input and body reaction as an output for a patient (dynamic system).
  • An example of a dynamic system model is the Subcutaneous Oral Glucose Minimal Model (SOGMM) having insulin and meal data for the inputs and glucose concentration as the output.
  • SOGMM Subcutaneous Oral Glucose Minimal Model
  • the example method has a broad range of potential applications for any dynamic system whose identification may suffer from the incomplete input data, such as human metabolic systems (e.g., the impact of a medication on vital signs), batch distillation, batch reactions, catalyst regeneration, pharmacokinetic and pharmacodynamic models, pricing models for markets, population growth models, resources distribution networks (such as a water distribution network), industrial processes, or aircraft dynamics.
  • human metabolic systems e.g., the impact of a medication on vital signs
  • batch distillation e.g., the impact of a medication on vital signs
  • batch distillation e.g., the impact of a medication on vital signs
  • batch reactions e.g., the impact of a medication on vital signs
  • catalyst regeneration e.g., the impact of a medication on vital signs
  • pharmacokinetic and pharmacodynamic models e.g., pricing models for markets, population growth models, resources distribution networks (such as a water distribution network), industrial processes, or aircraft dynamics.
  • behavioral and environmental inputs e.g., physical
  • the example method can be used in a variety of applications.
  • incomplete datasets and unknown measurements is a common condition, the estimation of the incomplete information can support decision-making and the detection of fault diagnosis.
  • the information on catastrophic events e.g., famine, wars, pandemics
  • the missing information on the user demands in a water distribution network may alter the inputoutput relationship.
  • a complete input information can aid the correct management of the water distribution without loss or waste of resources and can help enhance management of the aforementioned systems.
  • the individualized model with the optimized parameters and corrections may be used to improve the control strategy for insulin dosing.
  • controller parameters such as the penalties on the controller input and the system output may be tuned via the individualized model.
  • the predicted output of the individualized model may be output to display for analysis by a clinician or an insulin treatment parameter optimization program for the patient in relation to treatments for the patient.
  • the model identification in T1D is amenable to an approach using optimized sparse correction for incomplete input data while jointly identifying the model parameters of a dynamic system.
  • Numerous experiments with incomplete in-silico training data were conducted and show that models obtained from the example method are able to provide more accurate predictions on test data than models obtained from standard parameter identification methods.
  • Case studies on the performance of the proposed method for ten in- silico subjects from the United States Food and Drug Administration-accepted UVA/Padova T1D Metabolic Simulator and five real-life use-case scenarios demonstrate the effectiveness of the disclosed method in relation to the prediction of glucose levels.
  • FIG. 1 is a block diagram of an example insulin delivery system 100.
  • the system 100 allows determination of an individualized model and the application of the individualized model for treatment of a person with diabetes such as a patient 110.
  • the individualized model includes optimized parameters and input corrections specific to the patient 110.
  • the system 100 includes a glucose concentration measuring device such as a continuous glucose monitor (CGM) 120 that is worn by the patient 110.
  • CGM 120 outputs continuous measurements of the blood glucose levels of the patient 110 in this example.
  • An insulin pump 130 is also worn by the person 110. The insulin pump 130 supplies insulin in order to regulate blood glucose levels of the patient 110.
  • other insulin administration devices such as a connected insulin pen may be used.
  • a controller 140 provides instructions to the pump 130 to provide insulin boluses to the patient 110.
  • the controller 140 may include a control system that has one or more processors, memory and may include one or more control models, such as an individualized model 142, stored on a memory that process glucose data output from the CGM 120 and feedback insulin provided data from the pump 130, and other data to determine a bolus size of insulin that needs to be delivered to the patient 110 and sends the instructions to the pump 130.
  • other sensors may be used, such as an insulin level sensor that measures the insulin level in the patient 110. In the absence of the CGM 120, glucose levels may be determined by periodic testing of glucose by taking blood fluid samples.
  • the memory is a machine readable medium that stores machine executable code in the form of an application.
  • the controller 140 is a control system coupled to the memory.
  • the controller 140 may include one or more processors.
  • the application is executed by the controller 140 and executes machine executable code to collect individual data from the patient 110 to create the individualized model 142.
  • the controller 140 may be in communication with the pump 130 by a wired or wireless connection.
  • the controller 140 sends commands to the pump 130 to deliver the determined insulin dosage.
  • the CGM 120 may be any suitable sensor for continuous glucose monitoring and may be an under the skin sensor with a wireless connection to the controller 140. In other examples, it may be a non-invasive sensor and have a wired or wireless connection to the controller 140, for instance the Dexcom G6 CGM system manufactured by Dexcom, Inc.
  • the pump 130 may be any suitable insulin pump that is capable of receiving instructions from the controller 140 and delivering insulin boluses to the patient 110.
  • an artificial pancreas using a closed-loop system that includes an insulin pump may be used such as the Medtronic MiniMed 670G, the Tandem Control-IQ, the Omnipod 5, the CamAPS, the TidePool iLet Loop, and like products.
  • the pump 130 may also send a signal to the controller 140 indicating the amounts and times that insulin boluses are delivered to the patient 110.
  • a mobile computing device 150 such as a smartphone carried by the patient 110 includes the controller 140, that may be a CPU, and a memory 144. As will be explained, the mobile computing device 150 may execute an application that incorporates the individualized model 142 for the individual determination of bolus calculations and the control of the insulin pump 130. The application may also include an input interface to allow the patient 110 to enter carbohydrate intakes for inputs into the individualized model 142.
  • the mobile computing device 150 may include a communication interface such as a network transmitter receiver allowing communications to an external network such as the Internet or the Cloud 160 for upload of data related to the bolus calculation.
  • the bolus calculation described herein may be performed within applications running in the Cloud 160 in conjunction, or in place of the mobile computing device 150.
  • the availability of collected data as will be explained, may be accessed by other health care actors via the Cloud 160.
  • Other communication interfaces such as NFC and Bluetooth may be employed for communication with external devices such as the CGM 120 and the insulin pump 130.
  • the application may be preloaded on the mobile computer device 150 or may be offered as an application that may be downloaded to the mobile computer device 150 from a network device such as a server via a network. As explained above, parts or the entire application may be run on the Cloud 160.
  • the application may be executed by a controller on an insulin administration device such as the insulin pump 130.
  • the insulin pump 130 in this example would have a communication interface to receive data from the CGM 120.
  • the output of the application may be displayed to the user on either the pump 130 or on the mobile computer device 150.
  • the system in FIG. 1 is based on an example individualized model for T1D treatment and insulin-meal -glucose models.
  • the main treatment of T1D is external insulin injections, which can be administered manually or using automated insulin delivery systems - depending on the availability and personal choice.
  • the goal is to keep blood glucose levels in a desired range.
  • the two main elements of any insulin treatment regimen are (i) monitoring blood glucose levels, and (ii) timing and sizing the amount of insulin injections.
  • models are sought to provide predictions on glucose outcomes.
  • SOGMM Subcutaneous Oral Glucose Minimal Model
  • R a (t) fc ⁇ k abs ⁇ Q 2 (t)> (lh)
  • G (mg/dL) is the plasma glucose concentration
  • Gb (mg/dL) is the basal plasma glucose concentration
  • Xi (1/min) is the remote insulin compartment that corresponds to the portion of insulin that acts on glucose
  • Qi (mg) and Q2 (mg) are the two compartments representing the oral glucose transport
  • Ld and I SC 2 (mU) are the interstitial insulin transport compartments
  • Ip (mU) and I (mU/1) represent the insulin mass in plasma
  • R a (mg/min kg) represents the glucose rate of appearance.
  • BW(kg) is the person’s bodyweight V g (kg/dL) and Vi (L/kg) represent the distribution volume of glucose and insulin, respectively.
  • /?2, kd, kd are transport rates that act on the insulin-related compartments (i.e., Iscl. Isc2, Ip, Xi) while k T and kabs are the transport rates that act on meal related compartments (i.e., Qi, Q2).
  • S g (1/min) is the fractional glucose effectiveness
  • f c is the fractional meal absorption.
  • This model has two inputs: u m (mg/min) and m (mU/min) that are the carbohydrate size in an ingested meal and insulin infusion, respectively.
  • the parameters in the set of Equations (1) are generally assumed to be time-invariant.
  • the example process uses the SOGMM from the above description with an example model identification with input correction process to create the individualized model 142 in FIG. 1.
  • the majority of the parameters in the SOGMM are structurally identifiable.
  • the nominal values of the parameters Gb, h, and Si can be obtained based on average glucose of a patient over a period of time, such as three months, the prescribed basal insulin profile, and total daily insulin use, respectively.
  • plasma glucose concentration, Gb may be derived from data collected via a sensor such as the CGM 120 in FIG. 1.
  • the insulin at equilibrium h may be read from any suitable insulin recording device or application.
  • the insulin at equilibrium is read from a sensor on the insulin pump 130 in FIG. 1.
  • the meal transport rate, k T has a nominal value of 0.0018
  • the meal absorption rate, kabs has a nominal value of 0.0182
  • the insulin utilization rate, kd has a nominal value of 0.0100
  • the insulin clearance rate, kd has a nominal value of 0.1600.
  • the nominal values are obtained from existing literature for an average population.
  • a vector of insulin inputs u/ G R w , and a vector of meal inputs Um G T Equation (4) may be simulated forward to obtain the simulated outputs y G : in the following compact form: where the k th values of y, w, uTM are y(k), m(k), u m (k) in Equation (2), respectively. Note that/ 7 is linear in inputs w and u m but nonlinear in 0 P .
  • the initial state x(0) plays a negligible role in long-term simulations as its effect is dissipated after few hours into the simulation (warm-up period). In the case studies, the states are always initialized to their nominal values, which correspond to their equilibrium values derived from Equation (2) and the first 5 hours of data are used only for warm-up and then excluded from the performance evaluation.
  • Solving Equation (7) is the standard approach to find Op, but this solution is known to suffer from misreported/unreported input data.
  • the noise in y and u/ are relatively negligible, but the values in uTM are more prone to errors since the estimation can be subjective and recording is at individual discretion (e.g., a person may forget to report a meal during a busy time).
  • the goal of the present method is to identify the optimized parameters, 0* in a way that the optimized parameters capture the true dynamics regardless of inaccuracies/incompleteness in the data used for model identification.
  • the approach is to correct the carbohydrate concentrations, uTM by adding a correction signal, u c , such that resulting parameter values represent the dynamics of the individual’s metabolism with high fidelity as shown in FIG. 2.
  • T> as in Equation (6) integer y > 0, the optimal model parameters and the optimal corrections to meal input are found by minimizing the sum of the squared errors, while the number of non-zero meal correction instances does not exceed y: op,U c subject to
  • II Ih andll Ho are used to represent the Euclidean norm and the number of non-zero elements (zero semi-norm), respectively.
  • FIG. 2 shows a graph 200 of training data that includes a plot 202 representing a true output of glucose.
  • a second graph 210 represents the training input data with input inaccuracies.
  • a line 212 represents a recorded meal and a line 214 represents an unrecorded meal, which results in the training dataset to have inaccurate input.
  • a second set of graphs 220 and 230 relate to the parameter identification process if the model parameters are identified without correcting the training input.
  • the graph 220 represents uncorrected training input data.
  • the graph 220 only includes the recorded meal in a line 222. The unrecorded meal is not used in the training procedure.
  • the graph 230 represents the glucose output from the true measurement as shown in a plot 232 and the output predicted through a least squares model in Equation 7 based on the uncorrected training data as shown in a plot 234.
  • the graphs 240 and 250 represent the effect of the desired corrected data that is generated from the example method.
  • the graph 240 represents the recorded input data in a line 242.
  • the input correction obtained by the example method is shown as a line 244.
  • the correction obtained in a line 244 is similar to the unrecorded input data in the line 214.
  • the graph 250 represents the glucose output from the true measurement as shown in a plot 252 and that which is predicted through a least squares model with corrections such as that shown in the graph 240 as shown in a plot 254.
  • the outputs of the graph 250 that shows the glucose level predicted by the example optimized model is ideally reflective of the true glucose level due to application of the correction.
  • Two final graphs 270 and 280 show the learned dynamics of the model on the testing data.
  • the graph 270 is a graph of sample test data with accurate inputs as shown in a line 272 of an input meal.
  • the graph 280 is a set of outputs of the glucose level.
  • the graph 280 includes a true output data curve 282 and an identified output curve 284 generated with the model parameters obtained from the training phase with the example method.
  • the example method compensates for the missing data in the training phase that results in the output curve 284 being similar to the true output data curve 282.
  • the result of a model using the least square (LS) method with inaccurate data without corrections is shown as a curve 286.
  • Equation (8) An alternative to II udlo ⁇ y is adding lludl ⁇ to the cost in Equation (8) so that adding more and larger corrections would be penalized. Equation (8) is better suited to deal with incomplete inputs as corrections occur without a cost.
  • the only restriction is sparsity imposed by the cardinality constraint which suits the nature of typical errors in input data in T1D (e.g., incorrect meal report).
  • Equation (8) falls into the category of nonlinear optimization as T is nonlinear in 9p in Equation (1). Furthermore, the cardinality constraint HucllO ⁇ y is non-convex and makes the problem combinatorial. Therefore, solving Equation (8) is challenging as (i) gradient-based optimization is not suitable to handle the cardinality constraint and (ii) the problem size (the number of decision variables in Op and u c ) is too large for global search (using methods such as simulated annealing, genetic algorithm, etc.).
  • Equation (8) becomes a standard, smooth, nonlinear least squares problem, which can be solved to local optimality using standard gradient-based methods. Furthermore, HucllO ⁇ y can be cast as a set of mixed-integer constraints. While mixed integer optimization is computationally prohibitive, mature commercial solvers, such as Gurobi Optimization, are able to solve moderately sized mixed-integer linear/quadratic programs to (near) global optimality in reasonable times. However, such solvers cannot handle nonlinear constraints on Op.
  • the disclosed method is an alternating minimization approach to solve Equation (8) to obtain corrections for the model.
  • the correction signal, w c is fixed and the model parameters, 0 P is solved using a standard nonlinear least-squares solver.
  • the model parameters, 0 P is fixed and the correction signal, w c , is solved using Mixed- Integer Quadratic Programming (MIQP).
  • MIQP Mixed- Integer Quadratic Programming
  • the two steps are alternated until convergence.
  • convergence is defined as the latest solution cost of Equation (9) below being within 5% of the previous iteration’s cost while the maximum change in the parameters in 0 P remains within 10% range of their values from the previous iteration.
  • the pseudo-code of the described two-step process is shown in an Algorithm in FIG. 4. The details are explained in the following.
  • Equation (9) is a nonlinear program, it is smooth and can be solved to local optimality.
  • Equation (9) is warm started with the previous model parameters 0 P .
  • all elements of u c are set to zero and the population-level values for the model parameters 0 P are taken from known studies to provide a reasonable initial guess of the initial values for the model parameters 0 P .
  • the next phase is to choose the number of corrections (y) via an example algorithm such as the algorithm in FIG. 3.
  • a “significant difference” between the model prediction and the real output implies incomplete input within the few hours leading to the difference.
  • the significant difference is defined as an Absolute Relative Difference (ARD) value > 20/20%.
  • ARD Absolute Relative Difference
  • the selected ARD threshold is about twice the mean ARD of the FDA-approved CGM devices between 2016 and 2021, hence regular CGM noise does not increase the estimated number of corrections.
  • MIQP Mixed-Integer Quadratic Programming
  • a first scenario was complete information where all inputs are recorded accurately.
  • a second scenario is unreported snacks. This scenario reflects that persons tend to not report the snacks that they do not bolus for. This scenario implements this behavior for the bed-time snack on both days that are 40g and 30g, respectively.
  • a third scenario is miscalculated meals. Carbohydrate counting is a major challenge in daily T1D care and errors are frequent. For this scenario, a 30% overestimation on the biggest meal, lunch on Day 1, and a 30% underestimation error on the breakfast of Day 2 was implemented.
  • a fourth scenario is delays in eating the bolused meal. While this scenario is less frequent compared to the earlier scenarios, persons may have delays in eating the meal that they already bolused for.
  • this scenario involves both misreporting (a full meal is reported when no meal is actually eaten) and unreporting (a meal is eaten at a later time with no report) leading to a combined information inaccuracy in the input data.
  • This situation is implemented by introducing a one- hour delay to eating the lunch on Day 1 that was bolused for and reported at 2PM but eaten at 3PM.
  • a fifth scenario is pump occlusion where the misrecorded input is insulin. The input data suggests insulin being delivered while the insulin is not actually delivered, leading to high glucose values. It is difficult to detect pump occlusion in practice. Pump occlusion for the breakfast bolus and the following two-hours on Day 2 were implemented.
  • FIG. 5 A is a table 500 showing the references symbols in the graphs 510, 520, 530, 540, and 550 in FIG. 5B.
  • the table 500 includes a column 502 that lists whether the data correctly reported, misreported, or unreported.
  • a second column 504 lists bolus insulin input for both injected and not injected bolus. These show either insulin injections from the pump or other source and mealindependent insulin injections, i.e., basal insulin and correction bolus.
  • a third column 506 lists meal inputs.
  • the graph 510 in FIG. 5B shows the description of a complete information scenario.
  • the graph 510 includes a series of dots such as the dot 512 that represent times of injected insulin and a series of lines such as the line 514 that represents a meal input.
  • a line 516 represents meal independent injected insulin.
  • the graph 520 shows the results of an unreported snacks scenario.
  • the graph 520 thus includes two lines 522 that represent unreported snacks at different times.
  • the graph 530 shows the results of a miscalculated meals scenario.
  • the graph 530 includes two dashed lines 532 that represent miscalculated meals at different times.
  • the graph 540 shows the results of delays in eating a bolused meal.
  • the graph 540 shows a line 542 that represents a miscalculated meal and a dashed line 544 that represents a misreported meal.
  • the graph 550 shows the results of a pump occlusion scenario.
  • the graph 550 shows a cross symbol 552 that represents a misreported not injected bolus input.
  • Meal times, meal sizes, and unbolused snack times were different from the training scenario and were designed as follows: Day 1 (breakfast of 40g at 7AM, snack of 20g at 10AM, lunch of 70g at 2PM, dinner of 60g at 8PM), Day 2 (breakfast of 30g at 8AM, snack of 10g at 10AM, lunch of 50g at 1PM, dinner of 80g at 8PM).
  • FIG. 6 shows a table 600 that shows the RMSE values across 10 subjects for models obtained from each of the five scenarios and the three methods.
  • results are obtained for the glucose predictions using the nominal model, a model identified by the known least squares (LS) method, and a model identified by the example least square correction (LS-C) method.
  • LS-C model identified by the example least square correction
  • the example LS-C based model suggested the lowest number of corrections for the scenario of complete information and the highest number of corrections for the fifth scenario of CGM pump occlusion, when evaluated via median [IQR].
  • the median number of corrections estimated across in-silico subjects was 0 for Scenario I with the interquartile range of the number of corrections estimated across in -silico subjects being 0.
  • the respective medium number of corrections and interquartile ranges were 2.5 [1] for Scenario II, 3 [1] for Scenario III, 3 [2] for Scenario IV, and 6 [1.75] for Scenario V.
  • FIG. 7A shows a series of boxplots 700 showing the estimated number of corrections (y) required for each of the five scenarios.
  • FIG. 7B shows a series of boxplots 750 of the estimated number of corrections (y) required for different ARD thresholds.
  • the data for the boxplot 700 was collected across all subjects per scenario for the selected ARD threshold (20/20%).
  • the boxplot 700 illustrates the distribution of the corrections, y across all subjects per scenario.
  • a posteriori analysis of how estimated y changes for different thresholds of ARD across all scenarios is shown in the boxplot 750 in FIG. 7B.
  • the plot 750 shows data collected across all subjects and scenarios for different ARD threshold choices.
  • FIGs. 8 A shows a graph 810 of the training input data and a graph 820 of the model training for glucose outputs and for the fourth scenario of delays in eating the bolused meal.
  • the graph 810 shows a set of insulin injections symbolized by dots 812, a carbohydrate input symbolized by lines 814, and lines 816 representing optimal corrections generated by the example method.
  • the graph 820 shows the sample training output for Subject #6 when 2PM meal was digested one hour after it was recorded.
  • the graph 820 includes a plot 822 of the true glucose output, a plot 824 of the model based on the least square method, and a plot 826 of the model based on the example least squares with corrections (LS-C) method.
  • the graph 820 shows model fits for identification with standard approach described in Equation 7 in the plot 824 and the identification of the example LS-C method in the plot 826 (green line) compared with the true data in the plot 822.
  • FIG. 8B shows a graph 850 shows the testing outcomes obtained by the models identified on the training data shown in the graph 810 for the fourth scenario of delays in eating the bloused meal.
  • the graph 850 includes a plot 852 of the true glucose output, a plot 854 of the output of the model based on the least square (LS) method, and a plot 856 of the output of the model based on the example least squares with corrections (LS-C) method.
  • Sample test output for Subject #6 with the models identified with standard approach described in Equation (7) in the plot 854 and with the example least-squares with correction (LS-C) algorithm in the plot 856 are compared with the test data in the plot 852.
  • the graph 860 shows the meal and insulin input data that are reported with a series of lines 862 and dots 864, respectively.
  • the graphs in FIG. 8A-8B show how the presence of inaccuracies in the data, e.g., 2PM meal on day 1 digested one hour after it was recorded, can compromise the correct identification of the model parameters when using the standard LS approach from Equation (7).
  • the graph 850 in FIG. 8B shows how the example LS-C method allows to handle the inaccuracy in the identification phase achieving a superior glucose prediction performance on the test data compared to LS.
  • the example method described above may be used to identify metabolic parameters of 8 participants who participated in a supervised pregnancy study. The run-in portion of the study is used to identify model parameters and the supervised portion of the study is used to test the performance of the identified model.
  • the insulin sensitivity has a periodic nature with a 24-hour period and most of the time individuals have their own specific diurnal pattern - albeit non-periodic events such as physical exercise may cause temporary significant changes.
  • a truncated Fourier series is used to model the periodic changes of insulin sensitivity:
  • FIG. 9 shows a graph 910 that shows a comparison of the true glucose data in a plot 912, a plot 914 of glucose simulated from the least squares method, and a plot 916 of glucose simulated from the example least square with correction method over a time index.
  • a graph 920 shows the carbohydrate input data.
  • a set of lines 922 constitute announced carbohydrate data.
  • a set of lines 924 show the corrections generated from the example least squares with correction method that lead to the more accurate simulation in the plot 916.
  • a graph 930 shows a series of lines 932 of bolus insulin levels relative to the time index. The graph 930 also shows a series of line 934 of basal insulin levels relative to the time index.
  • a flow diagram 1000 in FIG. 10 is representative of example machine readable instructions for the process of generating an individualized model for an individual.
  • the machine readable instructions comprise an algorithm for execution by: (a) a processor; (b) a controller; and/or (c) one or more other suitable processing device(s).
  • the algorithm may be embodied in software stored on tangible media such as flash memory, CD- ROM, floppy disk, hard drive, digital video (versatile) disk (DVD), or other memory devices.
  • a set of training input data such as carbohydrate levels over a period of time and corresponding times and a set of training output data such as glucose concentrations over the period of time is collected from the individual patient (1010).
  • the initial values for model parameters are set (1012).
  • the estimated maximum number of corrections is then determined from the input data and output data (1014).
  • the model parameters are solved with a non-linear least squares solver while holding the correction constant (1016).
  • the correction is solved by executing Mixed-Integer Quadratic Programming while incorporating the solved model parameters in the dynamic model (1018).
  • the routine determines whether the solution cost of the least squares solver is within a predetermined percentage of the previous iteration (1020). If the solution cost is outside of the predetermined percentage, the routine loops back to 1016 for another iteration of determining parameters and corrections. If the solution cost is within the predetermined percentage, the routine determines whether the maximum change in the parameters is within a predetermined percentage range (1022). If the maximum change is outside of the predetermined percentage range, the routine loops back to 1016 for another iteration of determining parameters and corrections. If the maximum change is inside the predetermined percentage range, the model with the final optimized parameters and the final corrections are stored (1024).
  • the results show the importance of accurate input information in model identification and that the example least squares with correction method (LS-C), is effective in reconciling the unfavorable effects of imperfect input in model identification for T1D.
  • LS-C example least squares with correction method
  • misreported meals are the primary motivation to apply corrections through the meal input
  • other correction may be applied for phenomena such as different macronutrients, acute psychological stress, and physical activities are unreported and unmodeled, that may have profound effects on insulin-glucose dynamics.
  • the effects of these can also be approximately modeled by artificial meals. For instance, physical exercise is known to cause a sharp decrease in glucose values and hence can be roughly viewed as the person has taken negative meal.
  • the example least square with correction method can also handle inputs in the form of a step or a ramp instead of an impulse allowing the integration of the effect of combined inaccuracies.
  • Embodiments of the present disclosure may comprise or utilize a special purpose or general-purpose computer including computer hardware, such as, for example, one or more processors and system memory, as discussed in greater detail below.
  • Embodiments within the scope of the present disclosure also include physical and other computer-readable media for carrying or storing computer-executable instructions and/or data structures.
  • one or more of the processes described herein may be implemented at least in part as instructions embodied in a non-transitory computer-readable medium and executable by one or more computing devices (e.g., any of the media content access devices described herein).
  • a processor receives instructions, from a non-transitory computer- readable medium, (e.g., a memory, etc.), and executes those instructions, thereby performing one or more processes, including one or more of the processes described herein.
  • a non-transitory computer- readable medium e.g., a memory, etc.
  • Computer-readable media can be any available media that can be accessed by a general purpose or special purpose computer system.
  • Computer-readable media that store computer-executable instructions are non-transitory computer-readable storage media (devices).
  • Computer-readable media that carry computer-executable instructions are transmission media.
  • embodiments of the disclosure can comprise at least two distinctly different kinds of computer-readable media: non-transitory computer-readable storage media (devices) and transmission media.
  • Non-transitory computer-readable storage media includes RAM, ROM, EEPROM, CD-ROM, solid state drives (“SSDs”) (e.g., based on RAM), Flash memory, phasechange memory (“PCM”), other types of memory, other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store desired program code means in the form of computer-executable instructions or data structures and which can be accessed by a general purpose or special purpose computer.
  • SSDs solid state drives
  • PCM phasechange memory
  • a “network” is defined as one or more data links that enable the transport of electronic data between computer systems and/or modules and/or other electronic devices.
  • a network or another communications connection can include a network and/or data links which can be used to carry desired program code means in the form of computer-executable instructions or data structures and which can be accessed by a general purpose or special purpose computer. Combinations of the above should also be included within the scope of computer-readable media.
  • program code means in the form of computer-executable instructions or data structures can be transferred automatically from transmission media to non-transitory computer-readable storage media (devices) (or vice versa).
  • computer-executable instructions or data structures received over a network or data link can be buffered in RAM within a network interface module (e.g., a “NIC”), and then eventually transferred to computer system RAM and/or to less volatile computer storage media (devices) at a computer system.
  • a network interface module e.g., a “NIC”
  • non-transitory computer-readable storage media (devices) can be included in computer system components that also (or even primarily) utilize transmission media.
  • Computer-executable instructions comprise, for example, instructions and data which, when executed at a processor, cause a general purpose computer, special purpose computer, or special purpose processing device to perform a certain function or group of functions.
  • computer-executable instructions are executed on a general purpose computer to turn the general purpose computer into a special purpose computer implementing elements of the disclosure.
  • the computer executable instructions may be, for example, binaries, intermediate format instructions such as assembly language, or even source code.
  • the disclosure may be practiced in network computing environments with many types of computer system configurations, including, personal computers, desktop computers, laptop computers, message processors, hand-held devices, multi-processor systems, microprocessor-based or programmable consumer electronics, network PCs, minicomputers, mainframe computers, mobile telephones, PDAs, tablets, pagers, routers, switches, and the like.
  • the disclosure may also be practiced in distributed system environments where local and remote computer systems, which are linked (either by hardwired data links, wireless data links, or by a combination of hardwired and wireless data links) through a network, both perform tasks.
  • program modules may be located in both local and remote memory storage devices.
  • Embodiments of the present disclosure can also be implemented in cloud computing environments.
  • “cloud computing” is defined as an un-subscription model for enabling on-demand network access to a shared pool of configurable computing resources.
  • cloud computing can be employed in the marketplace to offer ubiquitous and convenient on-demand access to the shared pool of configurable computing resources.
  • the shared pool of configurable computing resources can be rapidly provisioned via virtualization and released with low management effort or service provider interaction, and then scaled accordingly.
  • a cloud-computing un-subscription model can be composed of various characteristics such as, for example, on-demand self-service, broad network access, resource pooling, rapid elasticity, measured service, and so forth.
  • a cloud-computing un-subscription model can also expose various service un-subscription models, such as, for example, Software as a Service (“SaaS”), a web service, Platform as a Service (“PaaS”), and Infrastructure as a Service (“laaS”).
  • SaaS Software as a Service
  • PaaS Platform as a Service
  • laaS Infrastructure as a Service
  • a cloud-computing un-subscription model can also be deployed using different deployment un-subscription models such as private cloud, community cloud, public cloud, hybrid cloud, and so forth.
  • a “cloud-computing environment” is an environment in which cloud computing is employed.
  • a computing device may be configured to perform one or more of the processes described above, the computing device can comprise a processor, a memory, a storage device, an VO interface, and a communication interface, which may be communicatively coupled by way of a communication infrastructure.
  • the computing device can include fewer or more components than those described above.
  • the processor includes hardware for executing instructions, such as those making up a computer program.
  • the processor may retrieve (or fetch) the instructions from an internal register, an internal cache, the memory, or the storage device and decode and execute them.
  • the memory may be a volatile or non-volatile memory used for storing data, metadata, and programs for execution by the processor(s).
  • the storage device includes storage, such as a hard disk, flash disk drive, or other digital storage device, for storing data or instructions related to object digitizing processes (e.g., digital scans, digital models).
  • the I/O interface allows a user to provide input to, receive output from, and otherwise transfer data to and receive data from computing device.
  • the I/O interface may include a mouse, a keypad or a keyboard, a touch screen, a camera, an optical scanner, network interface, modem, other known I/O devices or a combination of such I/O interfaces.
  • the I/O interface may include one or more devices for presenting output to a user, including, but not limited to, a graphics engine, a display (e.g., a display screen), one or more output drivers (e.g., display drivers), one or more audio speakers, and one or more audio drivers.
  • the I/O interface is configured to provide graphical data to a display for presentation to a user.
  • the graphical data may be representative of one or more graphical user interfaces and/or any other graphical content as may serve a particular implementation.
  • the communication interface can include hardware, software, or both.
  • the communication interface can provide one or more interfaces for communication (such as, for example, packet-based communication) between the computing device and one or more other computing devices or networks.
  • the communication interface may include a network interface controller (NIC) or network adapter for communicating with an Ethernet or other wire-based network or a wireless NIC (WNIC) or wireless adapter for communicating with a wireless network, such as a WI-FI.
  • NIC network interface controller
  • WNIC wireless NIC
  • the communication interface may facilitate communications with various types of wired or wireless networks.
  • the communication interface may also facilitate communications using various communication protocols.
  • the communication infrastructure may also include hardware, software, or both that couples components of the computing device to each other.
  • the communication interface may use one or more networks and/or protocols to enable a plurality of computing devices connected by a particular infrastructure to communicate with each other to perform one or more aspects of the digitizing processes described herein.
  • the image compression process can allow a plurality of devices (e.g., server devices for performing image processing tasks of a large number of images) to exchange information using various communication networks and protocols for exchanging information about a selected workflow and image data for a plurality of images.
  • the disclosure herein may be implemented with any type of hardware and/or software, and may be a pre-programmed general purpose computing device.
  • the system may be implemented using a server, a personal computer, a portable computer, a thin client, or any suitable device or devices.
  • the disclosure and/or components thereof may be a single device at a single location, or multiple devices at a single, or multiple, locations that are connected together using any appropriate communication protocols over any communication medium such as electric cable, fiber optic cable, or in a wireless manner.
  • modules which perform particular functions. It should be understood that these modules are merely schematically illustrated based on their function for clarity purposes only, and do not necessary represent specific hardware or software. In this regard, these modules may be hardware and/or software implemented to substantially perform the particular functions discussed. Moreover, the modules may be combined together within the disclosure, or divided into additional modules based on the particular function desired. Thus, the disclosure should not be construed to limit the present invention, but merely be understood to illustrate one example implementation thereof.
  • the computing system can include clients and servers.
  • a client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other.
  • a server transmits data (e.g., an HTML page) to a client device (e.g., for purposes of displaying data to and receiving user input from a user interacting with the client device).
  • client device e.g., for purposes of displaying data to and receiving user input from a user interacting with the client device.
  • Data generated at the client device e.g., a result of the user interaction
  • Implementations of the subject matter described in this specification can be implemented in a computing system that includes a back end component, e.g., as a data server, or that includes a middleware component, e.g., an application server, or that includes a front end component, e.g., a client computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the subject matter described in this specification, or any combination of one or more such back end, middleware, or front end components.
  • the components of the system can be interconnected by any form or medium of digital data communication, e.g., a communication network. Examples of communication networks include a local area network (“LAN”) and a wide area network (“WAN”), an internetwork (e.g., the Internet), and peer-to-peer networks (e.g., ad hoc peer to-peer networks).
  • LAN local area network
  • WAN wide area network
  • Internet internetwork
  • peer-to-peer networks e.g.,
  • Implementations of the subject matter and the operations described in this specification can be implemented in digital electronic circuitry, or in computer software, firmware, or hardware, including the structures disclosed in this specification and their structural equivalents, or in combinations of one or more of them.
  • Implementations of the subject matter described in this specification can be implemented as one or more computer programs, i.e., one or more modules of computer program instructions, encoded on computer storage medium for execution by, or to control the operation of, data processing apparatus.
  • the program instructions can be encoded on an artificially generated propagated signal, e.g., a machine-generated electrical, optical, or electromagnetic signal that is generated to encode information for transmission to suitable receiver apparatus for execution by a data processing apparatus.
  • a computer storage medium can be, or be included in, a computer-readable storage device, a computer-readable storage substrate, a random or serial access memory array or device, or a combination of one or more of them. Moreover, while a computer storage medium is not a propagated signal, a computer storage medium can be a source or destination of computer program instructions encoded in an artificially generated propagated signal. The computer storage medium can also be, or be included in, one or more separate physical components or media (e.g., multiple CDs, disks, or other storage devices).
  • control system encompasses all kinds of apparatus, devices, and machines for processing data, including by way of example a programmable processor, a computer, a system on a chip, or multiple ones, or combinations, of the foregoing.
  • the apparatus can include special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application specific integrated circuit).
  • the apparatus can also include, in addition to hardware, code that creates an execution environment for the computer program in question, e.g., code that constitutes processor firmware, a protocol stack, a database management system, an operating system, a cross-platform runtime environment, a virtual machine, or a combination of one or more of them.
  • the apparatus and execution environment can realize various different computing model infrastructures, such as web services, distributed computing and grid computing infrastructures.
  • a computer program (also known as a program, software, software application, script, or code) can be written in any form of programming language, including compiled or interpreted languages, declarative or procedural languages, and it can be deployed in any form, including as a stand alone program or as a module, component, subroutine, object, or other unit suitable for use in a computing environment.
  • a computer program may, but need not, correspond to a file in a file system.
  • a program can be stored in a portion of a file that holds other programs or data (e.g., one or more scripts stored in a markup language document), in a single file dedicated to the program in question, or in multiple coordinated files (e.g., files that store one or more modules, sub programs, or portions of code).
  • a computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network.
  • the processes and logic flows described in this specification can be performed by one or more programmable processors executing one or more computer programs to perform actions by operating on input data and generating output.
  • the processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application specific integrated circuit).
  • processors suitable for the execution of a computer program include, by way of example, both general and special purpose microprocessors, and any one or more processors of any kind of digital computer.
  • a processor will receive instructions and data from a read only memory or a random access memory or both.
  • the essential elements of a computer are a processor for performing actions in accordance with instructions and one or more memory devices for storing instructions and data.
  • a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto optical disks, or optical disks.
  • mass storage devices for storing data, e.g., magnetic, magneto optical disks, or optical disks.
  • a computer need not have such devices.
  • a computer can be embedded in another device, e.g., a mobile telephone, a personal digital assistant (PDA), a mobile audio or video player, a game console, a Global Positioning System (GPS) receiver, or a portable storage device (e.g., a universal serial bus (USB) flash drive), to name just a few.
  • Devices suitable for storing computer program instructions and data include all forms of non volatile memory, media and memory devices, including by way of example semiconductor memory devices, e.g., EPROM, EEPROM, and flash memory devices; magnetic disks, e.g., internal hard disks or removable disks; magneto optical disks; and CD ROM and DVD-ROM disks.
  • the processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Epidemiology (AREA)
  • Data Mining & Analysis (AREA)
  • Business, Economics & Management (AREA)
  • Pathology (AREA)
  • Databases & Information Systems (AREA)
  • General Business, Economics & Management (AREA)
  • Chemical & Material Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Medicinal Chemistry (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Infusion, Injection, And Reservoir Apparatuses (AREA)

Abstract

L'invention concerne un système et un procédé pour identifier un modèle de dynamique de glucose personnalisé pour gérer le glucose d'un patient. Le procédé consiste à collecter des données de sortie de concentration de glucose du patient pendant une période. Les données d'entrée de taille des glucides et de perfusion d'insuline sont collectées à partir du patient pendant la période. Des paramètres d'un modèle de dynamique de glucose sont optimisés par résolution d'un problème d'optimisation itératif tel qu'un problème des moindres carrés non linéaire sur la base des données collectées parallèlement à l'établissement d'un temps et d'une quantité constants d'une valeur de correction. Le temps et la quantité d'une valeur de correction modifiée sont déterminés sur la base de l'exécution d'un programme quadratique entier mixte sur les données collectées parallèlement à l'établissement d'une valeur constante pour les paramètres du modèle de dynamique de glucose. Une valeur de correction finale et des paramètres optimisés finaux du modèle sont fournis sur la base de l'itération des paramètres d'optimisation et de la détermination de la valeur de correction modifiée.
PCT/US2023/085495 2022-12-22 2023-12-21 Identification de modèle avec des données d'entrée incomplètes de diabète de type i Ceased WO2024138024A2 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US202263434833P 2022-12-22 2022-12-22
US63/434,833 2022-12-22

Publications (2)

Publication Number Publication Date
WO2024138024A2 true WO2024138024A2 (fr) 2024-06-27
WO2024138024A3 WO2024138024A3 (fr) 2024-07-25

Family

ID=91590129

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2023/085495 Ceased WO2024138024A2 (fr) 2022-12-22 2023-12-21 Identification de modèle avec des données d'entrée incomplètes de diabète de type i

Country Status (1)

Country Link
WO (1) WO2024138024A2 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN120044213A (zh) * 2025-04-24 2025-05-27 内蒙古黄河工贸集团千里山煤焦化有限责任公司 基于小焦炉的炼焦指标监测分析方法、系统、设备和介质

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9057606B2 (en) * 2009-09-10 2015-06-16 Nextnav, Llc Wide area positioning system
WO2012142502A2 (fr) * 2011-04-15 2012-10-18 Dexcom Inc. Étalonnage avancé de capteur d'échantillon à analyser et détection d'erreur avancée
WO2015021441A1 (fr) * 2013-08-09 2015-02-12 Alere San Diego, Inc. Dosage et méthode de détermination d'insuline active

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN120044213A (zh) * 2025-04-24 2025-05-27 内蒙古黄河工贸集团千里山煤焦化有限责任公司 基于小焦炉的炼焦指标监测分析方法、系统、设备和介质

Also Published As

Publication number Publication date
WO2024138024A3 (fr) 2024-07-25

Similar Documents

Publication Publication Date Title
US11574742B2 (en) Diabetes management therapy advisor
US11024429B2 (en) Method, system and computer readable medium for adaptive and advisory control of diabetes
Garcia-Tirado et al. Advanced closed-loop control system improves postprandial glycemic control compared with a hybrid closed-loop system following unannounced meal
JP7823121B2 (ja) 血糖機能障害の評価および可視化
Robertson et al. Blood glucose prediction using artificial neural networks trained with the AIDA diabetes simulator: a proof‐of‐concept pilot study
Harvey et al. Quest for the artificial pancreas: combining technology with treatment
EP3844782B1 (fr) Prédiction de la dose d'insuline basée sur un horizon de réponse
Pesl et al. An advanced bolus calculator for type 1 diabetes: system architecture and usability results
US20210193285A1 (en) A system and method for use in disease treatment management
US8768673B2 (en) Computer-implemented system and method for improving glucose management through cloud-based modeling of circadian profiles
JP7712438B2 (ja) 機械学習システムによって低血糖の率を予測すること
JP7776421B2 (ja) 予測データと対応する受信データの差異を評価する統合状態推定予測
US20150347707A1 (en) Computer-Implemented System And Method For Improving Glucose Management Through Cloud-Based Modeling Of Circadian Profiles
WO2024138024A2 (fr) Identification de modèle avec des données d'entrée incomplètes de diabète de type i
Reenberg et al. High-performance uncertainty quantification in large-scale virtual clinical trials of closed-loop diabetes treatment
US11931548B2 (en) Method and system for determining optimal and recommended therapy parameters for diabetic subject
CN116529831A (zh) 具有基于风险的推注剂量校正的推注剂量建议器、无需碳水化合物的推注剂量推荐器和进餐确认
Ozaslan et al. Model identification with incomplete input data in type 1 diabetes
Zhang Robust optimal control for medical treatment decisions
Alhijailan et al. InsulinSync: A Type 1 Diabetes management application with AI-Decision Support System for Insulin Dosing and Computer Vision-Based Nutrition Prediction
Tejera-Perez et al. Algorithm-guided insulin therapy in hospitalized patients: current evidence, implementation challenges, and future perspectives
JP2026505427A (ja) インスリンの用量の決定におけるｑ学習の使用および関連するシステム、方法、および装置
García-Jaramillo et al. Interval simulator of the glucose-insulin system
JP2025515419A (ja) 基礎インスリン用量を滴定するためのシステム及び方法
Martínez Biomedical Engineering Degree

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23908569

Country of ref document: EP

Kind code of ref document: A2

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 23908569

Country of ref document: EP

Kind code of ref document: A2