# Patent

## Title
Method and apparatus for providing analyte monitoring system calibration accuracy

## Abstract

Methods, system and devices for providing improved calibration accuracy of continuous glucose monitoring system based on insulin delivery information are provided.

## Background
RELATED APPLICATION
The present application claims priority under 35 U.S.C. §119(e) to U.S. provisional application No. 61/230,686 filed Jul. 31, 2009, entitled “Method and Apparatus for Providing Analyte Monitoring System Calibration Accuracy”, the disclosure of which is incorporated herein by reference for all purposes.


BACKGROUND
As is known, Type-1 diabetes mellitus condition exists when the beta cells ((β-cells) which produce insulin to counteract the rise in glucose levels in the blood stream) in the pancreas either die or are unable to produce a sufficient amount of insulin naturally in response to elevated glucose levels. It is increasingly common for patients diagnosed with diabetic conditions to monitor their blood glucose levels using commercially available continuous glucose monitoring systems to take timely corrective actions. Some monitoring systems use sensors that require periodic calibration using a reference glucose measurement (for example, using an in vitro test strip). The FreeStyle Navigator® Continuous Glucose Monitoring System available from Abbott Diabetes Care Inc., of Alameda, Calif. is a continuous glucose monitoring system that provides the user with real time glucose level information. Using the continuous glucose monitoring system, for example, diabetics are able to determine when insulin is needed to lower glucose levels or when additional glucose is needed to raise the level of glucose.
Further, typical treatment of Type-1 diabetes includes the use of insulin pumps that are programmed for continuous delivery of insulin to the body through an infusion set. The use of insulin pumps to treat Type-2 diabetes (where the beta cells in the pancreas do produce insulin, but an inadequate quantity) has also become more prevalent. Such insulin delivery devices are preprogrammed with delivery rates such as basal profiles which are tailored to each user, and configured to provide the needed insulin to the user. In addition, continuous glucose monitoring systems have been developed to allow real time monitoring of fluctuation in glucose levels.
When the insulin delivery system and the glucose monitoring system are used separately, used together, or integrated into a single system, for example, in a single semi-closed loop or closed loop therapy system, the administered insulin (as well as other parameters or conditions) may affect some functions associated with the glucose monitoring system.
SUMMARY
In view of the foregoing, in aspects of the present disclosure, there are provided methods and apparatus for improving accuracy of the continuous glucose monitoring system calibration based at least on part on the insulin delivery information, and parameters associated with the administration of insulin.
Also provided are systems and kits.
INCORPORATION BY REFERENCE
The following patents, applications and/or publications are incorporated herein by reference for all purposes: U.S. Pat. Nos. 4,545,382; 4,711,245; 5,262,035; 5,262,305; 5,264,104; 5,320,715; 5,356,786; 5,509,410; 5,543,326; 5,593,852; 5,601,435; 5,628,890; 5,820,551; 5,822,715; 5,899,855; 5,918,603; 6,071,391; 6,103,033; 6,120,676; 6,121,009; 6,134,461; 6,143,164; 6,144,837; 6,161,095; 6,175,752; 6,270,455; 6,284,478; 6,299,757; 6,338,790; 6,377,894; 6,461,496; 6,503,381; 6,514,460; 6,514,718; 6,540,891; 6,560,471; 6,579,690; 6,591,125; 6,592,745; 6,600,997; 6,605,200; 6,605,201; 6,616,819; 6,618,934; 6,650,471; 6,654,625; 6,676,816; 6,730,200; 6,736,957; 6,746,582; 6,749,740; 6,764,581; 6,773,671; 6,881,551; 6,893,545; 6,932,892; 6,932,894; 6,942,518; 7,041,468; 7,167,818; and 7,299,082; U.S. Published Application Nos. 2004/0186365; 2005/0182306; 2006/0025662; 2006/0091006; 2007/0056858; 2007/0068807; 2007/0095661; 2007/0108048; 2007/0199818; 2007/0227911; 2007/0233013; 2008/0066305; 2008/0081977; 2008/0102441; 2008/0148873; 2008/0161666; 2008/0267823; and 2009/0054748; U.S. patent application Ser. Nos. 11/461,725; 12/131,012; 12/242,823; 12/363,712; 12/495,709; 12/698,124; 12/698,129; 12/714,439; 12/794,721; and 12/842,013; U.S. Provisional Application Ser. No. 61/347,754.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a block diagram illustrating an overall system in accordance with one embodiment of the present disclosure;
FIG. 2 is a flowchart illustrating calibration accuracy improvement routine in one aspect of the present disclosure;
FIG. 3 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 4 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 5 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 6 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 7 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 8 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 9 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 10 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 11 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure;
FIG. 12 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure; and
FIG. 13 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure.



DETAILED DESCRIPTION
Before embodiments of the present disclosure are described, it is to be understood that this disclosure is not limited to particular embodiments described, as such may, of course, vary. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting, since the scope of the present disclosure will be limited only by the appended claims.
Where a range of values is provided, it is understood that each intervening value, to the tenth of the unit of the lower limit unless the context clearly dictates otherwise, between the upper and lower limit of that range and any other stated or intervening value in that stated range, is encompassed within the disclosure. The upper and lower limits of these smaller ranges may independently be included in the smaller ranges is also encompassed within the disclosure, subject to any specifically excluded limit in the stated range. Where the stated range includes one or both of the limits, ranges excluding either or both of those included limits are also included in the disclosure.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure belongs. Although any methods and materials similar or equivalent to those described herein can also be used in the practice or testing of the present disclosure, the preferred methods and materials are now described. All publications mentioned herein are incorporated herein by reference to disclose and describe the methods and/or materials in connection with which the publications are cited.
It must be noted that as used herein and in the appended claims, the singular forms “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise.
The publications discussed herein are provided solely for their disclosure prior to the filing date of the present application. Nothing herein is to be construed as an admission that the present disclosure is not entitled to antedate such publication by virtue of prior disclosure. Further, the dates of publication provided may be different from the actual publication dates which may need to be independently confirmed.
As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present disclosure.
The figures shown herein are not necessarily drawn to scale, with some components and features being exaggerated for clarity.
Generally, embodiments of the present disclosure relate to methods and system for providing improved analyte sensor calibration accuracy based at least in part on the insulin delivery information. In certain embodiments, the present disclosure relates to the continuous and/or automatic in vivo monitoring of the level of an analyte using an analyte sensor, and under one or more control algorithms, determines appropriate or suitable conditions for performing calibration of the analyte sensor in view of the scheduled delivery of insulin or administered insulin amount. While the calibration accuracy of the analyte sensor is discussed in conjunction with the insulin delivery information, one or more other parameters or conditions may be incorporated to improve the calibration accuracy including, for example but not limited to, the physiological model associated with the patient using the analyte sensor, meal information, exercise information, activity information, disease information, and historical physiological condition information.
Embodiments include medication delivery devices such as external infusion pumps, implantable infusion pumps, on-body patch pumps, or any other processor controlled medication delivery devices that are in communication with one or more control units which also control the operation of the analyte monitoring devices. The medication delivery devices may include one or more reservoirs or containers to hold the medication for delivery in fluid connection with an infusion set, for example, including an infusion tubing and/or cannula. The cannula may be positioned so that the medication is delivered to the user or patient at a desired location, such as, for example, in the subcutaneous tissue under the skin layer of the user.
Embodiments include analyte monitoring devices and systems that include an analyte sensor, at least a portion of which is positionable beneath the skin of the user, for the in vivo detection of an analyte, such as glucose, lactate, and the like, in a body fluid. Embodiments include wholly implantable analyte sensors and analyte sensors in which only a portion of the sensor is positioned under the skin and a portion of the sensor resides above the skin, e.g., for contact to a transmitter, receiver, transceiver, processor, etc.
A sensor (and/or a sensor insertion apparatus) may be, for example, configured to be positionable in a patient for the continuous or periodic monitoring of a level of an analyte in a patient's dermal fluid. For the purposes of this description, continuous monitoring and periodic monitoring will be used interchangeably, unless noted otherwise.
The analyte level may be correlated and/or converted to analyte levels in blood or other fluids. In certain embodiments, an analyte sensor may be configured to be positioned in contact with dermal fluid to detect the level of glucose, which detected glucose may be used to infer the glucose level in the patient's bloodstream. For example, analyte sensors may be insertable through the skin layer and into the dermal layer under the skin surface at a depth of approximately 3 mm under the skin surface and containing dermal fluid. Embodiments of the analyte sensors of the subject disclosure may be configured for monitoring the level of the analyte over a time period which may range from minutes, hours, days, weeks, months, or longer.
Of interest are analyte sensors, such as glucose sensors, that are capable of in vivo detection of an analyte for about one hour or more, e.g., about a few hours or more, e.g., about a few days of more, e.g., about three or more days, e.g., about five days or more, e.g., about seven days or more, e.g., about several weeks or at least one month. Future analyte levels may be predicted based on information obtained, e.g., the current analyte level at time, the rate of change of the analyte, etc. Predictive alarms may notify the control unit (and/or the user) of predicted analyte levels that may be of concern in advance of the analyte level reaching the future level. This enables the control unit to determine a priori a suitable corrective action and implement such corrective action.
FIG. 1 is a block diagram illustrating an overall system in accordance with one embodiment of the present disclosure. Referring to FIG. 1, in one aspect, the system 100 includes an insulin delivery unit 120 that is connected to a body 110 of a user or patient to establish a fluid path to deliver medication such as insulin. In one aspect, the insulin delivery unit 120 may include an infusion tubing fluidly connecting the reservoir of the delivery unit 120 to the body 110 using a cannula with a portion thereof positioned in the subcutaneous tissue of the body 110.
Referring to FIG. 1, the system 100 also includes an analyte monitoring unit 130 that is configured to monitor the analyte level in the body 110. As shown in FIG. 1, a control unit 140 is provided to control the operation of the insulin delivery unit 120 and the analyte monitoring unit 130. In one embodiment, the control unit 140 may be a processor based control unit having provided therein one or more control algorithms to control the operation of the analyte monitoring unit 130 and the delivery unit 120. In one aspect, the control unit 140, the analyte monitoring unit 130 and the delivery unit 120 may be integrated in a single housing. In other embodiments, the control unit 140 may be provided in the housing of the delivery unit 120 and configured for communication (wireless or wired) with the analyte monitoring unit 130. In an alternate embodiment, the control unit may be integrated in the housing of the analyte monitoring unit 130 and configured for communication (wireless or wired) with the delivery unit 120. In yet another embodiment, the control unit 140 may be a separate component of the overall system 100 and configured for communication (wireless or wired) with both the delivery unit 120 and the analyte monitoring unit 130.
Referring back to FIG. 1, the analyte monitoring unit 130 may include an analyte sensor that is transcutaneously positioned through a skin layer of the body 110, and is in signal communication with a compact data transmitter provided on the skin layer of the body 110 which is configured to transmit the monitored analyte level substantially in real time to the analyte monitoring unit 130 for processing and/or display. In another aspect, the analyte sensor may be wholly implantable in the body 110 with a data transmitter and configured to wirelessly transmit the monitored analyte level to the analyte monitoring unit 130.
Referring still to FIG. 1, also shown in the overall system 100 is a data processing device 150 in signal communication with the one or more of the control unit 140, delivery unit 120 and the analyte monitoring unit 130. In one aspect, the data processing device 150 may include an optional or supplemental device in the overall system 100 to provide user input/output functions, data storage and processing. Examples of the data processing device 150 include, but are not limited to mobile telephones, personal digital assistants (PDAs), in vitro blood glucose meters, smart phone devices including Blackberry® devices, iPhone® devices, and Palm® devices, data paging devices, and the like, each of which include an output unit such as one or more of a display, audible and/or vibratory output, and/or an input unit such as a keypad, keyboard, input buttons and the like, and which are configured for communication (wired or wireless) to receive and/or transmit data, and further, which include memory devices such as random access memory, read only memory, volatile and/or non-volatile memory that store data.
Also shown in the overall system 100 is a data processing terminal 160 which may include a personal computer, a server terminal, a laptop computer, a handheld computing device, or other similar computing devices that are configured for data communication (over the internet, local area network (LAN), cellular network and the like) with the one or more of the control unit 140, the delivery unit 120, the analyte monitoring unit 130, and the data processing device 150, to process, analyze, store, archive, and update information.
It is to be understood that the analyte monitoring unit 130 of FIG. 1 may be configured to monitor a variety of analytes at the same time or at different times. Analytes that may be monitored include, but are not limited to, acetyl choline, amylase, bilirubin, cholesterol, chorionic gonadotropin, creatine kinase (e.g., CK-MB), creatine, DNA, fructosamine, glucose, glutamine, growth hormones, hormones, ketones, lactate, peroxide, prostate-specific antigen, prothrombin, RNA, thyroid stimulating hormone, and troponin. The concentration of drugs, such as, for example, antibiotics (e.g., gentamicin, vancomycin, and the like), digitoxin, digoxin, drugs of abuse, theophylline, and warfarin, may also be monitored. In those embodiments that monitor more than one analyte, the analytes may be monitored at the same or different times.
Additional detailed descriptions of embodiments of the continuous analyte monitoring device and system, calibrations protocols, embodiments of its various components are provided in, among others, U.S. Pat. Nos. 6,175,752, 6,284,478, 7,299,082 and U.S. patent application Ser. No. 10/745,878 filed Dec. 26, 2003 entitled “Continuous Glucose Monitoring System and Methods of Use”, the disclosures of each of which are incorporated herein by reference in their entirety for all purposes. Additional detailed description of systems including medication delivery units and analyte monitoring devices, embodiments of the various components are provided in, among others, U.S. patent application Ser. No. 11/386,915, entitled “Method and System for Providing Integrated Medication Infusion and Analyte Monitoring System”, the disclosure of which is incorporated herein by reference for all purposes. Moreover, additional detailed description of medication delivery devices and components are provided in, among others, U.S. Pat. No. 6,916,159, the disclosure of which is incorporated herein by reference for all purposes.
Referring back to FIG. 1, each of the components shown in the system 100 may be configured to be uniquely identified by one or more of the other components in the system so that communication conflict may be readily resolved between the various components, for example, by exchanging or pre-storing and/or verifying unique device identifiers as part of communication between the devices, by using periodic keep alive signals, or configuration of one or more devices or units in the overall system as a master-slave arrangement with periodic bi-directional communication to confirm integrity of signal communication therebetween.
Further, data communication may be encrypted or encoded (and subsequently decoded by the device or unit receiving the data), or transmitted using public-private keys, to ensure integrity of data exchange. Also, error detection and/or correction using, for example, cyclic redundancy check (CRC) or techniques may be used to detect and/or correct for errors in signals received and/or transmitted between the devices or units in the system 100. In certain aspects, data communication may be responsive to a command or data request received from another device in the system 100, while some aspects of the overall system 100 may be configured to periodically transmit data without prompting, such as the data transmitter, for example, in the analyte monitoring unit 130 periodically transmitting analyte related signals.
In certain embodiments, the communication between the devices or units in the system 100 may include one or more of an RF communication protocol, an infrared communication protocol, a Bluetooth® enabled communication protocol, an 802.11x wireless communication protocol, internet connection over a data network or an equivalent wireless communication protocol which would allow secure, wireless communication of several units (for example, per HIPAA requirements) while avoiding potential data collision and interference.
In certain embodiments, data processing device 150, analyte monitoring unit 130 and/or delivery unit 120 may include blood glucose meter functions or capability to receive blood glucose measurements which may be used, for example to calibrate the analyte sensor. For example, the housing of these devices may include a strip port to receive a blood glucose test strip with blood sample to determine the blood glucose level. Alternatively, a user input device such as an input button or keypad may be provided to manually enter such information. Still further, upon completion of a blood glucose measurement, the result may be wirelessly and/or automatically transmitted to another device in the system 100. For example, it is desirable to maintain a certain level of water tight seal on the housing of the delivery unit 120 during continuous use by the patient or user. In such case, incorporating a strip port to receive a blood glucose test strip may be undesirable. As such, the blood glucose meter function including the strip port may be integrated in the housing of another one of the devices or units in the system (such as in the analyte monitoring unit 130 and/or data processing device 150). In this case, the result from the blood glucose test, upon completion may be wirelessly transmitted to the delivery unit 120 for storage and further processing.
Any suitable test strip may be employed, e.g., test strips that only require a very small amount (e.g., one microliter or less, e.g., 0.5 microliter or less, e.g., 0.1 microliter or less), of applied sample to the strip in order to obtain accurate glucose information, e.g. FreeStyle® or Precision® blood glucose test strips from Abbott Diabetes Care Inc. Glucose information obtained by the in vitro glucose testing device may be used for a variety of purposes, computations, etc. For example, the information may be used to calibrate the analyte sensor, confirm results of the sensor to increase the confidence in the accuracy level thereof (e.g., in instances in which information obtained by sensor is employed in therapy related decisions), determine suitable amount of bolus dosage for administration by the delivery unit 120.
In certain embodiments, a sensor may be calibrated using only one sample of body fluid per calibration event. For example, a user need only lance a body part one time to obtain sample for a calibration event (e.g., for a test strip), or may lance more than one time within a short period of time if an insufficient volume of sample is obtained firstly. Embodiments include obtaining and using multiple samples of body fluid for a given calibration event, where glucose values of each sample are substantially similar. Data obtained from a given calibration event may be used independently to calibrate or combined with data obtained from previous calibration events, e.g., averaged including weighted averaged, etc., to calibrate.
One or more devices or components of the system 100 may include an alarm system that, e.g., based on information from control unit 140, warns the patient of a potentially detrimental condition of the analyte. For example, if glucose is the analyte, an alarm system may warn a user of conditions such as hypoglycemia and/or hyperglycemia and/or impending hypoglycemia, and/or impending hyperglycemia. An alarm system may be triggered when analyte levels reach or exceed a threshold value. An alarm system may also, or alternatively, be activated when the rate of change or acceleration of the rate of change in analyte level increase or decrease reaches or exceeds a threshold rate of change or acceleration. For example, in the case of the glucose monitoring unit 130, an alarm system may be activated if the rate of change in glucose concentration exceeds a threshold value which might indicate that a hyperglycemic or hypoglycemic condition is likely to occur. In the case of the delivery unit 120, alarms may be associated with occlusion conditions, low reservoir conditions, malfunction or anomaly in the fluid delivery and the like. System alarms may also notify a user of system information such as battery condition, calibration, sensor dislodgment, sensor malfunction, etc. Alarms may be, for example, auditory and/or visual. Other sensory-stimulating alarm systems may be used including alarm systems which heat, cool, vibrate, or produce a mild electrical shock when activated.
Referring yet again to FIG. 1, the control unit 140 of the system 100 may include one or more processors such as microprocessors and/or application specific integrated circuits (ASIC), volatile and/or non-volatile memory devices, and additional components that are configured to store and execute one or more control algorithms to dynamically control the operation of the delivery unit 120 and the analyte monitoring unit 130. The one or more closed loop control algorithms may be stored as a set of instructions in the one or more memory devices and executed by the one or more processors to vary the insulin delivery level based on, for example, glucose level information received from the analyte sensor.
An exemplary model describing the blood-to-interstitial glucose dynamics taking into account of insulin information is described below. More specifically, the model described herein provides for specific elaboration of model-based improvements discussed below. The example provided herein is based on a particular blood-to-interstitial glucose model, and while other models may result in a different particular relationship and parameter set, the underlying concepts and related description remain equally applicable.
Provided below is a model of blood-to-interstitial glucose as described by Wilinska et al. (Wilinska, Bodenlenz, Chassin, Schaller, Schaupp, Pieber, and Hovorka, “Interstitial Glucose Kinetics in Subjects With Type 1 Diabetes Under Physiologic Conditions”, Metabolism, v.53 n.11, November 2004, pp 1484-1492, the disclosure of which is incorporated herein by reference), where interstitial glucose dynamics comprises of a zero order removal of glucose from interstitial fluid F02, a constant decay rate constant k02, a constant glucose transport coefficient k21, and an insulin dependent glucose transport coefficient ki.

ġ(t)=−k02gi(t)+[k21+[ki[I(t)−Ib]]]gb(t)−F02  (1)

where gi corresponds to interstitial glucose, gb corresponds to blood glucose, the dot corresponds to the rate of change operation, (t) refers to variables that change over time as opposed to relatively static aforementioned coefficients, I corresponds to insulin concentration at any given time, and Ib corresponds to the steady-state insulin concentration required to maintain a net hepatic glucose balance.

It should be noted that the blood-to-interstitial glucose model described above is affected by insulin and accordingly, factoring in the insulin information will provide improvement to the sensor sensitivity determination.
The determination of insulin concentration (I) and the steady state insulin concentration required to maintain a net hepatic glucose balance (Ib) as shown in Equation (1) above may be achieved using insulin dosing history and an insulin pharmacokinetic and pharmacodynamic model. For example, based on a three compartment model of subcutaneous insulin dynamics into plasma insulin I as described by Hovorka, et al. (Hovorka, Canonico, Chassin, Haueter, Massi-Benedetti, Federici, Pieber, Schaller, Schaupp, Vering and Wilinska, “Nonlinear model predictive control of glucose concentration in subjects with type 1 diabetes”, Physiological Measurement, v.25, 2004, pp 905-920, the disclosure of which is incorporated herein by reference):




  
    
      
        
          
            
              
                I
                .
              
              1
            
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              
                -
                
                  k
                  a
                
              
              ⁢
              
                
                  I
                  1
                
                ⁡
                
                  (
                  t
                  )
                
              
            
            +
            
              
                u
                sc
              
              ⁡
              
                (
                t
                )
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              
                I
                .
              
              2
            
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              
                -
                
                  k
                  a
                
              
              ⁢
              
                
                  I
                  2
                
                ⁡
                
                  (
                  t
                  )
                
              
            
            +
            
              
                k
                a
              
              ⁢
              
                
                  I
                  1
                
                ⁡
                
                  (
                  t
                  )
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              I
              .
            
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              
                -
                
                  k
                  e
                
              
              ⁢
              
                I
                ⁡
                
                  (
                  t
                  )
                
              
            
            +
            
              
                
                  k
                  a
                
                V
              
              ⁢
              
                
                  I
                  2
                
                ⁡
                
                  (
                  t
                  )
                
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




where I1 and I2 are internal insulin compartments that describe the pathway from subcutaneous insulin injection into the plasma insulin compartment I. Ib is calculated by taking the steady-state average of I over a finite window of past and present period. The coefficients ka and ke describe the various decay and transport rates of the compartments, and V is the plasma insulin volume. Insulin action time is related to the parameter ka. The input usc to this model is described in terms of subcutaneous insulin infusion rate. Insulin dose/bolus may be converted into its delivery rate equivalent by monitoring or estimating the actual amount of bolus amount/dose delivered after every regular intervals of time (e.g. by monitoring of the amount of bolus/dose delivered every minute for a given executed bolus dose delivery).

For analyte monitoring systems, an uncalibrated sensor measurement yCGM is related to the true interstitial glucose by the following equation:

yCGM(t)=S[gi(t)+vi(t)]  (3)

where S is the calibration sensitivity to be identified, and vi is sensor noise.

Further, reference blood glucose measurement yBG when available at certain times, such as when requested for calibration at time to, contaminated by measurement error vb may be expressed as follows:

yBG(to)=gb(to)+vb(to)  (4)

Accordingly, the models and functional relationships described above provide some exemplary system components for providing improvement to the calibration accuracy in analyte monitoring systems whether used as a standalone system, or in conjunction with a medication delivery system such as with an insulin pump.
Determination of the suitable or appropriate time period to perform sensor calibration routine may be accomplished in several manners within the scope of the present disclosure. In one aspect, the calibration schedule may be predetermined or preset based on the initial sensor insertion or positioning in the patient or alternatively, scheduled based on each prior successful calibration event on a relative time basis. In some aspects, calibration routines are delayed or cancelled during high rates of glucose fluctuation because physiological lag between interstitial glucose measured by the analyte sensor and the blood glucose measured by discrete in vitro test strips may result in an error in the sensor sensitivity estimation.
In one aspect, calibration routine or function may be prevented or rejected when the interstitial glucose absolute rate of change is determined to exceed a predetermined threshold level. As the interstitial glucose level generally lags blood glucose level, there may be time periods where the blood glucose may be changing rapidly while the measured interstitial glucose level may not report similar fluctuations—it would change rapidly at some later, lagged time period. In such a case, a lag error may be introduced to the sensitivity determination. Accordingly, in one aspect, the execution of the calibration routine may be delayed or postponed when a sensor calibration request is detected by the system 100 during a time period when an insulin dose of sufficient magnitude is delivered, which may cause the rapid change in blood glucose to occur without a rapid change of interstitial glucose at that instance.
Referring now to the Figures, FIG. 2 is a flowchart illustrating overall calibration accuracy improvement routine in one aspect of the present disclosure. Referring to FIG. 2, when calibration start event is detected (210), for example, based on a predetermined calibration schedule from sensor insertion, or in response to a user calibration function initiation or execution, it is determined whether the initiated calibration routine is to be executed based on, for example, insulin information (220). Thereafter, one or more data or information associated with the determination is used to generate an output (230) which may, in one aspect, be provided to the user and/or stored in the system 100 (FIG. 1).
FIG. 3 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. As shown, when the calibration start event is detected (310), it is determined whether an insulin dose (for example, a bolus amount such as a carbohydrate bolus, or a correction bolus dose) was delivered or administered to the patient (320). In one aspect, as part of determining whether the insulin dose was delivered, it may be also determined whether the insulin dose was delivered within a time period measured from the detected calibration start event (and further, optionally, whether the determined insulin dose delivered amount meets a predetermined threshold level of insulin).
Referring again to FIG. 3, if it is determined that the insulin dose was delivered, then the routine proceeds to step 340 where the initiated calibration routine is not executed, and the routine returns to the beginning and awaits for the detection of the next or subsequent calibration start event. On the other hand, if at step 320 it is determined that the insulin dose was not delivered, then at step 330, the initiated calibration routine is executed to determine, for example, the corresponding sensor sensitivity based on a contemporaneously determined reference measurement (e.g., blood glucose measurement from an in vitro test strip, or another sensor data point that may be used as reference measurement) to calibrate the sensor.
FIG. 4 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring to FIG. 4, in the embodiment shown, when the calibration start event is detected (410) it is determined whether the insulin on board (IOB) level exceeds a predetermined threshold level (420). That is, in one aspect, the control algorithm may be configured to determine, in response to the detection of a calibration routine initialization, the IOB level. In one aspect, if it is determined that the IOB level exceeds the predetermined threshold level, then the initiated calibration routine is not contemporaneously executed (440), but rather, the called routine may be delayed, postponed, or cancelled, and the routine returns to the beginning to detect the subsequent calibration start event.
Referring to FIG. 4, if on the other hand it is determined that the IOB level is not greater than the predetermined threshold level at step 420, then the initiated calibration routine is executed at step 430 (530 (FIG. 5)), as discussed above, for example, to determine the corresponding analyte sensor sensitivity based on one or more reference glucose measurements to calibrate the sensor data.
FIG. 5 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Compared to the embodiment described in conjunction with FIG. 4, in the embodiment shown in FIG. 5, when the IOB level is determined to exceed the predetermined threshold level (520), then again, the initiated calibration routine is not executed (540), but prior to returning to the beginning of the routine to detect the subsequent calibration start event (510), a user notification function is called to notify the user of a failed (or delayed/postponed) calibration event (550). Such notification may include one or more of a visual indication, an audible indication, a vibratory indication, or one or more combinations thereof.
FIG. 6 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring to FIG. 6, in a further aspect of the present disclosure, when a calibration start event such as the initialization of a scheduled calibration routine is detected (610), it is determined whether a predetermined or categorized event has been logged at step 620. In particular, the control algorithm may be configured to determine whether an event such as a meal event, an activity event, an exercise event, or any other suitable or classified event has been logged at step 620. As discussed above, in one aspect, the control algorithm may be configured to additionally determine the time period of when such event was logged, if any, to determine whether the determined time period falls within a relevant time period with respect to the initiated calibration routine.
For example, if the logged meal event occurred with sufficient temporal distance relative to the initiated calibration routine, that it likely will have minimal relevance, if any to the calibration accuracy associated with the analyte sensor, then such logged event may be ignored. Alternatively, with each retrieved logged event at step 620, the routine may be configured to determine whether the logged event occurred within a specified or predetermined time period, in which case, the routine proceeds to step 640 where the initiated calibration routine is not executed and/or postponed or delayed. As further shown in FIG. 6, the routine thereafter returns to the beginning and monitors the system to determine whether a subsequent calibration start event is detected.
Referring back to FIG. 6, if at step 620 there are no events logged which are classified or categorized as relevant or associated with a parameter that is considered to be relevant, or alternatively, the one or more logged events detected or retrieved fall outside of the predetermined relevant time period (for example, within one hour prior to the calibration start event detected), then the initiated calibration routine proceeds at step 630 and is executed to determine, for example, the sensitivity associated with the analyte sensor based, for example, on a received reference blood glucose measurement, to calibrate the sensor data.
In aspects of the present disclosure, the duration and/or threshold described may be determined based on parameters including, for example, but not limited to insulin sensitivity, insulin action time, time of day, analyte sensor measured glucose level, glucose rate of change, and the like.
Moreover, in aspects of the present disclosure, as discussed, if the condition described above is detected, rather instead of delaying or postponing the execution of the calibration routine, the sensitivity determination may be altered as described in further detail below. That is, in one aspect, a correction factor may be applied to the sensitivity determination based on the insulin dose amount, elapsed time since the administration of the insulin dose, insulin sensitivity and insulin action time, for example. In one aspect, the correction factor may be a predetermined value or parameter, for example, based in part on the model applied to the patient's physiological condition, or may be a factor that is configured to be dynamically updated in accordance with the variation in the monitored parameters such as those described above.
In a further aspect, a glucose model of a patient may be used to predict or determine future glucose (blood and/or interstitial) levels and to estimate present glucose levels (blood and/or interstitial). More specifically, in aspects of the present disclosure, the model applied may also be used to estimate a rate-of-change of these variables and higher order moments of these variables in addition to statistical error estimates (for example, uncertainty estimates).
As discussed, the insulin delivery information and the measured glucose data from the analyte sensor (e.g., multiple measurements of each in time) are two of many input parameters used in conjunction with the embodiments described herein. Accordingly, in one aspect, the calibration routine may be configured to use the predicted output(s) as a check or verification to determine if the calibration routine should be postponed or delayed. For example, if the rate of change of blood glucose is determined to exceed a predetermined threshold, the calibration routine may be postponed or delayed for a predetermined time period. Alternatively, in a further aspect, if it is determined that the uncertainty in the interstitial estimate exceeds a predetermined threshold, the calibration routine may be configured to be postponed or delayed for a predetermined time period. The predetermined time period for a delayed or postponed calibration routine may be a preset time period, or alternatively, dynamically modified based on, for example, but not limited to the level of determined uncertainly in the interstitial estimate, the level of the predetermined threshold, and/or any other relevant parameters or factors monitored or otherwise provided or programmed in the system 100 (FIG. 1).
Referring now again to the Figures, FIG. 7 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring to FIG. 7, in the embodiment shown, when the calibration start event is detected at step 710 the routine determines one or more physiological model outputs based on one or more present and/or past input parameters and values (720) including, for example, monitored sensor data, insulin delivery information, blood glucose estimates, blood glucose rate of change estimate values, and the like. Thereafter, at step 730, it is determined whether the rate of change of the estimated glucose level deviates from a predetermined threshold (for example, where the estimated rate exceeds a preset positive value, or the estimated rate falls below a preset negative value). If it is determined that the estimated glucose rate of change is not within the predetermined threshold at step 730, then at step 750, the routine discontinues the calibration function (or postpones or delays the initiated calibration routine). Thereafter, as shown in FIG. 7, the routine returns to the beginning to detect the subsequent calibration start event at step 710.
Referring still to FIG. 7, if at step 730 it is determined that the estimated glucose rate of change is within the predetermined threshold, then at step 740, the routine proceeds with the execution of the calibration routine to determine, for example, the sensitivity associated with the analyte sensor by prompting the user to input a reference blood glucose measurement value (for example, based on an in vitro blood glucose testing), or the system may be configured to retrieve an existing or contemporaneously received reference measurement data to determine the sensitivity value for calibrating the sensor data.
FIG. 8 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring now to FIG. 8, when the calibration start event is detected at step 810 and the model outputs are determined based on one or more present and/or past input parameters or values (820) as discussed above in conjunction with FIG. 7, in the embodiment shown in FIG. 8, the calibration routine is executed based, in part on the estimated glucose value and/or the determined rate of change of the glucose level (830). That is, in one embodiment, when the scheduled calibration routine is initiated, the routine determines the most suitable or accurate parameters or values that are available to proceed with the execution of the calibration routine (as compared to determining whether or not the calibration condition is appropriate).
FIG. 9 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring to FIG. 9, in one aspect, when the calibration start event is detected at step 910, it is thereafter determined when an insulin dose has been delivered at step 920. That is, when a scheduled calibration routine is called or initiated, the routine determines whether there has been insulin dose delivery that may impact the conditions associated with the calibration of the analyte sensor. For example, in one aspect, the routine may determine whether the insulin dose is delivered within a predetermined time period measured from the initiation of the calibration routine (step 910) such as, within the past 1-2 hours, for example. That is, the system may be configured such that insulin dose administered outside such predetermined time period may be considered not sufficiently significant to adversely affect the conditions related to the calibration of the analyte sensor, and therefore, ignored.
Referring again to FIG. 9, when it is determined that the insulin dose was delivered (920) for example, during the relevant predetermined time period, the scheduled calibration function is delayed for a predetermined or programmed time period. That is, the scheduled calibration function is executed after the programmed time period has expired at step 940 (such that any potentially adverse affect of the detected insulin dose delivery (at step 920) has dissipated sufficiently during the programmed time period to proceed with the calibration routine). On the other hand, if it is determined that there is no insulin dose delivery detected (920) or any detected insulin dose delivery falls outside the relevant time period, then at step 930, the initiated calibration routine is performed as described above. In this manner, in one aspect of the present disclosure, when insulin dose administration such as bolus dose administration is detected within a relevant time period during a scheduled or user initiated calibration routine, a time delay function is provided to dissipate the effects of the administered insulin dose, before calibration routine resumes.
FIG. 10 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring to FIG. 10, in the embodiment shown, upon detection of the calibration start event 1010, it is determined whether a medication dose (such as insulin dose) was delivered (1020) (for example, during a relevant time period as described above in conjunction with FIG. 9 above). If not, then the calibration routine is executed to completion at step 1050. On the other hand, if it is determined that the medication dose was delivered during the relevant time period (1020) (for example, within 1-2 hours of the detected calibration start event), at step 1030, the amount of delivered medication dose is compared against a threshold level to determine whether the delivered medication dose exceeds the threshold level. If not, then the calibration routine is executed or performed to completion as described above at step 1050.
If on the other hand it is determined that the delivered medication dose exceeds the threshold level, then at step 1040, the detected start of the calibration event is delayed or postponed for a preprogrammed time period. In one aspect, the preprogrammed time period may be dynamically adjusted based on the amount of the medication dose that exceeds that threshold level, or alternatively, the preprogrammed time period may be a fixed value. In this manner, in one aspect, when it is determined that medication dose was administered contemporaneous to a scheduled calibration event, the routine may be configured to determine the relevance of the delivered medication dose to modify the calibration timing accordingly (for example, to continue with the execution of the calibration routine or to delay the calibration routine to minimize any potential adverse effect of the delivered medication dose).
FIG. 11 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring to FIG. 11, in the embodiment shown, when the calibration start event is detected 1110, a model based on one or more output values is determined based on one or more present and/or past input parameters or values (1120) as discussed above in conjunction with FIGS. 7 and 8 above, for example. It is to be noted that the model based determination as described herein may include one or more physiological models determined to a particular individual, condition and/or the severity of the condition or customized for one or more specific applications.
Referring to FIG. 11, after the model based outputs are determined at step 1120, it its determined whether the determined outputs or estimates of the outputs are within a predetermined threshold level at step 1130. That is, output parameters or values are determined based on one or more predetermined model applications relevant to, for example, the glycemic profile of a patient or a type of patients, and thereafter, the determined or estimated output parameters are compared to the predetermined threshold level. When it is determined that the estimated outputs are not within the threshold level, then at step 1150, the initiated calibration routine is delayed or postponed for a predetermined time period before executing the calibration function to completion as described above.
On the other hand, as shown in FIG. 11, if it is determined that the estimated output parameters or values are within the predetermined threshold value, the at step 1140, the calibration routine is executed, for example, to determine the sensitivity associated with the analyte sensor based on available reference glucose data, and thereafter calibrating the sensor data.
As discussed, in aspects of the present disclosure, the calibration accuracy routines may include other parameters or data such as, for example, meal intake information. For example, an aspect of the calibration routine may include confirming or determining whether a meal event has occurred for example, within the last hour prior to the scheduled calibration event, and further postpone or delay calibration if it is determined that the consumed meal during the past hour was sufficiently large or greater than a set threshold amount (for example, based on carbohydrate estimate). In one aspect, the meal intake information parameter used in conjunction with the calibration routine may be performed in conjunction with the insulin dose information as described above, or alternatively, as a separate routine for determining or improving the timing of performing the calibration routines.
In another aspect, the insulin dose information and/or other appropriate or suitable exogenous data/information may be used to improve the sensor sensitivity determination. For example, in one aspect, a model may be used to account for blood glucose and interstitial glucose, and insulin measurement data is used to help compensate for the lag between the two. The model would produce a blood glucose estimate that could be related to the reference blood glucose estimate in order to determine the sensitivity. Alternatively, the sensitivity could be part of the model and estimated. Additional detailed description related to pump information to improve analyte sensor accuracy is provided in U.S. patent application Ser. No. 12/024,101 entitled “Method and System for Determining Analyte Levels”, the disclosure of which is incorporated by reference for all purposes.
More specifically, referring back to and based on an example of the blood-to-interstitial glucose dynamics model which accounts for insulin, an estimated sensitivity at time t0 that is a function of available reference blood glucose (BG) measurement, analyte sensor measurement, and insulin information can be described as below:




  
    
      
        
          
            S
            ^
          
          ⁡
          
            (
            
              t
              o
            
            )
          
        
        =
        
          
            
              [
              
                
                  
                    
                      y
                      .
                    
                    CGM
                  
                  ⁡
                  
                    (
                    
                      t
                      o
                    
                    )
                  
                
                +
                
                  [
                  
                    
                      k
                      02
                    
                    ⁢
                    
                      
                        y
                        CGM
                      
                      ⁡
                      
                        (
                        
                          t
                          o
                        
                        )
                      
                    
                  
                  ]
                
              
              ]
            
            +
            
              F
              02
            
          
          
            
              [
              
                
                  k
                  21
                
                +
                
                  [
                  
                    
                      k
                      i
                    
                    ⁡
                    
                      [
                      
                        
                          I
                          ⁡
                          
                            (
                            
                              t
                              o
                            
                            )
                          
                        
                        -
                        
                          I
                          b
                        
                      
                      ]
                    
                  
                  ]
                
              
              ]
            
            ⁢
            
              
                y
                BG
              
              ⁡
              
                (
                
                  t
                  o
                
                )
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




It is to be noted that if insulin information is not accounted for, as shown in Equation 5 above, the denominator will be smaller, resulting in the sensitivity estimate larger than the actual value.

In another aspect, a closed loop control system is contemplated where a portion of the control algorithm seeks not only to prevent glucose excursions outside the euglycemic range, but also to provide improved conditions for calibration. While two particular conditions are described as examples, within the scope of the present disclosure, other conditions may be contemplated that are suitable or appropriate, depending on the type of analyte sensor used and/or other factors, variables or parameters.
In some cases, two conditions or states generally provide better calibration performance (i.e., better accuracy in sensitivity estimate)—calibrating during higher glucose periods and during low glucose rates-of-change. Calibrating during high glucose episodes is favorable because some errors tend to be unrelated to glucose level and will contribute to the sensitivity calculation proportionally less when glucose is high. In addition, as discussed above, error induced due to lag between blood glucose and interstitial glucose is minimized when glucose rate-of-change is low.
FIG. 12 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring to FIG. 12, in the embodiment shown, when the calibration start event is detected (1210), it is determined whether the current or an anticipated or estimated glucose level is within a predetermined threshold level (1220). In certain embodiments, the threshold level is a higher than average glucose level. As described above, a higher than average glucose level may be favorable in certain embodiments for calibration because some errors may be proportionally less when the glucose level is high. In one aspect, if it is determined that the glucose level is not within the predetermined threshold, then the initiated calibration routine is not contemporaneously executed, but rather, the scheduled calibration function is delayed for a predetermined or programmed time period (1250).
Referring still to FIG. 12, if it is determined that the glucose level is within the predetermined threshold at step 1220, it is then determined whether the glucose rate-of-change is within a predetermined threshold (1230). In certain embodiments, as described above, performing calibration when the glucose level is fluctuating at a low rate-of-change may minimize errors, for example, due to lag between blood glucose and interstitial glucose levels. In one aspect, if it is determined that the glucose rate-of-change is not within the predetermined threshold, then the initiated calibration routine is not contemporaneously executed, but rather, the scheduled calibration function is delayed for a predetermined or programmed time period (1250). On the other hand, if the rate-of-change is within the threshold then the calibration routine is executed or performed to completion as described above at step 1240. In other embodiments, the calibration routine may be executed if only one of the glucose level and the rate-of-change of the glucose level are within the corresponding threshold levels.
FIG. 13 is a flowchart illustrating calibration accuracy improvement routine in another aspect of the present disclosure. Referring to FIG. 13, in one aspect of the present disclosure, the calibration routine may be configured to notify or inform the closed-loop control process or algorithm that calibration is required (or soon to be required) (1310). It should be noted that calibration routine may also be requested or initiated by the patient or the caregiver (e.g., health care provider (HCP)). Upon detection or determination of an impending calibration start event (1310), whether by user initiation or automatic initiation (i.e. at a predetermined time interval or in response to an event), the closed-loop control routine may be configured to modify the glucose control target to a higher value (1320) (balancing with a value that may be too high as to be detrimental to the patient).
Referring again to FIG. 13, once the calibration start event is detected (1330) the calibration routine, using the modified glucose control target, may be configured to determine if the current glucose level is within a target threshold, such as the target set by the modified glucose control target (1340) and only request or execute the calibration function (1350) if the glucose level is within the target threshold. If it is determined that the glucose level is not within the predetermined threshold, then the initiated calibration routine is not contemporaneously executed, but rather, the calibration function is delayed for a predetermined or programmed time period (1370). At this point, in certain embodiments, the routing may wait a predetermined amount of time and then the routine is restarted. Once the calibration function is executed (1350), the glucose control target may be reset back to normal glucose control target settings (1360).
In addition, the closed-loop control routine in one aspect may be configured to switch to a control target of maintaining a low rate of change of glucose, where the control target may be configured to incorporate the desired glucose threshold or range.
In one embodiment, control algorithm may be programmed or configured to maintain multiple control targets for optimal calibration glucose profile and euglycemic management. In one aspect, euglycemic management is configured as a higher priority over optimal calibration profile for the safety of the patient, in the control algorithm.
In the case where a model-based control algorithm is implemented, a vector of state estimates x(t) are provided that accounts for plasma insulin, plasma glucose, and other relevant states, the state observer may be realized in the form of a Kalman Filter or other types of state observers, and configured to use the analyte sensor data as its source of measurement, in addition to the insulin delivery or dosing information. One example of a model-based control algorithm includes a Linear Quadratic (LQ) controller, where the objective function governs the tradeoff between minimizing tracking error and maximizing control effort efficiency. Then, the relative weights under normal operation and when calibration is near can be appropriately adjusted or modified.
For example, consider the following truth model:




  
    
      
        
          
            
              
                
                  
                    I
                    .
                  
                  1
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    -
                    
                      k
                      a
                    
                  
                  ⁢
                  
                    
                      I
                      1
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  
                    u
                    sc
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
              
            
          
        
        
          
            
              
                
                  
                    I
                    .
                  
                  2
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    -
                    
                      k
                      a
                    
                  
                  ⁢
                  
                    
                      I
                      2
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  
                    k
                    a
                  
                  ⁢
                  
                    
                      I
                      1
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
              
            
          
        
        
          
            
              
                
                  I
                  .
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    -
                    
                      k
                      e
                    
                  
                  ⁢
                  
                    I
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  
                    
                      k
                      a
                    
                    V
                  
                  ⁢
                  
                    
                      I
                      2
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
              
            
          
        
        
          
            
              
                
                  
                    r
                    .
                  
                  1
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    -
                    
                      k
                      
                        b
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                    
                  
                  ⁢
                  
                    
                      r
                      1
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  
                    k
                    
                      a
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                  
                  ⁢
                  
                    I
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
              
            
          
        
        
          
            
              
                
                  
                    r
                    .
                  
                  2
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    -
                    
                      k
                      
                        b
                        ⁢
                        
                          
                        
                        ⁢
                        2
                      
                    
                  
                  ⁢
                  
                    
                      r
                      2
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  
                    k
                    
                      a
                      ⁢
                      
                        
                      
                      ⁢
                      2
                    
                  
                  ⁢
                  
                    I
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
              
            
          
        
        
          
            
              
                
                  
                    r
                    .
                  
                  3
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    -
                    
                      k
                      
                        b
                        ⁢
                        
                          
                        
                        ⁢
                        3
                      
                    
                  
                  ⁢
                  
                    
                      r
                      3
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  
                    k
                    
                      a
                      ⁢
                      
                        
                      
                      ⁢
                      3
                    
                  
                  ⁢
                  
                    I
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
              
            
          
        
        
          
            
              
                
                  
                    
                      
                        
                          g
                          .
                        
                        b
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                    =
                    
                    ⁢
                    
                      
                        
                          -
                          
                            [
                            
                              
                                
                                  r
                                  1
                                
                                ⁡
                                
                                  (
                                  t
                                  )
                                
                              
                              +
                              
                                k
                                31
                              
                            
                            ]
                          
                        
                        ⁢
                        
                          
                            g
                            b
                          
                          ⁡
                          
                            (
                            t
                            )
                          
                        
                      
                      -
                    
                  
                
              
              
                
                  
                    
                    ⁢
                    
                      
                        F
                        R
                      
                      +
                      
                        
                          k
                          12
                        
                        ⁢
                        
                          
                            g
                            i
                          
                          ⁡
                          
                            (
                            t
                            )
                          
                        
                      
                      +
                      
                        
                          k
                          13
                        
                        ⁢
                        
                          
                            g
                            2
                          
                          ⁡
                          
                            (
                            t
                            )
                          
                        
                      
                      +
                      
                        EGP
                        ⁡
                        
                          (
                          
                            r
                            3
                          
                          )
                        
                      
                      +
                      
                        
                          g
                          m
                        
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                    
                  
                
              
            
          
        
        
          
            
              
                
                  
                    g
                    .
                  
                  2
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    -
                    
                      [
                      
                        
                          
                            r
                            2
                          
                          ⁡
                          
                            (
                            t
                            )
                          
                        
                        +
                        
                          k
                          13
                        
                      
                      ]
                    
                  
                  ⁢
                  
                    g
                    2
                  
                
                +
                
                  
                    
                      r
                      1
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                  ⁢
                  
                    
                      g
                      b
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
              
            
          
        
        
          
            
              
                
                  
                    g
                    .
                  
                  i
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    -
                    
                      k
                      02
                    
                  
                  ⁢
                  
                    
                      g
                      i
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  
                    [
                    
                      
                        k
                        21
                      
                      +
                      
                        [
                        
                          
                            k
                            i
                          
                          ⁡
                          
                            [
                            
                              
                                I
                                ⁡
                                
                                  (
                                  t
                                  )
                                
                              
                              -
                              
                                I
                                b
                              
                            
                            ]
                          
                        
                        ]
                      
                    
                    ]
                  
                  ⁢
                  
                    
                      g
                      b
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                  ⁢
                  
                    F
                    02
                  
                
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




where, in addition to Equations 1 and 2 above, other glucose compartments gb and g2 as well as effective insulin compartments r1, r2, and r3 have been included. In the case where the model for the control algorithm is configured to perform a local linearization at every time step:





  
    
      
        
          
            
              
                
                  x
                  .
                
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    A
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                  ⁢
                  
                    x
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                +
                
                  B
                  ⁢
                  
                    
                  
                  ⁢
                  
                    u
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
              
            
          
        
        
          
            
              
                u
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  u
                  sc
                
                ⁡
                
                  (
                  t
                  )
                
              
            
          
        
        
          
            
              
                y
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                
                  
                    y
                    CGM
                  
                  ⁡
                  
                    (
                    t
                    )
                  
                
                =
                
                  S
                  ⁡
                  
                    [
                    
                      
                        
                          g
                          i
                        
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                      +
                      
                        
                          v
                          i
                        
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                    
                    ]
                  
                
              
            
          
        
        
          
            
              
                x
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
                [
                
                  
                    
                      
                        
                          I
                          1
                        
                        -
                        
                          I
                          
                            1
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          I
                          2
                        
                        -
                        
                          I
                          
                            2
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          I
                          3
                        
                        -
                        
                          I
                          
                            3
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          r
                          1
                        
                        -
                        
                          r
                          
                            1
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          r
                          2
                        
                        -
                        
                          r
                          
                            2
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          r
                          3
                        
                        -
                        
                          r
                          
                            3
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          g
                          b
                        
                        -
                        
                          g
                          
                            b
                            ⁢
                            
                              
                            
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          g
                          2
                        
                        -
                        
                          g
                          
                            2
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          g
                          i
                        
                        -
                        
                          g
                          
                            i
                            ⁢
                            
                              
                            
                            ⁢
                            t
                          
                        
                      
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




It is to be noted that the states have been defined as the difference between the physiologically meaningful states of the truth model and their corresponding targets.

Further, an LQ optimal control is determined such that the objective function J is minimized:




  
    
      
        
          J
          =
          
            
              ∫
              t
              
                t
                +
                
                  t
                  p
                
              
            
            ⁢
            
              
                [
                
                  
                    [
                    
                      
                        
                          x
                          T
                        
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                      ⁢
                      Q
                      ⁢
                      
                        
                      
                      ⁢
                      
                        x
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                    
                    ]
                  
                  +
                  
                    [
                    
                      
                        
                          u
                          T
                        
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                      ⁢
                      R
                      ⁢
                      
                        
                      
                      ⁢
                      
                        u
                        ⁡
                        
                          (
                          t
                          )
                        
                      
                    
                    ]
                  
                
                ]
              
              ⁢
              
                d
                t
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            Q
            =
            
              [
              
                
                  
                    
                      q
                      
                        1
                        ,
                        1
                      
                    
                  
                  
                    ...
                  
                  
                    
                      q
                      
                        1
                        ,
                        9
                      
                    
                  
                
                
                  
                    ⋮
                  
                  
                    ⋱
                  
                  
                    ⋮
                  
                
                
                  
                    
                      q
                      
                        9
                        ,
                        1
                      
                    
                  
                  
                    ...
                  
                  
                    
                      q
                      
                        9
                        ,
                        9
                      
                    
                  
                
              
              ]
            
          
          ,
          
            
          
          ⁢
          
            R
            =
            
              [
              
                r
                sc
              
              ]
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




where tp is a finite future horizon in which the controller must be optimized for, Q is a positive semidefinite matrix that penalizes linear combinations of the states x, and R is a positive definite matrix that penalizes the control action.

In particular, the distinction between controlling for optimal calibration and controlling for optimal glucose regulation, using this LQ framework as an example, is described below. In the case of controlling for optimal glucose regulation, for a given desired strict plasma glucose target of 100 mg/dL, the quantity gbt is set to 100 mg/dL, so that when the objective function in Equation 8 is evaluated, any deviation of gb from this value will contribute to an increase in J.
If other states do not need to be regulated at any specific level, then the corresponding targets I1t, I2t, and so on, can be set to any arbitrary real value (such as zero), and Q must be tuned such that only q7,7 (which corresponds to the penalty for gb) be left nonzero. The relative magnitude between q7,7 and rsc then determines aggressive target tracking and conservative control action.
In the case of controlling for optimal calibration, a combination of strict plasma glucose target and zero glucose rate is obtained, which, in one aspect may be approximated by setting the rate of change of the glucose rates to zero. As a result, the corresponding targets for the glucose compartments can be estimated as follows:




  
    
      
        
          [
          
            
              
                
                  g
                  
                    b
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                
              
            
            
              
                
                  g
                  
                    2
                    ⁢
                    t
                  
                
              
            
            
              
                
                  g
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                
              
            
          
          ]
        
        =
        
          
            inv
            ⁡
            
              (
              
                [
                
                  
                    
                      
                        -
                        
                          [
                          
                            
                              r
                              1
                            
                            +
                            
                              k
                              31
                            
                          
                          ]
                        
                      
                    
                    
                      
                        k
                        13
                      
                    
                    
                      
                        k
                        12
                      
                    
                  
                  
                    
                      
                        r
                        1
                      
                    
                    
                      
                        -
                        
                          [
                          
                            
                              r
                              2
                            
                            +
                            
                              k
                              13
                            
                          
                          ]
                        
                      
                    
                    
                      0
                    
                  
                  
                    
                      
                        
                          k
                          21
                        
                        +
                        
                          [
                          
                            
                              k
                              i
                            
                            ⁡
                            
                              [
                              
                                I
                                -
                                
                                  I
                                  b
                                
                              
                              ]
                            
                          
                          ]
                        
                      
                    
                    
                      0
                    
                    
                      
                        k
                        02
                      
                    
                  
                
                ]
              
              )
            
          
          ⁡
          
            [
            
              
                
                  
                    
                      F
                      R
                    
                    -
                    
                      EGP
                      ⁡
                      
                        (
                        
                          r
                          3
                        
                        )
                      
                    
                    -
                    
                      g
                      m
                    
                  
                
              
              
                
                  0
                
              
              
                
                  
                    F
                    02
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        9
        )
      
    
  




The above targets can be assigned to the glucose compartments, and as in the optimal glucose regulation case, other targets can be set to zero. The proper state weighting matrix Q must be set such that the glucose states track the established targets.

If calibration favors not only steady glucose but also a particular blood glucose value, then the target for blood glucose may be set explicitly (e.g. gbt=100 mg/dL), and the other glucose targets can be derived such that the following is satisfied:




  
    
      
        
          
            [
            
              
                
                  
                    k
                    13
                  
                
                
                  
                    k
                    12
                  
                
              
              
                
                  
                    -
                    
                      [
                      
                        
                          r
                          2
                        
                        +
                        
                          k
                          13
                        
                      
                      ]
                    
                  
                
                
                  0
                
              
              
                
                  0
                
                
                  
                    -
                    
                      k
                      02
                    
                  
                
              
            
            ]
          
          ⁡
          
            [
            
              
                
                  
                    g
                    2
                  
                
              
              
                
                  
                    g
                    i
                  
                
              
            
            ]
          
        
        =
        
          [
          
            
              
                
                  
                    F
                    R
                  
                  -
                  
                    EGP
                    ⁡
                    
                      (
                      
                        r
                        3
                      
                      )
                    
                  
                  -
                  
                    g
                    m
                  
                  +
                  
                    
                      [
                      
                        
                          r
                          1
                        
                        +
                        
                          k
                          31
                        
                      
                      ]
                    
                    ⁢
                    
                      g
                      
                        b
                        ⁢
                        
                          
                        
                        ⁢
                        t
                      
                    
                  
                
              
            
            
              
                
                  
                    -
                    
                      r
                      1
                    
                  
                  ⁢
                  
                    g
                    
                      b
                      ⁢
                      
                        
                      
                      ⁢
                      t
                    
                  
                
              
            
            
              
                
                  
                    F
                    02
                  
                  -
                  
                    
                      [
                      
                        
                          k
                          21
                        
                        +
                        
                          [
                          
                            
                              k
                              i
                            
                            ⁡
                            
                              [
                              
                                I
                                -
                                
                                  I
                                  b
                                
                              
                              ]
                            
                          
                          ]
                        
                      
                      ]
                    
                    ⁢
                    
                      g
                      
                        b
                        ⁢
                        
                          
                        
                        ⁢
                        t
                      
                    
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        10
        )
      
    
  




The targets for g2 and gi can then be computed using the least-squares error approximation shown:





  
    
      
        
          [
          
            
              
                
                  g
                  
                    2
                    ⁢
                    t
                  
                
              
            
            
              
                
                  g
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                
              
            
          
          ]
        
        =
        
          
            inv
            ⁡
            
              (
              
                
                  [
                  
                    
                      
                        
                          k
                          13
                        
                      
                      
                        
                          -
                          
                            [
                            
                              
                                r
                                2
                              
                              +
                              
                                k
                                13
                              
                            
                            ]
                          
                        
                      
                      
                        0
                      
                    
                    
                      
                        
                          k
                          12
                        
                      
                      
                        0
                      
                      
                        
                          -
                          
                            k
                            02
                          
                        
                      
                    
                  
                  ]
                
                ⁡
                
                  [
                  
                    
                      
                        
                          k
                          13
                        
                      
                      
                        
                          k
                          12
                        
                      
                    
                    
                      
                        
                          -
                          
                            [
                            
                              
                                r
                                2
                              
                              +
                              
                                k
                                13
                              
                            
                            ]
                          
                        
                      
                      
                        0
                      
                    
                    
                      
                        0
                      
                      
                        
                          -
                          
                            k
                            02
                          
                        
                      
                    
                  
                  ]
                
              
              )
            
          
          ⁢
          
             
            
              
                [
                
                  
                    
                      
                        k
                        13
                      
                    
                    
                      
                        -
                        
                          [
                          
                            
                              r
                              2
                            
                            +
                            
                              k
                              13
                            
                          
                          ]
                        
                      
                    
                    
                      0
                    
                  
                  
                    
                      
                        k
                        12
                      
                    
                    
                      0
                    
                    
                      
                        -
                        
                          k
                          02
                        
                      
                    
                  
                
                ]
              
              ⁡
              
                [
                
                  
                    
                      
                        
                          F
                          R
                        
                        -
                        
                          EGP
                          ⁡
                          
                            (
                            
                              r
                              3
                            
                            )
                          
                        
                        -
                        
                          g
                          m
                        
                        +
                        
                          
                            [
                            
                              
                                r
                                1
                              
                              +
                              
                                k
                                31
                              
                            
                            ]
                          
                          ⁢
                          
                            g
                            
                              b
                              ⁢
                              
                                
                              
                              ⁢
                              t
                            
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          -
                          
                            r
                            1
                          
                        
                        ⁢
                        
                          g
                          
                            b
                            ⁢
                            
                              
                            
                            ⁢
                            t
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          F
                          02
                        
                        -
                        
                          
                            [
                            
                              
                                k
                                21
                              
                              +
                              
                                [
                                
                                  
                                    k
                                    i
                                  
                                  ⁡
                                  
                                    [
                                    
                                      I
                                      -
                                      
                                        I
                                        b
                                      
                                    
                                    ]
                                  
                                
                                ]
                              
                            
                            ]
                          
                          ⁢
                          
                            g
                            
                              b
                              ⁢
                              
                                
                              
                              ⁢
                              t
                            
                          
                        
                      
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




In the manner described above, in accordance with aspects of the present disclosure, one or more parameters or information of events that may impact the level of blood glucose or glucose measurements, if available during the analyte sensor calibration process, may be factored in to improve the sensor calibration accuracy, for example, by improving the accuracy of the sensor sensitivity determination. Events or conditions referred to herein include, but not limited to exercise information, meal intake information, patient health information, medication information, disease information, physiological profile information, and insulin delivery information. While the various embodiments described above in conjunction with the improvement of the sensor calibration accuracy include insulin delivery information, within the scope of the present disclosure, any exogenous information that are available to the and during the calibration process or routine that may have an impact on the level of glucose may be considered.
In one aspect, the user or the patient may provide this information into one or more components of the system 100 (FIG. 1) which includes a user interface for entering events and/or data. Alternatively, this information may be entered manually into another device and transferred electronically to the processor(s) performing the calibration process/routine. Finally, this information may be recorded by either the device(s) that perform the calibration process/routine, or by a separate device that transfers the information electronically to the device(s) that perform the calibration process/routine.
In one embodiment, the medication delivery device is configured to deliver appropriate medication based on one or more delivery profiles stored therein, and in addition, configured to record the amount of medication delivered with delivery time association in an electronic log or database. The medication delivery device may be configured to periodically (automatically, or in response to one or more commands from the controller/another device) transfer medication delivery data/information to the controller (or another device) electronic log or database. In this manner, the analyte monitoring device including the receiver/controller unit may be provided with software programming that can be executed to perform the sensor calibration routine and provided with access to all relevant information received from the medication delivery unit, the analyte sensor/transmitter, user input information, as well as previously stored information.
In this manner, in one aspect of the present disclosure, the accuracy of the sensor sensitivity determination may be improved based on the insulin delivery information which provides additional data to determine or anticipate future glucose values, and may help to compensate for potential error in the sensor readings or measurements due to lag, in particular, when the level of glucose is undergoing a rapid fluctuation. In addition, the insulin information may be used to adjust or determine the suitable or appropriate time to perform the sensor calibration routine. For example, this information may be used to determine or anticipate periods of high rates of glucose change which would not be an ideal condition for determining sensor sensitivity for performing sensor calibration.
Within the scope of the present disclosure, the programming, instructions or software for performing the calibration routine, user interaction, data processing and/or communication may be incorporated in the analyte monitoring device, the medication delivery device, the control unit, or any other component of the overall system 100 shown in FIG. 1, and further, may also be provided in multiple devices or components to provide redundancy. Additionally, embodiments described herein may also be integrated in a closed loop control system which is programmed to control insulin delivery so as to provide, in part, conditions that are suitable for performing sensor calibration in the closed loop control system.
In one embodiment, a method may include detecting an analyte sensor calibration start event, determining one or more parameters associated with a calibration routine corresponding to the detected calibration start event, and executing the calibration routine based on the one or more determined parameters, wherein the one or more determined parameters includes a medication delivery information.
Detecting the calibration start event may include monitoring an elapsed time period from initial analyte sensor placement.
Detecting the calibration start event may be based at least in part on a predetermined schedule.
The predetermined schedule may include approximately once every twenty four hours.
The determined one or more parameters may include an amount of insulin dose delivered, a time period of the delivered insulin dose, an insulin sensitivity parameter, an insulin on board information, an exercise information, a meal intake information, an activity information, or one or more combinations thereof.
The medication delivery information may include an insulin delivery amount and time information relative to the detected calibration start event.
Executing the calibration routine may include delaying the calibration routine by a predetermined time period.
The predetermined time period may include approximately 1-2 hours.
The calibration routine may not be executed when one of the one or more determined parameters deviates from a predetermined threshold level.
The predetermined threshold level may be dynamically modified based on a variation in the corresponding one or more determined parameters.
The predetermined threshold level may be user defined.
Executing the calibration routine may include determining a reference measurement value.
Determining the reference measurement value may include prompting for a blood glucose measurement, and receiving data corresponding to the measured blood glucose level.
Executing the calibration routine may include determining a sensitivity value associated with the analyte sensor.
Executing the calibration routine may include calibrating the analyte sensor.
In another embodiment, a device may include one or more processors and a memory operatively coupled to the one or more processors, the memory for storing instructions which, when executed by the one or more processors, causes the one or more processors to detect an analyte sensor calibration start event, to determine one or more parameters associated with a calibration routine corresponding to the detected calibration start event, and to execute the calibration routine based on the one or more determined parameters, wherein the one or more determined parameters includes a medication delivery information.
The analyte sensor may include a glucose sensor.
The medication delivery information may include information associated with insulin dose administered.
Furthermore, an output unit may be operatively coupled to the one or more processors for outputting one or more data or signals associated with the calibration start event or the calibration routine.
In yet another embodiment, a method may include initializing an analyte sensor, receiving a data stream from the initialized analyte sensor, detecting a calibration start event associated with the initialized analyte sensor, determining one or more parameters associated with insulin dose administration, and executing a calibration routine based on the one or more determined parameters.
In yet another embodiment, a method may include detecting an impending glucose sensor calibration start event, modifying a medical treatment profile to a higher than average target glucose level upon detection of the impending glucose sensor calibration start event, determining one or more parameters associated with a calibration routine corresponding to the detected impending calibration start event, wherein the one or more determined parameters includes a current glucose level, executing the calibration routine based on the one or more determined parameters, and resetting the medical treatment profile to an average target glucose level.
The calibration routine may be executed only if the current glucose level is above a predetermined threshold.
The predetermined threshold may be higher than the average glucose level.
In one aspect, the method may include delaying execution of the calibration routine until the current glucose level is above the predetermined threshold.
In another aspect, the method may include outputting one or more data or signals associated with the calibration routine.
The medical treatment profile may include insulin dose administration information.
Various other modifications and alterations in the structure and method of operation of this disclosure will be apparent to those skilled in the art without departing from the scope and spirit of the present disclosure. Although the present disclosure has been described in connection with specific embodiments, it should be understood that the present disclosure as claimed should not be unduly limited to such specific embodiments. It is intended that the following claims define the scope of the present disclosure and that structures and methods within the scope of these claims and their equivalents be covered thereby.

## Claims

1. A method, comprising:
detecting, using one or more processors, an analyte sensor calibration start event;
determining, using the one or more processors, one or more parameters associated with a calibration routine corresponding to the detected calibration start event; and
executing, using the one or more processors, the calibration routine based on the one or more determined parameters;
wherein the one or more determined parameters include a medication delivery information.

2. The method of claim 1 wherein detecting the calibration start event includes monitoring an elapsed time period from initial analyte sensor placement.

3. The method of claim 1 wherein detecting the calibration start event is based at least in part on a predetermined schedule.

4. The method of claim 3 wherein the predetermined schedule includes approximately once every twenty four hours.

5. The method of claim 1 wherein the determined one or more parameters include an amount of insulin dose delivered, a time period of the delivered insulin dose, an insulin sensitivity parameter, an insulin on board information, an exercise information, a meal intake information, an activity information, or one or more combinations thereof.

6. The method of claim 1 wherein the medication delivery information includes an insulin delivery amount and time information relative to the detected calibration start event.

7. The method of claim 1 wherein executing the calibration routine includes delaying the calibration routine by a predetermined time period.

8. The method of claim 7 wherein the predetermined time period includes approximately 1-2 hours.

9. The method of claim 1 wherein the calibration routine is not executed when one of the one or more determined parameters deviates from a predetermined threshold level.

10. The method of claim 9 wherein the predetermined threshold level is dynamically modified based on a variation in the one or more determined parameters.

11. The method of claim 9 wherein the predetermined threshold level is user defined.

12. The method of claim 1 wherein executing the calibration routine includes determining a reference measurement value.

13. The method of claim 12 wherein determining the reference measurement value includes:
prompting for a blood glucose measurement; and
receiving data corresponding to a measured blood glucose level.

14. The method of claim 1 wherein executing the calibration routine includes determining a sensitivity value associated with an analyte sensor.

15. The method of claim 1 wherein executing the calibration routine includes calibrating an analyte sensor.

16. A device, comprising:
one or more processors; and
a memory operatively coupled to the one or more processors, the memory for storing instructions which, when executed by the one or more processors, causes the one or more processors to detect an analyte sensor calibration start event, to determine one or more parameters associated with a calibration routine corresponding to the detected calibration start event, and to execute the calibration routine based on the one or more determined parameters, wherein the one or more determined parameters include a medication delivery information.

17. The device of claim 16 wherein the analyte sensor calibration start event includes a glucose sensor calibration start event.

18. The device of claim 16 wherein the medication delivery information includes information associated with an administered insulin dose.

19. The device of claim 16 further including an output unit operatively coupled to the one or more processors for outputting one or more data or signals associated with the calibration start event or the calibration routine.

20. A method, comprising:
initializing, using one or more processors, an analyte sensor;
receiving, using the one or more processors, a data stream from the initialized analyte sensor;
detecting, using the one or more processors, a calibration start event associated with the initialized analyte sensor;
determining, using the one or more processors, one or more parameters associated with insulin dose administration; and
executing, using the one or more processors, a calibration routine based on the one or more determined parameters.

21. A method, comprising:
detecting, using one or more processors, an impending glucose sensor calibration start event;
modifying, using the one or more processors, a medical treatment profile to a higher than average target glucose level upon detection of the impending glucose sensor calibration start event;
determining, using the one or more processors, one or more parameters associated with a calibration routine corresponding to the detected impending calibration start event, wherein the one or more determined parameters include a current glucose level;
executing, using the one or more processors, the calibration routine based on the one or more determined parameters; and
resetting, using the one or more processors, the medical treatment profile to an average target glucose level.

22. The method of claim 21 wherein the calibration routine is executed only if the current glucose level is above a predetermined threshold.

23. The method of claim 22 wherein the predetermined threshold is higher than the average glucose level.

24. The method of claim 22 further comprising delaying execution of the calibration routine until the current glucose level is above the predetermined threshold.

25. The method of claim 21 further comprising outputting one or more data or signals associated with the calibration routine.

26. The method of claim 21 wherein the medical treatment profile includes insulin dose administration information.