Assessment and visualization of glycemic dysfunction

Abstract

To provide evaluation and visualization of glycemic dysfunction.SOLUTION: An amount of glycemic dysfunction associated with mis-timing (e.g., delay) of meal boluses based on replay analysis is determined. The amount of dysfunction of historical or estimated bolusing compared to an optimally timed bolus based on the replay analysis is quantified and visualized. Inferences may be made about diabetes meal management regarding inputs from a patient.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] INCORPORATION-BY-REFERENCE TO RELATED APPLICATIONS This application claims priority to U.S. Provisional Patent Application No. 62 / 786,149, entitled "EVALUATION AND VISUALIZATION OF GLYCEMIC DYSFUNCTION," filed December 28, 2018. The foregoing application is incorporated by reference in its entirety and expressly made a part hereof. [Background technology]

[0002] With the widespread adoption of CGM (continuous glucose monitoring) and connected devices, the availability and reliability of glucose time-series data has increased in recent years. However, despite the availability of reliable glucose data, accurately tracking insulin and meal data and optimally and effectively timing mealtime insulin boluses continues to be a problem for many people with diabetes, resulting in poor glucose control. Summary of the Invention [Means for solving the problem]

[0003] Based on the replay analysis, the amount of glycemic dysfunction associated with mistimed (e.g., delayed) meal boluses is determined. The amount of dysfunction of past or estimated boluses compared to optimally timed boluses based on the replay analysis is quantified and visualized. Inferences can be made about diabetes dietary management regarding input from the patient.

[0004] According to one embodiment, the estimated time, user-announced time, adjusted estimated time, or adjusted user-announced time for a meal is compared to when the subject administered insulin, and the benefits (e.g., risk reduction) of altering the insulin to an optimal dose and optimal pre-meal timing are evaluated.

[0005] In an implementation, a method for determining an amount of glycemic dysfunction is provided, the method including receiving continuous glucose monitoring (CGM) and insulin data for a subject, the CGM and insulin data including insulin bolus amount and timing of the insulin bolus amount, performing a replay analysis of the CGM and insulin data, quantifying the amount of glycemic dysfunction using the replay analysis, and providing an output representative of the amount of glycemic dysfunction.

[0006] In an implementation, a system for determining an amount of glycemic dysfunction is provided, the system including a processor that receives continuous glucose monitoring (CGM) and insulin data for a subject, the CGM and insulin data including insulin bolus amounts and timing of the insulin bolus amounts, a replay analyzer that generates a replay analysis using the CGM and insulin data, a quantifier that uses the replay analysis to quantify the amount of glycemic dysfunction, and an output device that provides an output representative of the amount of glycemic dysfunction.

[0007] In an implementation, a method for determining an amount of glycemic dysfunction is provided, the method including receiving historical blood glucose (BG) data and historical insulin data with a replay analyzer, performing a plurality of replay analyses of past boluses and estimated meals using the historical BG data and the historical insulin data, performing at least one residual insulin (IOB) data analysis using at least one of the plurality of replay analyses, performing at least one BG data analysis using at least one of the plurality of replay analyses, performing at least one discrepancy analysis and at least one comparison using the at least one IOB data analysis and the at least one BG data analysis, and providing an output based on the at least one discrepancy analysis and the at least one comparison.

[0008] In an implementation, a system for determining an amount of glycemic dysfunction includes an input device configured to receive blood glucose (BG) data, insulin data, and a plurality of functions, a replay analyzer configured to receive the BG data, the insulin data, and the plurality of functions and to perform a plurality of replay analyses on past boluses and on estimated meals using the BG data, the insulin data, and the plurality of functions, a residual insulin (IOB) assessor configured to perform a plurality of IOB assessments using at least one of the insulin data or the plurality of replay analyses, a BG risk profiler configured to perform a plurality of BG risk profiles using the BG data or at least one of the plurality of replay analyses, and a replay analyzer configured to perform a plurality of IOB assessments using the BG data or the plurality of replay analyses. an IOB profiler configured to perform a plurality of IOB profiles using at least one of the IOB risk profiles; a BG risk comparator configured to determine at least one of a risk profile plot or a BG risk score using the output of the BG risk profiler; a discrepancy analyzer configured to perform a plurality of discrepancy analyses using the IOB profiles and the IOB assessment to determine at least one of an IOB profile plot, an IOB discrepancy plot, or an IOB trace similarity score; and an output device configured to output at least one of the risk profile plot, the BG risk score, the IOB profile plot, the IOB discrepancy plot, or the IOB trace similarity score.

[0009] In one implementation, a method for providing a recommendation based on glycemic risk is provided, the method including: receiving at least one of a subject's past blood glucose (BG) data or past insulin data, the BG data and insulin data including an insulin bolus amount and a timing of the insulin bolus amount; determining a first risk score for glycemic dysfunction for the subject based on the at least one of the past BG data or the past insulin data; receiving at least one of optimal BG data or optimal insulin data for the subject; determining a second risk score for glycemic dysfunction for the subject based on the at least one of the past BG data or the past insulin data; comparing the first risk score with the second risk score; and outputting a recommendation based on the comparison of the first risk score and the second risk score.

[0010] In an implementation, a method for assessing glycemic risk is provided, the method including: determining a historical risk score for glycemic dysfunction based on at least one of historical blood glucose (BG) data or historical insulin data, where the historical BG data and the historical insulin data include an insulin bolus amount and a timing of the insulin bolus amount; determining an optimal risk score for glycemic dysfunction based on at least one of optimal BG data or optimal insulin data; determining a difference between the historical risk score and the optimal risk score; and providing a recommendation based on the difference and based on a threshold value.

[0011] In an implementation, the method includes calculating a first IOB trace based on the actual bolus recording, calculating a second IOB trace of an optimal therapy, comparing the first IOB trace to the second IOB trace to determine an IOB similarity score, and using the IOB similarity score to modify at least one of smart alert behavior, user recommendations, reporting functions, and visualizations.

[0012] In an implementation, the system includes a discrepancy analyzer configured to determine an IOB similarity score using a comparison of a first IOB trace based on the actual bolus recording and a second IOB trace of the optimal treatment, and an output device configured to output the IOB similarity score.

[0013] This Summary is provided to introduce some concepts in a simplified form that are further described below in the Detailed Description. This Summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used to limit the scope of the claimed subject matter.

[0014] The foregoing summary, as well as the following detailed description of exemplary embodiments, will be better understood when read in conjunction with the appended drawings. For the purpose of illustrating the embodiments, there are shown in the drawings exemplary configurations of the embodiments; however, the embodiments are not limited to the specific methods and instrumentalities disclosed. The drawings are as follows: [Brief explanation of the drawings]

[0015] [Figure 1] FIG. 1 is a high-level functional block diagram of one embodiment of the present invention. [Figure 2] 1 is an illustration of an exemplary environment for assessing and visualizing glycemic dysfunction. [Figure 3A] 1 shows a schematic block diagram of an exemplary system implementation with related methods for assessing and visualizing glycemic dysfunction. [Figure 3B] 1 shows a schematic block diagram of an exemplary system implementation with related methods for assessing and visualizing glycemic dysfunction. [Figure 4] 1A-1C are diagrams of various risk profile plots, according to an embodiment. [Figure 5] FIG. 10 is a diagram of risk zones, according to an embodiment. [Figure 6] FIG. 10 is a diagram of zones that may be significantly improved by altering bolus dosage or bolus timing, according to an example. [Figure 7A] 1A-1C are diagrams of various residual insulin (IOB) profile plots according to an example. [Figure 7B] 1A-1C are diagrams of various residual insulin (IOB) profile plots according to an example. [Figure 8] FIG. 10 is an IOB discrepancy plot of a bolus dose at a past time point, according to an example. [Figure 9] FIG. 1 is an illustration of an IOB mismatch zone plot, according to an example. [Figure 10] FIG. 10 is an IOB discrepancy plot for optimal bolus administration at estimated mealtime, according to an example. [Figure 11] FIG. 10 is an IOB discrepancy error zone plot for optimal bolus administration at estimated mealtime, according to an example. [Figure 12] FIG. 10 is an exemplary summary report. [Figure 13] 10A-10C are diagrams of various risk profile plots according to another embodiment. [Figure 14A] 10A-10C are diagrams of various IOB profile plots according to another embodiment. [Figure 14B] 10A-10C are diagrams of various IOB profile plots according to another embodiment. [Figure 15A] 10A-10C are diagrams of various IOB mismatch plots according to another embodiment. [Figure 15B] 10A-10C are diagrams of various IOB mismatch plots according to another embodiment. [Figure 16] 1 is an operational flow of an implementation of a method for determining the amount of glycemic dysfunction. [Figure 17] 1 is an operational flow of another implementation of a method for determining the amount of glycemic dysfunction. [Figure 18] 1 is an operational flow of another implementation of a method for determining the amount of glycemic dysfunction. [Figure 19] 1 is an operational flow of another implementation of a method for determining the amount of glycemic dysfunction. [Figure 20]1 is an operational flow of another implementation of a method for determining the amount of glycemic dysfunction. [Figure 21] 1 illustrates an exemplary computing environment in which exemplary embodiments and aspects may be implemented. DETAILED DESCRIPTION OF THE INVENTION

[0016] The claimed subject matter is described with reference to the drawings, wherein like reference numerals are used to refer to like elements throughout. In the following description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the claimed subject matter. It may be apparent, however, that the claimed subject matter can be practiced without these specific details. In other instances, structures and devices are shown in block diagram form to facilitate description of the claimed subject matter.

[0017] The present invention is intended for replay use, meaning that the impact of an inefficient or suboptimal bolus in a subject's own historical data can be explained in terms of glycemic risk or insulin delivery profile relative to optimal bolus dosing in the same dataset. Replay bolus timing effects are determined. Pre-meal bolus compliance is determined. Implementations of the present invention aim to quantify the impact of meal misinformation (both in terms of timing and configuration) in a specific way. Insulin benefits (e.g., risk reduction) are assessed with respect to optimal doses and with respect to changes in optimal pre-meal timing. A meal can be of any length, depending on the implementation. A meal can be an event having a certain length or an event spanning a time window. A time window can be of any duration or length, depending on the implementation. In some implementations, the length of time or time window can vary from meal to meal or can be a fixed duration. Those skilled in the art will understand that a meal is a time window, and that time window can be defined in various ways.

[0018] As further described herein, the amount of glycemic dysfunction associated with a subject's mistimed (e.g., delayed) meal bolus may be determined based on the replay analysis, and the amount of dysfunction of the past or estimated bolus administration compared to an optimally timed bolus based on the replay analysis may then be quantified and output in a visual form for the subject or user. Inferences about diabetes dietary management may be made using input provided by the subject or user. Adjustments to past and replay meal boluses may be used to distinguish between two categories for improving diabetes management: mistimed boluses and missed meal boluses.

[0019] The embodiments described herein address the problem of identifying patterns of inefficient and mistimed bolus dosing and its impact on blood glucose levels in diabetic patients. Complicating this problem is the fact that a patient's perception of bolus efficiency and / or timing is determined by whether they consistently follow prescribed insulin self-treatment guidelines. For example, under physician guidance to administer boluses at or immediately before meals, patients may generally believe that they are already maximizing the benefits of their insulin therapy by implementing this strategy for only two of their three daily meals. However, delaying bolus administration for the third meal may actually contribute disproportionately to glycemic dysfunction in the form of hypoglycemic or hyperglycemic exposure. Vague memories of "compliant behavior" often fall short of accurate assessment, even with thorough analysis of retrospective data. As a result, conventional methods involving subjective recall or partially complete data records do not clearly convey to patients and their caregivers the existence or impact of inefficient insulin therapy.

[0020] Implementations described herein provide a technically non-trivial solution to (i) deriving a mechanism by thorough analysis of blood glucose records, associated records of insulin data, and possibly other metabolic-related information, (ii) a series of meal events (including carbohydrate, fat, and protein content) from meals published or estimated in conjunction with blood glucose and insulin records, and (iii) deriving a mechanism by which the glycemic impact of alternative dosing can be estimated from past records, (iv) systematically applying this mechanism to determine the impact on blood glucose levels of numerically optimal dosing at past bolus times and at validated meal event times, including performing replay simulations at past boluses or before estimated meals, (v) comprehensively quantifying glycemic dysfunction associated with inefficient insulin delivery, including the impact in terms of risk of exposure to hypoglycemia and hyperglycemia, and (vi) considering systematic discrepancies in insulin delivery between past boluses and numerically optimal boluses, including an overall similarity score between optimal bolus dosing and past insulin dosing patterns, and identifying times of day when, on average, past insulin dosing deviates from numerically optimal insulin dosing.

[0021] Advantages are offered over existing self-assessment or clinical tools. The systems and methods described herein are systematic, taking into account all available blood glucose and insulin data, far surpassing "pencil and paper" analysis and even "spreadsheet" analysis using standard statistical assessments. The systems and methods are comprehensive in identifying correlations in the complex relationships between blood glucose, insulin, and food records through the lens of patient-specific mathematical models. In this way, patterns of inefficient bowling and patterns of glycemic abnormalities are revealed through statistical analysis of these patient-specific interactions. The systems and methods described herein are not subject to cognitive biases introduced (i) by patients when recalling insulin delivery and / or meal events or (ii) by physicians when formulating questions about the patient's experience, which can lead to misinterpretation of the patient's treatment routine.

[0022] Depending on the implementation, various calculated metrics can be used, including, but not limited to, for example, hg_min, the minimum level of historical glycemic dysfunction required (which can be measured as BG risk, time in range, average, etc. and can be broken down by time of day in hours, minutes, or other useful increments), gdiff_min, the minimum level of improvement in glycemic outcome measured as the difference between historical and simulated glucose-derived metrics, and ifdiff_min, the minimum level of normalized change in insulin treatment required for improved outcome (which can be frequency-based and / or magnitude-based).

[0023] FIG. 1 is a high-level functional block diagram 100 of an embodiment of the present invention. A processor 130 communicates with an insulin device 110 and a glucose monitor 120. The insulin device 110 and the glucose monitor 120 communicate with a subject 140 to deliver insulin to the subject 140 and monitor the subject's 140 glucose level, respectively. The processor 130 is configured to perform calculations described further herein. The insulin device 110 and the glucose monitor 120 may be implemented as separate devices, as a single device, within a single device, or across multiple devices. The processor 130 may be implemented locally to the insulin device 110, the glucose monitor 120, or as a standalone device, indicated by dashed lines as standalone device 132 (or in any combination of two or more of the insulin device 110, the glucose monitor 120, or a standalone device). The processor 130 or portions of the system may be located remotely, such as in a server or cloud-based system.

[0024] Examples of insulin devices, such as insulin device 110, include insulin syringes, external pumps, and patch pumps that deliver insulin to a subject, typically into subcutaneous tissue. Insulin device 110 also includes devices that deliver insulin by different means, such as insulin inhalers, insulin jet injectors, intravenous infusion pumps, and implantable insulin pumps. In some embodiments, a subject will use a combination of two or more insulin delivery devices, for example, injecting long-acting insulin with a syringe and using inhaled insulin before meals. In other embodiments, these devices can deliver other drugs that help control glucose levels, such as glucagon, pramlintide, or glucose-like peptide-1 (GLP1).

[0025] Examples of glucose monitors, such as glucose monitor 120, include continuous glucose monitors that record glucose values ​​at regular intervals, such as every 1, 5, or 10 minutes. These continuous glucose monitors can be, for example, transcutaneously inserted, fully implanted, or can use electrochemical or optical sensors that measure tissue non-invasively. Examples of glucose monitors, such as glucose monitor 120, also include devices that periodically draw blood or other fluids to measure glucose, such as intravenous blood glucose monitors, microperfusion sampling, or periodic fingersticks. In some embodiments, glucose measurements are provided in near real time. In other embodiments, glucose measurements determined by the glucose monitor can be stored within the glucose monitor itself for subsequent retrieval.

[0026] The insulin device 110, glucose monitor 120, and standalone device 132 can be implemented using a variety of computing devices, such as smartphones, desktop computers, laptop computers, and tablets. Other types of computing devices may be supported. A suitable computing device is shown in FIG. 21 as computing device 2100 and cloud-based applications.

[0027] The insulin device 110, glucose monitor 120, and standalone device 132 can communicate over a network. The network can be a variety of network types, including the public switched telephone network (PSTN), a cellular network, and a packet-switched network (e.g., the Internet). While only one insulin device 110, one glucose monitor 120, and one processor 130 (in one standalone device 132) are shown in FIG. 1, there is no limit to the number of insulin devices, glucose monitors, processors, and standalone devices that can be supported. An activity monitor 150 and / or smartphone 160 can also be used to collect dietary and / or activity data from or associated with the subject 140 and provide the dietary and / or activity data to the processor 130.

[0028] The processor 130 may execute an operating system and one or more applications. The operating system may control which applications are executed by the insulin device 110, the glucose monitor 120, and the standalone device 132 and how the applications interact with one or more sensors, services, or other resources of the insulin device 110, the glucose monitor 120, and the standalone device 132.

[0029] The processor 130 receives data from the insulin device 110 and glucose monitor 120, and in some implementations, from the subject 140, and performs a replay analysis to evaluate the patient's current meal bolus dosing strategy and compare it with alternative bolus dosing times and insulin amounts.

[0030] More specifically, Figure 2 is a diagram of an exemplary environment 200 for assessing and visualizing glycemic dysfunction. Computing device 201 includes one or more processors, such as processor 130, that can implement a replay analyzer 220, an residual insulin (IOB) assessor 225, a blood glucose (BG) profiler 230, an IOB profiler 235, a BG risk comparator 240, a discrepancy analyzer 245, and an output device 250. Computing device 201 can be implemented in the computing device shown in Figure 21 as computing device 2100, and can also be implemented using a processor located on a microcontroller or within an insulin device 110, glucose monitor 120, or other insulin delivery device, or as a standalone component (or a combination of two or more of insulin device 110, glucose monitor 120, or other insulin delivery device, or standalone component).

[0031] The computing device 201 is in communication with the target computing device 205, for example, via a network and / or a communication interface, and is configured to receive data 210 and functions 215 via the network and / or communication interface from the insulin device 110 and the glucose monitor 120, as well as one or more devices or processors maintained by or associated with a physician, caregiver, healthcare provider, hospital, clinic, university, etc.

[0032] A subject computing device 205 (e.g., activity monitor 150, smartphone 160, or other type of computing device, such as computing device 2100 described with respect to FIG. 21 ) may be used by a subject 140 to input data (subject-provided data 207) regarding the subject's meals and bolus doses (e.g., timing, amount, etc.). In some implementations, this data may be considered unreliable because it is subject to human error, bias, and / or behavioral inconsistencies or anomalies. Data 210 may include data that is not subject-provided data 207, such as past blood glucose levels, past insulin doses, dietary information, the subject's initial (e.g., pre-prandial) state, insulin dosing strategy, etc. In some implementations, this data is considered reliable because it is machine-provided, machine-verified, and determined to have a confidence value above a predetermined threshold, and / or is free from, resistant to, and / or unaffected by human error, bias, or behavioral inconsistencies or anomalies. The functions 215 may include one or more metabolic modeling functions, such as a diet compiler (e.g., diet compiler 303) and a replay compiler (e.g., replay compiler 304) that generates a replay prediction function (e.g., replay prediction function 305), as well as one or more risk, profiling, and assessment functions as further described herein.

[0033] The replay analyzer 220 uses the subject-provided data 207, data 210, and functions 215 to perform one or more replay analyses, such as replay with optimal boluses for past bolus times and replay with optimal boluses for assessed meal times, as further described herein. Depending on the implementation, one or more of the outputs from the replay analyzer 220 may be provided to one or more of the IOB assessor 225, the BG risk profiler 230, and the IOB profiler 235. Thus, in implementations, the replay analyzer 220 can function as a data processor that receives CGM and insulin data associated with the subject. In some implementations, alternatively or additionally, the replay analyzer 220 can function as a data generator that generates CGM and insulin data associated with the subject.

[0034] IOB assessor 225 can perform one or more IOB assessments, such as an IOB assessment of historical data, an IOB assessment of optimal bolus at past times, and an IOB assessment of optimal bolus at estimated meal times, etc. Depending on the implementation, one or more of the outputs from IOB assessor 225 can be provided to one or more of IOB profiler 235 and discrepancy analyzer 245.

[0035] The BG risk profiler 230 may generate one or more BG risk profiles, such as a BG risk profile for optimal bolus at a past time and a BG risk profile for optimal bolus at an estimated mealtime. Depending on the implementation, one or more of the outputs from the BG risk profiler 230 may be provided to one or more of the BG risk comparator 240 and the discrepancy analyzer 245. Thus, in an implementation, the BG risk profiler 230 functions as a quantifier that quantifies the amount of glycemic dysfunction using replay analysis.

[0036] The IOB profiler 235 may generate one or more IOB profiles, such as an IOB profiling for an optimal bolus at a past time and an IOB profiling for an optimal bolus at an estimated mealtime, etc. Depending on the implementation, one or more of the outputs from the IOB profiler 235 may be provided to a discrepancy analyzer 245.

[0037] The BG risk comparator 240 performs a BG risk profile comparison function to generate a risk profile plot and an overall BG risk score, which may be stored, displayed, and / or transmitted by the output device 250 .

[0038] The discrepancy analyzer 245 may perform one or more IOB discrepancy analyses and provide one or more plots and scores to the output device 250 for storage, display, and / or transmission.

[0039] 3A and 3B show a schematic block diagram of an exemplary implementation of a system 300 in a related method for assessing and visualizing glycemic dysfunction, directed to a replay analysis tool that evaluates a subject's current meal bolus strategy and compares it with alternative bolus times and insulin amounts.

[0040] The output of system 300 can be used to identify behavioral problems, for example, by isolating the effects of bolus timing, carbohydrate count, and carbohydrate ratio relative to a subject's meal. Replay simulation analysis can be used to assess the impact of a numerically optimal bolus (removing inaccurate carbohydrate counts or improper carbohydrate ratios) on past bolus times, and to assess the impact of a numerically optimal bolus on estimated past meal times.

[0041] In an implementation, once the replay simulator determines it is time to calculate a numerically optimal bolus, it generates a set of candidate boluses, simulates future BG values ​​associated with each candidate bolus, "scores" the simulated BG trajectory for each candidate bolus using risk analysis, such as rate-dependent risk analysis, selects and implements the best bolus, and continues replaying the simulation until it is time for the next bolus. In addition to or instead of replaying candidate boluses at their exact times, replays can be made at times within a target window near adjusted meal times or past medication times.

[0042] The system 300 uses accurate assessment of both meal timing and amount in historical data. Inputs, including data 210 (e.g., received from a storage device) and / or subject-provided data 207 from the subject computing device 205, are directed to CGM (continuous glucose monitoring) data and insulin data, including insulin bolus amounts and timing. In an implementation, historical BG data 301 ({bg(t)} t=0、1、 ...) and past insulin 302 ({i(t)} t=0、1、 …) is a time series and is used as input. Note that this is the original data.

[0043] Undocumented mealtimes are estimated using CGM data (e.g., historical BG data 301) and insulin data (e.g., historical insulin data 302), and possibly user acknowledgement of meals, to estimate the time and amount of all carbohydrate events. Other inputs that may be used in other implementations include exercise (e.g., to determine how to optimize therapy for exercise), rescue carbohydrates (e.g., to determine the amount needed), other meal bolus strategies (e.g., split, delayed, extended bolus), and Type 2 treatments including oral medications, etc. This data may be reliable or unreliable.

[0044] In an implementation, the CGM and insulin data are "discretized" in time (e.g., snapped to the nearest 5-minute sample) and aggregated appropriately. For example, if there are two or more boluses in a 5-minute interval, the insulin value for the corresponding sample will be the sum of the boluses (and the basal effect). This CGM and insulin data form the basic input for replay state estimation, which is then the step for building a replay simulator for that dataset.

[0045] The historical BG data 301 and historical insulin data 302 are provided to a meal compiler 303 and a replay compiler 304, which then output a replay prediction function 305 and an estimated meal assessed sequence 306. The replay prediction function 305 provides the ability to predict, from any point in the historical record, the blood glucose outcomes that would be associated with an alternative insulin delivery (i.e., deviating from the historical insulin delivery). Specifically, the alternative insulin time series 302 ({i alt (t)} t=0,1,… ) is used to predict the time series of blood glucose data ({bg alt (t)} t=0,1,… (which can be expressed as

[0046] The meal compiler 303 analyzes historical BG data and insulin delivery data, and possibly other sources of historical data, including meal events recorded and / or recalled directly from the patient. The output sequence of estimated meals 306 represents the reconciliation between (i) meal events that can be detected / estimated from correlations between historical blood glucose and insulin delivery data and (ii) meal events inferred or directly recalled from other sources of data. The detected / estimated meals from historical blood glucose and insulin data can be derived from various mechanisms, including deconvolution or normalized deconvolution (see, e.g., Patek SD et al., “Empirical Representation of Blood Glucose Variability in a Compartmental Model,” Kirchsteiger H. et al. (eds) Prediction Methods for Blood Glucose Concentration (Lecture Notes in Bioengineering), 2016).

[0047] Meal times can also be estimated from the same data set using other approaches known to those skilled in the art of metabolic modeling. Metabolic models describe whole-body and local glucose metabolism (see, for example, Cobelli C. et al., "Modeling Glucose Metabolism in Man: Theory and Practice," Hormone and Metabolic Research, Supplemental Series, 1990, 24:1-10, https: / / www.ncbi.nlm.nih.gov / pubmed / 2272612). These models describe the levels and distribution of glucose and other metabolites in different compartments of the body, the levels and distribution of insulin and other hormones and drugs in different compartments of the body, and the kinetics of glucose transport between compartments, representing consumption and production. These models can be used in clinical settings, where samples (e.g., venous blood or interstitial fluid) are taken and clinical analyzers are used to track compartmental inputs and outputs. Tracer studies also exist that track the distribution of metabolites. In these studies, inputs (e.g., diet, insulin, exercise) and outputs are carefully controlled and measured. The resulting metabolic model can provide a patient-specific model of glucose excursions from a meal, which can be compared with the patient's glucose monitor to estimate the characteristics of the meal.

[0048] There are also methods for detecting the timing and magnitude of meal excursions from patterns and features in glucose monitor traces. Examples include detecting accelerations or upward changes in CGM data or fitting meal excursion models (e.g., minimal models) to the time series. Meal information can also be obtained directly from clinical studies or patient records. Reconciliation between detected / estimated meals and sources of uncertain or unreliable meal data can be achieved, for example, by clustering detected and reported meals and averaging meal content. An assessed sequence 306 of estimated meals is determined and output from the compilers 303, 304.

[0049] In implementations, the meal compiler 303 estimates meal times and / or effective sizes from past blood glucose measurements (e.g., CGM), insulin dosages (e.g., pump or pen), and optionally meal announcements.

[0050] The replay compiler 304 also analyzes other sources of historical data, including historical BG and insulin delivery data, and potentially meal events recorded directly by the patient and / or recalled. This generates a replay prediction function 305, a mathematical function that can be expressed in a variety of ways (including computer algorithms) to provide the ability to predict, from any point in the historical record, the blood glucose outcomes that would be associated with alternative insulin deliveries (i.e., deviations from past insulin deliveries). The replay prediction function 305 evaluates whether alternative insulin dosing strategies for the same meal would result in better outcomes (e.g., lower glycemic risk). The replay prediction function 305 takes as input the patient's initial state (pre-meal), one or more meal events, and insulin dosing strategy, and maps them to the resulting glucose excursion for that meal / day. Replay prediction function functionality can be derived from blood glucose data by inserting validated sequences of estimated meals into the mathematical models from which they are derived (see, e.g., Patek SD et al., "Empirical Representation of Blood Glucose Variability in a Compartmental Model," Kirchsteiger H. et al. (eds) Prediction Methods for Blood Glucose Concentration (Lecture Notes in Bioengineering), 2016). Alternatively, other blood glucose simulation methodologies can be applied to validated sequences of estimated meals, including simulators derived from compartmental models or data-driven methods.

[0051] The past BG 301 may also be provided to a BG risk profiler (such as BG risk profiler 230) for use in determining BG risk profiling for an optimal bolus at a past time 322. The past insulin 302 may also be provided to an IOB assessor (such as IOB assessor 225) for use in determining an IOB assessment of the past data 316, and may also be provided to a replay analyzer (such as replay analyzer 220) for use in replaying an optimal bolus at the time of the past bolus 312.

[0052] Other functions that the system 300 can use include a BG risk function 307 used to identify a numerically optimal bolus to insert into the replay prediction function 305 at a past bolus time or at a verified meal event time, a BG risk profiling function 308 used to derive an aggregate assessment of the risk of exposure to hypoglycemia and hyperglycemia in either past blood glucose records or replayed blood glucose records from the replay prediction function, and an IOB assessment function 310.

[0053] The IOB assessment function 310 is a calculation that tracks how much insulin remains in the body from the previous dose, including the time course over which the injected insulin acts in the body, including how the insulin is absorbed by the body, how it begins to affect glucose levels, and how it clears the body. For example, a rapid-acting insulin may reach peak action in one hour, with a total duration of action of four hours. The IOB function can take various mathematical forms to represent the time it takes for insulin to reach its peak value and the subsequent decay curve. These IOB functions are used in bolus calculators (insulin pumps and phone apps) that help diabetic subjects estimate bolus insulin dosages (see, e.g., Zisser H. et al., "Bolus Calculator: A Review of Four 'Smart' Insulin Pumps," Diabetes Technology & Therapeutics, December 2008, 10(6):441-4, doi:10.1089 / dia.2007.0284, https: / / www.ncbi.nlm.nih.gov / pubmed / 19049372). In some embodiments, the IOB function uses parameters provided by the patient, their healthcare provider, or the insulin manufacturer. In other embodiments, the IOB assessment function can be optimized to generate a patient-specific insulin action curve.

[0054] In an implementation, the regenerative prediction function 305 maps the original data (e.g., meals and insulin) to the observations (e.g., CGM) with sufficient accuracy, maps alternative dosing strategies to the resulting glucose excursions in a way that models the patient's physiology (e.g., insulin activity and carbohydrate sensitivity), and provides a reliable and interpretable solution to the problem, such as an estimate of IOB that is a stable, smooth function.

[0055] In an implementation, known inputs and user-declared meals are used to generate a function that can replay alternative insulin strategies under the same input conditions. This is a structure where a function takes input data, extracts knowledge, and creates a new function.

[0056] The replay prediction function 305 can also be constructed from the same data set using other approaches known to those skilled in the art of metabolic modeling. One exemplary approach is to fit the parameters of a metabolic model to subject data or clinical studies. Another possible approach is to train a neural network on similar data to predict glucose excursions.

[0057] For system 300 with respect to elements 312-342, taking into account behavioral issues such as missed meals or erroneous user input, the system determines and quantifies glycemic dysfunction associated with improper timing of boluses relative to meals based on CGM and insulin data, which can be used to identify behavioral issues and to isolate the effects of bolus timing relative to meals, as well as carbohydrate counting and carbohydrate ratios, as described further herein, e.g., with respect to Figures 4-15B.

[0058] At 312, a replay of the optimal bolus at the time of the past bolus is performed (e.g., by the replay analyzer 220) using as input the replay prediction function 305 and the past insulin 302. The output of 312 includes the replayed CGM for the optimal bolus at the past time and the replayed insulin for the optimal bolus at the past time.

[0059] At 314, a replay of the optimal bolus at the estimated meal time is performed (e.g., by the replay analyzer 220) using the replay prediction function 305 and the estimated meal data 306 as inputs. The output of 314 includes a replayed CGM for the optimal bolus at the estimated meal time and a replayed insulin for the optimal bolus at the estimated meal time. Note that the output of 312 and 314 is simulated data using replay. This simulates how the risk, IOB, BG, and related data and outcomes would differ if the subject were to receive a bolus at these times.

[0060] At 316, IOB assessment of past data is performed (e.g., by IOB assessor 225) using past insulin 302 as input. Any known technique for performing IOB assessment can be used. The output of 316 is used in IOB profiling of optimal boluses at 324 and at past times at 338 and 342, which are further described below.

[0061] At 318, IOB assessment of the optimal bolus at past time is performed (e.g., by IOB assessor 225) using as input the replayed insulin of the optimal bolus at past time from 312. The output of 318 is used in IOB profiling of the optimal bolus at past time at 328 and at 338, which are further described below.

[0062] At 320, an IOB assessment of the optimal bolus at the estimated mealtime is performed (e.g., by IOB assessor 225) using the replayed insulin of the optimal bolus at the estimated mealtime as input. The output of 320 is used in IOB profiling of the optimal bolus at the estimated mealtime at 332 and at 342, which are further described below.

[0063] At 322, BG risk profiling of optimal bolus at past times is performed (e.g., by BG risk profiler 230) using past BG 301 as input. The output is used by BG risk profile comparison function 334, which is further described below.

[0064] As described above, IOB profiling of the optimal bolus in the past is performed (e.g., by IOB profiler 235) at 324. The output of 324 is used at 336 and 340, which are further described below.

[0065] At 326, BG risk profiling of the optimal bolus at past time is performed (e.g., by BG risk profiler 230) using the replayed CGM of the optimal bolus at past time generated at 312. The output of 326 is used at 334, which is further described below.

[0066] At 328, IOB profiling of the optimal bolus at past times is performed (e.g., by IOB profiler 235) using as input the output of the IOB assessment of the optimal bolus at past times 318. The output of 328 is used at 336, which is further described below.

[0067] At 330, BG risk profiling for the estimated optimal bolus at mealtimes is performed (e.g., by BG risk profiler 230) using the replayed CGM of the estimated optimal bolus at mealtimes generated at 314. The output of 330 is used at 334, which is further described below.

[0068] At 332, IOB profiling for optimal bolus at estimated mealtime is performed (e.g., by IOB profiler 235) using as input the output of IOB assessment for optimal bolus at estimated mealtime 320. The output of 332 is used at 340, which is further described below.

[0069] With respect to output, charts and / or reports can be generated, transmitted, and / or displayed to show that a better outcome for the subject would have been achieved if the patient had not had the behavioral problem. In some implementations, the methods and systems contemplated herein determine whether to display certain recommendations, make certain graphical inferences, or perform other actions based on thresholds as further described herein.

[0070] The comparison of BG risk profiles is performed (e.g., by BG risk comparator 240) using a BG risk profile comparison function 334 that receives as input the outputs of the BG risk profile for optimal bolus at past time 322 and 326 and the BG risk profile for optimal bolus at estimated mealtime 330. The output of the BG risk profile comparison function 334 may include a risk profile plot and an overall BG risk score, each having data and / or indicators directed to the past information, optimal bolus at past bolus time information, and optimal bolus at estimated mealtime information, e.g., as further described herein with respect to Figures 4-15B.

[0071] The difference in the comparative BG risk profiles can be quantified to generate numerical outputs for various applications and / or feedback loops, including historical CGM, replay-simulated CGM with numerically optimal boluses at historical bolus times, and replay-simulated CGM with numerically optimal boluses inserted at estimated meal times.

[0072] Differences in comparative BG risk profiles can be quantified in various ways, such as using plots and visualizations, using risk indices indicating high blood glucose risk, low glucose risk, total glycemic risk, and showing the number and extent of excursions into out-of-range BG levels, for example, as described with respect to FIG. 4.

[0073] At 336, an IOB profile discrepancy analysis is performed (e.g., by discrepancy analyzer 245) to compare past insulin with optimal bolus at past times. The input to 336 includes the IOB profiling outputs of optimal bolus at past times of 324 and 328, and the output of 336 includes data and / or indicators directed to the past information and an IOB profile plot with optimal bolus at past bolus time information.

[0074] At 338, an IOB trace discrepancy analysis is performed (e.g., by discrepancy analyzer 245) to compare past insulin to the optimal bolus at the past time. The input to 338 includes the output of the IOB assessment of past data at 316 and the output of the IOB assessment of the optimal bolus at the past time at 318. The output of 338 includes an IOB discrepancy plot and an IOB trace similarity score with data and / or indicators directed to the optimal bolus at the past bolus time versus the past bolus information.

[0075] At 340, an IOB profile discrepancy analysis is performed (e.g., by discrepancy analyzer 245) comparing past insulin to the optimal bolus at the estimated mealtime. Inputs to 340 include the output of IOB profiling for optimal bolus at past time 324 and the output of IOB profiling for optimal bolus at estimated mealtime 332. Outputs of 340 include an IOB profile plot with data and / or indicators directed to the optimal bolus at past information and estimated mealtime information.

[0076] At 342, an IOB trace discrepancy analysis is performed (e.g., by discrepancy analyzer 245) to compare past insulin to the optimal bolus at the estimated mealtime. The input to 342 includes the output of the IOB assessment of past data at 316 and the output of the IOB assessment of the optimal bolus at the estimated mealtime at 320. The output of 342 includes an IOB discrepancy plot and an IOB trace similarity score with data and / or indicators directed to the optimal bolus at the estimated mealtime against the past bolus information.

[0077] Thus, the output of the replay simulation shows glycemic dysfunction in the past. For example, a historical / estimated risk profile on a graph can show differences in glycemic exposure. A historical CGM can be displayed against a replay simulation CGM with an optimal bolus at a past bolus time, and against a replay simulation CGM with an optimal bolus at an estimated meal time and past bolus time.

[0078] Risk indicators such as high blood glucose risk, low blood glucose risk, and total glycemic risk may be displayed. These risks may be displayed as representations on a chart including numbers, colors, quantified values, etc. The number and range of excursions into out-of-range BG levels may also be identified and shown. Numerical representations could be used to indicate glycemic function / good or bad, such as comparing two numbers to one another.

[0079] Another output may be an indication of compliance with the pre-meal bolus requirement. The subject's ideal pre-meal bolus is determined, and then the subject's deviation from this ideal pre-meal bolus is determined. A score may be generated and output, for example, numerically or in a chart. In this embodiment, quantifying glycemic dysfunction may be considered as quantifying compliance (or deviation) with the ideal pre-meal bolus timing. In implementations, visualizations (showing the effect of compliance) or recommendations (e.g., "bolus 20 minutes before this type of meal in the future") may be generated and output.

[0080] The various outputs described herein can be used as risk stratification tools to help healthcare professionals identify subjects with pre-meal bolus compliance and potential benefit from treatment adjustments. The outputs can be used to highlight the most impactful potential changes in clarity reports (e.g., changing insulin timing relative to bolus volume or highlighting the best meals). The outputs can be used to trigger coaching comments and discussion points via coaching calls or chatbot prompts.

[0081] 4 is a diagram of various risk profile plots 400, according to an embodiment. In this example, for the subject, the total risk CGM is 0.99, the total risk with optimal bolus at past bolus time is 0.46, the total risk with optimal bolus at estimated mealtime is 0.39, the IOB trace similarity score for optimal bolus at past bolus time is 0.69, and the IOB trace similarity score for optimal bolus at estimated mealtime is 0.59. The number of days in the sample was 13.

[0082] The total risk CGM is a value that captures both hyperglycemia and hypoglycemia risk over a time segment as a single value. Any known technique for determining the value can be used. In implementations, a positive value with a low value indicates good glycemic control. A total risk of 1 or less is a good outcome. A score of 2 or greater indicates larger or persistent excursions outside the range. The CGM risk in these examples is based on a target range of 95 to 195 mg / dL, so a score of 1 occurs if the subject is consistently within one of these limits. This analysis could be based on other target ranges or risk functions (e.g., 70 to 180 mg / dL), but these targets are consistent with medical concerns observed in clinical studies. Any threshold can be selected for the target zone. The CGM risk metric penalizes glucose excursions outside these limits (e.g., a parabolic function) that drive up the score, especially on the hypoglycemic side.

[0083] Regarding the total risk associated with optimal bolus dosing over a past time, Replay Bolus identifies the optimal amount of insulin at the same bolus time (e.g., pump or pen bolus time) as the recorded data. The optimal bolus amount halves the CGM risk. A reduction from 0.99 to 0.46 indicates an opportunity for improvement through therapy / bolus calculator changes. If Replay analysis indicates that the optimal strategy can significantly reduce the total risk (e.g., 25%, 50%, etc.), the subject or user may be prompted with a treatment suggestion (e.g., visualization or message) that changing insulin dosing can reduce their exposure to glycemic risk. There are numerous ways in which these potential risk reductions can be integrated into subsequent reports, analyses, and calculations, so the scope is not limited to displaying the results in this way.

[0084] Regarding the total risk associated with optimal bolus delivery at estimated mealtimes, adjusting the timing of the bolus can further reduce and improve the total risk. In this example, the subject's risk reduction is 0.46-0.39 compared to the optimal dose replay described above. Because this incremental risk reduction from timing changes is not as large as the previous step, it may indicate that the historical data included insulin boluses closer to mealtimes. In other words, the subject's risk is a matter of dose, not timing. In some embodiments, replay analysis (comparing historical to mealtime bolus delivery) can identify subject behavior that could be improved by dosing earlier or later than historical values. This is a comparison of recorded insulin bolus times to estimated / adjusted mealtimes. In addition to the pump, connected pens record when boluses are delivered, but there may be questions about priming and whether there are any ad hoc / additional doses from other pens or syringes.

[0085] The IOB trace's similarity score to the optimal bolus at the previous bolus time is calculated based on the actual bolus recording, and the optimal treatment IOB trace is calculated. These traces are compared as a dot product between the unit magnitude traces, which can be assigned a value ranging from 0 to 1, where 1 indicates that the actual and optimal dosing shapes are most similar. A value of 1 indicates the same pattern in the traces, even if the overall magnitude is different. In this example, a value of 0.69 indicates that the optimal bolus amount at the previous bolus time resulted in an IOB profile shape that differs from the current result. This similarity score can be thresholded or compared to current behavior, and the output of this comparison can change smart alert behavior, user recommendations, reporting, and visualization. Depending on the implementation, the information can be conveyed, for example, as a qualitative score of amount and timing, visualization, or interpretation (e.g., "poor," "good," "excellent," etc.). At a high level, the amount or timing of insulin dosing can help compare current behavior to optimal behavior, which can be converted into recommendations or visualizations (e.g., score rankings, or prioritization of needs or benefits). Trace similarity can also examine how close a subject is to achieving ideal insulin levels, and interpret the results to make recommendations (replay reports and / or real-time).

[0086] The score can be compared to a threshold to determine whether there is a risk. If there is a risk, a warning or other indicator can be generated and output, e.g., provided to the subject or a physician. The higher the score, the greater the risk of glycemic dysfunction.

[0087] Comparison of the historical score with the estimated score (and similarly comparison of the risk) shows how much improvement (risk reduction) can be achieved with optimal, or at least better, bolus administration.

[0088] Area 410 shows the subject's historical hyperglycemic risk profile (above the x-axis) and hypoglycemic risk profile (below the x-axis) and identifies times when the subject is consistently experiencing glycemic dysfunction. The peak in the hyperglycemic risk profile identified as 410 indicates that the subject is consistently experiencing out-of-range BG between 1500 and 1600 hours.

[0089] Area 420 shows the hyperglycemia and hyperglycemia risk profile calculated from a replay of the subject's data, with past insulin boluses replaced by numerically optimal bolus amounts. The fact that the peak in the past hyperglycemia risk profile between 15:00 and 16:00 (area 410) is no longer present in area 420 indicates that the subject can significantly reduce their exposure to afternoon hyperglycemia by injecting / infusing a different amount of insulin without taking the additional step of selecting a different bolus time. Similarly, hyperglycemia risk reduction is possible (i) by injecting / infusing a different amount before 22:00 h (by administering a more accurate bolus at dinner time) or (ii) by injecting / infusing a different amount between 23:00 h and 02:00 h (by administering a bedtime snack bolus).

[0090] Region 430 shows the hyperglycemia and hyperglycemia risk profile calculated from a replay of the patient's data, where past insulin boluses have been replaced with the numerically optimal bolus amount for the estimated mealtime. The fact that hyperglycemia risk is further reduced in the 1500-1600 time frame (compared to region 420) indicates that the subject may be able to reduce their exposure to afternoon hyperglycemia by injecting / infusing insulin closer to meals, although the improvement is modest based on this data for this subject. More consistent optimal bolus administration closer to meals / snacks between 2300 and 0600 hours may result in a more significant improvement in hyperglycemia risk.

[0091] FIG. 5 is a diagram 500 of risk zones 510, 520, 530, and 540, according to an embodiment. As shown, the risk profile plot 400 can be converted into daily zones, which in this example indicate missed medications at lunchtime and dinnertime. The risk zones 510, 520, 530, and 540 indicate intervals of significant past hyperglycemic risk that can be completely eliminated by administering boluses at different times or in different amounts. The diagram may show the effect of time of day across multiple days. In an implementation, the diagram 500 may be generated by the BG risk comparator 240 and output by the output device 250 based on the results of the BG risk profile comparison function 334. Past hyperglycemic (above the x-axis) and hypoglycemic (below the x-axis) zones may be identified through the application of thresholds applied to the past BG risk profile generated by the profiling function 322.

[0092] FIG. 6 is a diagram 600 of zones of potential significant improvement (e.g., with respect to risk profile plot 400) by changing bolus dosage or bolus timing, according to an embodiment. Interval 610 indicates an interval in which BG results can be improved through changes in bolus amount (in the past). Interval 620 indicates an interval in which BG results can be improved through bolus administration closer to mealtimes (e.g., earlier in the day). Similar to the generation and output of diagram 500, in an implementation, diagram 600 can be generated by BG risk comparator 240 and output by output device 250 based on the results of BG risk profile comparison function 334. For numerically optimal boluses at past bolus times and numerically optimal boluses at verified estimated mealtimes, historical hyperglycemic (above the x-axis) and hypoglycemic (below the x-axis) zones can be identified by application of thresholds applied to the historical BG risk profiles generated by profiling functions 326 and 330, respectively.

[0093] 7A and 7B are diagrams of various residual insulin (IOB) profile plots 700, 750, according to examples. Figure 7A shows (i) the mean IOB of past bolus doses (line 740) + / - 1 standard deviation (area 720), and (ii) the mean IOB of the numerically optimal bolus dose at past bolus times (line 730) + / - 1 standard deviation (area 710). While the examples herein use the mean and standard deviation, other statistics such as the median, confidence interval, standard error, coefficient of variation, etc., can be used depending on the implementation.

[0094] Similarly, Figure 7B shows (i) the mean IOB of past bolus doses (line 790) + / - 1 standard deviation (area 770), and (ii) the mean IOB of the numerically optimal bolus dose at the estimated bolus time (line 780) + / - 1 standard deviation (area 760).

[0095] Note that in both Figures 7A and 7B, the difference between lines 730, 780 and lines 740, 790 corresponds to different amounts of insulin recently infused / injected. For example, the gap between lines 730 and 740 in Figure 7A between 1200 and 1600 hours corresponds to the fact that a numerically optimal bolus (at a past bolus time) delivers (and "remains") more insulin than previous boluses administered by the patient. A similar gap can be seen between lines 780 and 790 in Figure 7B, even at times with zero or very low IOB in the past, corresponding to the difference between the optimal bolus at the estimated mealtime and the optimal bolus at a past bolus time.

[0096] Area 710 shows the variability (+ / - 1 standard deviation) of IOB about the mean (line 730) associated with the numerically optimal bolus delivery at the past bolus time.

[0097] Area 720 shows the variability (+ / - 1 standard deviation) of IOB about the mean (line 740) associated with the numerically optimal bolus delivery at the estimated bolus time.

[0098] Line 730 shows the average IOB associated with the numerically optimal bolus delivery at the past bolus time, and line 740 shows the average IOB associated with the numerically optimal bolus delivery at the estimated mealtime.

[0099] FIG. 8 is an illustration of an IOB discrepancy plot 800 showing the average difference between (i) the IOB associated with the numerically optimal bolus at the past bolus time and (ii) the IOB associated with past boluses. Line 810 corresponds to the mean IOB difference, with lines 820 and 830 corresponding to + / - 1 standard deviation, respectively. Positive mean values ​​correspond to times when the numerically optimal bolus delivery results in a higher IOB (on average) compared to past boluses. Smaller standard deviation values ​​correspond to more consistent discrepancies between the numerically optimal and past IOB-related boluses. Positive values ​​on plot 800 typically indicate more insulin in the optimal dose than in past doses. Thus, an error profile of bolus delivery over time is provided.

[0100] Line 810 shows the average difference between the IOB associated with the numerically optimal bolus delivery at the past bolus time and the past bolus itself.

[0101] Line 820 shows the mean difference plus one standard deviation, providing a visual indication of the degree of uncertainty about the difference between the means. Note that in some implementations, a confidence interval may be used.

[0102] Line 830 shows the mean difference minus one standard deviation, providing a visual indication of the degree of uncertainty about the difference between the means. Note that in some implementations, a confidence interval may be used.

[0103] 9 is an illustration of an IOB mismatch zone plot 900 for illustrating times when the IOB associated with a numerically optimal bolus delivery is substantially higher than the IOB associated with a past bolus. The bars above the x-axis correspond to times at past bolus times when the subject should consider taking a larger bolus. Although not shown in this example, the bars below the x-axis correspond to times at past bolus times when the subject should consider taking a smaller bolus.

[0104] FIG. 10 is an illustration of an IOB discrepancy plot 1000 for optimal bolus delivery at an estimated mealtime, according to an embodiment. Similar to FIG. 8, this plot shows the mean difference between the IOB associated with the numerically optimal bolus and previous boluses, here along with an envelope of + / - standard deviation for the difference between the numerically optimal bolus and previous boluses at the estimated mealtime. Positive (negative) mean values ​​correspond to times when the numerically optimal bolus delivery at the estimated mealtime resulted in, on average, more (less, respectively) active insulin than previous boluses. Here, the negative discrepancy occurring at 0900 corresponds to the fact that the subject is taking a bolus at a different time than the estimated meal. Thus, an error profile for optimal bolus delivery at the estimated mealtime is provided.

[0105] 11 is a diagram of IOB mismatch zones 1100 for optimal bolus delivery at estimated mealtimes, according to an embodiment. The zone above the x-axis is a time when a subject may benefit from more insulin by delivering a bolus larger than in the past, or by delivering a bolus not delivered previously. Here, the zone below the x-axis at 0900 corresponds to times when a subject should not bolus; instead, a bolus closer to actual mealtimes is desirable.

[0106] FIG. 12 is a diagram of an exemplary summary report 1200 showing glycemic risk based on the determined glycemic function. Summary report 1200 is merely one example of a report that may be generated and provided for a subject or caregiver and may be generated by overlaying the plots and charts of FIGS. 5, 6, 9, and 11. Panel 1210 shows time segments of the profile with glycemic risk. Panel 1220 identifies time segments that would have had significantly lower risk in replay with a change in bolus strategy. Panel 1230 identifies time segments (e.g., meals) where a new bolus amount (more insulin in this example) would act to significantly change IOB. Panel 1240 identifies time segments (e.g., meals) where different bolus timing would act to significantly change IOB. This allows prioritization of actions and therapy changes, identifying periods of the day where a larger bolus could reduce risk. Suboptimal times to boluses may also be identified in panel 1240, and better alternatives may be suggested.

[0107] 13 is a diagram of various risk profile plots 1300 according to another embodiment. In this example, for the subject, the total risk CGM is 1.75, the total risk with optimal bolus at past bolus time is 1.21, the total risk with optimal bolus at estimated mealtime is 0.55, the IOB trace similarity score for optimal bolus at past bolus time is 0.60, and the IOB trace similarity score for optimal bolus at estimated mealtime is 0.47. The sample duration was 79 days.

[0108] Region 1310 indicates a strong tendency to experience hypoglycemia at 1200 hours in past records, with a significant risk of hypoglycemia more generally from 1000 hours to 1600 hours, and a moderate risk of hypoglycemia at 2000 hours.

[0109] Region 1320 shows that the risk of moderate hypoglycemia around 2000 hours can be addressed by administering a different bolus at a previous bolus time, but if the bolus time is not changed, the reduction in evening hypoglycemia comes at the cost of a substantial new risk of hyperglycemia around 2200 hours.

[0110] Region 1330 shows that the risk of severe hypoglycemia around 1200 hours can be addressed by administering a different bolus at earlier bolus times, but if the bolus time is not changed, the reduction in daytime hypoglycemia comes at the cost of a substantial new risk of hyperglycemia between 1000 and 1800 hours.

[0111] Region 1340 shows that administering a numerically optimal bolus at the estimated mealtime effectively avoids the subject being exposed to both hypoglycemia and hyperglycemia.

[0112] 14A and 14B are diagrams of various IOB profile plots 1400, 1450, respectively, according to another embodiment. Line 1410 shows the larger mean IOB (with greater variability) at 2100 hours obtained from past bolus administrations compared to the smaller mean IOB (with less variability) resulting from numerically optimal boluses at past bolus times.

[0113] Line 1420 shows high variability in the mean IOB resulting from a numerically optimal bolus delivery at a past bolus time around 0730 hours. Here, the mean IOB from the numerically optimal bolus delivery is itself very small because there were not many past boluses at that time. Note in Figure 14B that at the same time, the mean IOB associated with the bolus at the estimated mealtime is greater than the mean IOB associated with the past bolus, indicating that the subject may be consistently delaying the bolus for his first meal of the day.

[0114] Line 1430 shows the larger mean IOB (with greater variability) at 1700 hours obtained from past bolus administrations compared to the smaller mean IOB (with less variability) resulting from numerically optimal boluses at past bolus times.

[0115] Line 1440 shows that the IOB at 0830 with the numerically optimal bolus is lower than the average IOB associated with previous bolus records. This is also true for Figure 14B at the same time. Thus, regardless of the timing of the morning bolus, the subject is taking too much insulin, which is associated with a stronger tendency to experience morning hypoglycemia.

[0116] 15A and 15B are diagrams of various IOB discrepancy plots 1500, 1550, respectively, according to another embodiment. In FIG. 15A, the fact that the discrepancy (shown by solid line 1510) is negative throughout the day (e.g., from about 07:30 to about 23:00) indicates that the subject should reduce their bolus if constrained to bolus administration at past bolus times. In FIG. 15B, the fact that the mean discrepancy (shown by solid line 1560) fluctuates between positive and negative values ​​indicates that the subject has a substantially different IOB throughout the day (and that this is beneficial).

[0117] The solid line 1510 shows the average difference between the IOB associated with the numerically optimal bolus delivery at the past bolus time and the past bolus itself.

[0118] Dashed line 1520 indicates the mean difference above plus one standard deviation, providing a visual indication of the degree of uncertainty about the difference between the means.

[0119] Dashed line 1530 indicates the mean difference above plus one standard deviation, providing a visual indication of the degree of uncertainty about the difference between the means.

[0120] 16 is an operational flow of an implementation of a method 1600 for determining the amount of glycemic dysfunction. At 1610, CGM data and insulin data associated with a subject (such as subject 140) are received at a replay analyzer, such as replay analyzer 220. The CGM data and insulin data may include subject-provided data 207, data 210, and / or function 215. The CGM data and insulin data may be received from a storage device, such as a database, containing the subject's 140 CGM and insulin records. Alternatively or additionally, the data may be received from the subject 140 (e.g., via the subject computing device 205), from a CGM and / or insulin device or monitor (e.g., insulin device 110, glucose monitor 120, activity monitor 150, and / or smartphone 160), or from other equipment or entities. In implementations used for replay analysis, historical data can be used. Accurate estimation of both the timing and amount of meals in historical data is useful.

[0121] At 1620, a replay analysis is performed on the CGM data and insulin data, for example, by replay analyzer 220. This analysis step includes (i) the use of a meal estimation function (e.g., meal compiler 303) that generates estimates (e.g., estimated meal data 306) of meal events (particularly carbohydrate intake) in past records based on an analysis of trends and patterns in the time series data of the CGM data and insulin data, and (ii) the use of a replay prediction function 305 (e.g., the output of replay compiler 304) that uses the estimated meal time series to simulate past data starting from any point in the past time series from any particular metabolic initial state, allowing for the specification of either a past or optimized insulin bolus.

[0122] In 1630, the amount of glycemic dysfunction is quantified based on replay analysis, as further described herein, using, for example, one or more of IOB assessor 225, BG risk profiler 230, IOB profiler 235, BG risk comparator, and discrepancy analyzer 245. The analysis includes (i) calculating a hyperglycemia and hypoglycemia risk profile (such as the profiling performed in 322) based on the entirety of the historical CGM data, where the profile value at any point in a 24-hour day describes the clinical significance of past exposure to hyperglycemia (or hypoglycemia) at that time in the day, and where the assignment of hyperglycemia (respectively, hypoglycemia) risk thresholds is normalized to mark intervals of significant glycemic risk within the day, and (ii) simulating the entire record using a replay prediction function (as described with respect to element 312), replacing past boluses with numerically optimal boluses based on the predicted responses from the replay function. (iii) simulating (replaying) the entire recording canceling past boluses and inserting numerically optimal boluses at estimated meal times using a replay prediction function (as described with respect to element 314); (iv) calculating hyperglycemia and hypoglycemia risk profiles for the replayed time series using the numerically optimal boluses at the times of past boluses (as described with respect to element 326); and (v) calculating hyperglycemia and hypoglycemia risk profiles for the replayed time series using the numerically optimal boluses at the times of estimated meals (as described with respect to element 330).

[0123] The BG risk profiler 230 includes calculations for converting one or more blood glucose values ​​into a relative risk score consistent with the risk of diabetic complications or the degree of glycemic control. Risk scores can be assigned a single value or summarized over a time segment, such as a daily average risk score. In one example, the metric is a positive value, with lower values ​​indicating better glycemic control, such as staying within the target zone or closer to the target value. For a target blood glucose range of 95-195 mg / dL, an exemplary risk function maps a reading at 95 mg / dL to a risk of 1.0, a reading at 120 mg / dL to a risk of 0.2, and a reading at 195 mg / dL to a risk of 1.0. As a result, subjects who stay within the target range have a total risk less than 1, while subjects who are tightly controlled have a risk close to zero.

[0124] Continuing with this example, blood glucose values ​​in the hyperglycemic and hypoglycemic ranges are given increasingly higher risk scores as glucose deviates from the target value. The risk function maps a reading at 65 mg / dL to a risk of 2.5 and a reading at 225 mg / dL to a risk of 2.0. Thus, a score of 2 or higher indicates a larger or persistent excursion outside the target range. In this case, the penalty function is not symmetrical because low blood glucose values ​​are heavily penalized. The shape of the risk function can be based on mathematical functions such as a quadratic or lognormal curve, or a set of linear functions. One known metric is the average daily risk range (see, for example, Kovatchev BP et al., "Evaluation of a New Measure of Blood Glucose Variability in Diabetes," Diabetes Care, November 2006, 29(11):2433-8, https: / / www.ncbi.nlm.nih.gov / pubmed / 17065680). In other cases, it can be empirically derived from historical data of the subject or further tailored for the therapeutic goal.

[0125] The risk function can include and emphasize other factors related to diabetes risk. One example is a risk metric that explicitly considers the rate of change of glucose as a patient threat factor (e.g., hypoglycemia and hyperglycemia risk levels are amplified in decreasing and increasing glucose trends, respectively) (see, e.g., U.S. Patent No. 9,439,602, entitled "Alert System for Hypo and Hyperglycemia Prevention Based on Clinical Risk").

[0126] The three sets of risk profiles are evaluated with a BG risk profile comparison function (e.g., BG risk profile comparison function 334), which calculates a total blood glucose score for each of the past data, the replayed data with the numerically optimal bolus at the past bolus time, and the replayed data with the numerically optimal bolus at the estimated meal time. Differences between these respective scores indicate impairments associated with suboptimal boluses at past bolus times and / or boluses that do not occur at the estimated meal times (corresponding to the replay bolus timing effect).

[0127] Additionally, to calculate a measure of pre-meal bolus compliance, time series of IOB are calculated for (i) the past time series of insulin delivery (as described with respect to element 316), (ii) the replayed time series of the numerically optimal bolus at past bolus times (as described with respect to element 318), and (iii) the replayed time series of the numerically optimal bolus at the estimated meal time (as described with respect to element 320). From these IOB time series, a 24-hour profile of IOB can be calculated, where the value of the profile at any time corresponds to a first-order statistic about the amount of insulin (e.g., mean insulin) remaining for the patient at that time in the past, or at the numerically optimal bolus at the past bolus time, or at the numerically optimal bolus at the estimated meal time. These respective IOB profiles allow for assessment of IOB discrepancies between (i) past bolus administrations and numerically optimal boluses at past bolus times (as described with respect to element 336) and (ii) past bolus administrations and numerically optimal boluses at estimated meal times (as described with respect to element 340). Additionally, the three sets of IOB time series allow for the calculation of numerical similarity measures between (i) past bolus administrations and numerically optimal boluses at past bolus times (as described with respect to element 338) and (ii) past bolus administrations and numerically optimal boluses at estimated meal times (as described with respect to element 342).

[0128] In 1640, an output representing the amount of glycemic dysfunction is provided via an output device, such as output device 250, to, for example, a computing device associated with the subject, a physician, or another entity. The output may include plots or other information and recommendations, as further described herein (e.g., a summary report, such as report 1200 of FIG. 12). Such output may provide indicators of replay bolus timing effectiveness and pre-meal bolus compliance. In one or more implementations, the impact of meal misinformation (both in terms of timing and configuration) is quantified. Thus, insulin benefit (e.g., risk reduction) is assessed with respect to optimal doses and with respect to changes in optimal pre-meal timing.

[0129] 17 is an operational flow of another implementation of a method 1700 for determining the amount of glycemic dysfunction. At 1710, input is received at a computing device, such as computing device 201. In an implementation, the input may include subject-provided data 207, data 210, and / or functions 215 and may be received at a replay analyzer 220. The input may include CGM data and insulin data similar to that described with respect to 1610 (e.g., past BG 301 and past insulin 302).

[0130] At 1720, a replay analysis is performed (e.g., by replay analyzer 220) on past boluses and estimated meals, for example, using the embodiments described with respect to 312 and 314. In this way, it can be determined what would have happened if the subject had taken a bolus at past bolus times, and what would have happened if the subject had taken a bolus at estimated meal times.

[0131] At 1730, IOB data analysis is performed, such as described with respect to Figures 3A and 3B. In an implementation, IOB data analysis may be performed by IOB assessor 225 and / or IOB profiler 235. At 1740, BG data analysis is performed, such as described with respect to Figures 3A and 3B. In an implementation, BG data analysis may be performed by BG risk profiler 230 and / or BG risk comparator 240. Such analysis may include various assessments and profiling, as further described above.

[0132] At 1750, discrepancies are determined along with the comparison, as described in detail above with respect to Figures 3A and 3B. In implementations, discrepancies may be determined by discrepancy analyzer 245. At 1760, plots and scores as further described herein may be generated and output, for example, via output device 250.

[0133] 18 is an operational flow of another implementation of a method 1800 for determining the amount of glycemic dysfunction. The method may be performed by a computing device 201 in some implementations.

[0134] At 1810, past BG data and past insulin data are received at a replay analyzer, such as replay analyzer 220. Multiple replay analyses are performed at 1820 by replay analyzer 220 on past boluses and on meals estimated using the past BG data and past insulin data.

[0135] At 1830, at least one IOB data analysis is performed using at least one of a plurality of replay analyses. The IOB data analysis may be performed using, for example, one or more of IOB assessor 225 or IOB profiler 235. At 1840, at least one BG data analysis is performed using at least one of a plurality of replay analyses. The BG data analysis may be performed using, for example, one or more of BG risk profiler 230 or BG risk comparator 240.

[0136] At least one discrepancy analysis is performed and at least one comparison is performed using at least one IOB data analysis and at least one BG data analysis in 1850. The discrepancy analysis and / or comparison may be performed by, for example, discrepancy analyzer 245.

[0137] At 1860, an output is provided as further described herein. The output may be provided by output device 250 and may be based on at least one discrepancy analysis and at least one comparison.

[0138] 19 is an operational flow of another implementation of a method 1900 for determining the amount of glycemic dysfunction. The method may be performed by a computing device 201 in some implementations.

[0139] In 1910, at least one of the subject's historical BG data or historical insulin data is received in a replay analyzer, such as replay analyzer 220. In 1920, a first risk score is determined for the subject's glycemic dysfunction based on at least one of the historical BG data or historical insulin data. The risk score may be determined, for example, by replay analyzer 220. In one example, the risk score may be a value from BG risk profiler 230. In another example, the output from BG risk profiler 230 is converted into a risk score having a set of thresholds or zones, such that a risk of less than 1 is assigned as a low risk score, a risk of 1-2 is assigned as a moderate risk, and a risk of greater than 2 is assigned as a high risk. These thresholds and zones may be preset or customized for each patient based on historical values ​​or current goals. The risk score may also be used to summarize BG risk profiles associated with different times of day, meals, and / or other relevant events.

[0140] At 1930, at least one of the subject's optimal BG data or optimal insulin data is received in a replay analyzer, such as replay analyzer 220. At 1940, a second risk score is determined for the subject's glycemic dysfunction based on at least one of the optimal BG data or optimal insulin data. The risk score may be determined, for example, by replay analyzer 220. This risk score may use the same attributes as the risk score calculated in 1910 so that potential improvement can be assessed. For example, a risk score for breakfast is compared between historical data and optimal data.

[0141] In 1950, the first risk score is compared to the second risk score, e.g., by the replay analyzer 220. A recommendation may be generated and output in 1960 based on the comparison in 1940, e.g., in some implementations, by the replay analyzer 220 and / or the output device 250. The recommendation may be optimal BG data or optimal insulin data for the subject or other user, regarding instructions or other suggestions for improving or maintaining blood glucose information (e.g., in some implementations, such as that shown in FIG. 12 , for reducing glycemic dysfunction). If the first risk score is close to the second risk score (e.g., less than a predetermined percentage difference, within a score range or difference, less than a threshold difference, etc.), the subject's behavior is identified and encouraged as the insulin amount and timing are near ideal. If the first risk score is significantly greater than the second risk score (e.g., by more than a predetermined percentage difference, outside a score range or difference, more than a threshold difference, etc.), there is an opportunity to optimize therapy by adjusting the amount or timing of insulin boluses for future similar meals.

[0142] The output recommendations in 1960 can be generated in a variety of formats to suit the subject's needs. In one embodiment, the recommendations can be qualitative, for example, identifying that dosing early before breakfast will help avoid high glucose levels later in the day. For example, in some implementations, the recommendations can be an on-screen prompt or overlaid on a plot showing glucose excursions. Alternatively or additionally, the recommendations can be a summary report spanning a period of time, such as a week or a month. The summary report can identify periods of high blood glucose levels throughout the day and assign possible causes. For example, the report might identify low afternoon glucose levels as the result of taking too much insulin at lunch combined with taking the medication after eating a meal.

[0143] FIG. 20 is an operational flow of another implementation of a method 2000 for determining the amount of glycemic dysfunction. In some implementations, the method may be executed by a computing device 201. Operation 2010 involves determining a past risk score based on at least one of past blood glucose data or past insulin data, which, in the implementation of FIGS. 3A and 3B, involves elements 322 (for calculating a risk profile) and 334 (for calculating a risk score). Operation 2020 involves calculating an optimal risk score for glycemic dysfunction based on at least one of optimal BG data or optimal insulin, which, in the implementation of FIGS. 3A and 3B, involves elements 312 (for a numerically optimal bolus delivery at a past bolus time) or 314 (for a numerically optimal bolus delivery at an estimated mealtime), 326 or 330 (for an optimal risk profile, respectively), and 334 (for calculating a risk score, respectively). Operation 2030 involves determining the difference between the historical risk score and the optimal risk score, which, for the implementation of Figures 3A and 3B, involves element 334. Operation 2040 provides a recommendation based on the difference and based on a threshold.

[0144] At 2010, a historical risk score for glycemic dysfunction is determined based on at least one of the subject's historical BG data or historical insulin data, which may have been received, for example, at replay analyzer 220.

[0145] At 2020, an optimal risk score for glycemic dysfunction is determined based on at least one of the subject's optimal BG data or optimal insulin data, which may have been received, for example, at replay analyzer 220.

[0146] At 2030, the difference between the historical risk score and the optimal risk score is determined, for example, by replay analyzer 220.

[0147] At 2040, a recommendation may be generated and output based on the determined difference, for example, in some implementations, by the replay analyzer 220 and / or the output device 250.

[0148] 21 illustrates an exemplary computing environment in which exemplary embodiments and aspects may be implemented. The computing device environment is only one example of a suitable computing environment and is not intended to suggest any limitation as to the scope of use or functionality.

[0149] Many other general-purpose or special-purpose computing device environments or configurations can be used. Examples of well-known computing devices, environments, and / or configurations that may be suitable for use include, but are not limited to, personal computers, server computers, handheld or laptop devices, multiprocessor systems, microprocessor-based systems, networked personal computers (PCs), minicomputers, mainframe computers, embedded systems, distributed computing environments that include any of the above systems or devices, etc.

[0150] Computer-executable instructions, such as program modules, may be used that are executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Distributed computing environments may be used in which tasks are performed by remote processing devices that are linked through a communications network or other data transmission medium. In a distributed computing environment, program modules and other data may be located in both local and remote computer storage media, including memory storage devices.

[0151] 21 , an exemplary system for implementing aspects described herein includes a computing device, such as computing device 2100. In its most basic configuration, computing device 2100 typically includes at least one processing unit 2102 and memory 2104. Depending on the exact configuration and type of computing device, memory 2104 may be volatile (such as random access memory (RAM)), non-volatile (such as read-only memory (ROM), flash memory), or some combination of the two. This most basic configuration is indicated in FIG. 21 by dashed line 2106.

[0152] Computing device 2100 may have additional features / functionality. For example, computing device 2100 may include additional storage (removable and / or non-removable) including, but not limited to, magnetic or optical disks or tape. Such additional storage is illustrated in FIG. 21 by removable storage 2108 and non-removable storage 2110.

[0153] Computing device 2100 typically includes a variety of computer-readable media, which can be any available media that can be accessed by device 2100 and includes both volatile and nonvolatile media, removable and non-removable media.

[0154] Computer storage media include volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information, such as computer-readable instructions, data structures, program modules, and other data. Memory 2104, removable storage 2108, and non-removable storage 2110 are all examples of computer storage media. Computer storage media include, but are not limited to, RAM, ROM, electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, CD-ROM, digital versatile disks (DVDs) or other optical storage devices, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to store the desired information and that can be accessed by computing device 2100. Any such computer storage media may be part of computing device 2100.

[0155] Computing device 2100 may contain communications connections 2112 that allow the device to communicate with other devices. Computing device 2100 may also have input devices 2114, such as a keyboard, mouse, pen, voice input device, touch input device, etc. Output devices 2116, such as a display, speakers, printer, etc. All of these devices are well known in the art and need not be described at length here.

[0156] In an implementation, a method for determining an amount of glycemic dysfunction is provided, the method including receiving continuous glucose monitoring (CGM) and insulin data for a subject, the CGM and insulin data including insulin bolus amount and timing of the insulin bolus amount, performing a replay analysis of the CGM and insulin data, quantifying the amount of glycemic dysfunction using the replay analysis, and providing an output representative of the amount of glycemic dysfunction.

[0157] Implementations may include some or all of the following features: A first portion of CGM and insulin data is estimated and a second portion of CGM and insulin data is reported; The CGM and insulin data includes at least one of estimated timing of the subject's meals or estimated composition of the subject's meals; The CGM and insulin data includes misinformation regarding at least one of the timing of the subject's meals or composition of the subject's meals; A portion of the CGM and insulin data is received from the subject; The CGM and insulin data includes an estimated metabolic state in time series form, an adjusted meal history, and a delivered insulin history; The CGM and insulin data includes an estimated metabolic state, adjusted estimated metabolic inputs, and known metabolic inputs; The CGM and insulin data is discretized in time.

[0158] Implementations may also include some or all of the following features: Performing replay analysis on the CGM and insulin data includes isolating the effects of estimated bolus timing, carbohydrate counting, and carbohydrate ratios relative to meals; Performing replay analysis on the CGM and insulin data includes assessing the impact of numerically optimal boluses at past bolus times by eliminating inaccurate carbohydrate counting and improper carbohydrate ratios using replay simulation analysis; Performing replay analysis on the CGM and insulin data includes assessing the impact of numerically optimal boluses at estimated past meal times; Performing replay analysis on the CGM and insulin data includes running a replay simulation using at least one of past CGM and insulin data, numerically optimal boluses at past bolus times, or numerically optimal boluses at estimated meal times. Performing replay analysis on the CGM and insulin data includes generating a set of candidate boluses, simulating future blood glucose (BG) values ​​associated with each candidate bolus, scoring the simulated BG trajectory for each candidate bolus using risk analysis, selecting and implementing the best bolus based on the score for each candidate bolus, and continuing to replay the simulation until the next bolus time. The risk analysis is a rate-dependent risk analysis.

[0159] Implementations may also include some or all of the following features: Quantifying the amount of glycemic dysfunction includes quantifying the amount of dysfunction of an estimated bolus dose compared to an optimally timed bolus dose; Quantifying the amount of glycemic dysfunction includes quantifying the amount of dysfunction of a past bolus dose compared to an optimally timed bolus dose; Quantifying the amount of glycemic dysfunction includes quantifying patient compliance with ideal pre-meal bolus timing.

[0160] Implementations may also include some or all of the following features: The output includes at least one of a plot or visualization; The output includes a risk index including at least one of high blood glucose risk, low blood glucose risk, or total glycemic risk; The output includes the number and severity of excursions to out-of-range blood glucose; The output indicates glycemic dysfunction over a past time; Providing the output includes providing a visualization showing at least one of the behavioral impact of glycemic dysfunction, past CGM and insulin data versus a simulated CGM replayed with an optimal bolus at a past bolus time, or past CGM and insulin data versus a simulated CGM replayed with an optimal bolus at an estimated meal time; The output includes a visualization of patient compliance with ideal pre-meal bolus timing; Meal management information is provided to the subject based on the amount of glycemic dysfunction; The meal management information includes the patient's pre-bolus timing information; The pre-bolus timing information includes a recommendation to the patient regarding a bolus administration before a certain type of future meal.

[0161] In an implementation, a method for determining an amount of glycemic dysfunction is provided, the system including a processor that receives continuous glucose monitoring (CGM) and insulin data for a subject, the CGM and insulin data including insulin bolus amounts and timing of the insulin bolus amounts, a replay analyzer that generates a replay analysis using the CGM and insulin data, a quantifier that uses the replay analysis to quantify the amount of glycemic dysfunction, and an output device that provides an output representative of the amount of glycemic dysfunction.

[0162] Implementations may include some or all of the following features: A first portion of CGM and insulin data is estimated and a second portion of CGM and insulin data is reported; The CGM and insulin data includes at least one of estimated timing of the subject's meals or estimated composition of the subject's meals; The CGM and insulin data includes misinformation regarding at least one of the timing of the subject's meals or composition of the subject's meals; A portion of the CGM and insulin data is received from the subject's computing device; The CGM and insulin data includes an estimated metabolic state in time series form, an adjusted meal history, and a delivered insulin history; The CGM and insulin data includes an estimated metabolic state, adjusted estimated metabolic inputs, and known metabolic inputs; The CGM and insulin data is discretized in time.

[0163] Implementations may also include some or all of the following features: A replay analyzer separates the effects of estimated bolus timing, carbohydrate counting, and carbohydrate ratios relative to meals. The replay analyzer assesses the impact of numerically optimal boluses at past bolus times by eliminating inaccurate carbohydrate counts and improper carbohydrate ratios using replay simulation analysis. The replay analyzer assesses the impact of numerically optimal boluses at estimated past meal times. The replay analyzer performs a replay simulation using at least one of past CGM and insulin data, numerically optimal boluses at past bolus times, or numerically optimal boluses at estimated meal times. The replay analyzer is configured to generate a set of candidate boluses, simulate future blood glucose (BG) values ​​associated with each candidate bolus, score the simulated BG trajectory for each candidate bolus using risk analysis, such as rate-dependent risk analysis, select and implement the best bolus based on the score of each candidate bolus, and continue replaying the simulation until the next bolus time.

[0164] Implementations may also include some or all of the following features: A quantifier quantifies the amount of impairment of an estimated bolus dose compared to an optimally timed bolus dose; A quantifier quantifies the amount of impairment of a past bolus dose compared to an optimally timed bolus dose; A quantifier quantifies patient compliance with ideal pre-meal bolus timing.

[0165] Implementations may also include some or all of the following features: The output includes at least one of a plot or visualization; The output includes a risk index including at least one of hyperglycemic risk, hypoglycemic risk, or total glycemic risk; The output includes the number and severity of excursions to out-of-range blood glucose; The output indicates glycemic dysfunction over a past time; The output device provides a visualization showing at least one of the behavioral impact of glycemic dysfunction, past CGM and insulin data versus a simulated CGM replayed with an optimal bolus at a past bolus time, or past CGM and insulin data versus a simulated CGM replayed with an optimal bolus at an estimated meal time; A compliance engine configured to provide a visualization of a patient's compliance with ideal pre-meal bolus timing is provided; The output device provides meal management information to the subject based on the amount of glycemic dysfunction; The meal management information includes the patient's pre-bolus timing information; The pre-bolus timing information includes a recommendation to the patient regarding a bolus administration before a certain type of future meal.

[0166] In an implementation, a method for determining an amount of glycemic dysfunction is provided, the method including receiving historical blood glucose (BG) data and historical insulin data with a replay analyzer, performing a plurality of replay analyses of past boluses and estimated meals using the historical BG data and the historical insulin data, performing at least one residual insulin (IOB) data analysis using at least one of the plurality of replay analyses, performing at least one BG data analysis using at least one of the plurality of replay analyses, performing at least one discrepancy analysis and at least one comparison using the at least one IOB data analysis and the at least one BG data analysis, and providing an output based on the at least one discrepancy analysis and the at least one comparison.

[0167] Implementations may include some or all of the following features: The BG data and insulin data include an insulin bolus amount and timing of the insulin bolus amount; A first portion of the BG data and insulin data is estimated, and a second portion of the BG data and insulin data is reported; The BG data and insulin data include at least one of estimated timing of the subject's meals or estimated composition of the subject's meals; The BG data and insulin data include misinformation regarding at least one of the timing of the subject's meals or composition of the subject's meals; A portion of the BG data and insulin data is received from the subject; The BG data and insulin data includes an estimated metabolic state, an adjusted meal history, and a delivered insulin history in time series form; The BG data and insulin data includes an estimated metabolic state, adjusted estimated metabolic inputs, and known metabolic inputs; The BG data and insulin data are discretized in time.

[0168] Implementations may also include some or all of the following features: Performing replay analysis includes isolating the effects of estimated bolus timing, carbohydrate count, and carbohydrate ratio relative to a meal; Performing replay analysis includes assessing the impact of a numerically optimal bolus at a past bolus time by eliminating inaccurate carbohydrate counts and improper carbohydrate ratios using replay simulation analysis; Performing replay analysis includes assessing the impact of a numerically optimal bolus at an estimated past meal time; Performing replay analysis includes running a replay simulation at a past bolus time or at a time before the estimated meal; Performing replay analysis includes running a replay simulation using at least one of BG data and insulin data, a numerically optimal bolus at a past bolus time, or a numerically optimal bolus at an estimated meal time. Performing the replay analysis includes generating a set of candidate boluses, simulating future blood glucose (BG) values ​​associated with each candidate bolus, scoring the simulated BG trajectory for each candidate bolus using a risk analysis, selecting and implementing a best bolus based on the score for each candidate bolus, and continuing to replay the simulation until the next bolus time. The risk analysis is a rate-dependent risk analysis. Performing at least one discrepancy analysis includes quantifying an amount of impairment of an estimated bolus administration compared to an optimally timed bolus administration. Performing at least one discrepancy analysis includes quantifying an amount of impairment of a past bolus administration compared to an optimally timed bolus administration. Performing at least one discrepancy analysis includes quantifying patient compliance with ideal pre-meal bolus timing.

[0169] Implementations may also include some or all of the following features: The output includes at least one of a plot or visualization; The output includes a risk index including at least one of high blood glucose risk, low blood glucose risk, or total glycemic risk; The output includes the number and severity of excursions to out-of-range blood glucose; The output indicates glycemic dysfunction over a past time; Providing the output provides a visualization showing at least one of the behavioral impact of glycemic dysfunction, BG data and insulin data versus BG data replayed with an optimal bolus at a past bolus time, or BG data and insulin data versus BG data replayed with an optimal bolus at an estimated meal time; The output includes a visualization of patient compliance with ideal pre-meal bolus timing; Providing meal management information to a subject based on the amount of glycemic dysfunction; The meal management information includes pre-bolus timing information for the patient. The pre-bolus timing information includes a recommendation to the patient regarding a bolus administration before a certain type of future meal.

[0170] In an implementation, a system for determining an amount of glycemic dysfunction is provided, the system including: an input device configured to receive blood glucose (BG) data, insulin data, and a plurality of functions; a replay analyzer configured to receive the BG data, the insulin data, and the plurality of functions and to perform a plurality of replay analyses on past boluses and on estimated meals using the BG data, the insulin data, and the plurality of functions; a residual insulin (IOB) assessor configured to perform a plurality of IOB assessments using at least one of the insulin data or the plurality of replay analyses; a BG risk profiler configured to perform a plurality of BG risk profiles using at least one of the BG data or the plurality of replay analyses; and a replay analyzer configured to perform a plurality of IOB assessments using at least one of the BG data or the plurality of replay analyses. an IOB profiler configured to perform a plurality of IOB profiles using at least one of the IOB risk profiles; a BG risk comparator configured to determine at least one of a risk profile plot or a BG risk score using the output of the BG risk profiler; a discrepancy analyzer configured to perform a plurality of discrepancy analyses using the IOB profiles and the IOB assessment to determine at least one of an IOB profile plot, an IOB discrepancy plot, or an IOB trace similarity score; and an output device configured to output at least one of the risk profile plot, the BG risk score, the IOB profile plot, the IOB discrepancy plot, or the IOB trace similarity score.

[0171] Implementations may include some or all of the following features: The BG data and insulin data include an insulin bolus amount and timing of the insulin bolus amount; A first portion of the BG data and insulin data is estimated, and a second portion of the BG data and insulin data is reported; The BG data and insulin data include at least one of estimated timing of the subject's meals or estimated composition of the subject's meals; The BG data and insulin data include misinformation regarding at least one of the timing of the subject's meals or composition of the subject's meals; A portion of the BG data and insulin data is received from the subject's computing device; The BG data and insulin data include an estimated metabolic state, an adjusted meal history, and a delivered insulin history in time series form; The BG data and insulin data include an estimated metabolic state, adjusted estimated metabolic inputs, and known metabolic inputs; The BG data and insulin data are discretized in time.

[0172] Implementations may also include some or all of the following features: A replay analyzer separates the effects of estimated bolus timing, carbohydrate counting, and carbohydrate ratio relative to a meal. The replay analyzer assesses the impact of a numerically optimal bolus at a past bolus time by eliminating inaccurate carbohydrate counting and improper carbohydrate ratios using replay simulation analysis. The replay analyzer assesses the impact of a numerically optimal bolus at an estimated past meal time. The replay analyzer runs a replay simulation at a past bolus time or at a time before an estimated meal. The replay analyzer runs a replay simulation using at least one of BG and insulin data, a numerically optimal bolus at a past bolus time, or a numerically optimal bolus at an estimated meal time. The replay analyzer is configured to generate a set of candidate boluses, simulate future blood glucose (BG) values ​​associated with each candidate bolus, score the simulated BG trajectory for each candidate bolus using a risk analysis, select and implement a best bolus based on the score for each candidate bolus, and continue replaying the simulation until the next bolus time. The risk analysis is a rate-dependent risk analysis. The discrepancy analyzer quantifies the amount of impairment of an estimated bolus delivery compared to an optimally timed bolus delivery. The discrepancy analyzer quantifies the amount of impairment of a past bolus delivery compared to an optimally timed bolus delivery. The discrepancy analyzer quantifies patient compliance with ideal pre-meal bolus timing.

[0173] Implementations may also include some or all of the following features: The output includes at least one of a plot or visualization; The output includes a risk index including at least one of high blood glucose risk, low blood glucose risk, or total glycemic risk; The output includes the number and severity of excursions to out-of-range blood glucose; The output indicates glycemic dysfunction over a past time; The output device provides a visualization showing at least one of the behavioral impact of glycemic dysfunction, BG data and insulin data versus BG data replayed with an optimal bolus at a past bolus time, or CGM replayed with an optimal bolus at an estimated meal time; The output device is configured to provide a visualization of the patient's compliance with ideal pre-meal bolus timing; The output device provides meal management information to the subject based on the amount of glycemic dysfunction; The meal management information includes the patient's pre-bolus timing information; The pre-bolus timing information includes a recommendation to the patient regarding a bolus administration before a certain type of future meal.

[0174] In one implementation, a method for providing a recommendation based on glycemic risk is provided, the method including receiving at least one of historical blood glucose (BG) data or historical insulin data of a subject, the BG data and insulin data including an insulin bolus amount and a timing of the insulin bolus amount, determining a first risk score of glycemic dysfunction for the subject based on the at least one of the historical BG data or the historical insulin data, receiving at least one of optimal BG data or optimal insulin data for the subject, determining a second risk score of glycemic dysfunction for the subject based on the at least one of the historical BG data or the historical insulin data, comparing the first risk score to the second risk score, and outputting a recommendation based on the comparison of the first risk score and the second risk score.

[0175] Implementations may include some or all of the following features: A first portion of BG data and insulin data is estimated, and a second portion of BG data and insulin data is reported; The BG data and insulin data includes at least one of estimated timing of the subject's meals or estimated composition of the subject's meals; The BG data and insulin data includes misinformation regarding at least one of the timing of the subject's meals or composition of the subject's meals; A portion of the BG data and insulin data is received from the subject; The BG data and insulin data includes an estimated metabolic state in time series form, an adjusted meal history, and a delivered insulin history; The BG data and insulin data includes an estimated metabolic state, adjusted estimated metabolic inputs, and known metabolic inputs; The BG data and insulin data are discretized in time.

[0176] Implementations may also include some or all of the following features: comparing the first risk score to the second risk score includes quantifying an amount of impairment of an estimated bolus dose compared to an optimally timed bolus dose; comparing the first risk score to the second risk score includes quantifying an amount of impairment of a past bolus dose compared to an optimally timed bolus dose; and comparing the first risk score to the second risk score includes quantifying patient compliance with ideal pre-meal bolus timing.

[0177] Implementations may also include some or all of the following features: The recommendation includes at least one of a plot or visualization; The recommendation includes a risk index including at least one of high blood glucose risk, low blood glucose risk, or total glycemic risk; The recommendation includes the number and severity of excursions to out-of-range blood glucose; The output indicates glycemic dysfunction in a past time; Outputting the recommendation includes providing a visualization showing at least one of the behavioral impact of glycemic dysfunction, BG data and insulin data versus BG data replayed with an optimal bolus at a past bolus time, or BG data and insulin data versus BG data replayed with an optimal bolus at an estimated meal time; The recommendation includes visualization of patient compliance with ideal pre-meal bolus timing.

[0178] Implementations may also include some or all of the following features: The method further includes providing dietary management information to the subject based on the amount of glycemic dysfunction; The dietary management information includes pre-bolus timing information for the patient; The pre-bolus timing information includes a recommendation to the patient regarding a bolus administration before a certain type of future meal; The method further includes converting one or more blood glucose values ​​into a relative risk score consistent with the risk of diabetic complications or the degree of glycemic control; Each of the first risk score and the second risk score is assigned a single value or summarized over a time segment, such as a daily average risk score; The method further includes identifying and encouraging a subject action if the first risk score is close to the second risk score because the amount and timing of insulin are near ideal; and identifying and encouraging a subject action if the comparison between the first risk score and the second risk score is at least one of less than a predetermined percentage difference, within a score range or difference, or less than a threshold difference because the amount and timing of insulin are near ideal. The method further includes adjusting the amount or timing of an insulin bolus for a future similar meal if the first risk score is significantly greater than the second risk score, and adjusting the amount or timing of an insulin bolus for a future similar meal if the comparison of the first risk score to the second risk score is at least one of a predetermined percentage difference or greater, a score range or difference outside, or a threshold difference or greater.

[0179] In an implementation, a method for assessing glycemic risk is provided, the method including: determining a historical risk score for glycemic dysfunction based on at least one of historical blood glucose (BG) data or historical insulin data, where the historical BG data and the historical insulin data include an insulin bolus amount and a timing of the insulin bolus amount; determining an optimal risk score for glycemic dysfunction based on at least one of optimal BG data or optimal insulin data; determining a difference between the historical risk score and the optimal risk score; and providing a recommendation based on the difference and based on a threshold value.

[0180] Implementations may include some or all of the following features: A first portion of BG data and insulin data is estimated, and a second portion of BG data and insulin data is reported. The BG data and insulin data include at least one of estimated timing of the subject's meals or estimated composition of the subject's meals. The BG data and insulin data include misinformation regarding at least one of the timing of the subject's meals or composition of the subject's meals. A portion of the BG data and insulin data is received from the subject. The BG data and insulin data include an estimated metabolic state in time series form, an adjusted meal history, and a delivered insulin history. The BG data and insulin data include an estimated metabolic state, adjusted estimated metabolic inputs, and known metabolic inputs. The BG data and insulin data are discretized in time. Determining a difference between the historical risk score and an optimal risk score includes quantifying an amount of impairment of the estimated bolus dose compared to an optimally timed bolus dose. Determining a difference between the historical risk score and an optimal risk score includes quantifying an amount of impairment of the historical bolus dose compared to an optimally timed bolus dose. Determining the difference between the historical risk score and the optimal risk score involves quantifying the patient's compliance with ideal pre-meal bolus timing.

[0181] Implementations may also include some or all of the following features: The recommendation includes at least one of a plot or visualization; The recommendation includes a risk index including at least one of high blood glucose risk, low blood glucose risk, or total glycemic risk; The recommendation includes the number and severity of excursions to out-of-range blood glucose; The output indicates glycemic dysfunction over a past time; Providing the recommendation includes providing a visualization showing at least one of the behavioral impact of glycemic dysfunction, BG data and insulin data versus BG data replayed with an optimal bolus at a past bolus time, or BG data and insulin data versus BG data replayed with an optimal bolus at an estimated meal time; The recommendation includes visualization of patient compliance with ideal pre-meal bolus timing.

[0182] Implementations may also include some or all of the following features: The method further includes providing dietary management information to the subject based on the amount of glycemic dysfunction; The dietary management information includes pre-bolus timing information for the patient; The pre-bolus timing information includes a recommendation to the patient regarding a bolus administration before a certain type of future meal; The method further includes converting one or more blood glucose values ​​into a relative risk score consistent with the risk of diabetic complications or the degree of glycemic control; The historical risk score and the optimal risk score are each assigned a single value or summarized over a time segment, such as a daily average risk score; If the historical risk score is close to the optimal risk score, the method further includes identifying and encouraging the subject's behavior because the insulin amount and timing are near ideal; If the difference between the historical risk score and the optimal risk score is at least one of less than a predetermined percentage difference, within a score range or difference, or less than a threshold difference, the method further includes identifying and encouraging the subject's behavior because the insulin amount and timing are near ideal. The method further includes adjusting the amount or timing of an insulin bolus for a future similar meal if the past risk score is significantly greater than the optimal risk score, and adjusting the amount or timing of an insulin bolus for a future similar meal if the difference between the past risk score and the optimal risk score is at least one of a predetermined percentage difference or greater, a score range or difference outside, or a threshold difference or greater.

[0183] In an implementation, the method includes calculating a first IOB trace based on the actual bolus recording, calculating a second IOB trace of an optimal therapy, comparing the first IOB trace to the second IOB trace to determine an IOB similarity score, and using the IOB similarity score to modify at least one of smart alert behavior, user recommendations, reporting functions, and visualizations.

[0184] Implementations may include some or all of the following features: Comparing the traces includes determining a dot product between the unit magnitude traces to provide a value ranging from 0 to 1, with 1 indicating the actual and optimal dosing shapes are most similar; The method further includes thresholding or comparing the IOB similarity score against current behavior and using the thresholding or comparison to modify smart alert behaviors, user recommendations, reporting functions, and visualizations.

[0185] In an implementation, the system includes a discrepancy analyzer configured to determine an IOB similarity score using a comparison of a first IOB trace based on an actual bolus recording and a second IOB trace of an optimal treatment, and an output device configured to output the IOB similarity score.

[0186] Implementations may include some or all of the following features: The system further comprises an input device configured to calculate a first IOB trace based on the recording of the actual bolus and to calculate a second IOB trace of an optimal therapy.

[0187] It should be understood that the various techniques described herein may be implemented in connection with hardware or software components or, where appropriate, with a combination of both. Illustrative types of hardware components that may be used include field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems-on-chips (SOCs), complex programmable logic devices (CPLDs), etc. The methods and apparatus of the presently disclosed subject matter, or certain aspects or portions thereof, may take the form of program code (i.e., instructions) embodied in a tangible medium, such as a floppy disk, a CD-ROM, a hard drive, or any other machine-readable storage medium such that, when the program code is loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the presently disclosed subject matter.

[0188] While example implementations may refer to utilizing aspects of the presently disclosed subject matter in the context of one or more stand-alone computer systems, the subject matter is not so limited, but rather may be implemented in connection with any computing environment, such as a network or distributed computing environment. Still further, aspects of the presently disclosed subject matter may be implemented within or across multiple processing chips or devices, and storage devices may similarly be implemented across multiple devices. Such devices may include, for example, personal computers, network servers, and handheld devices.

[0189] Although the subject matter has been described in language specific to structural features and / or methodological acts, it is to be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or acts described above. Rather, the specific features and acts described above are disclosed as example forms of implementing the claims.

Claims

1. 1. A system for determining an amount of glycemic dysfunction, comprising: a data processor that receives continuous glucose monitoring (CGM) and insulin data for a subject, the insulin data including an insulin bolus amount and a timing of the insulin bolus amount; a replay analyzer that uses the CGM and insulin data to generate a replay analysis to determine the impact of previously delivered inefficient or suboptimal boluses and assess the impact on blood glucose levels of a numerically optimal bolus at an estimated past mealtime; a quantifier that quantifies the amount of glycemic dysfunction using the replay analysis; and an output device that provides an output representative of the amount of glycemic dysfunction for use in improving diabetes management.

2. The system of claim 1 , wherein a first portion of the CGM and insulin data is estimated and a second portion of the CGM and insulin data is reported.

3. 10. The system of claim 1, wherein the CGM and insulin data includes at least one of estimated timing of the subject's meals or estimated composition of the subject's meals.

4. 10. The system of claim 1, wherein the CGM and insulin data includes misinformation regarding at least one of the timing of the subject's meals or the composition of the subject's meals.

5. The system of claim 1 , wherein a portion of the CGM and insulin data is received from the subject's computing device.

6. The system of claim 1 , wherein the CGM and insulin data includes an estimated metabolic state and delivered insulin history in time series form.

7. The system of claim 1 , wherein the CGM and insulin data includes an estimated metabolic state, an adjusted estimated metabolic input, and a known metabolic input.

8. The system of claim 1 , wherein the CGM and insulin data are discretized in time.

9. 10. The system of claim 1, wherein the replay analyzer separates the effects of estimated bolus timing, carbohydrate count, and carbohydrate ratio relative to the meal.

10. 10. The system of claim 1, wherein the replay analyzer uses replay simulation analysis to assess the impact of a numerically optimal bolus on past bolus times by eliminating inaccurate carbohydrate counts and inappropriate carbohydrate ratios.

11. The system of claim 1 , wherein the replay analyzer performs a replay simulation at a past bolus time or at a time before an estimated meal.

12. 10. The system of claim 1, wherein the replay analyzer performs a replay simulation using at least one of past CGM and insulin data, a numerically optimal bolus at past bolus times, or a numerically optimal bolus at estimated meal times.

13. the replay analyzer: generating a set of candidate boluses; simulating future blood glucose (BG) values ​​associated with each candidate bolus; scoring the simulated BG trajectory for each candidate bolus using risk analysis; selecting and implementing the best bolus based on the score of each candidate bolus; and continuing to replay the simulation until the next bolus time.

14. The system of claim 13 , wherein the risk analysis is a rate-dependent risk analysis.

15. The system of claim 1 , wherein the quantifier quantifies the amount of dysfunction of the estimated bolus administration compared to an optimally timed bolus administration.

16. The system of claim 1 , wherein the quantifier quantifies the amount of dysfunction of past bolus doses compared to optimally timed bolus doses.

17. The system of claim 1 , wherein the quantifier quantifies a degree of patient compliance with ideal pre-meal bolus timing.

18. The system of claim 1 , wherein the output comprises at least one of a plot or a visualization.

19. The system of claim 1 , wherein the output comprises a risk index comprising at least one of high blood glucose risk, low blood glucose risk, or total glycemic risk.

20. The system of claim 1 , wherein the output includes the number and degree of excursions into out-of-range blood glucose levels.

21. The system of claim 1 , wherein the output indicates glycemic dysfunction at a past time point.

22. 10. The system of claim 1, wherein the output device provides visualizations showing at least one of the behavioral impact of the glycemic dysfunction, past CGM and insulin data versus CGM replay simulated with optimal bolus at past bolus times, or past CGM and insulin data versus CGM replay simulated with optimal bolus at estimated meal times.

23. 10. The system of claim 1, further comprising a compliance engine configured to provide a visualization of a patient's compliance with ideal pre-meal bolus timing.

24. The system of claim 1 , wherein the output device provides dietary management information to the subject based on the amount of glycemic dysfunction.

25. 25. The system of claim 24, wherein the meal management information includes patient pre-bolus timing information.

26. 26. The system of claim 25, wherein the pre-bolus timing information includes a recommendation to the patient regarding a bolus administration before a certain type of future meal.

Images