Method for determining glucose changes in a subject - Patent Application 20070122997

Abstract

To provide a computer-implemented method and system for determining a glucose change in a subject.SOLUTION: A method includes: receiving an actual glucose measured value of a subject; receiving past subject model parameters; estimating subject model parameters of a state-based model of the subject on the basis of, the actual glucose measured value and the past subject model parameters; determining, by using a Kalman filter, an innovation parameter and an innovation covariance parameter on the basis of the subject model parameters and a previous state of the subject; calculating a test statistic on the basis of the determined innovation parameter and innovation covariance parameter; comparing the calculated test statistic with a given threshold; and in response to the calculated test statistic being above the given threshold, outputting an indication of the glucose change.SELECTED DRAWING: Figure 4

Description

[Technical Field]

[0001] CROSS-REFERENCE TO RELATED APPLICATIONS This application claims the benefit of U.S. Provisional Patent Application No. 62 / 871,931, filed July 9, 2019, which is incorporated herein by reference in its entirety.

[0002] The present technology relates generally to drug monitoring systems, and more particularly to methods and systems for determining whether a diabetic subject's glucose changes are abnormal or indicative of a problem, e.g., whether a meal consumed by a subject in an artificial pancreas system is not being recorded. [Background technology]

[0003] In healthy individuals, plasma glucose concentrations are tightly regulated by the action of hormones secreted by the endocrine pancreas, primarily insulin and glucagon. Insulin is secreted by pancreatic beta cells and signals the organ to absorb glucose, while glucagon is secreted by pancreatic alpha cells and signals the liver to produce glucose. In type 1 diabetes, insulin secretion is lost due to autoimmune destruction of beta cells.

[0004] Type 1 diabetes is currently treated with lifelong insulin replacement therapy administered using multiple injections (MDI) or continuous subcutaneous (subcutaneous) insulin infusion (CSII) via a portable pump. Both treatments follow a basal-bolus insulin injection pattern aimed at mimicking the physiological plasma insulin secretion seen in healthy individuals. Basal insulin represents the insulin required to maintain constant blood glucose levels in the fasting state, while insulin bolus is the dose of insulin typically given to cover the expected increase in glucose from a consumed meal.

[0005] Strict glucose control is important for patients with type 1 diabetes. Persistently elevated blood glucose levels (hyperglycemia) can lead to long-term complications, such as heart disease, blindness, kidney failure, and lower-limb amputation. Low blood glucose levels (hypoglycemia) are the limiting factor in glycemic control, as nonsevere hypoglycemia can lead to anxiety, nausea, confusion, blurred vision, and difficulty speaking, while severe hypoglycemia can cause coma or seizures and require assistance. A target HbA1c (a biomarker correlated with average blood glucose levels over a three-month period) of less than 7.0% is recommended for most patients with type 1 diabetes.

[0006] Despite advances in insulin analogs, insulin pumps, and continuous glucose sensors, most patients do not achieve acceptable glucose targets. Advances in glucose sensors have spurred the development of the artificial pancreas (AP), a closed-loop insulin delivery system that automatically regulates blood glucose levels in patients with type 1 diabetes. In the artificial pancreas, a control procedure adjusts the pump insulin infusion rate based on continuous glucose sensor readings. Artificial pancreas systems are considered the most promising treatment for type 1 diabetes. While attempts at fully automated closed-loop insulin delivery systems have been investigated, the most popular artificial pancreas systems still rely on user prompts to deliver insulin boluses with meals.

[0007] With conventional insulin therapy, the leading cause of poor glucose control in adolescents is missed mealtime insulin bolus delivery. Mealtime bolus failure was observed at least once per week in 65% of adolescents, and this was associated with significantly higher HbA1c levels compared with adolescents who missed boluses less than once per week (8.8% vs. 8.0%, respectively). Another study observed that more than one-third of adolescents missed 15% or more of their required boluses. Similar to conventional insulin therapy, the performance of closed-loop (CL) insulin delivery can also be affected after missed boluses. Adding a meal detection module to an artificial pancreas system that detects missing meals and signals the delivery of more insulin may improve the performance of the artificial pancreas.

[0008] In artificial pancreas systems, when an unknown meal is consumed, a closed-loop feedback mechanism responds to changes in blood glucose levels by modifying the pump's insulin basal rate. Generally, a significant amount of insulin is required to cover the glucose increase from the meal, sometimes up to 20% of the patient's total daily insulin dose. Consequently, without delivering an insulin bolus, the artificial pancreas cannot provide the required amount of insulin in a short period of time. Therefore, hyperglycemic events with undesirably high blood glucose levels are inevitable. Furthermore, if the feedback mechanism responds aggressively by injecting large amounts of insulin to prevent further increases in glucose, delayed postprandial hypoglycemia may occur due to slow absorption of the insulin delivery (as the delivered insulin continues to act beyond the absorption of the meal). A clear strategy is needed to mitigate hyperglycemia and hypoglycemia after a bolus failure. Summary of the Invention

[0009] It is an object of the present technology to ameliorate at least some of the inconveniences present in the prior art. One or more embodiments of the present technology may provide and / or extend the scope of the present technology's objectives and approaches and / or methods for achieving the objectives.

[0010] One or more embodiments of the present technology were developed based on the inventors' assessment that glucose control deteriorates significantly after a meal bolus miss. If the control algorithm is enhanced with meal detection technology, the performance of the closed-loop delivery system after a bolus miss can be improved.

[0011] The inventors appreciate that automatically detecting meals (for which a bolus was not delivered) and notifying diabetic subjects can improve the quality of life and health of diabetic users. In one non-limiting example, the system can notify the user, who can take action, such as delivering the forgotten insulin to themselves. In another non-limiting example, users of traditional pump therapy or multiple injection therapy can be reminded if they forget to eat a meal and deliver a bolus.

[0012] Such a system can be used to detect disorders that cause elevated glucose levels, such as a malfunctioning infusion set or a skipped meal.

[0013] The inventors also appreciate that such techniques can be used online or offline to analyze and model data, validate algorithm performance, and identify unannounced meals and hypoglycemic treatments, as non-limiting examples.

[0014] Accordingly, one or more embodiments of the present technology are directed to methods and systems for detecting glucose changes in a subject.

[0015] In accordance with a broad aspect of the present technology, there is provided a computer-implemented method, executable by an electronic device, for determining glucose change in a subject, the method including receiving subject model parameters of a state-based model of the subject, determining an innovation parameter and an innovation covariance parameter using a Kalman filter based on the subject model parameters and a prior state of the subject, calculating a test statistic based on the determined innovation parameter and innovation covariance parameter, comparing the calculated test statistic to a given threshold, and outputting an indication of glucose change in response to the calculated test statistic exceeding the given threshold.

[0016] In one or more embodiments of this method, the method further includes, prior to receiving the subject model parameters, receiving, by the electronic device, actual glucose measurements for the subject and receiving past subject model parameters, wherein receiving the subject model parameters for the state-based model for the subject includes estimating the subject model parameters based on the actual glucose measurements and the past subject model parameters.

[0017] In one or more embodiments of this method, the method further includes transmitting the indication to at least one of a display interface of the electronic device and an artificial pancreas system of the subject.

[0018] In one or more embodiments of this method, a test statistic above a given threshold indicates that the Kalman filter is inconsistent.

[0019] In one or more embodiments of the method, estimating said subject model parameters includes using maximum a posteriori (MAP) estimation.

[0020] In one or more embodiments of this method, estimating said subject model parameters is further based on previous glucose measurements, previous insulin measurements, and previous meals consumed.

[0021] In one or more embodiments of this method, a test statistic above a given threshold indicates that the innovation parameters are not independent and not identically distributed with a zero-mean Gaussian distribution with a covariance corresponding to the covariance of the innovation parameters.

[0022] In one or more embodiments of this method, the glucose change is indicative of an unknown meal, the unknown meal not being recorded by the subject.

[0023] In one or more embodiments of this method, the given threshold is based on a predetermined number of false positives.

[0024] In one or more embodiments of this method, the method further includes initializing the historical subject model parameters based on the subject's total daily dose, basal insulin, and carbohydrate ratio prior to receiving the historical subject model parameters.

[0025] In one or more embodiments of this method, the actual glucose measurement is received from a glucose sensor connected to the electronic device.

[0026] In one or more embodiments of this method, the method further includes determining an insulin bolus for an unknown meal not recorded by the given user based on the remaining meal, the patient's carbohydrate ratio, and blood glucose level prior to transmitting the indicator to at least one of the display interface of the electronic device and the subject's artificial pancreas system, wherein transmitting the indicator includes transmitting the insulin bolus.

[0027] In one or more embodiments of this method, the method further includes determining an unknown meal amount and an unknown meal time based on the innovation parameter and the innovation covariance parameter prior to determining the insulin bolus.

[0028] In one or more embodiments of this method, the calculated test statistic represents a cumulative sum of the correlation between the innovation parameter and glucose change based on unknown meal amounts and unknown meal times weighted by the innovation covariance parameter.

[0029] In one or more embodiments of this method, the given threshold is determined based on a given false positive rate of a random variable having a zero-mean Gaussian distribution and a covariance proportional to the square of the most likely glucose increase due to the most likely meal size and meal time weighted by an innovation covariance parameter.

[0030] In accordance with broad aspects of the present technology, there is provided a computer-implemented method for detecting a meal consumed by a patient, the method being executed by a processor, the method including: determining a discrepancy between an actual glucose measurement and a predicted glucose measurement; determining a probability that a meal was consumed based at least in part on the determined discrepancy; and determining a medication bolus in response to the determined probability.

[0031] In one or more embodiments of the method, determining the probability that the meal has been consumed is based, at least in part, on the actual blood glucose level, the target blood glucose level, and the remaining insulin.

[0032] In one or more embodiments of this method, the method further comprises estimating meal size and meal consumption time.

[0033] In one or more embodiments of this method, determining said medication bolus is based, at least in part, on at least one of an estimated meal size and an estimated time of consumption of the meal.

[0034] In one or more embodiments of the method, determining that said meal has been consumed is responsive to the determined probability exceeding a threshold.

[0035] In accordance with a broad aspect of the present technology, there is provided a system for determining glucose change in a subject, the system including a processor and a non-transitory storage medium operatively connected to the processor, the storage medium including computer-readable instructions, the processor, when executed, being configured to: receive subject model parameters of a state-based model for the subject, determine an innovation parameter and an innovation covariance parameter based on the subject model parameters and a prior state of the subject using a Kalman filter, calculate a test statistic based on the determined innovation parameter and innovation covariance parameter, compare the calculated test statistic to a given threshold, and output an indication of glucose change in response to the calculated test statistic being above the given threshold.

[0036] In one or more embodiments of this system, the processor is further configured to receive actual glucose measurements for the subject and receive past subject model parameters before receiving the subject model parameters, and receiving the subject model parameters for the state-based model for the subject includes estimating the subject model parameters based on the actual glucose measurements and the past subject model parameters.

[0037] In one or more embodiments of this system, the processor is further configured to send instructions to at least one of a display interface operably connected to the processor and the subject's artificial pancreas system.

[0038] In one or more embodiments of the system, a test statistic above a given threshold indicates an inconsistent Kalman filter.

[0039] In one or more embodiments of the system, estimating the subject model parameters includes using maximum a posteriori (MAP) estimation.

[0040] In one or more embodiments of the system, the estimation is further based on previous glucose measurements, previous insulin measurements, and previously consumed meals.

[0041] In one or more embodiments of this system, a test statistic above a given threshold indicates that the innovation parameters are not independent and not identically distributed with a zero-mean Gaussian distribution with a covariance corresponding to the covariance of the innovation parameters.

[0042] In one or more embodiments of the system, the glucose change is indicative of an unknown meal, where the unknown meal was not recorded by the subject.

[0043] In one or more embodiments of the system, the given threshold is based on a predetermined number of false positives.

[0044] In one or more embodiments of this system, the processor is further configured to, prior to receiving the aforementioned historical subject model parameters, initialize the historical subject model parameters based on the subject's total daily dose, basal insulin, and carbohydrate ratio.

[0045] In one or more embodiments of the system, the actual glucose measurement is received from a glucose sensor connected to the processor.

[0046] In one or more embodiments of this system, the processor is further configured to determine an insulin bolus for an unknown meal not recorded by the given user based on the remaining meal, the patient's carbohydrate ratio, and blood glucose level prior to transmitting the indicator to at least one of the display interface operably connected to the processor and the subject's artificial pancreas system, wherein transmitting the indicator includes transmitting the insulin bolus.

[0047] In one or more embodiments of this system, the processor is further configured to determine an unknown meal amount and an unknown meal time based on the innovation parameter and the innovation covariance parameter prior to determining the insulin bolus.

[0048] In one or more embodiments of this system, the test statistic represents a cumulative sum of the correlation between the innovation parameter and the glucose change based on unknown meal amounts and unknown meal times weighted by the innovation covariance parameter.

[0049] In one or more embodiments of this system, the given threshold is determined based on a given false positive rate of a random variable having a zero-mean Gaussian distribution and a covariance proportional to the square of the most likely glucose increase due to the most likely meal size and meal time weighted by an innovation covariance parameter.

[0050] According to another broad aspect, a computer-implemented method for detecting a meal consumed by a patient is provided. The method includes determining a discrepancy between an actual glucose measurement and a predicted glucose measurement. Based at least in part on the determined discrepancy, the method includes determining a probability that a meal was consumed. In response to the determined probability, the method includes determining a medication bolus.

[0051] In one embodiment of this method, the probability that a meal was consumed is based, at least in part, on the actual blood glucose level, the target blood glucose level, and the active insulin.

[0052] In one embodiment of this method, the method further comprises estimating meal size and meal consumption time.

[0053] In one embodiment of this method, the amount of the medication bolus is based, at least in part, on the estimated size of the meal and / or the estimated time of consumption of the meal.

[0054] In one embodiment of this method, the method further includes determining that the meal has been consumed in response to the determined probability exceeding a threshold.

[0055] According to another broad aspect, there is provided a system for detecting a meal consumed by a patient, the system including a processor and a non-transitory storage medium operatively connected to the processor, the storage medium including computer-readable instructions, the processor, when executed, configured to: determine a discrepancy between an actual glucose measurement and a predicted glucose measurement; determine a probability that a meal was consumed based at least in part on the determined discrepancy; and determine a medication bolus in response to the determined probability.

[0056] In one or more embodiments of the system, determining the probability that the aforementioned meal has been consumed is based, at least in part, on the actual blood glucose level, the target blood glucose level, and the remaining insulin.

[0057] In one or more embodiments of the system, the method further includes estimating meal size and meal consumption time.

[0058] In one or more embodiments of the system, determining said medication bolus is based, at least in part, on at least one of an estimated meal size and an estimated time of consumption of the meal. In one or more embodiments of the system, determining that the aforementioned meal has been consumed is responsive to the determined probability exceeding a threshold.

[0059] In the context of this specification, an "electronic device" is any computing device or computer hardware capable of executing software appropriate for the relevant task at hand. Accordingly, some (non-limiting) examples of electronic devices include general-purpose personal computers (desktops, laptops, netbooks, etc.), mobile computing devices, smartphones, tablets, and network equipment such as routers, switches, and gateways. Note that an electronic device in this context does not exclude acting as a server to other electronic devices. The use of the phrase "electronic device" does not exclude multiple electronic devices from being used to receive / send, perform, or cause to be performed tasks or requests, or the results of tasks or requests, or steps of methods described herein. In the context of this specification, a "client device" refers to any of a range of end-user client electronic devices associated with a user, such as a personal computer, tablet, smartphone, etc.

[0060] In the context of this specification, the expression "computer-readable storage medium" (also referred to as "storage medium" and "storage device") is intended to include non-transitory media of all natures and types, including but not limited to RAM, ROM, disks (CD-ROM, DVD, floppy disk, hard drive, etc.), USB keys, solid state drives, tape drives, etc. Multiple components may be combined to form a computer information storage medium that includes two or more media components of the same type and / or two or more media components of different types.

[0061] In the context of this specification, a "database" is a structured collection of data, regardless of its particular structure, database management software, or computer hardware on which the data is stored, implemented, or otherwise made available. A database may reside on the same hardware as the processes that store or utilize the information stored in the database, or it may reside on separate hardware, such as a dedicated server or multiple servers.

[0062] In the context of this specification, the expression "information" includes information of any nature or type that can be stored in a database, including, but not limited to, audiovisual works (images, videos, sound recordings, presentations, etc.), data (location data, numerical data, etc.), text (opinions, comments, questions, messages, etc.), documents, spreadsheets, lists of words, etc.

[0063] In the context of this specification, unless otherwise specified, an “indication” of an information element can be the information element itself or a pointer, reference, link, or other indirect mechanism that allows a recipient of the indication to locate a network, memory, database, or other computer-readable medium where the information element can be obtained. For example, an indication of a document might include the document itself (i.e., its contents), a unique document descriptor that identifies the file with respect to a particular file system, or some other means that directs the recipient of the indication to a network location, memory address, database table, or other place where the file can be accessed. As those skilled in the art will recognize, the degree of precision required for such an indication depends on the degree of prior understanding of the interpretation to be given to the information exchanged between the sender and recipient of the indication. For example, if, prior to communication between the sender and recipient, it is understood that an indication of an information element takes the form of a database key for an entry in a particular table of a predetermined database that contains the information element, transmitting the database key is all that is necessary to effectively convey the information element to the recipient, even if the information element itself was not transmitted between the sender and recipient of the indication.

[0064] In the context of this specification, the expression "communications network" is intended to include telecommunications networks such as computer networks, the Internet, telephone networks, telex networks, TCP / IP data networks (e.g., WAN networks, LAN networks, etc.), etc. The term "communications network" includes wired networks or direct-wired connections, as well as wireless media such as acoustic, radio frequency (RF), infrared and other wireless media, and combinations of any of the above.

[0065] In the context of this specification, words such as “first,” “second,” and “third” are used as adjectives solely to distinguish the nouns they modify, and not to describe any particular relationship between those nouns. Thus, for example, use of the terms “first server” and “third server” is not intended to imply a particular order, type, chronology, hierarchy, or ranking (for example) of / among servers, nor is their use (in and of itself) intended to imply that a “second server” must necessarily be present in a given situation. Furthermore, in other contexts, as discussed herein, reference to a “first” element and a “second” element does not exclude the two elements from being the same actual, real-world element. Thus, for example, in some cases, the “first” server and the “second” server may be the same software and / or hardware, while in other cases, they may be different software and / or hardware.

[0066] Each implementation of the present technology will have at least one, but not necessarily all, of the above-described objects and / or aspects, and it will be understood that some aspects of the present technology that arise from seeking to achieve the above-described object may not meet that object and / or may meet other objects not specifically set forth herein.

[0067] Additional and / or alternative features, aspects, and advantages of implementations of the present technology will become apparent from the following description, the accompanying drawings, and the appended claims. [Brief explanation of the drawings]

[0068] For a better understanding of the present technology, as well as other aspects and further features thereof, reference is made to the following description taken in conjunction with the accompanying drawings.

[0069] [Figure 1] 1 shows a schematic diagram of an electronic device in accordance with a non-limiting embodiment of the present technology; [Figure 2] 1 shows a schematic diagram of a system in accordance with a non-limiting embodiment of the present technology; [Figure 3] 1 shows a schematic diagram of a procedure for detecting missing meals according to a non-limiting embodiment of the present technology; [Figure 4] 1 shows a block diagram of a flowchart of a method for determining glucose change in a subject, the method being performed in accordance with a non-limiting embodiment of the present technology. [Figure 5A] An exemplary plot of the results of a sample simulation is shown, where the meal detection procedure detects a predicted meal and provides a 2 U bolus. Model variability often causes blood glucose levels to rise and fall for no apparent reason, making the meal detection procedure challenging. [Figure 5B] An example plot of a simulation is shown in which a false positive (FP) occurs, where a meal is flagged at 15:30, 3.5 hours after eating lunch, the algorithm delivers a 1.8 U bolus, and no hypoglycemia is observed for the next 4.5 hours. [Figure 6] Figure 1 shows exemplary plots of the percentage of time spent in hypoglycemia and hyperglycemia (compared to the 8 h after lunch) for three performed experiments (n=1536), where CL+B corresponds to no meal detection, lunch was announced and bloated, CL+MD corresponds to the use of the meal detection procedure with lunch not announced, and CL corresponds to no meal detection with lunch not announced. [Figure 7] 1 shows an exemplary plot of clinical data illustrating the performance of the meal detection procedure, in which a 60 g unknown meal was consumed at 13:00, the meal was detected at 13:40, and a 0.9 U bolus was delivered. [Figure 8] 1 shows exemplary plots of incremental glucose after consuming a meal without a bolus for four patients using conventional pump therapy, closed loop, or closed loop with meal detection, where diamonds indicate when a correction bolus was delivered for safety or automatically by the meal detection procedure. DETAILED DESCRIPTION OF THE INVENTION

[0070] The examples and conditional language described herein are intended primarily to aid the reader in understanding the principles of the technology and are not intended to limit its scope to such specifically described examples and conditions. Those skilled in the art will appreciate that various arrangements may be devised that, while not expressly described or shown herein, nevertheless embody the principles of the technology and are within its spirit and scope.

[0071] Furthermore, to aid in understanding, the following description may describe a relatively simplified implementation of the technology. As those skilled in the art will appreciate, various implementations of the technology may be more complex.

[0072] In some cases, what are believed to be useful examples of modifications to the technology may also be described. This is done merely to aid in understanding, and again, not to define the scope or describe the limits of the technology. These modifications are not an exhaustive list, and one of ordinary skill in the art may make other modifications while still remaining within the scope of the technology. Furthermore, if examples of modifications are not described, it should not be construed that the modification is not possible and / or that what is described is the only way to implement that element of the technology.

[0073] Furthermore, all statements herein describing principles, aspects, and implementations of the present technology, as well as specific examples thereof, are intended to encompass both structural and functional equivalents, whether now known or developed in the future. Thus, for example, those skilled in the art will understand that any block diagrams herein represent conceptual diagrams of illustrative circuitry embodying the principles of the present technology. Similarly, any flowcharts, flow diagrams, state transition diagrams, pseudocode, and the like, whether explicitly depicted as a computer or processor, will be understood to represent various processes substantially embodied in a computer-readable medium and that may be so executed by such a computer or processor.

[0074] The functionality of the various elements illustrated in the figures, including functional blocks labeled "processor" or "graphics processing unit," can be provided using dedicated hardware and hardware capable of executing software in association with appropriate software. When provided by a processor, the functionality may be provided by a single dedicated processor, a single shared processor, or multiple individual processors, some of which may be shared. In some non-limiting embodiments of the present technology, the processor may be a general-purpose processor such as a central processing unit (CPU) or a processor dedicated to a specific purpose, such as a graphics processing unit (GPU). Furthermore, the explicit use of the terms "processor" or "controller" should not be construed to refer only to hardware capable of executing software, but may implicitly include, but is not limited to, digital signal processor (DSP) hardware, network processors, application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), read-only memory (ROM), random access memory (RAM), and non-volatile storage for storing software. Other hardware, conventional and / or custom, may also be included.

[0075] Software modules, or modules that are simply implied to be software, may be represented herein as any combination of flowchart elements or other elements that illustrate the execution of process steps and / or textual representations. Such modules may be executed by explicitly or implicitly shown hardware.

[0076] With these foundations in place, some non-limiting examples will be considered to illustrate various implementations of aspects of the present technology.

[0077] Electronic Devices Referring to FIG. 1, an electronic device 100 suitable for use with some implementations of the present technology is shown, which includes various hardware components including one or more single-core or multi-core processors collectively represented by processor 110, a graphics processing unit (GPU) 111, a solid-state drive 120, random access memory 130, a display interface 140, and an input / output interface 150.

[0078] Communication between the various components of electronic device 100 may be enabled by one or more internal and / or external buses 160 (e.g., PCI bus, Universal Serial Bus, IEEE 1394 "Firewire" bus, SCSI bus, Serial ATA bus, etc.) to which the various hardware components are electronically coupled.

[0079] Input / output interface 150 may be coupled to touchscreen 190 and / or one or more internal and / or external buses 160. Touchscreen 190 may be part of a display. In some embodiments, touchscreen 190 is a display. Touchscreen 190 may also be referred to as screen 190. In the embodiment shown in FIG. 1 , touchscreen 190 includes touch hardware 194 (e.g., pressure-sensitive cells embedded in a layer of the display that enable detection of physical interaction between a user and the display) and touch input / output controller 192, enabling communication with display interface 140 and / or one or more internal and / or external buses 160. In some embodiments, input / output interface 150 may be connected to a keyboard (not shown), a mouse (not shown), or a trackpad (not shown), allowing a user to interact with electronic device 100 in addition to or instead of touchscreen 190.

[0080] According to an implementation of the present technology, solid-state drive 120 stores program instructions suitable for loading into random access memory 130 and executing by processor 110 and / or GPU 111 to determine whether a diabetic subject has consumed a meal. For example, the program instructions may be part of a library or an application.

[0081] The electronic device 100 may be a server, a desktop computer, a laptop computer, a tablet, a smartphone, a personal digital assistant, or any device that may be configured to implement the present technology, as can be understood by one skilled in the art.

[0082] system Referring to FIG. 2 , a schematic diagram of a system 200 is shown, which is suitable for implementing a non-limiting embodiment of the present technology. It should be clearly understood that the illustrated system 200 is merely an exemplary implementation of the present technology. Accordingly, the following description thereof is intended only as an explanation of an example of the use of the present technology. This description is not intended to define the scope or describe limits of the present technology. In some cases, what are believed to be useful examples of modifications to the system 200 may be described below. This is done merely to facilitate understanding, and again, it is not intended to define the scope or describe limits of the present technology. These modifications are not an exhaustive list, and other modifications may be possible, as will be understood by those skilled in the art. Furthermore, if this is not done (i.e., if an example of a modification is not described), it should not be construed that the modification is not possible and / or that what is described is the only way to implement that element of the present technology. As will be understood by those skilled in the art, this is not the case. Furthermore, it should be understood that the system 200 may, in some cases, provide a simple implementation of the present technology, and in such cases, they are presented in this manner to facilitate understanding. As one skilled in the art will appreciate, various implementations of this technology may be more complex.

[0083] The system 200 includes, among other things, the electronic device 100, a database 250, and an artificial pancreas system 220.

[0084] The system 200 is associated with a diabetic subject 205 or diabetic user 205 .

[0085] Electronic device 100 is associated with diabetic user 205. As a non-limiting example, electronic device 100 may be diabetic user 205's smartphone. Diabetic user 205 may input their own health and diabetes-related information into electronic device 100, which stores the information in database 250. In one embodiment, electronic device 100 may be part of an artificial pancreas system (e.g., a component of artificial pancreas system 220). In an alternative embodiment, electronic device 100 may be diabetic user 205's desktop computer.

[0086] Electronic device 100 is configured, among other things, to (i) model the glucose regulation system of diabetic user 205, (ii) predict glucose measurements, (iii) determine, based on pre-screened measurements, whether an insulin bolus was missed because a meal consumed by user 205 was not recorded in electronic device 100, and (iv) send information to artificial pancreas system 220 for insulin delivery to user 205. How electronic device 100 is configured to achieve its objectives is described in more detail below.

[0087] The artificial pancreas system 220, also known as a closed-loop system, an automatic insulin delivery system, or an autonomous system for glycemic control, is configured to mimic the glucose-regulating functions of a healthy pancreas. The artificial pancreas system 220 is operably connected to and associated with a diabetic user 205.

[0088] The artificial pancreas system 220 includes a continuous glucose monitoring system (CGM) 230, an insulin infusion pump 240, and a control process 245.

[0089] The CGM system 230 provides a steady stream of information reflecting the blood glucose level of the user 205. The CGM 230 includes a sensor placed subcutaneously under the patient's skin (not shown) that measures glucose in the fluid surrounding the cells (interstitial fluid), which is associated with blood glucose levels. The CGM system 230 may have a user interface, such as a screen or touchscreen (not shown), and / or may transmit the glucose-related information to the electronic device 100 of the user 205 or another electronic device (not shown) via a communications link (not numbered) over a communications network (not shown).

[0090] In one embodiment, glucose monitoring system 230 transmits information reflecting user's 205 blood glucose levels for storage in database 250 .

[0091] In one embodiment, electronic device 100 executes control procedures 245 that receive information from CGM 230 and perform a series of mathematical calculations. Based on these calculations, electronic device 100 sends medication instructions to the infusion pump. In alternative embodiments, control procedures 245 can execute on any number of devices, including insulin infusion pump 240, such as, but not limited to, a desktop computer, a remote server, and a smartphone.

[0092] The control procedure 245 includes a meal detection procedure 300, which is described in more detail below.

[0093] The insulin infusion pump 240 adjusts insulin delivery based on instructions received from the control procedure 245 .

[0094] In one embodiment, database 250 is configured to store a user-specific set of parameters 260 for user 205. User-specific set of parameters 260 may be used to model the glucoregulatory system of user 205. User-specific set of parameters 260 may include one or more of the patient's age, the patient's weight, endogenous glucose production, insulin-independent glucose flux, activation rate of insulin remote action, the patient's insulin sensitivity (e.g., insulin sensitivity of glucose transport, insulin sensitivity of glucose disposal, insulin sensitivity of suppression of EGP), insulin absorption rate, insulin elimination rate, time to maximum CHO absorption, insulin volume of distribution, the patient's total daily dose, the patient's basal insulin, the patient's carbohydrate ratio, the patient's meal, and glucose volume of distribution.

[0095] The database 250 is configured to store glucose measurements 262 for the user 205. In one embodiment, the glucose measurements 262 are received from the CGM 230. The glucose measurements 262 include, by way of non-limiting example, interstitial glucose concentrations. By way of non-limiting example, a glucose appearance rate from a meal can be calculated based on the glucose measurements 262.

[0096] The database 250 is configured to store insulin measurements 264 for the user 205. In one embodiment, the delivered insulin measurements 264 are received from the insulin infusion pump 240. The delivered insulin measurements 264 include the amount of subcutaneous insulin delivered and one or more of the amount of insulin pending to be delivered (i.e., held by request but not yet delivered), the amount of subcutaneous insulin not delivered, residual insulin, and an insulin pump malfunction or error.

[0097] The database 250 is configured to store consumed meal information 266 for the user 205. The user 205 may record metrics of consumed meals on the user's electronic device 100, which may transmit the metrics of consumed meals to the database 250. The consumed meal information 266 may include one or more of the following: meal composition, meal weight, meal composition, meal type, amount of protein in the meal, amount of fiber in the meal, amount of carbohydrates in the meal, or estimates thereof.

[0098] The database 250 is configured to store, for a given period of time, a set of model parameters 270 for the user 205. Generally speaking, the set of model parameters 270 are parameters that describe the glucoregulatory system of the user 205. The set of model parameters 270 typically change over time to adapt to the user 205. How the set of model parameters 270 is determined is described in more detail below.

[0099] Database 250 is configured to store a state estimate 280 for user 205. Generally speaking, state estimate 280 represents the state of diabetic user 205 at a given moment in time. Determining state estimate 280 is described in more detail below.

[0100] Meal detection procedure Turning now to FIG. 3, a schematic diagram of an unknown meal detection procedure 300 is shown, in accordance with a non-limiting embodiment of the present technology.

[0101] The unknown meal detection procedure 300 is performed by an electronic device including a processor, such as electronic device 100. In one embodiment, the unknown meal detection procedure 300 may be performed by the artificial pancreas system 220 or by another electronic device (not shown). It is contemplated that the unknown meal detection procedure 300 may be performed by different devices in a distributed manner.

[0102] In one embodiment, the unknown meal detection procedure 300 is part of the control procedure 245 .

[0103] The unknown meal detection procedure 300 is adapted to generate a glucoregulatory system model for the user 205 based on historical data for the user 205, predict glucose measurements using the glucoregulatory system model for the user 205, and compare the predicted glucose measurements with current glucose measurements to determine if the user 205 has missed a meal. In one embodiment, the unknown meal detection procedure 300 sends an indication of bolus noncompliance to the artificial pancreas system 220, which can cause the artificial pancreas system 220 to deliver an insulin bolus. In one embodiment, the indication of bolus noncompliance includes a recommendation of a bolus to be delivered. The unknown meal detection procedure 300 uses a state-space representation of the glucoregulatory system for the user 205.

[0104] The unknown meal detection procedure 300 includes a state space modeling procedure 320 , a probabilistic detection procedure 360 ​​, and an insulin bolus determination procedure 380 .

[0105] State Space Modeling Procedure The purpose of the state space modeling procedure 320 is to model the glucose regulatory system of the user 205. The state space modeling procedure 320 generates a mathematical model that describes one or more of the absorption of insulin from the subcutaneous tissue, the absorption of carbohydrates from a consumed meal, the change in glucose due to insulin action, and the change in glucose due to absorbed carbohydrates. The state space modeling procedure 320 predicts glucose measurements using Kalman filtering.

[0106] In one embodiment, a model of the user's 205 glucose regulatory system can be represented by a set of differential equations. In one embodiment, the user's 205 glucose regulatory system is described using a linear time-invariant model. As a non-limiting example, the Bergman model can be linearized to describe the user's 205 glucose regulatory system.

[0107] In one embodiment, the internal state of the model can be represented by: · Amount of subcutaneous insulin delivered; ·Plasma insulin concentration; · Amount of food digested; - Rate of glucose appearance from meals; glucose plasma concentration; and Interstitial glucose concentration

[0108] In one embodiment, the state space modeling procedure 320 includes a variable p n The set of model parameters 270 allows for representing the observed glucose measurements 262. The state space modeling procedure 320 uses the state space representation to determine state estimates 280. The state estimates 280 of the state space representation are values ​​that change over time in a way that depends on the values ​​they have at any time and also on the externally imposed values ​​of the input variables. The values ​​of the output variables depend on the values ​​of the state estimates.

[0109] Kalman filtering is then used to determine whether the glucose measurements are explained by the set of model parameters 270 and the delivered insulin measurements 264 and consumed meal information 266 .

[0110] A Kalman filter, also known as a linear quadratic estimation (LQE), is an algorithm that uses a series of measurements over time, which may contain noise and / or imprecision, to generate estimates of unknown variables that may be more accurate than those based on a single measurement. In other words, it is a set of equations that implement a predictor-correcting estimator to minimize the estimated covariance when conditions are respected, and the equations are executed recursively by an electronic device such as electronic device 100.

[0111] The state space modeling procedure 320 is configured to receive actual glucose measurements. In one embodiment, the state space modeling procedure 320 receives actual glucose measurements from the artificial pancreas system 220.

[0112] The state space modeling procedure 320 is configured to receive glucose measurements 262 from the CGM 230 and / or database 250. The glucose measurements 262 include N previous glucose measurements, z n ={z n-N+1 ,…,z n},

[0113] The state space modeling procedure 320 is configured to receive the amount of insulin delivered 264 from the insulin infusion pump 240 and / or the database 250. The amount of insulin delivered 264 includes the amount of insulin for time nN.

[0114] The state space modeling procedure 320 is configured to receive consumed meal information 266 from the database 250. The consumed meal information 266 includes consumed meals recorded by the user during time nN.

[0115] The amount of insulin delivered 264 and the meal consumed information 266 are both n ={U n-N ,…,U n-1} Note that other inputs can be added that may affect the user's 205 blood glucose level, such as, but not limited to, exercise and heart rate.

[0116] In one embodiment, the state X of the user 205 at time n n About p n The set of model parameters 270, denoted by , varies with the state according to the state space model.

[0117] X n =A(p n )X n-1+ B(p n )U n、 (1a)

[0118] y n =C(p n )X n、 (1b)

[0119] where U n are all inputs to the system, including the amount of insulin delivered 264 and the meal consumed information 266 at time n, as well as (A(p n ),B(p n ),C(p n )) is a set of parameters 270p n is a set of state, input, and output matrices.

[0120] In one embodiment, a standard linear Kalman filter is represented by the following equation:

[0121]

number

[0122] P n|n-1 =AP n-1 A T +Q, (2b)

[0123] S n =CP n|n-1 C T +R, (2c)

[0124]

number

[0125]

number

[0126] P n =P n|n-1 -K n CP n|n-1 , (2f)

[0127] During the ceremony,

number

[0128] Q is the process noise covariance matrix, R is the measurement noise covariance matrix, and K Nは This is the Kalman gain.

[0129] v n =z n -y n is the actual glucose reading Z n is an innovation parameter that indicates the discrepancy between the predicted measurements by the state space modeling procedure 320 and

number

[0130] In one embodiment, the innovation parameter can be considered to be or include an innovation covariance parameter. In one embodiment, the innovation parameter can be considered to be or include a test statistic.

[0131] The innovation parameter is the actual glucose measurement z n and the predicted measurement value y n In one embodiment, the innovation parameter quantifies how much the actual glucose measurement z n and the predicted measurement value y n Therefore, the higher the innovation parameter value, the greater the difference between the actual glucose reading z n and the predicted measurement value y n Conversely, the lower the innovation parameter value, the greater the discrepancy between the actual glucose reading Z n and the predicted measurement Y N The discrepancy between

[0132] Actual glucose reading Z n and the predicted measurement value y nIt will be appreciated that the innovation parameter indicative of discrepancy (or lack thereof) with the actual glucose measurement z may be determined in a variety of ways, and a correction factor or threshold may be used to determine the innovation parameter. n and the predicted measurement y n The value of the innovation parameter indicative of the discrepancy between the actual glucose measurement z n and the predicted measurement y n If the difference between the actual glucose measurement z is above (or below) a threshold, the value of the innovation parameter may be rounded to another value. n and the predicted measurement y n can be considered "equal" if they are within a given range.

[0133] S n is the innovation parameter v n is the covariance of

[0134] The Kalman filter calculates the true state X n The probability distribution function of

number

[0135] The state space modeling procedure 320 calculates p(p) of a set of model parameters 270 based on specific characteristics of the patient, e.g., a set of user-specific parameters 260, and general knowledge such as the total daily insulin dose, e.g., the amount of insulin delivered per day 264. n ) may be determined or received.

[0136] In the context of the present technology, the state space modeling procedure 320 adjusts or updates the set of model parameters 270 to fit the latest glucose trends, i.e., glucose measurements 262 received from the CGM 230 and / or database 250, insulin measurements 264 received from the insulin infusion pump 240 and / or database 250, and meal information 266 received from the user 205.

[0137] In one embodiment, X n-N is a known state at time nN, the state space modeling procedure 320 calculates the set of model parameters 270p n , state matrix

number

number

number

[0138] In one embodiment, a set of parameters 270p describing the last N glucose measurements. n The maximum likelihood estimator of

number

[0139]

number

[0140] It is contemplated that other methods may be used, such as, for example, recursive least squares.

[0141] Set parameters 270p n The maximum a posteriori probability estimator (MAP) of is given by:

[0142]

number

[0143] Assuming that the measurements are conditionally independent of each other when conditioned on the corresponding states and inputs, the distribution of glucose measurements 262 given the state, amount of insulin delivered 264, and set of parameters 270 can be expressed as follows:

[0144]

number

[0145] constant covariance r 2 Assuming zero-mean Gaussian measurement noise with ∑ k∈[nN,n], the distribution can be expressed as:

[0146]

number

[0147] The state space modeling procedure 320 uses maximum a posteriori estimation to adjust the set of model parameters 270. Next, the glucose measurements 262

number

number

[0148] In one embodiment, the state space modeling procedure 320 is configured to perform the following: At time k, the state space modeling procedure 320 calculates p n Glucose measurement Z based on a set of patient parameters 270 represented by n The Kalman state estimate corresponding to

number

number

number

number

number

[0149] The state space modeling procedure 320 stores in the database 250 a set of model parameters 270 and state estimates 280 for each iteration.

[0150] The state space modeling procedure 320 stores the actual glucose measurements z in the database 250. n and predicted glucose reading y n The innovation parameter v indicates the discrepancy between n , predicted glucose reading y n , the innovation parameter v n Covariance S n , and the Kalman gain K n In one embodiment, the value may be obtained from the artificial pancreas system 200.

[0151] In one embodiment, the state space modeling procedure may be performed by the artificial pancreas 220 and the output may be forwarded to the unknown meal detection procedure 320 performed by the electronic device 100.

[0152] Probabilistic Detection Procedures A probabilistic detection procedure 360 ​​is performed based on the state space modeling procedure 320 to determine if the user 205 is not logging meals via the electronic device 100, which causes a change in the glucose reading.

[0153] The innovation parameter, which indicates a discrepancy between the glucose measurement and the predicted glucose measurement, may have a large value (ie, relative to the other values ​​of the state parameters), which may be caused by a disturbance to the system.

[0154] Because disturbances may be due to other factors, the probabilistic meal detection procedure 360 ​​may use hypothesis testing to determine whether the disturbance is caused by a meal not recorded by the user 205. Two hypotheses are possible: H0: No unknown meals were consumed in the last M iterations (the Kalman filter is consistent). H1: A meal of size m is consumed without informing the system at time p∈[nM, n] (Kalman filter inconsistency).

[0155] For complex hypotheses that depend on unknown parameters θ (in this case θ = (p, m) unknown meal times and sizes by users 205), a generalized likelihood ratio test (GLRT) can be used. If Θ is the parameter space of θ, then the two hypotheses must satisfy the following condition:

[0156]

number

[0157] where Θ is a discrete set Θ={(p,m)|p∈[nM,n],m∈[m min ,m min +Δm,…,m max ]}, where m min , m max are the minimum and maximum detectable unknown diets, and Δm is the minimum detectable difference in the unknown diets.

number

[0158] The generalized likelihood ratio test (GLRT) is used. The GLRT statistic is written as follows:

[0159]

number

[0160] where V θ is a random variable with a probability distribution function that depends on θ. In this case, V θ is the Kalman filter innovation {v n-M ,…,v n} is a random variable that represents the process.

[0161] The null hypothesis that the Kalman filter is consistent is P(Vθ|H θ ) can be expressed as follows:

[0162]

number

[0163] Under the alternative hypothesis, P(V θ |H1) is written as follows for θ = (p, m):

[0164]

number

[0165] And for k∈[p+1,n],

[0166]

number

[0167] where U m is a column vector with zeros and the value of the meal input channel m, and I is the identity matrix.

[0168] If a meal of size m is consumed at time p, the correct state prediction of the Kalman filter hypothesis is

number

number

number

[0169] therefore,

number

number

number

number

number

number

[0170] By recursion, for k∈[p+1,n],

[0171]

number

[0172] Therefore, the true innovation parameter is

number

[0173]

number

[0174]

number

number

[0175] therefore:

number

[0176] In one embodiment, θ *= (p * ,m * )∈arg max P(V θ=(p,m) |H1) is defined in the probabilistic detection procedure 360 ​​as the most likely time and size of a hypothetical unknown meal.

[0177]

number

number

[0178]

number

[0179] Under the null hypothesis, λ follows a zero-mean Gaussian distribution with covariance

number

[0180] Therefore, the probabilistic detection procedure 360 ​​determines whether the parameter θ is greater than the threshold η when λ is less than the threshold η that satisfies P(λ>η|H0)<α. * In one embodiment, α=0.05. Other values ​​of α are contemplated as possible.

[0181] The probabilistic detection procedure 360 ​​sends information to the insulin bolus determination procedure 380 .

[0182] Insulin bolus determination procedure The insulin bolus determination procedure 380 receives an indication from the probabilistic detection procedure 360 ​​that a meal failure is likely.

[0183] Once a meal is detected by the probabilistic meal detection procedure 360, the insulin bolus determination procedure 380 determines the meal size m * and time p * of θ *= (p * ,m * )∈arg max P(V θ=(p,m) |H1). It is judged as follows.

[0184] The insulin bolus determination procedure 38 is *The system is configured to run another Kalman filter routine using new information about the patient's condition, resulting in a new state that includes a better estimate of the patient's condition.

number

number

[0185] The insulin bolus determination procedure 380 determines an insulin bolus, where the insulin bolus u is determined based on the remaining meal, the patient's carbohydrate ratio CR, blood glucose level G, and glucose target G. target , a patient-specific correction factor CF, and the remaining residual insulin (IOB). The insulin bolus can be expressed as:

[0186]

number

[0187] The insulin bolus determination procedure 380 sends an indication of the insulin bolus to the insulin infusion pump 240, which then injects the insulin bolus u. In one embodiment, the insulin bolus determination procedure 380 sends an indication of the insulin bolus for display to the user (e.g., as a notification on the electronic device 100) who can take appropriate action.

[0188] Show me how FIG. 4 shows a flowchart of a method 400 for determining glucose change in a subject, according to a non-limiting embodiment of the present technology.

[0189] In one embodiment, method 400 is performed by an electronic device including a processor operably connected to a non-transitory storage medium, such as electronic device 100.

[0190] In one embodiment, solid-state drive 120 stores computer-readable instructions suitable for being loaded into random access memory 130 and executed by processor 110 and / or GPU 111 of electronic device 100. Processor 110 is configured or operable to perform method 400 upon execution of the computer-readable instructions.

[0191] The method 400 begins at step 402 .

[0192] In step 402, electronic device 100 receives an actual glucose measurement of a subject, i.e., diabetic user 205. In one embodiment, the actual glucose measurement is received from CGM 230. In other embodiments, the actual glucose measurement may be stored in another non-transitory storage medium or may be received from another electronic device (not shown).

[0193] In step 404, the processor 110 receives the historical subject model parameters. In one embodiment, the historical subject model parameters are a set of model parameters 270 that are parameters representative of the user's 205 glucoregulatory system.

[0194] In step 406, processor 110 estimates subject model parameters of the state-based model for the subject based on the actual glucose measurements and the past subject model parameters. In one embodiment, electronic device 100 determines predicted glucose measurements based on the estimated subject model parameters. In another embodiment, steps 402 through 406 can be replaced by a single step of receiving the subject model parameters, which may have been determined by another electronic device (not shown).

[0195] In step 408, processor 110 uses a Kalman filter to determine an innovation parameter and an innovation covariance parameter based on the subject model parameters and the subject's previous state. In one embodiment, the innovation parameter is determined by multiplying the actual glucose measurement z n and predicted measurements

number

[0196] In step 410, processor 110 calculates a test statistic based on the determined innovation parameter and innovation covariance parameter. In one embodiment, the test statistic is calculated using equation (18).

[0197] In step 412, processor 110 compares the calculated test statistic to a given threshold value. In one embodiment, the given threshold value is predetermined based on a number of false positives.

[0198] In step 414, processor 110 outputs an indication that the meal was consumed by the subject in response to the calculated test statistic exceeding a given threshold. In one embodiment, electronic device 100 calculates a bolus value based on the calculated test statistic and transmits the bolus value to the artificial pancreas system.

[0199] The method 400 ends.

[0200] Turning now to Figures 5 through 8, several plots of simulations and clinical data experiments are shown.

[0201] Simulation Verification The simulation experiment was carried out with the following objectives: Calculate the sensitivity of the meal detection procedure, i.e., the number of unknown meals detected relative to the total number of unknown meals. Calculate the false alarm rate, i.e., the number of times the algorithm detected a meal when no meal was consumed. ·Evaluating the impact on overall glycemic control of introducing a meal detection procedure alongside a conventional closed-loop insulin dosing algorithm.

[0202] Simulation Setup The glucose regulatory system in T1D patients is nonlinear and time-varying. To simulate intra- and inter-patient variability, a simulation model with time-varying parameters, as presented by Wilinska et al., was implemented. To account for inter-patient variability, model parameters are randomly sampled from a prior distribution. Furthermore, intra-individual variability is considered by periodically varying some parameters (with random frequency and phase) (Table I). The simulation is augmented with correlated noise in the glucose measurements (coefficient of variation 7% and correlation 80%). [Table 1]

[0203] A simulation experiment, designated "CL+MD," is conducted using 512 virtual patients randomly sampled from the distribution in Table I. The meal detection procedure is implemented with a closed loop using a model predictive controller (MPC). The simulation experiment (Figure 5A) consists of a 13-hour simulation in which the virtual patient consumes a 40g carbohydrate (CHO) breakfast at 7:00 AM and a noon lunch consisting of either 40g, 60g, or 80g CHO.

[0204] Morning breakfast is entered into the dosing algorithm and a meal-accompanying bolus is given at breakfast. Lunch is given to the hypothetical patient but is not predicted by the insulin dosing algorithm. No meal is consumed after lunch, as the effect of unknown meals and a given bolus on the meal detection procedure is investigated. If plasma glucose is below 2.7 mmol / L, 15 g of rescue CHO is given to the hypothetical patient.

[0205] A total of 1,536 simulations (3 meal sizes x 512 hypothetical patients) were conducted in which the lunch meal was not flagged by the dosing algorithm. A true positive (TP) was counted when the meal detection procedure successfully flagged a meal within 120 minutes of the lunch meal. A false negative (FN) was counted when a meal was not flagged by the algorithm within 120 minutes of the lunch meal. Sensitivity is the ratio of TPs to the total number of unknown meals. The sensitivity of the meal detection procedure for all meals combined (40g, 60g, and 80g) is 93.23%. Other statistics can be found in Table II. Because the detection procedure is driven by increasing glucose, it is expected that the sensitivity of the algorithm will decrease with meal size (the lowest sensitivity is for a 40g meal). For an unknown medium-sized meal of 60g CHO, it is detected 96.29% of the time. On average, the algorithm detected a meal after a glucose jump above the threshold of 2.6 ± 1.2 mmol / L, with a detection time of approximately 40 minutes for unknown meals. These values ​​seem reasonable for identifying the effect of a meal on glucose increases. Similar detection time values ​​have been observed in other studies.

[0206] [Table 2]

[0207] A false positive (FP) is when a meal is detected in the absence of an unknown meal. Over 19,968 hours of simulation (13 hours x 1,536 simulations), 64 FPs occurred, representing a 4.17% FP rate per simulation. The relatively high rate of FPs after the 40 g meal (34 out of 64 false positives) is mostly due to delayed detection of the unknown meal (after the 120 min threshold) due to a slight increase in glucose. If 180 min were considered instead, the FP count would be 18 (instead of 34). Figure 5B shows a case where FP detection occurred after a slow glucose increase. The delivered bolus was safe and did not cause hypoglycemia.

[0208] Effect on blood sugar control Because classification algorithms are prone to flagging FPs, it is important to evaluate the impact of such events. It is also necessary to investigate the benefits on glucose control of adding a meal detection procedure to a closed-loop system. Therefore, two other simulation experiments were conducted to answer these two questions. Both experiments had the same structure as the CL+MD experiment. 1,536 simulations (3 meal sizes × 512 virtual patients) were conducted in which virtual patients used the closed-loop algorithm and consumed two meals: breakfast and lunch. However, in both experiments, the closed-loop algorithm consisted only of MPC without the meal detection procedure.

[0209] The first experiment, designated "CL+B," simulates a scenario in which lunch is predicted and a bolus is administered. The second experiment, designated "CL," simulates a scenario in which lunch is not predicted and the MPC only responds to changes in blood glucose. The two experiments serve to establish baseline values ​​for the expected time spent in hypoglycemia and time spent in hyperglycemia.

[0210] Figure 6 shows that adding the meal detection procedure to the closed-loop algorithm significantly improved the time spent in hyperglycemia from 34.9% to 30.4%, validating the effectiveness of the proposed meal detection procedure. Table III compares in more detail the increment in area under the curve (AUC) across three experiments with various meals. On average, the AUC improves by 19% from CL to CL+MD (the baseline is CL+B). [Table 3]

[0211] The meal detection procedure (CL+MD) is safe, as no increase in hypoglycemia was observed compared to when an accurate bolus was delivered (CL+B) (Figure 6). To further investigate the safety of the meal detection procedure when an FP is flagged and an unwanted bolus is delivered, the time spent in hyperglycemia was compared between simulations in which an FP was flagged (n=64) and simulations in which there was no FP (n=1472). It was found that the time spent in hypoglycemia when an FP was flagged (1.1 ± 0.35%) was not significantly different (p=0.38) from the time spent in hypoglycemia when there was no FP (0.76 ± 0.08%). This suggests that there is no obvious correlation between the detection of an FP and the cause of hypoglycemia according to the developed algorithm. The safety of the algorithm after an FP is attributed to the way the insulin bolus delivered after a meal was flagged was calculated. The calculated bolus is a term used to return blood glucose to the target value ((G - G target ) / CF-IOB), and a combination of terms covering the detected meals consumed.

number

number

[0212] Clinical validation Experiment Description Preliminary results from an ongoing clinical study evaluating the safety and efficacy of closed-loop insulin delivery with or without a meal detection module after bolus failure versus conventional pump therapy in hospitalized adolescents with type 1 disease (T1D) are presented. The study consisted of three randomized interventions per patient. Each patient consumed breakfast with an insulin bolus. Four hours after breakfast, patients were given a 60g lunch without a bolus. Depending on the intervention, insulin dosing was based on either the closed-loop algorithm, the closed-loop algorithm with a meal detection module, or the patient's conventional pump therapy. The intervention ended six hours after lunch. Figure 7 shows data from the intervention in which the meal detection procedure was used.

[0213] For patient safety, a correction bolus was delivered if blood glucose levels remained above 18 mmol / L. If this occurred, blood glucose levels were considered to remain constant until the end of the intervention. Figure 8 shows the incremental AUC for four patients who completed all interventions. A trend was observed indicating that the meal detection procedure may reduce the increase in AUC after bolus failure. In fact, the AUC decreased by 39% with the meal detection procedure compared to 16% without meal detection (baseline being conventional insulin therapy).

[0214] To further investigate the meal detection procedure, we performed the meal detection procedure offline using 108 hours of clinical data (4 patients x 3 visits x 9 hours). All 12 unknown meals were successfully detected and none were flagged as FP. The meal detection time was 35 minutes. The glucose increase at the time of meal detection was 2.89 ± 1.72 mmol / L, and the glucose increase 10 minutes before meal detection was 0.45 ± 0.73 mmol / L.

[0215] While the technology has been described in the context of an artificial pancreas system, it is contemplated that the technology could be used to notify a user of a missed dose of insulin and recommend a specific dose. The user could then take action, such as delivering the missed dose to themselves. In another application, a user of traditional pump or multiple injection therapy could be reminded if they forgot to eat and deliver a bolus.

[0216] The technology can also be used to detect disturbances that elevate glucose levels, such as infusion set failure or skipped meals. The technology can be used online or offline to analyze and model data, validate algorithm performance, and identify missing meals and hypoglycemic treatments, as non-limiting examples.

[0217] It should be clearly understood that not all technical advantages referred to herein need to be enjoyed in every embodiment of the present technology. For example, embodiments of the present technology may be implemented without users enjoying some of these advantages, while other non-limiting embodiments may be implemented without users enjoying other advantages or none at all.

[0218] Some of these steps and signal transmission and reception are well known in the art and, therefore, have been omitted in certain portions of this description for simplicity. Signals can be transmitted and received using optical means (such as fiber optic connections), electronic means (such as wired or wireless connections), and mechanical means (such as pressure-based, temperature-based, or other suitable physical parameter-based).

[0219] Modifications and improvements to the above-described implementations of the technology may become apparent to those skilled in the art. The foregoing description is intended to be illustrative, not limiting.

Claims

1. 1. A computer-implemented method for determining whether an unannounced meal has been consumed by a subject, the method being executable by an electronic device, the method comprising: receiving subject model parameters for a state-based model of the subject; For each set of hypothetical meals (each hypothetical meal having a meal size and a meal duration), calculating a test statistic according to a) and b) below; a) determining an innovation parameter and an innovation covariance parameter based on the subject model parameters, the subject's previous state, and the meal size and meal duration of the hypothetical meal using a linear Kalman filter, wherein the innovation parameter indicates a difference between an actual glucose measurement and a predicted glucose measurement; b) calculating the test statistic based on the determined innovation parameter and the innovation covariance parameter; and comparing each calculated test statistic to a given threshold. and in response to the calculated test statistic exceeding the given threshold, determining whether an unannounced meal was consumed by the subject and outputting an indication of the unannounced meal.

2. prior to receiving the subject model parameters; receiving, by the electronic device, an actual glucose measurement for the subject; receiving, by the electronic device, past subject model parameters; 2. The method of claim 1 , wherein receiving the subject model parameters of the state-based model for the subject comprises estimating the subject model parameters based on the actual glucose measurements and the past subject model parameters.

3. 3. The method of claim 1 or 2, wherein outputting the indication of the unannounced meal further comprises transmitting an indication of an insulin bolus to at least one of a display interface of the electronic device and an insulin delivery system of the subject.

4. The method of any one of claims 1 to 3, wherein the test statistic exceeding the given threshold indicates that the Kalman filter is inconsistent.

5. The method of claim 2 , wherein estimating the subject model parameters comprises using maximum a posteriori (MAP) estimation.

6. The method of claim 2 or 5, wherein estimating the subject model parameters is further based on previous glucose measurements, previous insulin measurements, and previously consumed meals.

7. 7. The method of any one of claims 1 to 6, wherein the test statistic exceeding the given threshold indicates that the innovation parameters are not independent and not distributed similarly to a zero-mean Gaussian distribution with a covariance corresponding to the innovation covariance parameter.

8. The method according to any one of claims 1 to 7, wherein the subject model parameters represent the glucoregulatory system of the subject.

9. The method of any one of claims 1 to 8, wherein the given threshold is based on a predetermined number of false positives.

10. prior to receiving the past subject model parameters; 7. The method of any one of claims 2, 5 or 6, further comprising initializing the past subject model parameters based on the subject's total daily dose, basal insulin, and carbohydrate ratio.

11. 11. The method of claim 2, 5, 6 or 10, wherein the actual glucose measurement is received from a glucose sensor connected to the electronic device.

12. prior to transmitting the indicator to the at least one of the display interface of the electronic device and the insulin delivery system of the subject; 4. The method of claim 3, further comprising determining an insulin bolus for the unannounced meal based on remaining meal size, the patient's carbohydrate ratio, and blood glucose level.

13. prior to determining said insulin bolus, The method of claim 12 , further comprising determining meal sizes and meal times for the unannounced meal based on the innovation parameter and the innovation covariance parameter.

14. 14. The method of claim 13, wherein the calculated test statistic represents a cumulative sum of correlations between the innovation parameter and glucose change based on the meal size and the meal time weighted by the innovation covariance parameter.

15. 15. The method of claim 14, wherein the given threshold is determined based on a given false positive rate of a random variable having a zero-mean Gaussian distribution and a covariance proportional to the square of the most likely glucose increase due to the most likely meal size and meal time weighted by the innovation covariance parameter.

16. 1. A system for determining an unannounced meal consumed by a subject, comprising: a processor; a non-transitory storage medium operably connected to the processor, the storage medium containing computer-readable instructions; When the processor executes the computer-readable instructions, receiving subject model parameters for a state-based model of the subject; For each set of hypothetical meals (each hypothetical meal having a meal size and a meal duration), calculating a test statistic according to a) and b) below; a) determining an innovation parameter and an innovation covariance parameter based on the subject model parameters, the subject's previous state, and the meal size and meal duration of the hypothetical meal using a linear Kalman filter, wherein the innovation parameter indicates a difference between an actual glucose measurement and a predicted glucose measurement; b) calculating the test statistic based on the determined innovation parameter and the innovation covariance parameter; and comparing each of the calculated test statistics to a given threshold value; In response to the calculated test statistic exceeding the given threshold, the system is configured to determine whether an unannounced meal was consumed by the subject and output an indication of the unannounced meal.

17. Prior to the processor receiving the subject model parameters, receiving an actual glucose measurement for the subject; receiving historical subject model parameters; 17. The system of claim 16, wherein receiving the subject model parameters of the state-based model for the subject comprises estimating the subject model parameters based on the actual glucose measurements and the past subject model parameters.

18. a display interface operatively connected to the processor; an insulin delivery system operably connected to the processor; 18. The system of claim 16 or 17, wherein the processor is further configured to transmit to at least one of the display interface and the subject's insulin delivery system an indication of an insulin bolus based on the estimated meal size of the unannounced meal.

19. The system of any one of claims 16 to 18, wherein the test statistic exceeding the given threshold indicates that the Kalman filter is inconsistent.

20. 20. The system of claim 17, wherein estimating the subject model parameters comprises using maximum a posteriori (MAP) estimation.

21. 21. The system of claim 17 or 20, wherein the estimating is further based on previous glucose measurements, previous insulin measurements, and previously consumed meals.

22. 22. The system of any one of claims 16 to 21, wherein the test statistic exceeding the given threshold indicates that the innovation parameters are not independent and not distributed similarly to a zero-mean Gaussian distribution with a covariance corresponding to the innovation covariance parameter.

23. The system according to any one of claims 16 to 22, wherein the subject model parameters represent the subject's glucoregulatory system.

24. The system of any one of claims 16 to 23, wherein the given threshold is based on a predetermined number of false positives.

25. Prior to the processor receiving the past subject model parameters, 22. The system of claim 17, 20, or 21, further configured to initialize the past subject model parameters based on the subject's total daily dose, basal insulin, and carbohydrate ratio.

26. further comprising a glucose sensor operably connected to the processor; 20. The system of claim 17, wherein the actual glucose measurement is received from the glucose sensor.

27. prior to the processor transmitting the indicator to at least one of the display interface operatively connected to the processor and the insulin delivery system of the subject; 20. The system of claim 18, further configured to determine an insulin bolus for the unannounced meal based on remaining meal size, the patient's carbohydrate ratio, and blood glucose level.

28. Prior to the processor determining the insulin bolus, 28. The system of claim 27, further configured to determine meal sizes and meal times for the unannounced meal based on the innovation parameter and the innovation covariance parameter.

29. 29. The system of claim 28, wherein the test statistic represents a cumulative sum of correlations between the innovation parameter and glucose change based on the meal size and the meal time weighted by the innovation covariance parameter.

30. 30. The system of claim 29, wherein the given threshold is determined based on a given false positive rate of a random variable having a zero-mean Gaussian distribution and a covariance proportional to the square of the most likely glucose increase due to the most likely meal size and meal time weighted by the innovation covariance parameter.

Images