WO2002024065A1 - Method and apparatus for real-time estimation and control of pysiological parameters - Google Patents

Method and apparatus for real-time estimation and control of pysiological parameters Download PDF

Info

Publication number
WO2002024065A1
WO2002024065A1 PCT/US2001/029649 US0129649W WO0224065A1 WO 2002024065 A1 WO2002024065 A1 WO 2002024065A1 US 0129649 W US0129649 W US 0129649W WO 0224065 A1 WO0224065 A1 WO 0224065A1
Authority
WO
WIPO (PCT)
Prior art keywords
glucose
estimator
time
real
measurement
Prior art date
Application number
PCT/US2001/029649
Other languages
French (fr)
Inventor
Edward J. Knobbe
Wah L. Lim
Bruce A. Buckingham
Original Assignee
Knobbe, Lim & Buckingham
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Knobbe, Lim & Buckingham filed Critical Knobbe, Lim & Buckingham
Priority to AU2001291189A priority Critical patent/AU2001291189A1/en
Publication of WO2002024065A1 publication Critical patent/WO2002024065A1/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0002Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14532Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/40ICT specially adapted for the handling or processing of patient-related medical or healthcare data for data related to laboratory analysis, e.g. patient specimen analysis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
    • G16H20/13ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients delivered from dispensers
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2560/00Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
    • A61B2560/02Operational features
    • A61B2560/0242Operational features adapted to measure environmental factors, e.g. temperature, pollution
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M2230/00Measuring parameters of the user
    • A61M2230/20Blood composition characteristics
    • A61M2230/201Glucose concentration
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M5/00Devices for bringing media into the body in a subcutaneous, intra-vascular or intramuscular way; Accessories therefor, e.g. filling or cleaning devices, arm-rests
    • A61M5/14Infusion devices, e.g. infusing by gravity; Blood infusion; Accessories therefor
    • A61M5/142Pressure infusion, e.g. using pumps
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/30ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to physical therapies or activities, e.g. physiotherapy, acupressure or exercising
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S128/00Surgery
    • Y10S128/92Computer assisted medical diagnostics
    • Y10S128/923Computer assisted medical diagnostics by comparison of patient data to other data
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S128/00Surgery
    • Y10S128/92Computer assisted medical diagnostics
    • Y10S128/923Computer assisted medical diagnostics by comparison of patient data to other data
    • Y10S128/924Computer assisted medical diagnostics by comparison of patient data to other data using artificial intelligence

Definitions

  • This present invention relates generally to a method and apparatus for estimating and controlling physiological parameters and more particularly to an optimal estimator and an optimal controller for estimating and controlling glucose levels in a patient.
  • Glucose monitoring is typically performed by people with diabetes mellitus which is a medical condition involving a body's inability to produce the quantity or quality of insulin needed to maintain a normal circulating blood glucose. Frequent monitoring of glucose is generally necessary to provide effective treatment and to prevent long term complications of diabetes (e.g., blindness, kidney failure, heart failure, etc.). New methods of monitoring glucose are fast, painless and convenient alternatives to the typical capillary blood glucose (CBG) measurements which involve finger pricks that are painful, inconvenient and difficult to perform for long term.
  • CBG capillary blood glucose
  • Optical measurement of glucose is performed by focusing a beam of light onto the body.
  • Optical sensors determine glucose concentration by analyzing optical signal changes in wavelength, polarization or intensity of light.
  • sensor characteristics e.g., aging
  • environmental variations e.g., changes in temperature, humidity, skin hydration, pH, etc.
  • physiological variations e.g., changes in tissue fluid due to activity, diet, medication or hormone fluctuations
  • One method utilizes calibration models developed by initially measuring known glucose concentrations to correct subsequent sensor measurements. The calibration models become inaccurate over time due to dynamic changes in physiological processes.
  • Another method e.g., adaptive noise canceling
  • the present invention solves these and other problems by providing a method and apparatus for making optimal estimates of a physiological parameter (e.g., glucose level), assessing reliability of the optimal estimates, and/or providing optimal control of the physiological parameter in real time using one or more sensor measurements at each measurement time epoch (or interval).
  • the sensor measurements can be time-based (e.g., every five minutes) to provide continuous monitoring and/or regulation of the physiological parameter.
  • the sensor measurements are a function of the physiological parameter within specified uncertainties.
  • An optimal estimator provides an accurate estimate of glucose level in real time using a sensor with at least one output.
  • the optimal estimator is integrated with the sensor and an output display to be a compact glucose monitoring device which can be worn by a patient for continuous monitoring and real-time display of glucose level.
  • the optimal estimator is a separate unit which can interface with different types of sensors and provide one or more outputs for display, further processing by another device, or storage on a memory device.
  • the optimal estimator employs a priori deterministic dynamic models developed with stochastic variables and uncertain parameters to make estimates of glucose level.
  • glucose level is defined as one of the stochastic (or random) variables.
  • Dynamic mathematical models define process propagation (i.e., how physiological and sensor parameters change in time) and measurement relationship (i.e., how physiological and sensor parameters relate to environmental conditions).
  • Environmental conditions e.g., temperature, humidity, pH, patient activity, etc.
  • the optimal estimator uses dynamic models to propagate estimates of respective stochastic variables, error variances, and error covariances forward in time.
  • the optimal estimator At each measurement time epoch, the optimal estimator generates real-time estimates of the stochastic variables using one or more sensor outputs and any ancillary input related to environmental conditions.
  • the optimal estimator employs a linearized Kalman filter to perform optimal estimation of the stochastic variables (e.g., glucose level).
  • an extended Kalman filter is used to accommodate nonlinear stochastic models.
  • the optimal estimator is initialized by providing initial values for the stochastic variables, error variances, and error covariances. For example, a CBG measurement or another direct glucose measurement is performed at initialization to provide a starting value for the stochastic variable corresponding to glucose level.
  • the optimal estimator provides one or more outputs to a patient health monitor which is capable of optimized real-time decisions and displays.
  • the patient health monitor evaluates system performance by assessing the performance of the sensor and/or optimal estimator in real time. For example, the patient health monitor applies statistical testing to determine the reliability of the real-time estimates of the stochastic variables by the optimal estimator. The statistical testing is performed in real time on residual errors of the optimal estimator to establish performance measures.
  • the patient health monitor acts as an input/output interface between the patient or medical staff (e.g., a doctor, nurse, or other healthcare provider) and the optimal estimator.
  • the patient or medical staff e.g., a doctor, nurse, or other healthcare provider
  • environmental conditions can be provided to the patient health monitor for forwarding to the optimal estimator.
  • Optimal estimator outputs can be provided to the patient health monitor for display or forwarding to an external device
  • a computer e.g., a computer or a data storage device.
  • the optimal estimator provides one or more outputs to an optimal controller which can regulate in real time the physiological parameter being monitored.
  • an optimal controller responds to realtime optimal estimator outputs and provides an output to operate an actuator.
  • the actuator can be a dispenser or a pump which secretes insulin to correct a relatively high glucose level and glucagon to correct a relatively low glucose level.
  • the optimal controller takes advantage of a priori information regarding the statistical characteristics of the actuator and is able to control the output of the actuator to be within specified uncertainties.
  • the optimal estimator and the optimal controller form an optimal closed-loop system.
  • a glucose sensor, an optimal estimator, an optimal controller, and an insulin/glucagon dispenser work together as an artificial pancreas to continuously regulate glucose level.
  • the glucose sensor can be internal or external to a patient's body.
  • the optimal controller provides a control feedback to the optimal estimator to account for delivery of the insulin/glucagon.
  • the optimal closed-loop system is effective in a variety of biomedical applications.
  • cardiovascular functions can be continuously regulated by using sensors to detect blood pressure, blood oxygen level, physical activity and the like, an optimal estimator to process the sensor measurements and make real-time estimates of heart function parameters, and an optimal controller to control operations of an artificial device (e.g., a pacemaker) in real time based on the real-time estimates from the optimal estimator to achieve a set of desired heart function parameter values.
  • an artificial device e.g., a pacemaker
  • Other artificial devices e.g., artificial limbs, bionic ears, and bionic eyes
  • the optimal closed-loop system is also effective in optimal treatment of chronic illnesses (e.g., HIV). Some medications for treatment of chronic illnesses are relatively toxic to the body. Over delivery of medication generally has adverse effects on the patient.
  • the optimal closed-loop system is capable of providing effective and safe treatment for the patient. For example, an optimal estimator provides real-time estimates of key physiological parameters using one or more sensors, and an optimal controller controls a slow infusion of medication in real time based on the real-time estimates from the optimal estimator to obtain desirable values for the key physiological parameters.
  • the optimal estimator, patient health monitor, and optimal controller are software algorithms which can be implemented using respective microprocessors. New information regarding process propagation or measurement relationship can be easily incorporated by modifying, reconfiguring, and/or adding to the software algorithms.
  • the optimal estimator, patient health monitor, and optimal controller can be implemented as one joint algorithm or separate respective algorithms which function together to provide an optimal closed-loop system.
  • Figure 1 A is a block diagram of one embodiment of an estimator.
  • Figure 1 B is a block diagram of another embodiment of a glucose-monitoring device.
  • Figure 2 illustrates one embodiment of modeling physiological processes in a linearized Kalman filter application.
  • Figure 3A describes one embodiment of an estimation function which depicts a linearized Kalman filter formulation.
  • Figure 3B illustrates an initialization in one embodiment of an optimal estimator.
  • Figure 3C illustrates a time-update cycle in accordance with one embodiment of an optimal estimator.
  • Figure 3D illustrates a measurement-update cycle in accordance with one embodiment of an optimal estimator.
  • Figures 4A and 4B illustrate a first set of time history plots of optimal glucose estimates and CBG measurements with respect to time.
  • Figures 5A and 5B illustrate a second set of time history plots of optimal glucose estimates and CBG measurements with respect to time.
  • Figure 6A is a block diagram of one embodiment of a patient health monitor.
  • Figure 6B illustrates one embodiment of a residual test data process in the patient health monitor.
  • Figure 6C illustrates one embodiment of a statistical test process in the patient health monitor.
  • Figure 6D illustrates one embodiment of an input/output interface in the patient health monitor.
  • Figures 7A and 7B illustrate time history plots of residual test data with respect to measurement time.
  • Figure 8 is a functional diagram of one embodiment of an artificial pancreas.
  • Figure 9A illustrates one embodiment of an optimal controller for a closed loop system.
  • Figure 9B illustrates one embodiment of a control model for a controller in accordance with the present invention.
  • Figure 9C illustrates one embodiment of a control algorithm.
  • the present invention involves application of real-time optimal estimation, optimized real-time decisions and displays (e.g., a patient health monitor), or optimal real-time control to physiological processes.
  • the real-time optimal estimation, the optimized real-time decisions and displays, and the optimal real-time control are implemented as separate modules which can be combined functionally.
  • the real-time optimal estimation, the optimized real-time decisions and displays, and the optimal real-time control are implemented as one joint algorithm.
  • input to an optimal estimator is provided by a physical sensor (or a plurality of sensors) which measures some arbitrary, but known, function (or functions) of variables or parameters to be estimated to within specified uncertainties and whose statistical characteristics are known.
  • an output of a real-time controller is provided to a physical controllable dispenser, or actuator, whose output is some known function of parameters to be controlled within specified uncertainties and whose statistical characteristics are known.
  • a decision and display function utilizes statistical testing of estimator residual errors using internally computed, and updated, estimator variances and covariances.
  • the present invention is implemented as a software algorithm.
  • the present invention uses models (e.g., dynamic process and measurement models). For best performance, the models should reflect the latest and most complete information available. As new and more complete information is developed, performance can be improved through incorporation of this information by simply modifying the software algorithm of the present invention.
  • Figure 1A is a block diagram of one embodiment of an estimator.
  • the estimator uses a linearized Kalman filter.
  • the linearized Kalman filter accommodates nonlinear process models and/or nonlinear measurement models.
  • the linearized Kalman filter is a discrete extended Kalman filter which is linearized after each update using best estimates.
  • a general formulation of a continuous-discrete extended Kalman filter is provided in Table 1.
  • the linearized Kalman filter is linearized about a nominal set for which a general formulation is given in
  • the estimator computes an estimator gain based on time updated and measurement updated error variable variances and covariances.
  • the estimator is implemented using discrete formulations.
  • the estimator is implemented using continuous formulations.
  • a Covariance formulation or an Information formulation depending on initialization uncertainty considerations.
  • the use of Bierman factorization techniques (UDITT) in the implementation leads to numerically stable algorithms which are excellent for operation over very long time periods.
  • the estimator of Figure 1A is applied to the problem of monitoring patient glucose levels. Any type of physical sensor which measures some function of glucose can be used.
  • one or more capillary blood glucose (CBG) measurements are obtained on initialization of the estimator and when estimated glucose deviations exceed computed variance levels.
  • CBG capillary blood glucose
  • additional capillary glucose values are obtained on a periodic basis to assess sensor function.
  • CBG measurements may be obtained once or twice a day pre and one hour post prandial or when the estimator determines that the glucose values are out of range of predetermined limits.
  • the estimator is contained in a small portable package.
  • the estimator can be operated by a patient.
  • the estimator can be operated by medical staff in a hospital or clinic.
  • the estimator e.g., Kalman estimator
  • the Kalman estimator can also use additional information to more closely predict time propagation of glucose levels, such as exercise, food intake, insulin administration, or other factors which influence glucose levels or the function of the sensor (e.g., local pH, temperature or oxygen tension near the physical sensor).
  • FIG. 1B is a block diagram of one embodiment of a glucose-monitoring device 100.
  • the glucose-monitoring device 100 includes a glucose sensor 108, an ancillary sensor 110, a glucose estimator 114, and a patient health monitor 114.
  • the glucose-monitoring device 100 is used to provide real-time glucose estimates of a patient 102.
  • the glucose-monitoring device 100 is an integrated unit which is portable by the patient 102 to provide continuous glucose monitoring and real-time displays.
  • the glucose sensor 108 (e.g., a probe, patch, infrared, or laser sensor) is coupled to the patient 102 and outputs a measurement f(g) which is a function of the glucose level of the patient 102.
  • the glucose sensor 108 makes measurements within specified uncertainties and has known statistical characteristics.
  • the glucose sensor 108 provides the measurement f(g) to the glucose estimator 112.
  • the glucose sensor 108 makes measurements periodically.
  • the glucose sensor 108 makes measurements intermittently or upon command.
  • the ancillary sensor 110 is coupled to the patient 102, the glucose sensor 108, and/or surroundings of the patient 102 and/or glucose sensor 108 to provide information regarding environmental and/or glucose sensor conditions (e.g., temperature, humidity, local pH, etc.) which affect the measurement.
  • environmental and/or glucose sensor conditions e.g., temperature, humidity, local pH, etc.
  • the ancillary sensor 110 are provided to the glucose estimator 112.
  • the ancillary sensor 110 provides outputs to the patient health monitor 114 which can process the information for display or for forwarding to the glucose estimator 112.
  • the ancillary sensor 110 is not a part of the glucose-monitoring device 100.
  • the patient 102 and/or a medical staff 104 can provide information on the environmental and glucose sensor conditions as well as other information affecting the measurement.
  • the patient i.e., a user
  • the patient health monitor 114 inputs information (e.g., exercise activity, food intake, insulin administration, etc.) using the patient health monitor 114.
  • the patient health monitor 114 acts as an input/output interface or a means for the user to configure the glucose estimator 112.
  • the patient health monitor 114 forwards the information to the glucose estimator 112.
  • the patient health monitor 114 is a display device.
  • an output of the glucose estimator 1 12 e.g., an optimal real-time estimate of glucose
  • the patient health monitor 114 can also contemporaneously display information provided by the glucose sensor 108, the ancillary sensor 1 10, the patient 102 and/or the medical staff 104 which affects the measurement used to make the optimal glucose estimate.
  • the patient health monitor 1 14 is a status indicator.
  • the glucose estimator 1 12 provides outputs (e.g., residuals and variances) to the patient health monitor 1 14 which applies statistical testing to determine the reliability of the glucose sensor 108 and/or the glucose estimator 1 12.
  • the patient health monitor 114 provides a warning when poor performance is detected.
  • the glucose estimator 112 is a linearized Kalman filter (e.g., a discrete extended Kalman filter) to account for a nonlinear process model and/or measurement model.
  • the glucose estimator (or Kalman estimator) 1 12 provides real-time estimates of the glucose level in the patient 102.
  • An initialization measurement 106 e.g., a CBG measurement
  • Figure 2 is a flow chart illustrating one embodiment of a modeling process for physiological processes.
  • physiological processes are described by nonlinear stochastic models.
  • the flow chart of Figure 2 illustrates one method of developing a dynamic model for an optimal estimator. The method includes steps for state vector development, nominal dynamic/measurement model development, linearized model development, and uncertainty modeling and nominal model verification.
  • Figure 2 illustrates modeling of physiological processes in a linearized Kalman filter application.
  • glucose is an estimation variable.
  • the time rate of change of glucose might become another estimation variable.
  • An uncertain parameter might be the glucose sensor scale factor, and additionally, the rate of change of scale factor over time could be another.
  • the way in which variables and parameters nominally propagate over time may change with conditions.
  • inputs can be used to identify patient related activities: eating; exercising; sleeping; insulin injection; etc. With these identified patient related activities, additional state variables can be identified and modeled.
  • dynamics can be implemented within the estimator model which will propagate (or extrapolate) a rise in the glucose level over that time interval.
  • This rise may be modeled by appropriate functions whose variables contain uncertainties, which may be added as elements of the state vector.
  • decreases in glucose levels e.g., due to insulin injections
  • These models may be general in nature, or they may be patient specific. Consequently, patient related activities that have a significant impact on glucose levels, or the rate of change of glucose levels, can be accounted for within the dynamic process model.
  • infrared sensor measurement bias errors vary with, among other things, temperature. If this variation with temperature can be modeled, and included in the process model, then a temperature measurement will improve estimator performance.
  • a particular physical glucose sensor may have a scale factor which has a characteristic decay in sensitivity over time as discussed in further detail below. There may be other variables that can be measured which will affect the physical function of a sensor such as local pH, oxygen tension, etc.
  • the nominal modeling is comprised of three types: 1) predictable characteristics of a physical sensor function over time (e.g., a fixed rate in decline of sensor output), 2) other measurable physical variables which may affect sensor function (e.g., local temperature, pH, etc.), and 3) predictable changes in the model which occur with patient related activities (e.g., exercise, eating or insulin administration). Changes to the dynamic process model may add variables and/or uncertain parameters to the state vector and changes to measurement models. As a result, the activities indicated in blocks 1 and 2 of Figure 2 constitute an iterative process.
  • the relationship between variables, parameters, and measurements determines which parts of the processes and measurements are nonlinear and are linearized. Further, the relationship between variables, parameters, and measurements determines if the variables and parameters are observable and can be estimated. In certain cases, observability can be enhanced through introduction of additional modeling information. For the glucose example, the ability to estimate sensor scale factor and/or detect sensor failure may be improved by modeling glucose propagation changes due to insulin injections or ingestion of sugar. For example, tracking known changes enhances scale factor observability through estimator generated correlations. Further, if it is known that glucose levels vary and the sensor measurement does not change accordingly within prescribed levels of uncertainty, a sensor problem is indicated.
  • the development of a relatively large database is used to empirically verify and/or modify the nominal nonlinear dynamic process/measurement models and derive uncertainty levels associated with the variables, parameters, and measurements.
  • the empirical data can be fitted to nonlinear functions using a nonlinear regression package contained in a commercially available software application program such as "Mathematica" from Wolfram Research, Inc. Analytical functions may be added or modified using the test database.
  • the repeatability of the fit over nominal ranges of the patient environment determines the uncertain parameters and the variations establish uncertainty levels.
  • the more accurate the dynamic process and measurement models the more the uncertainties are reduced, and the greater the estimator performance.
  • the modeling becomes better defined through iterations between blocks 4 and 2. Certain portions of the models may be developed on an individual basis (e.g., for a specific patient).
  • a database is used to empirically develop and verify models.
  • An embodiment discussed below uses the database to develop two separate dynamic process models. For example, after processing a number of data sets from a physical sensor that was used to monitor the glucose level of various patients, it was observed that the sensor scale factor was equally likely to move up or down over the first fifteen to twenty-five hours. However, the scale factor tended to decay for the remaining life of the sensor after this period of time, usually three or four days. These observations are incorporated in an embodiment of the estimator discussed in further detail below.
  • Figure 3A describes an estimation function which depicts a linearized Kalman filter formulation.
  • the particular form is that of a discrete extended Kalman filter which linearizes after each update using best estimates.
  • a vector whose elements comprise variables and/or parameters (with uncertainties) for making estimates defines the state of a system.
  • variables e.g., random variables
  • parameters with uncertainty e.g., parameters with uncertainty.
  • Variables such as glucose and rate of change of glucose, are estimated and can be controlled (if control is implemented).
  • Parameters with uncertainty are part of the model structure not known precisely and are estimated and updated (like variables) but not controlled, e.g., a glucose sensor scale factor or insulin dispensing controller scale factor.
  • real-time variable and parameter estimates are used to re-linearize the model following each update.
  • Inputs to the estimator can consist of any measurement which is related to, or can be correlated with, any element in the state vector.
  • other inputs can consist of dynamic process configuration control based on patient related activities or other circumstances.
  • the time update and measurement update cycles form a recursive loop.
  • the time update period is the time interval between the receipt of measurements. This is a function of the sensor and of acceptable latency in the estimates.
  • Figures 3B, 3C and 3D illustrate, respectively, Initialization and the recursive Time Update and Measurement Update cycles according to one embodiment of the present invention.
  • Table 3 defines, in more detail, symbols used in these figures.
  • dynamic process and measurement modeling is contained in software algorithms with parameter and structure updates in real time.
  • a state vector estimate, Xeo contains the initial estimates of the process variables and model parameters with uncertainty, while the covariance matrix, Po, contains the initial variances and covariances associated with the Xeo elements.
  • Figure 3B shows an example embodiment with a two-element state vector which is based on empirical observations described above. As an example, the point in time when the scale factor begins to decay was chosen as
  • alpha is not modeled as an uncertain parameter; a value of 0.012 was chosen as a nominal value, over the database, for the five-minute cycle time of this physical sensor.
  • the initial estimate of glucose is 150 mg/dl (state vector element) with an initial uncertainty variance of 100 mg/dl squared (a covariance matrix element).
  • the glucose element was initialized by setting it substantially equal to the first capillary blood glucose measurement.
  • the nominal initial scale factor value, for the physical sensor is 0.25 nano-amps/(mg/dl) with a parameter uncertainty (variance) of 0.1 nano-amps/(mg/dl) squared.
  • the initial covariance between the glucose variable and the scale factor parameter is zero. Correlation between glucose and scale factor will develop as the estimator processes the sensor measurements. Each type of sensor will have its own model and characteristics. The measurement uncertainty is 5 nano-amps squared, an element of the R matrix. In one embodiment, this is a scalar measurement and the R matrix contains a single element.
  • a second measurement in this example embodiment is an occasional direct measurement, such as a capillary blood glucose measurement (CBG), with a unity scale factor, and a measurement error uncertainty of 15 mg/dl squared (a single element in a second R matrix). This is probably better modeled as 15% of the measured glucose value.
  • CBG capillary blood glucose measurement
  • the growth in uncertainty of glucose and scale factor from measurement to measurement is, respectively, 20 mg/dl squared and 0.002 nano-amps squared (elements of the process noise matrix, Q).
  • Figure 3C illustrates a time-update process of the recursive process in one embodiment of an optimal estimator.
  • a negative superscript indicates a time update while a positive superscript indicates a measurement update.
  • a letter "i" indicates time at the ith interval and (i-1) indicates time at the previous time interval.
  • a state vector is first updated since these elements are used to update matrices and to bring the time epoch of the estimated measurement to be consistent with that of the next measurement to be received.
  • the dynamic process is linear.
  • the first derivative of the measured physiologic variable is zero, corresponding to the case when the level (on average) is as likely to either go up or to go down. No additional a priori information is assumed about the time propagation of glucose.
  • the solution to the scale factor propagation after the first 20 hours is defined by the exponential shown in the second column. Consequently, for this embodiment, the dynamic process function (f), is linear, is not a function of state vector elements, and, from linear system theory, the transition matrix
  • (A) is the 2 by 2 identity matrix for the first 20 hours and thereafter is defined by the 2 by 2 matrix in the second column.
  • Figure 3D illustrates a measurement-update process of one embodiment of an optimal estimator.
  • the measurement update sequence begins with the computation of the gain matrix, K(i).
  • the covariance matrix is then measurement updated, reflecting the level of uncertainty in the estimates following the processing of a measurement.
  • a second measurement such as a
  • the measurement sequence is again cycled through, starting with a new gain computation, and using the " appropriate measurement matrix and new best estimate of the next measurement.
  • the updated state vector and covariance matrix are then used to start the time update for the (i+ 1)th time epoch which begins the next cycle.
  • Table 4 is an algorithm programmed in the MATLAB language (from Math Works). This printout defines a working program and has been used to process a significant number of data sets.
  • the estimation results from the process of two data files, both gathered from the same patient, and taken about a month apart, are shown in Figures 4A, 4B, 5A, and 5B.
  • sensor inputs are provided and processed every 5 minutes.
  • Occasional CBGs are also provided. For example, two CBGs per day were processed by an estimator; and additional CBG values were used to judge estimator performance by comparing glucose estimates with actual CBG values not used by the estimator.
  • Figures 4A and 4B illustrate a first set of time history plots of optimal glucose estimates and CBG measurements with respect to time.
  • Figure 4A shows a time history of the real-time estimates of glucose (e.g., every 5 minutes) along with all available discrete CBGs.
  • Figure 4B shows a time history of the glucose estimates along with CBGs that were processed. In this figure, the estimated glucose value took several rapid swings between approximately 48 hours and 58 hours.
  • Figure 4B shows a time history of the glucose estimates along with CBGs that were processed. In this figure, the estimated glucose value took several rapid swings between approximately 48 hours and 58 hours.
  • Figures 5A and 5B illustrate a second set of time history plots of optimal glucose estimates and CBG measurements with respect to time.
  • the time duration for these runs is about 4 days without patient inputs. Over this time period, the patient ate, slept, exercised, and took insulin injections. Dynamic models to account for these activities, much like the decaying scale factor, could be implemented and called into use upon command.
  • Figure 6A is a block diagram of one embodiment of a patient health monitor.
  • an algorithm for the patient health monitor is in Table 4 which generates statistical test data based on Kalman filter residuals as well as test displays.
  • the patient health monitor generates real-time decisions and displays which are integrated with a Kalman filter. The patient health monitor allows the patient or medical staff to interact with a
  • Kalman estimator and/or Kalman controller described herein.
  • the patient health monitor provides insight into how well a Kalman filter is working through the filter residual that is the difference between the estimate of the measurement at the time the measurement is received and the actual measurement.
  • other checks are used from time to time, such as the CBGs in the glucose example described above.
  • An example of another check on the status of a sensor is through the use of patient inputs or signals indicating that something is changing in a prescribed way and then noting whether or not the sensor is observing this change within prescribed uncertainties.
  • the real-time displays and decisions of the patient health monitor uses some occasional outside checks but relies substantially on results of statistical testing performed on filter residuals. If estimates of the measurements are, on average, good (e.g., residuals are small and unbiased), then the estimator is generally working well, and vice versa. More specifically, elements of the covariance matrix can be used to construct statistical test applications.
  • real-time displays of glucose estimates along with real-time displays of estimator performance test results can be important visual inputs to the medical staff. Requests for additional CBG measurements or the sounding of an alarm in the event glucose estimates exceed critical limits may also prove useful. In one embodiment, a reduced number of outputs is provided in relatively small estimators for individual use.
  • Figure 6B illustrates one embodiment of a residual test data process in the patient health monitor.
  • a time history plot of the residuals, Sy(i) in Figure 6B should appear random, zero mean, and unbaised. If they are summed over time, the deviation of the sum from zero should, on average, grow as the square root of time, as should its absolute value, ASy(i).
  • E[y(i)*y(j)] 0, for i not equal to j, also means that the sum of the variances, SV(i) in Figure 6B, which is easily computed, is equal to the variance of the sum of residuals for a properly performing estimator.
  • the standard deviation of the sum of residuals, StdSV(i) also grows as the square root of time.
  • Figure 6C illustrates one embodiment of a statistical test process in the patient health monitor.
  • Figure 6C defines statistical tests which can be constructed based on filter residual test data. One of these compares the absolute value of the sum of the residuals with the standard deviation of the sum of residuals. On average, the ASy(i) should be bounded by the StdSV(i). If not, this indicates that the deviations of Sy(i) are growing faster than that of a white noise sequence, implying degraded estimator performance.
  • Real-time displays of ASy(i) and StdSV(i) histories can provide a visual picture of estimator performance.
  • Figures 7A and 7B illustrate time history plots of residual test data with respect to measurement time.
  • Figure 6D illustrates one embodiment of an input/output interface in the patient health monitor.
  • the patient health monitor provides a capability for the patient and/or staff to communicate and interact in real time with estimation and control processes using simple commands, visual displays, and audio outputs described above. Confidence levels in estimator performance can be established and a genuine interface established whereby the estimation and control processes could request additional information to check and insure confidence in estimates, physical sensors, and physical controllers.
  • the staff can provide useful real-time inputs to augment this process.
  • Figure 8 is a block diagram of one embodiment of a Kalman optimal stochastic control solution as applied to physiological processes.
  • the application includes both optimal stochastic regulator and optimal stochastic tracking control solutions.
  • Tracking control involves a controlled variable following a reference value, constant or dynamic, as closely as possible.
  • Control is applied to physiological processes, wherein a control gain is computed based upon optimization criteria which minimizes controlled variable errors while minimizing application of control based on cost weightings.
  • the Kalman estimator gain minimizes estimation error variances.
  • the Kalman control gain minimizes variances of error between the controlled variables and the reference variables while minimizing the level of control applied.
  • the optimal control function includes specification of the controlled variables and their associated costs as well as costs associated with the amount of control to be applied. For example, smaller control variable error costs and larger control application costs will allow the controlled variable to deviate farther from the reference, but with reduced application of control.
  • linearization techniques described above in association with the Kalman estimator is applied to physiological nonlinear stochastic processes. For example, linearization about nominal values or about best estimates are provided by the Kalman estimator. Uncertain parameters associated with a controllable dispenser, or actuator, are included in an estimator state vector.
  • the optimal stochastic controller can use linear or nonlinear formulations and discrete or continuous time formulations.
  • the optimal stochastic controller is used with an optimal estimator described herein and/or an optimized decision and display function also described herein to form a closed loop system.
  • the closed loop system works as an artificial pancreas when applied to a glucose problem in one embodiment.
  • Figure 8 is a functional diagram showing a controllable dispenser, or pump, with the capability to secrete insulin and glucagon to control high and low glucose levels, respectively.
  • an estimator and a controller share a state vector wherein the estimator estimates it and the controller controls designated elements of it. Unlike the estimator, the time varying gain computations for the controller is computationally intensive and may not be used in all applications.
  • the controller includes time varying or steady state gain formulations.
  • FIG 9A illustrates one embodiment of an optimal controller for the artificial pancreas.
  • the glucose control problem is a tracking problem since the glucose level is controlled to a desired level which may either be constant or a function of time.
  • the dynamic process model includes a quantitative description of how glucose levels propagate, in time, as a function of insulin/glucagon secreted by the controllable pump.
  • the dynamic process model is described by a system of first order differential or difference equations. Excluding the pump, much of the modeling for the glucose control is available through the estimator modeling development.
  • Figure 9B illustrates one embodiment of a control model for a controller.
  • the dynamic process model is described by three first order differential equations that are forced by a control variable, u.
  • the dynamics of this process include a first order time lag and a scale factor associated with the pump, and a first order time lag and a scale factor associated with the physiological process.
  • the model for the glucose sensors is the same as in the example glucose estimator embodiment.
  • alpha g, s, and d are the inverse of the first order time lags associated with glucose, glucose sensor scale factor, and pump, respectively.
  • Beta g is the glucose scale factor multiplied by alpha g.
  • Delta t is the time interval between measurement/control application epochs.
  • the controlled error is the difference between the estimate of glucose (provided by the estimator) and the glucose control point input.
  • Beta d is a pump scale factor multiplied by alpha d.
  • the measurements are the same as for the example glucose estimator, except that the measurements contain an additional zero since the measurements are not functions of the variable insulin/glucagon.
  • the cost function which is minimized by the optimal stochastic control, contains costs associated with the glucose error and the application of control. Choosing a value of Cu that is much larger than Cg will result in a relative gentle application of control.
  • Another example embodiment would utilize a higher order model in which the first derivative of glucose would be included in the state vector, and included as a control variable. If a relatively large cost is associated with the first derivative relative to the glucose control point error, then the control will be very active when rapid changes in patient glucose occur while relatively gentle otherwise.
  • Other example embodiments would include pump scale factor as a state vector element (uncertain parameter) as well as insulin/glucagon measurements to the estimator.
  • Figure 9C illustrates one embodiment of a control algorithm based on the model defined in Figure 9B.
  • control is applied at each epoch based upon a difference between an estimate of a patient glucose level and a glucose control point.
  • a control, u is applied to a pump over each time interval, and a state vector is time updated as indicated in this figure using the control variable, u, as a forcing function.
  • Other time updates are performed in accordance with embodiments discussed above in association with estimator equations.
  • a potential problem in the application of closed loop control to that of physiological processes is due to potentially long time delays that may be de-stabilizing. These delays can be in the form of transport lags.
  • a transport lag is the time between when control is applied and when the process action begins.
  • X, Xe, Xe(0) denote, respectively, the state vector, the best estimate (or optimal estimate) of the state vector, and the initial best estimate of the state vector. They contain both variables and parameters as defined above and, with regard to estimation, there are no distinctions.
  • P, P(0) Denotes, respectively, the matrix of variances and co-variances associated with the errors of each of the state variable estimates (within the estimated state vector) and their initial values,
  • A, A(0) Denotes the transition matrix, which is used to propagate the covariance matrix forward it time (time update), and it's initial value, respectively; in the linear case it is also used to time update the state vector estimate.
  • Q Denotes the matrix of variances and co-variances associated with error growth uncertainty accumulated in the state variable estimates since the last measurement update
  • R Denotes the matrix of variances and co-variances associated with measurement error uncertainties.
  • K Denotes the Kalman gain matrix which, when multiplied by the difference between the actual measurement, and the best estimate of the measurement, yields the estimated state correction. This estimated state correction, when added to the old best estimate, becomes the new best estimate.
  • figure (1) plot (t, xhistory ( : ,2) ) ; xlabel ( ' Measurement Time (hrs)'); ylabel ( ' Scale Factor Estimate; (ISIG Units/ (mg/dl) ) * ) ; title ('Scale Factor Estimate vs Tir ⁇ e'); figure (2) plot (t, gmhistory ( : ) , fc, xhistor ( :,1) ) ; xlabel ( ' Time (hrs ) ' ) ;
  • ⁇ k) c + ' - ) + w 2 , (2)
  • Y is an m x 1 observation vector
  • X is an n x 1 random state vector whose initial uncertainty is uncorrelated with W, and W 2 and with initial covariance Q 0
  • U is the control input vector.
  • the system representation defined by Eqs. (1) and (2) can be cast in stochastic regulator form by augmenting the state vector with the time- delayed states, that is,
  • h max(p, q).
  • Eq. (3) the control remains the same as in Eq. ( 1) because, physically, past states cannot be changed or controlled.
  • Eq. (A) the observation Y and the observation noise W 2 also remain unchanged.
  • the dimensions of B and Wj , in Eq. (3) must also be consistent with the integer h. Therefore, without loss of generality, the system defined by Eqs. ( 1 ) and (2) can be represented by the augmented system model described by Eqs. (3) and (4).
  • the optimal linear stochastic control of the augmented system is given by a deterministic, optimal linear controller with state input X (or estimated state feedback), which is provided by an optimal one-step predictor using the augmented model [4]. That is,
  • X(k + 1) A(k)X(k) + B(k)V(k) + K(k)[Y(k) - C(k)X(k) ⁇ , (9) the estimator gain K satisfies
  • the final value of P used to "initialize” Eq. (8) (which is solved backward in time) is the final value defined in the quadratic cost function of Eq. (5), that is, P(k ⁇ ) — P,.
  • the above solution is the optimal solution without qualification; it not, then it is the optima] linear control solution.
  • the expected system performance is determined by analyzing the augmented system as a linear, stochastic regulator problem.
  • control gain Eqs. (7) and (8) and estimator gain Eqs. (10) and (11) it can be shown that the required dimensions of the controller and estimator are not, in general, equal.
  • the optimal control is given by
  • V ⁇ k - £ F t (k)X(k + 1 - 0, (12)
  • X(k) is the one-step predicted value of the original system state vector
  • X(k — 1) is the filtered value
  • X(k — 2) is the one-step smoothed value
  • X (k + 1 — ft) is the (ft — 2)th smoothed value.
  • the real-time computational requirements associated with the implementation of time-varying optimal control are always stressing due to the "backward-in-time" recursion that is required to obtain solutions P of Eq. (8).
  • the A, B, R t , and R 2 matrices of the augmented system and cost function can be treated as time invariant over the time intervals of interest.
  • the control gain F will converge to a unique value such that the steady- state optimal control law is time invariant, asymptotically stable, and minimizes the quadratic cost function of Eq. (5) as &, ⁇ oo.
  • the steady-state gain matrix F B can be computed off-line and stored for realtime use; the real-time computations required to implement this steady-state control law are negligible.
  • the steady-state control gains and time-varying gains are such that the initial value of the time-varying gain is equal to the steady-state gain but is less than or equal to the steady-state gain as time progresses, that is,
  • the steady-slate gain tends to maintain the controlled state closer to the estimated state, but at the cost of control energy. If accuracy is the significant criterion, then the steady-state gain is not only easy to implement but also provides essentially equivalent or better accuracy.
  • the estimator gain is usually significantly larger than the steady-state gain, to account for initial uncertainties in the knowledge of the system state. As time progresses, the time-varying gain converges to the steady-state gain, that is,
  • the estimator Riccati equation can be separated into a time update and a measurement update, where the time update for the augmented system becomes primarily one of data transfer. If Q, from Eq. (11), is defined as
  • the iteration interval ⁇ t (implied in the discrete-time system representation) or some integer number of ⁇ f's should be set equal to the time delay.
  • it may be advantageous to use (a) variable iteration intervals for the estimator, (b) fixed iteration intervals for the controller, and (c) variable-time updates of the state estimate to time synch the estimator output with the controller input
  • the filter equation becomes a partial differential equation with a boundary condition
  • the covariance equation becomes a partial differential matrix equation with three boundary conditions (see [2]). Consequently, the most practical control solution is obtained by discretizing the continuous-time system representation and then applying the approach of Section III.
  • One technique for discretizing the continuous-time representation is by using the Z transform method; see, for example, [5]. Using a scalar differential equation with one delayed state as an example, we have
  • A-(.) ⁇ fl.-X(f - i ⁇ .) + l/( . (18)

Abstract

A real-time glucose estimator uses a linearized Kalman filter to determine a best estimate of glucose level in real time. The real-time glucose estimator receives at least one measurement corresponding to glucose level. The measurement can be obtained with one or more sensors and is provided to the linearized Kalman filter in real time. The linearized Kalman filter has dynamic models and executes a recursive routine to determine the best estimate of glucose level based upon the measurement. Additional information can be provided to the linearized Kalman filter for initialization, configuration, and the like. Outputs of the linearized Kalman filter can be provided to a patient health monitor for display or for statistical testing to determine status of the real-time glucose estimator. The real-time glucose estimator can be implemented using a software algorithm. A real-time controller operating as an artificial pancreas uses a Kalman control algorithm to control glucose level of a patient in real time. The real-time controller receives an estimate of the patient glucose level and a reference glucose level. The estimate of the patient glucose level can be provided by an optimal estimator implemented using a linearized Kalman filter. The estimated glucose level and the reference glucose level are processed by the Kalman control algorithm to determine a control command in real time. The Kalman control algorithm has a dynamic process forced by the control command a cost function determining a relative level of control. The control command is provided to a dispenser which secretes insulin or glucagon in response to the control command to correct a relatively high glucose level or a relatively low glucose level.

Description

METHOD AND APPARATUS FOR REAL-TIME ESTIMATION AND CONTROL OF PHYSIOLOGICAL PARAMETERS
Copyright Rights A portion of the disclosure of this patent document contains material that is subject to copyright protection. The copyright owner has no objection to the facsimile reproduction by anyone of the patent document or of the patent disclosure as it appears in the Patent and Trademark Office patent files or records, but otherwise reserves all copyright rights whatsoever.
Background of the Invention Field of the Invention
This present invention relates generally to a method and apparatus for estimating and controlling physiological parameters and more particularly to an optimal estimator and an optimal controller for estimating and controlling glucose levels in a patient.
Description of the Related Art
Different types of sensors (e.g., optical sensors) are available for monitoring of physiological parameters (e.g., glucose concentration). Glucose monitoring is typically performed by people with diabetes mellitus which is a medical condition involving a body's inability to produce the quantity or quality of insulin needed to maintain a normal circulating blood glucose. Frequent monitoring of glucose is generally necessary to provide effective treatment and to prevent long term complications of diabetes (e.g., blindness, kidney failure, heart failure, etc.). New methods of monitoring glucose are fast, painless and convenient alternatives to the typical capillary blood glucose (CBG) measurements which involve finger pricks that are painful, inconvenient and difficult to perform for long term. Optical measurement of glucose is performed by focusing a beam of light onto the body. Optical sensors determine glucose concentration by analyzing optical signal changes in wavelength, polarization or intensity of light. However, many factors other than glucose concentration also contribute to the optical signal changes. For example, sensor characteristics (e.g., aging), environmental variations (e.g., changes in temperature, humidity, skin hydration, pH, etc.), and physiological variations (e.g., changes in tissue fluid due to activity, diet, medication or hormone fluctuations) affect sensor measurements.
Various methods are used to improve the accuracy of the sensor measurements. One method (e.g., multivariate spectral analysis) utilizes calibration models developed by initially measuring known glucose concentrations to correct subsequent sensor measurements. The calibration models become inaccurate over time due to dynamic changes in physiological processes. Another method (e.g., adaptive noise canceling) utilizes signal processing to cancel portions of the sensor measurements unrelated to glucose concentration. For example, two substantially simultaneous sensor measurements at different wavelengths make up a composite signal which can be processed to cancel its unknown and erratic portions. However, many sensors do not provide substantially simultaneous measurements at two different wavelengths.
Summary of the Invention
The present invention solves these and other problems by providing a method and apparatus for making optimal estimates of a physiological parameter (e.g., glucose level), assessing reliability of the optimal estimates, and/or providing optimal control of the physiological parameter in real time using one or more sensor measurements at each measurement time epoch (or interval). The sensor measurements can be time-based (e.g., every five minutes) to provide continuous monitoring and/or regulation of the physiological parameter. The sensor measurements are a function of the physiological parameter within specified uncertainties.
An optimal estimator provides an accurate estimate of glucose level in real time using a sensor with at least one output. In one embodiment, the optimal estimator is integrated with the sensor and an output display to be a compact glucose monitoring device which can be worn by a patient for continuous monitoring and real-time display of glucose level. In an alternate embodiment, the optimal estimator is a separate unit which can interface with different types of sensors and provide one or more outputs for display, further processing by another device, or storage on a memory device.
In one embodiment, the optimal estimator employs a priori deterministic dynamic models developed with stochastic variables and uncertain parameters to make estimates of glucose level. For example, glucose level is defined as one of the stochastic (or random) variables. Dynamic mathematical models define process propagation (i.e., how physiological and sensor parameters change in time) and measurement relationship (i.e., how physiological and sensor parameters relate to environmental conditions). Environmental conditions (e.g., temperature, humidity, pH, patient activity, etc.) can be provided to the optimal estimator intermittently or periodically via environment sensors and/or data entries by a patient or a doctor. The optimal estimator uses dynamic models to propagate estimates of respective stochastic variables, error variances, and error covariances forward in time. At each measurement time epoch, the optimal estimator generates real-time estimates of the stochastic variables using one or more sensor outputs and any ancillary input related to environmental conditions. In one embodiment, the optimal estimator employs a linearized Kalman filter to perform optimal estimation of the stochastic variables (e.g., glucose level). In particular, an extended Kalman filter is used to accommodate nonlinear stochastic models.
Before making real-time estimates, the optimal estimator is initialized by providing initial values for the stochastic variables, error variances, and error covariances. For example, a CBG measurement or another direct glucose measurement is performed at initialization to provide a starting value for the stochastic variable corresponding to glucose level. In one embodiment, the optimal estimator provides one or more outputs to a patient health monitor which is capable of optimized real-time decisions and displays. The patient health monitor evaluates system performance by assessing the performance of the sensor and/or optimal estimator in real time. For example, the patient health monitor applies statistical testing to determine the reliability of the real-time estimates of the stochastic variables by the optimal estimator. The statistical testing is performed in real time on residual errors of the optimal estimator to establish performance measures.
In one embodiment, the patient health monitor acts as an input/output interface between the patient or medical staff (e.g., a doctor, nurse, or other healthcare provider) and the optimal estimator. For example, environmental conditions can be provided to the patient health monitor for forwarding to the optimal estimator. Optimal estimator outputs can be provided to the patient health monitor for display or forwarding to an external device
(e.g., a computer or a data storage device).
In one embodiment, the optimal estimator provides one or more outputs to an optimal controller which can regulate in real time the physiological parameter being monitored. For example, an optimal controller responds to realtime optimal estimator outputs and provides an output to operate an actuator. In the case of glucose control, the actuator can be a dispenser or a pump which secretes insulin to correct a relatively high glucose level and glucagon to correct a relatively low glucose level. The optimal controller takes advantage of a priori information regarding the statistical characteristics of the actuator and is able to control the output of the actuator to be within specified uncertainties.
In one embodiment, the optimal estimator and the optimal controller form an optimal closed-loop system. For example, a glucose sensor, an optimal estimator, an optimal controller, and an insulin/glucagon dispenser work together as an artificial pancreas to continuously regulate glucose level. The glucose sensor can be internal or external to a patient's body. The optimal controller provides a control feedback to the optimal estimator to account for delivery of the insulin/glucagon.
The optimal closed-loop system is effective in a variety of biomedical applications. For example, cardiovascular functions can be continuously regulated by using sensors to detect blood pressure, blood oxygen level, physical activity and the like, an optimal estimator to process the sensor measurements and make real-time estimates of heart function parameters, and an optimal controller to control operations of an artificial device (e.g., a pacemaker) in real time based on the real-time estimates from the optimal estimator to achieve a set of desired heart function parameter values. Other artificial devices (e.g., artificial limbs, bionic ears, and bionic eyes) can be part of similar optimal closed-loop systems with sensors detecting nerve signals or other appropriate signals.
The optimal closed-loop system is also effective in optimal treatment of chronic illnesses (e.g., HIV). Some medications for treatment of chronic illnesses are relatively toxic to the body. Over delivery of medication generally has adverse effects on the patient. The optimal closed-loop system is capable of providing effective and safe treatment for the patient. For example, an optimal estimator provides real-time estimates of key physiological parameters using one or more sensors, and an optimal controller controls a slow infusion of medication in real time based on the real-time estimates from the optimal estimator to obtain desirable values for the key physiological parameters.
In one embodiment, the optimal estimator, patient health monitor, and optimal controller are software algorithms which can be implemented using respective microprocessors. New information regarding process propagation or measurement relationship can be easily incorporated by modifying, reconfiguring, and/or adding to the software algorithms. The optimal estimator, patient health monitor, and optimal controller can be implemented as one joint algorithm or separate respective algorithms which function together to provide an optimal closed-loop system.
Brief Description of the Drawings Figure 1 A is a block diagram of one embodiment of an estimator.
Figure 1 B is a block diagram of another embodiment of a glucose-monitoring device. Figure 2 illustrates one embodiment of modeling physiological processes in a linearized Kalman filter application.
Figure 3A describes one embodiment of an estimation function which depicts a linearized Kalman filter formulation.
Figure 3B illustrates an initialization in one embodiment of an optimal estimator. Figure 3C illustrates a time-update cycle in accordance with one embodiment of an optimal estimator. Figure 3D illustrates a measurement-update cycle in accordance with one embodiment of an optimal estimator. Figures 4A and 4B illustrate a first set of time history plots of optimal glucose estimates and CBG measurements with respect to time.
Figures 5A and 5B illustrate a second set of time history plots of optimal glucose estimates and CBG measurements with respect to time.
Figure 6A is a block diagram of one embodiment of a patient health monitor. Figure 6B illustrates one embodiment of a residual test data process in the patient health monitor.
Figure 6C illustrates one embodiment of a statistical test process in the patient health monitor. Figure 6D illustrates one embodiment of an input/output interface in the patient health monitor. Figures 7A and 7B illustrate time history plots of residual test data with respect to measurement time. Figure 8 is a functional diagram of one embodiment of an artificial pancreas. Figure 9A illustrates one embodiment of an optimal controller for a closed loop system.
Figure 9B illustrates one embodiment of a control model for a controller in accordance with the present invention.
Figure 9C illustrates one embodiment of a control algorithm.
Detailed Description of the Invention The present invention involves application of real-time optimal estimation, optimized real-time decisions and displays (e.g., a patient health monitor), or optimal real-time control to physiological processes. In one embodiment, the real-time optimal estimation, the optimized real-time decisions and displays, and the optimal real-time control are implemented as separate modules which can be combined functionally. In an alternate embodiment, the real-time optimal estimation, the optimized real-time decisions and displays, and the optimal real-time control are implemented as one joint algorithm.
In one embodiment, input to an optimal estimator is provided by a physical sensor (or a plurality of sensors) which measures some arbitrary, but known, function (or functions) of variables or parameters to be estimated to within specified uncertainties and whose statistical characteristics are known. In one embodiment, an output of a real-time controller is provided to a physical controllable dispenser, or actuator, whose output is some known function of parameters to be controlled within specified uncertainties and whose statistical characteristics are known. In one embodiment, a decision and display function utilizes statistical testing of estimator residual errors using internally computed, and updated, estimator variances and covariances.
In one embodiment, the present invention is implemented as a software algorithm. The present invention uses models (e.g., dynamic process and measurement models). For best performance, the models should reflect the latest and most complete information available. As new and more complete information is developed, performance can be improved through incorporation of this information by simply modifying the software algorithm of the present invention.
Embodiments of the present invention will be described hereinafter with reference to the drawings. Figure 1A is a block diagram of one embodiment of an estimator. The estimator uses a linearized Kalman filter. The linearized Kalman filter accommodates nonlinear process models and/or nonlinear measurement models. In one embodiment, the linearized Kalman filter is a discrete extended Kalman filter which is linearized after each update using best estimates.
A general formulation of a continuous-discrete extended Kalman filter is provided in Table 1. In an alternate embodiment, the linearized Kalman filter is linearized about a nominal set for which a general formulation is given in
Table 2.
The estimator computes an estimator gain based on time updated and measurement updated error variable variances and covariances. In one embodiment, the estimator is implemented using discrete formulations. In an alternate embodiment, the estimator is implemented using continuous formulations. In the actual development of algorithms, one can choose a Covariance formulation or an Information formulation depending on initialization uncertainty considerations. Further, the use of Bierman factorization techniques (UDITT) in the implementation leads to numerically stable algorithms which are excellent for operation over very long time periods.
In one embodiment, the estimator of Figure 1A is applied to the problem of monitoring patient glucose levels. Any type of physical sensor which measures some function of glucose can be used. In one embodiment, one or more capillary blood glucose (CBG) measurements are obtained on initialization of the estimator and when estimated glucose deviations exceed computed variance levels.
In another embodiment, additional capillary glucose values are obtained on a periodic basis to assess sensor function. As an example, CBG measurements may be obtained once or twice a day pre and one hour post prandial or when the estimator determines that the glucose values are out of range of predetermined limits.
In one embodiment, the estimator is contained in a small portable package. The estimator can be operated by a patient. Alternatively, the estimator can be operated by medical staff in a hospital or clinic. The estimator (e.g., Kalman estimator) uses a Kalman filter to provided real-time estimates of patient glucose levels based upon glucose sensor measurements. The Kalman estimator can also use additional information to more closely predict time propagation of glucose levels, such as exercise, food intake, insulin administration, or other factors which influence glucose levels or the function of the sensor (e.g., local pH, temperature or oxygen tension near the physical sensor).
Figure 1B is a block diagram of one embodiment of a glucose-monitoring device 100. The glucose-monitoring device 100 includes a glucose sensor 108, an ancillary sensor 110, a glucose estimator 114, and a patient health monitor 114. The glucose-monitoring device 100 is used to provide real-time glucose estimates of a patient 102. In one embodiment, the glucose-monitoring device 100 is an integrated unit which is portable by the patient 102 to provide continuous glucose monitoring and real-time displays.
The glucose sensor 108 (e.g., a probe, patch, infrared, or laser sensor) is coupled to the patient 102 and outputs a measurement f(g) which is a function of the glucose level of the patient 102. The glucose sensor 108 makes measurements within specified uncertainties and has known statistical characteristics. The glucose sensor 108 provides the measurement f(g) to the glucose estimator 112. In one embodiment, the glucose sensor 108 makes measurements periodically. In an alternate embodiment, the glucose sensor 108 makes measurements intermittently or upon command.
The ancillary sensor 110 is coupled to the patient 102, the glucose sensor 108, and/or surroundings of the patient 102 and/or glucose sensor 108 to provide information regarding environmental and/or glucose sensor conditions (e.g., temperature, humidity, local pH, etc.) which affect the measurement. One or more outputs of the ancillary sensor
110 are provided to the glucose estimator 112. In one embodiment, the ancillary sensor 110 provides outputs to the patient health monitor 114 which can process the information for display or for forwarding to the glucose estimator 112. In an alternate embodiment, the ancillary sensor 110 is not a part of the glucose-monitoring device 100.
The patient 102 and/or a medical staff 104 (i.e., a user) can provide information on the environmental and glucose sensor conditions as well as other information affecting the measurement. In one embodiment, the patient
102 and/or medical staff 104 inputs information (e.g., exercise activity, food intake, insulin administration, etc.) using the patient health monitor 114. The patient health monitor 114 acts as an input/output interface or a means for the user to configure the glucose estimator 112. The patient health monitor 114 forwards the information to the glucose estimator 112. In one embodiment, the patient health monitor 114 is a display device. For example, an output of the glucose estimator 1 12 (e.g., an optimal real-time estimate of glucose) is provided to the patient health monitor 114 which displays the information in a comprehensible format for the patient 102 and/or the medical staff 104. The patient health monitor 114 can also contemporaneously display information provided by the glucose sensor 108, the ancillary sensor 1 10, the patient 102 and/or the medical staff 104 which affects the measurement used to make the optimal glucose estimate.
In another embodiment, the patient health monitor 1 14 is a status indicator. For example, the glucose estimator 1 12 provides outputs (e.g., residuals and variances) to the patient health monitor 1 14 which applies statistical testing to determine the reliability of the glucose sensor 108 and/or the glucose estimator 1 12. The patient health monitor 114 provides a warning when poor performance is detected.
In one embodiment, the glucose estimator 112 is a linearized Kalman filter (e.g., a discrete extended Kalman filter) to account for a nonlinear process model and/or measurement model. The glucose estimator (or Kalman estimator) 1 12 provides real-time estimates of the glucose level in the patient 102. An initialization measurement 106 (e.g., a CBG measurement) is obtained from the patient 102 and provided to the Kalman estimator 112 to initialize the Kalman estimator 112.
Figure 2 is a flow chart illustrating one embodiment of a modeling process for physiological processes. In one embodiment, physiological processes are described by nonlinear stochastic models. In one embodiment, the flow chart of Figure 2 illustrates one method of developing a dynamic model for an optimal estimator. The method includes steps for state vector development, nominal dynamic/measurement model development, linearized model development, and uncertainty modeling and nominal model verification.
Figure 2 illustrates modeling of physiological processes in a linearized Kalman filter application. As an example in glucose estimation, in a first block of Figure 2, glucose is an estimation variable. The time rate of change of glucose might become another estimation variable. An uncertain parameter might be the glucose sensor scale factor, and additionally, the rate of change of scale factor over time could be another. In one embodiment of a second block of Figure 2, the way in which variables and parameters nominally propagate over time may change with conditions. In the example of glucose estimation, inputs can be used to identify patient related activities: eating; exercising; sleeping; insulin injection; etc. With these identified patient related activities, additional state variables can be identified and modeled. For example, if patient eating can be related to a change in glucose level over some specified time interval, dynamics can be implemented within the estimator model which will propagate (or extrapolate) a rise in the glucose level over that time interval. This rise may be modeled by appropriate functions whose variables contain uncertainties, which may be added as elements of the state vector. In an analogous way, decreases in glucose levels (e.g., due to insulin injections) can be modeled. These models may be general in nature, or they may be patient specific. Consequently, patient related activities that have a significant impact on glucose levels, or the rate of change of glucose levels, can be accounted for within the dynamic process model. Insofar as sensor modeling is concerned, experience shows that, for example, infrared sensor measurement bias errors vary with, among other things, temperature. If this variation with temperature can be modeled, and included in the process model, then a temperature measurement will improve estimator performance. A particular physical glucose sensor may have a scale factor which has a characteristic decay in sensitivity over time as discussed in further detail below. There may be other variables that can be measured which will affect the physical function of a sensor such as local pH, oxygen tension, etc.
In one embodiment, the nominal modeling is comprised of three types: 1) predictable characteristics of a physical sensor function over time (e.g., a fixed rate in decline of sensor output), 2) other measurable physical variables which may affect sensor function (e.g., local temperature, pH, etc.), and 3) predictable changes in the model which occur with patient related activities (e.g., exercise, eating or insulin administration). Changes to the dynamic process model may add variables and/or uncertain parameters to the state vector and changes to measurement models. As a result, the activities indicated in blocks 1 and 2 of Figure 2 constitute an iterative process.
In one embodiment of block 3 of Figure 2, the relationship between variables, parameters, and measurements determines which parts of the processes and measurements are nonlinear and are linearized. Further, the relationship between variables, parameters, and measurements determines if the variables and parameters are observable and can be estimated. In certain cases, observability can be enhanced through introduction of additional modeling information. For the glucose example, the ability to estimate sensor scale factor and/or detect sensor failure may be improved by modeling glucose propagation changes due to insulin injections or ingestion of sugar. For example, tracking known changes enhances scale factor observability through estimator generated correlations. Further, if it is known that glucose levels vary and the sensor measurement does not change accordingly within prescribed levels of uncertainty, a sensor problem is indicated.
In one embodiment of block 4 of Figure 2, the development of a relatively large database is used to empirically verify and/or modify the nominal nonlinear dynamic process/measurement models and derive uncertainty levels associated with the variables, parameters, and measurements. In one embodiment, the empirical data can be fitted to nonlinear functions using a nonlinear regression package contained in a commercially available software application program such as "Mathematica" from Wolfram Research, Inc. Analytical functions may be added or modified using the test database. The repeatability of the fit over nominal ranges of the patient environment determines the uncertain parameters and the variations establish uncertainty levels. The more accurate the dynamic process and measurement models, the more the uncertainties are reduced, and the greater the estimator performance. As the process evolves, the modeling becomes better defined through iterations between blocks 4 and 2. Certain portions of the models may be developed on an individual basis (e.g., for a specific patient).
A database is used to empirically develop and verify models. An embodiment discussed below uses the database to develop two separate dynamic process models. For example, after processing a number of data sets from a physical sensor that was used to monitor the glucose level of various patients, it was observed that the sensor scale factor was equally likely to move up or down over the first fifteen to twenty-five hours. However, the scale factor tended to decay for the remaining life of the sensor after this period of time, usually three or four days. These observations are incorporated in an embodiment of the estimator discussed in further detail below.
Figure 3A describes an estimation function which depicts a linearized Kalman filter formulation. In one embodiment, the particular form is that of a discrete extended Kalman filter which linearizes after each update using best estimates. In the formulation, a vector whose elements comprise variables and/or parameters (with uncertainties) for making estimates defines the state of a system.
There is a distinction between variables (e.g., random variables) and parameters with uncertainty. Variables, such as glucose and rate of change of glucose, are estimated and can be controlled (if control is implemented). Parameters with uncertainty are part of the model structure not known precisely and are estimated and updated (like variables) but not controlled, e.g., a glucose sensor scale factor or insulin dispensing controller scale factor.
In one embodiment, real-time variable and parameter estimates are used to re-linearize the model following each update. Inputs to the estimator can consist of any measurement which is related to, or can be correlated with, any element in the state vector. In the case of glucose estimation, other inputs can consist of dynamic process configuration control based on patient related activities or other circumstances. Following initialization, the time update and measurement update cycles form a recursive loop. The time update period is the time interval between the receipt of measurements. This is a function of the sensor and of acceptable latency in the estimates.
Figures 3B, 3C and 3D illustrate, respectively, Initialization and the recursive Time Update and Measurement Update cycles according to one embodiment of the present invention. Table 3 defines, in more detail, symbols used in these figures. In one embodiment, dynamic process and measurement modeling is contained in software algorithms with parameter and structure updates in real time.
One embodiment of an initialization of the estimation process is described in Figure 3B. A state vector estimate, Xeo, contains the initial estimates of the process variables and model parameters with uncertainty, while the covariance matrix, Po, contains the initial variances and covariances associated with the Xeo elements.
Figure 3B shows an example embodiment with a two-element state vector which is based on empirical observations described above. As an example, the point in time when the scale factor begins to decay was chosen as
20 hours, a nominal value over the database. An exponential decay was chosen to model this decay rate and is consistent with the first derivative of scale factor equal to a parameter, alpha, multiplied by the variable scale factor.
In this case, alpha is not modeled as an uncertain parameter; a value of 0.012 was chosen as a nominal value, over the database, for the five-minute cycle time of this physical sensor. The initial estimate of glucose is 150 mg/dl (state vector element) with an initial uncertainty variance of 100 mg/dl squared (a covariance matrix element). When processing the data, the glucose element was initialized by setting it substantially equal to the first capillary blood glucose measurement. The nominal initial scale factor value, for the physical sensor, is 0.25 nano-amps/(mg/dl) with a parameter uncertainty (variance) of 0.1 nano-amps/(mg/dl) squared.
The initial covariance between the glucose variable and the scale factor parameter is zero. Correlation between glucose and scale factor will develop as the estimator processes the sensor measurements. Each type of sensor will have its own model and characteristics. The measurement uncertainty is 5 nano-amps squared, an element of the R matrix. In one embodiment, this is a scalar measurement and the R matrix contains a single element.
A second measurement in this example embodiment is an occasional direct measurement, such as a capillary blood glucose measurement (CBG), with a unity scale factor, and a measurement error uncertainty of 15 mg/dl squared (a single element in a second R matrix). This is probably better modeled as 15% of the measured glucose value. The growth in uncertainty of glucose and scale factor from measurement to measurement is, respectively, 20 mg/dl squared and 0.002 nano-amps squared (elements of the process noise matrix, Q).
Figure 3C illustrates a time-update process of the recursive process in one embodiment of an optimal estimator. In this figure, a negative superscript indicates a time update while a positive superscript indicates a measurement update. In the brackets, a letter "i" indicates time at the ith interval and (i-1) indicates time at the previous time interval. In an update process, a state vector is first updated since these elements are used to update matrices and to bring the time epoch of the estimated measurement to be consistent with that of the next measurement to be received.
In an example embodiment with a two-element state vector, the dynamic process is linear. With no patient inputs, the first derivative of the measured physiologic variable is zero, corresponding to the case when the level (on average) is as likely to either go up or to go down. No additional a priori information is assumed about the time propagation of glucose.
In one embodiment of estimating glucose, the solution to the scale factor propagation after the first 20 hours is defined by the exponential shown in the second column. Consequently, for this embodiment, the dynamic process function (f), is linear, is not a function of state vector elements, and, from linear system theory, the transition matrix
(A) is the 2 by 2 identity matrix for the first 20 hours and thereafter is defined by the 2 by 2 matrix in the second column.
In one embodiment, the measurement function (h) for the sensor measurement is non-linear in the state vector elements. If Ge and Se are used to denote glucose and scale factor estimates, respectively, then Ye=Se*Ge. Definitions of the above terms are provided in Table 3. When linearized using best estimates, the linearized measurement matrix H - [Se Ge] and is of the same functional form both before and after 20 hours.
Figure 3D illustrates a measurement-update process of one embodiment of an optimal estimator. The measurement update sequence begins with the computation of the gain matrix, K(i). The difference between the actual sensor measurement and the best estimate of the measurement is computed: y(i) = Ym(i) - Se(i)*Ge(i). This difference, or residual, when multiplied by the gain matrix and added to the time-updated estimate produces the measurement-updated estimate.
The covariance matrix is then measurement updated, reflecting the level of uncertainty in the estimates following the processing of a measurement. In the glucose example, if a second measurement is available, such as a
CBG, then the measurement sequence is again cycled through, starting with a new gain computation, and using the" appropriate measurement matrix and new best estimate of the next measurement. Following the processing of all available measurements at the ith time epoch, the updated state vector and covariance matrix are then used to start the time update for the (i+ 1)th time epoch which begins the next cycle.
One embodiment of an estimation algorithm illustrating the initialization, the time update process, and the measurement update process discussed above is provided in Table 4. Table 4 is an algorithm programmed in the MATLAB language (from Math Works). This printout defines a working program and has been used to process a significant number of data sets. In the glucose example, the estimation results from the process of two data files, both gathered from the same patient, and taken about a month apart, are shown in Figures 4A, 4B, 5A, and 5B.
In one embodiment, sensor inputs are provided and processed every 5 minutes. Occasional CBGs are also provided. For example, two CBGs per day were processed by an estimator; and additional CBG values were used to judge estimator performance by comparing glucose estimates with actual CBG values not used by the estimator.
Figures 4A and 4B illustrate a first set of time history plots of optimal glucose estimates and CBG measurements with respect to time. Figure 4A shows a time history of the real-time estimates of glucose (e.g., every 5 minutes) along with all available discrete CBGs.
Figure 4B shows a time history of the glucose estimates along with CBGs that were processed. In this figure, the estimated glucose value took several rapid swings between approximately 48 hours and 58 hours. Figure
4B indicates that only sensor measurements were processed during that interval and no CBGs were processed.
However, the CBGs plotted in Figure 4A indicate that the glucose estimates did tend to follow the excursions of the patient glucose levels.
Figures 5A and 5B illustrate a second set of time history plots of optimal glucose estimates and CBG measurements with respect to time. The time duration for these runs is about 4 days without patient inputs. Over this time period, the patient ate, slept, exercised, and took insulin injections. Dynamic models to account for these activities, much like the decaying scale factor, could be implemented and called into use upon command.
Figure 6A is a block diagram of one embodiment of a patient health monitor. On example of an algorithm for the patient health monitor is in Table 4 which generates statistical test data based on Kalman filter residuals as well as test displays. In one embodiment, the patient health monitor generates real-time decisions and displays which are integrated with a Kalman filter. The patient health monitor allows the patient or medical staff to interact with a
Kalman estimator and/or Kalman controller described herein.
In one embodiment, the patient health monitor provides insight into how well a Kalman filter is working through the filter residual that is the difference between the estimate of the measurement at the time the measurement is received and the actual measurement. In another embodiment, other checks are used from time to time, such as the CBGs in the glucose example described above. An example of another check on the status of a sensor is through the use of patient inputs or signals indicating that something is changing in a prescribed way and then noting whether or not the sensor is observing this change within prescribed uncertainties.
In one embodiment, the real-time displays and decisions of the patient health monitor uses some occasional outside checks but relies substantially on results of statistical testing performed on filter residuals. If estimates of the measurements are, on average, good (e.g., residuals are small and unbiased), then the estimator is generally working well, and vice versa. More specifically, elements of the covariance matrix can be used to construct statistical test applications.
For hospital applications, real-time displays of glucose estimates along with real-time displays of estimator performance test results can be important visual inputs to the medical staff. Requests for additional CBG measurements or the sounding of an alarm in the event glucose estimates exceed critical limits may also prove useful. In one embodiment, a reduced number of outputs is provided in relatively small estimators for individual use.
Figure 6B illustrates one embodiment of a residual test data process in the patient health monitor. For example, the residual (y), covariance matrix (P), measurement matrix (H), and measurement noise matrix (R) are available from an estimator algorithm at each time epoch (i). If the estimator is operating properly, the sequence of residuals has the property of zero mean, white noise sequence, i.e., any two residuals taken at different times are uncorrelated (E[y(i)*y(j)] = 0, for all j not equal to i). This condition provides a unique means for constructing statistical tests.
Visually, a time history plot of the residuals, Sy(i) in Figure 6B, should appear random, zero mean, and unbaised. If they are summed over time, the deviation of the sum from zero should, on average, grow as the square root of time, as should its absolute value, ASy(i). The fact that E[y(i)*y(j)] = 0, for i not equal to j, also means that the sum of the variances, SV(i) in Figure 6B, which is easily computed, is equal to the variance of the sum of residuals for a properly performing estimator. The standard deviation of the sum of residuals, StdSV(i) also grows as the square root of time. Figure 6C illustrates one embodiment of a statistical test process in the patient health monitor. Figure 6C defines statistical tests which can be constructed based on filter residual test data. One of these compares the absolute value of the sum of the residuals with the standard deviation of the sum of residuals. On average, the ASy(i) should be bounded by the StdSV(i). If not, this indicates that the deviations of Sy(i) are growing faster than that of a white noise sequence, implying degraded estimator performance. Real-time displays of ASy(i) and StdSV(i) histories can provide a visual picture of estimator performance. An example is provided in Figures 7A and 7B which illustrate time history plots of residual test data with respect to measurement time. These two plots are one form of residual tests for the two sets of glucose estimation results provided in Figures 4A, 4B, 5A, and 5B respectively. In both cases, the sum of residuals are well behaved and was bounded by the Standard Deviation, Std, of the sum of residuals. The example embodiment for the glucose application of test data, test generation, and test result display described above is implemented by the algorithm in Table 4. Other real-time quantitative tests can also be constructed using these data. Tests on individual residuals can be performed using individual variances. If a measurement is received which causes the residual to exceed a four signal level, for example, the action might be to emit a warning and request an immediate CBG measurement. Other tests are identified in Figure 6C and will be evaluated as the system develops. Figure 6D illustrates one embodiment of an input/output interface in the patient health monitor. The patient health monitor provides a capability for the patient and/or staff to communicate and interact in real time with estimation and control processes using simple commands, visual displays, and audio outputs described above. Confidence levels in estimator performance can be established and a genuine interface established whereby the estimation and control processes could request additional information to check and insure confidence in estimates, physical sensors, and physical controllers. The staff can provide useful real-time inputs to augment this process.
Figure 8 is a block diagram of one embodiment of a Kalman optimal stochastic control solution as applied to physiological processes. In one embodiment, the application includes both optimal stochastic regulator and optimal stochastic tracking control solutions. Tracking control involves a controlled variable following a reference value, constant or dynamic, as closely as possible. Control is applied to physiological processes, wherein a control gain is computed based upon optimization criteria which minimizes controlled variable errors while minimizing application of control based on cost weightings.
There is a duality between the computation of a Kalman estimator gain and a Kalman control gain. The Kalman estimator gain minimizes estimation error variances. The Kalman control gain minimizes variances of error between the controlled variables and the reference variables while minimizing the level of control applied. As a result, the optimal control function includes specification of the controlled variables and their associated costs as well as costs associated with the amount of control to be applied. For example, smaller control variable error costs and larger control application costs will allow the controlled variable to deviate farther from the reference, but with reduced application of control. In one embodiment, linearization techniques described above in association with the Kalman estimator is applied to physiological nonlinear stochastic processes. For example, linearization about nominal values or about best estimates are provided by the Kalman estimator. Uncertain parameters associated with a controllable dispenser, or actuator, are included in an estimator state vector.
The optimal stochastic controller can use linear or nonlinear formulations and discrete or continuous time formulations. In one embodiment, the optimal stochastic controller is used with an optimal estimator described herein and/or an optimized decision and display function also described herein to form a closed loop system. The closed loop system works as an artificial pancreas when applied to a glucose problem in one embodiment.
Figure 8 is a functional diagram showing a controllable dispenser, or pump, with the capability to secrete insulin and glucagon to control high and low glucose levels, respectively. In one embodiment of a closed loop system, an estimator and a controller share a state vector wherein the estimator estimates it and the controller controls designated elements of it. Unlike the estimator, the time varying gain computations for the controller is computationally intensive and may not be used in all applications. The controller includes time varying or steady state gain formulations.
Figure 9A illustrates one embodiment of an optimal controller for the artificial pancreas. The glucose control problem is a tracking problem since the glucose level is controlled to a desired level which may either be constant or a function of time. In one embodiment, the dynamic process model includes a quantitative description of how glucose levels propagate, in time, as a function of insulin/glucagon secreted by the controllable pump. The dynamic process model is described by a system of first order differential or difference equations. Excluding the pump, much of the modeling for the glucose control is available through the estimator modeling development. Figure 9B illustrates one embodiment of a control model for a controller. In an example embodiment, the dynamic process model is described by three first order differential equations that are forced by a control variable, u. The dynamics of this process include a first order time lag and a scale factor associated with the pump, and a first order time lag and a scale factor associated with the physiological process. The model for the glucose sensors is the same as in the example glucose estimator embodiment. In the definition of the elements of the transition matrix, alpha g, s, and d are the inverse of the first order time lags associated with glucose, glucose sensor scale factor, and pump, respectively. Beta g is the glucose scale factor multiplied by alpha g. Delta t is the time interval between measurement/control application epochs. The controlled error is the difference between the estimate of glucose (provided by the estimator) and the glucose control point input.
In a controller command. Beta d is a pump scale factor multiplied by alpha d. The measurements are the same as for the example glucose estimator, except that the measurements contain an additional zero since the measurements are not functions of the variable insulin/glucagon.
The cost function, which is minimized by the optimal stochastic control, contains costs associated with the glucose error and the application of control. Choosing a value of Cu that is much larger than Cg will result in a relative gentle application of control. Another example embodiment would utilize a higher order model in which the first derivative of glucose would be included in the state vector, and included as a control variable. If a relatively large cost is associated with the first derivative relative to the glucose control point error, then the control will be very active when rapid changes in patient glucose occur while relatively gentle otherwise. Other example embodiments would include pump scale factor as a state vector element (uncertain parameter) as well as insulin/glucagon measurements to the estimator. Figure 9C illustrates one embodiment of a control algorithm based on the model defined in Figure 9B. Using either a time varying gain or a pre-computed steady state gain, control is applied at each epoch based upon a difference between an estimate of a patient glucose level and a glucose control point. A control, u, is applied to a pump over each time interval, and a state vector is time updated as indicated in this figure using the control variable, u, as a forcing function. Other time updates are performed in accordance with embodiments discussed above in association with estimator equations.
A potential problem in the application of closed loop control to that of physiological processes is due to potentially long time delays that may be de-stabilizing. These delays can be in the form of transport lags. A transport lag is the time between when control is applied and when the process action begins. An optimal technique for control using delay states is discussed in "Optimal Control of Linear Stochastic Systems with Process and Observation Time Delays" by E. J. Knobbe (Academic Press, Inc., 1989) and is hereby incorporated herein in its entirety by reference thereto. A discussion of principles developed for optimal control with process and observation time delays is provided in Table 5.
Although described above in connection with particular embodiments of the present invention, it should be understood that the descriptions of the embodiments are illustrative of the invention and are not intended to be limiting. Various modifications and applications may occur to those skilled in the art without departing from the true spirit and scope of the invention.
TABLE 1
SUMMARY OF roNTJNUOUSΦISCRFTE EXTENDED KALMAN FIWER
Figure imgf000018_0001
TABLE 2
SUMMARY OF C0 ΪINUOUS-D1SCRETB LINEARIZED KALMAN FILTER
Figure imgf000019_0001
TABLE 3
State Vector '
X, Xe, Xe(0) : These symbols denote, respectively, the state vector, the best estimate (or optimal estimate) of the state vector, and the initial best estimate of the state vector. They contain both variables and parameters as defined above and, with regard to estimation, there are no distinctions.
Dynamic Process Model dX= f(X,t) , dXe=f(Xe,t) : These symbols denote the time derivative of the state vector and estimated state vector, respectively. They represent a system of first order differential equations (or difference equations) which describe the manner in which the state elements propagate in time. The function, f, may be linear or nonlinear.
Measurement Model
Ym=h(X) + e, Ye=h(Xe) : These symbols denote the actual measurement and the best estimate of the measurement, respectively. The function, h, defines the arbitrary, but known, way in which the state vector elements are related to the measurement, and e represents the sensor measurement error.
Covariance Matrix
P, P(0) : Denotes, respectively, the matrix of variances and co-variances associated with the errors of each of the state variable estimates (within the estimated state vector) and their initial values,
Transition Matrix
A, A(0) : Denotes the transition matrix, which is used to propagate the covariance matrix forward it time (time update), and it's initial value, respectively; in the linear case it is also used to time update the state vector estimate.
Process Noise Matrix
Q : Denotes the matrix of variances and co-variances associated with error growth uncertainty accumulated in the state variable estimates since the last measurement update,
Measurement Matrix
H, H(0) : Denotes the measurement matrix and it's initial value, respectively. This matrix defines the linear functional relationship between the measurement and the state vector elements. If, the measurement model is linear, then H = h; if h is non-linear, then H is defined by linearizing h using partial derivatives or perturbation techniques.
Measurement Noise Matrix
R : Denotes the matrix of variances and co-variances associated with measurement error uncertainties.
Kalman Gain Matrix
K : Denotes the Kalman gain matrix which, when multiplied by the difference between the actual measurement, and the best estimate of the measurement, yields the estimated state correction. This estimated state correction, when added to the old best estimate, becomes the new best estimate.
Estimated Measurement Error (Residual) y = Ym-Ye : Denotes the difference between the actual measurement and the best estimate of the measurement. This difference multiplied by the Kalman gain yields the correction to the previous state vector estimate. TABLE 4 % Optimal Glucose Estimator, Extended Kalman Filter; Computes Statistical Test Data for Patient Health Monitor %Ed_JTP_Estimator; Two-step scale-factor model, Exp = .999 for t>20; r=15Λ2 for y<10 ;
x=[150;.25J; hh=[l 0] ,-p=fl00Λ2 0 ; 0 .lΛ2];a=U 0 ; 0 1] ; q=[20Λ2 0 ; 0 .002Λ2J; rr= 15Λ2; r= 5Λ2; I=[l 0;0 1] ; load EdJYS28.prn,
[m,n]=size(Ed Y528) , t=Ed Y528(:,l)/12; y=EdJY528(:,2) ; g=EdJY528(:,3) ; grα=EdJY528(:,4); x=[gm(l);.25]; sz=0.0;
SVZ = p[l,l)*x{2)Λ2 + p{2,2)*x{l)Λ2; for i = l:m; . - - r = 5Λ2;... zz = 0; ... if y(i) < 10 r = 15Λ2;... end; if t(i) > 20 a=[l 0 ; 0 .999];... end; p=a*p*a'+q; ... x=a*x; ... if gm(i} > 0 k=p*hh't{hhVhh'+rr);, zz=(grα(i)-x(l)); x=x+k*zz; p=(I-k*hh)*p;... end;
.
Figure imgf000021_0001
. az=abs (z) ; ... sz=sz+z; ... asz = abs(sz);...
VZ = p{l,l)*x(2)Λ2 + p(2,2)*x(l)Λ2 + 2* (1,2) *x(l) *x(2)+ r sigZ = sqrt (VZ) ; ...
SVZ = SVZ + VZ; sigSZ = sqr (SVZ);
VZZ = p(l,l) + rr;... sigZZ = sqrt (VZZ); azz = abs (zz) ; . . . x istory (i,l) = x(l);... xhistory (i,2) = x(2);... gmhistory =gm; ... ghistory =g; ... zhistory(i) =z; ... sigsz_history (i) = sigSZ;... sigz_history (i) = 2*sigZ;... aszhistor (i) = asz; ... azhistory(i) = az; ... zzhistory(i) = zz; ... sigzz_history (i) = 2*sigZZ;... azzhistory(i) = azz;... sig_histor (i, l)=sqr (p (1,1) ) ; • • • sig_histor (i,2)=sqrt (p{2,2) ) ; ... yhistory=y; ... end; figure (1) plot (t, xhistory ( : ,2) ) ; xlabel ( ' Measurement Time (hrs)'); ylabel ( ' Scale Factor Estimate; (ISIG Units/ (mg/dl) ) * ) ; title ('Scale Factor Estimate vs Tirαe'); figure (2) plot (t, gmhistory ( : ) , fc, xhistor ( :,1) ) ; xlabel ( ' Time (hrs ) ' ) ;
Ylabel ( 'Glucose Meas & Estimated Glucose; (mg/dl)'); title ( 'All Glucose Measurements & Estimated Glucose vs Time'); figure (3) plot (t,yhistory( : ) , t, zhistory ( : ) ) ; xlabel ('Time (hrs)');
Ylabel ('ISIG Meas & z Residual; (ISIG Units) '); title ('ISIG Measurement & z residual vs Time'); figure (4) plot (t,xhistory ( : ,1) , t, ghistory( : ,1) ) ; xlabel ( ' Time (hrs)'); ylabel (' Glucose & Estimated Glucose; (mg/dl)'); title ( 'Glucose Estimate & All Glucose CBGs vs Time'); figure (5) plot (t, sigzz_history( : ) , t,azzhistory ( : ) ) ; xlabel ( ' Time (hrs) ' ' ylabel ( 'Abs (zz) Residual S two sigma(zz); (mg/dl)'); title ( 'Residual Test: 2 Sigma (zz) & Abs(zz) Residual vs Time ' ) ; figure (6) plot (t,azhistor ( : ) , t, sigz_history ( : ) ) ; xlabel ( ' Measurement Time (hrs) ' ) ; ylabel ('z Residual & Standard Deviation; (ISIG Units) title ( 'Residual Test : Abs(z) & 2 Sigma (z) vs Time'); figure (7) plot (t, aszhistory ( : ) ,t, sigsz_history( : ) ) ; xlabel ( ' Measurement Time (hrs)'); ylabel ( ' sz Residual & Sigma; (ISIG Units) '); title ( 'Residual Test : Abs(sz) S Sigma(sz) vs Time'); TABLE 5
PROBLEM DEFINITION
The general, discrete-time system representation with explicit process and observation time delays is given by
X(k + 1) = f, k)X k + ! - .) + B(k)U(k) + W,(k), (1)
γ k) = c + ' - ) + w2 , (2) where p and q are integers > 1; Wγ and W2 are zero-mean white-noise sequences such that
'I.JΛ'
Figure imgf000023_0001
Y is an m x 1 observation vector; X is an n x 1 random state vector whose initial uncertainty is uncorrelated with W, and W2 and with initial covariance Q0; and U is the control input vector. The objective is to find the control function (functional) U(k) for k = 1,2,... that minimizes an expected quadratic cost function for the linear stochastic regulator defined by Eqs. (1) and (2). Because linear stochastic tracking problems can be formulated as linear stochastic regulator problems by combining the reference and plant models in an augmented system [3], this system representation applies equally to the tracking problem.
PROBLEM FORMULATION AND SOLUTION
The system representation defined by Eqs. (1) and (2) can be cast in stochastic regulator form by augmenting the state vector with the time- delayed states, that is,
X(k + 1) = Mk)X(k) + B{k)V(k) + Wt(k), (3) and
Y(k) = C(k)X(k) + W2(k) (4)
0 / 0 0
X(k - h + 1) 0 X(k - h + 2)
A(k) X(k)
0 (* - 1)
0 ••• 0 / X(k)
Ah(k)
Figure imgf000024_0001
and the dimensionality of the system is defined by h, where h = max(p, q).
In Eq. (3), the control remains the same as in Eq. ( 1) because, physically, past states cannot be changed or controlled. In Eq. (A) the observation Y and the observation noise W2 also remain unchanged.
In Eqs. (1) and (2), p and q are not, in general, equal. For example, if the observation is not a linear function of all the delayed states contained in X, then h = q > p, and the appropriate (q — p) submatrices in C. of Eq. (4), arc set equal to zero. If the process evolution is not a linear function of all the delayed states contained in X, then h — p > q, and the appropriate submatrices in A, of Eq. (3), are set equal to zero. The dimensions of B and Wj , in Eq. (3), must also be consistent with the integer h. Therefore, without loss of generality, the system defined by Eqs. ( 1 ) and (2) can be represented by the augmented system model described by Eqs. (3) and (4).
Given that the control is to minimize the expected value of a quadratic cost function of the form,
Figure imgf000024_0002
then the separation principle applies [3]. Also, the optimal linear stochastic control of the augmented system is given by a deterministic, optimal linear controller with state input X (or estimated state feedback), which is provided by an optimal one-step predictor using the augmented model [4]. That is,
U(k) = -F(k)X(k); k = k0, k0 + !,...,*„ (6) where the control gain F satisfies
F(k) = {R2(k) + EFWRtlk + 1) + P(k + l]B(A)}-lBT(k)[R(k + 1)
+ P(k + \)}A(k) (?) the matrix P satisfies the recursive matrix Riccati equation,
P(k) = ^ k^R^k + I) + P(k + l][ ϊ(fc) - B(k)F(k)l, (8) the one-step predictor output X is defined by
X(k + 1) = A(k)X(k) + B(k)V(k) + K(k)[Y(k) - C(k)X(k)ϊ, (9) the estimator gain K satisfies
K(k) = A(k)Q(k)C k)ZC(k)Q(k)Cτ(k) + V2(k)T *; (10) and the state estimation error covariance matrix Q satisfies the recursive matrix Riccati equation,
Q(k + 1) = lA ) - R(k)C(mQ(k)F(k) + VM (11)
The final value of P used to "initialize" Eq. (8) (which is solved backward in time) is the final value defined in the quadratic cost function of Eq. (5), that is, P(kι) — P,. The initial value of Q used to initialize Eq. (11) is the error covariance of the initial estimate of X, that is, Q(0) = Q0.
If the system statistics defined in Eqs, (1) and (2) are gaussian, then the above solution is the optimal solution without qualification; it not, then it is the optima] linear control solution. The expected system performance is determined by analyzing the augmented system as a linear, stochastic regulator problem.
PRACTICAL APPLICATIONS
From control gain Eqs. (7) and (8) and estimator gain Eqs. (10) and (11), it can be shown that the required dimensions of the controller and estimator are not, in general, equal. The dimensions of the control matrix Riccati equation (8) is determined by the number of delayed states in the process evolution, that is, P has the dimensions of n • q x n • q. Consequently, for the special case of measurement delays only (q = 1), the controller implementation is unaffected. However, the dimensions of the estimator matrix Riccati equation (11) are determined by the maximum of p and q. Since h = max(p, q), then Q has the dimensions n • ft x n • h. For the general case, the optimal control is given by
V{k) = - £ Ft(k)X(k + 1 - 0, (12) where
FW & ZFJk) F,.^) ... F2(k) F.(*)], (13) and Xτ(k) A [X T(fe + i - ft) ... Jrτ(fc + l - «) ••• Xτ(k - Ϊ) Xτ(k)
(14)
We note that, in Eq. (14), X(k) is the one-step predicted value of the original system state vector; X(k — 1) is the filtered value; X(k — 2) is the one-step smoothed value; and, finally, X (k + 1 — ft) is the (ft — 2)th smoothed value.
Typically, the real-time computational requirements associated with the implementation of time-varying optimal control are always stressing due to the "backward-in-time" recursion that is required to obtain solutions P of Eq. (8). However, in many high-accuracy applications, the A, B, Rt, and R2 matrices of the augmented system and cost function can be treated as time invariant over the time intervals of interest. Further, if the augmented system satisfies the relatively minor requirements of stabilizability and detectability, then the control gain F will converge to a unique value such that the steady- state optimal control law is time invariant, asymptotically stable, and minimizes the quadratic cost function of Eq. (5) as &, → oo. For this case, the steady-state gain matrix FB can be computed off-line and stored for realtime use; the real-time computations required to implement this steady-state control law are negligible.
Usually, for tracking and regulator problems, the steady-state control gains and time-varying gains are such that the initial value of the time-varying gain is equal to the steady-state gain but is less than or equal to the steady-state gain as time progresses, that is,
F„ ≥ F(k) for k0 < k ≤ kλ.
As a consequence, the steady-slate gain tends to maintain the controlled state closer to the estimated state, but at the cost of control energy. If accuracy is the significant criterion, then the steady-state gain is not only easy to implement but also provides essentially equivalent or better accuracy.
This, however, is not the case for the estimator gain, even though in most cases where linear steady-state optimal stochastic control is implemented, both the steady-state control and the estimator gains arc used. For the estimator, and the same general conditions as before, the initial time-varying estimator gain is usually significantly larger than the steady-state gain, to account for initial uncertainties in the knowledge of the system state. As time progresses, the time-varying gain converges to the steady-state gain, that is,
Kα ≤ K{k) for k0 ≤ k ≤ kι. Consequently, initial system performance is significantly degraded if „ is used. In fact, if applied to a "linearized" system, the estimator can actually diverge given the initial small steady-state filter gains.
Hence a good compromise between performance and computational complexity is to choose the steady-state controller with the time-varying estimator. Further, because the estimator matrix Riccati equation is solved forward in time, the computations associated with the time-varying filter gain are orders of magnitude less than with the time-varying control gain and can usually be implemented in real time.
Although the computations associated with the augmented system estimator are significantly increased because of the increased dimensions, some simplifications can be made. The estimator Riccati equation can be separated into a time update and a measurement update, where the time update for the augmented system becomes primarily one of data transfer. If Q, from Eq. (11), is defined as
Figure imgf000027_0001
then, for the special case where q = 1 and p > 2, the time update becomes β i,j = l,.-.,p - 1 β 2,, P(k + l) = Qi+i.p(k)A](k) (16) Qp k + 1) = At(k)Q,Jik)A](k) A- V,(k)
For the worst case, where q > p > 1, the time update is given by
Qt k + 1) = β,+ w+,(*λ U = l,..-,9 - J.
Qi k + 1) = t Quj k)A + l)^ kl i = },..., q - I,
(17)
.8. β* + - ∑ ∑ W+ l)-l(*)βt/(*). T + »_,(*) + K,(fc)
1= 1 j= l
In practice, another computational simplification results because rarely, if ever, is the actual system measurement a function of every time-delayed state element, and likewise for the actual system process model. Hence the submatrices At and (for i > 1) in A and C, respectively, are usually of significantly reduced dimension. This, in turn, significantly reduces the dimensions of the augmented system model.
We note that the iteration interval Δt (implied in the discrete-time system representation) or some integer number of Δf's should be set equal to the time delay. For variable time delays, it may be advantageous to use (a) variable iteration intervals for the estimator, (b) fixed iteration intervals for the controller, and (c) variable-time updates of the state estimate to time synch the estimator output with the controller input
Finally, for the analogous continuous-time problem, the filter equation becomes a partial differential equation with a boundary condition, and the covariance equation becomes a partial differential matrix equation with three boundary conditions (see [2]). Consequently, the most practical control solution is obtained by discretizing the continuous-time system representation and then applying the approach of Section III. One technique for discretizing the continuous-time representation is by using the Z transform method; see, for example, [5]. Using a scalar differential equation with one delayed state as an example, we have
A-(.) = ∑ fl.-X(f - iΔ.) + l/( . (18)
If we assume that the output solution X(t) is sampled in discrete time, and if the system process can be reasonably approximated by a continuous system in which X is driven by the output of a high-frequency sampling of the right side of Eq. (18), then the discrete-time process model is given by
X(k + !) = (! + a0άt)X(k) + (α, At)X(k - At) + (Δ.)U(fc). (19)
The approximate discrete-time solution to Eq. (18), which is given by Eq. (19), is now in a form consistent with Eq. (1) of Section II, and the approach of Section HI can be applied.

Claims

WHAT IS CLAIMED IS:
1. A method for providing a best estimate of glucose level in real time comprising the acts of: obtaining a measurement which is a function of glucose level, wherein noise associated with the measurement is within limits of a prede ined measurement uncertainty; supplying the measurement to an extended Kalman filter in real time, wherein the extended Kalman filter has a dynamic process model, a dynamic measurement model, a state vector with at least one element corresponding to glucose level, and an error covariance matrix of the state vector; and determining the best estimate of glucose level in real time using the extended Kalman filter.
2. The method of Claim 1, wherein the extended Kalman filter is implemented using a software algorithm.
3. The method of Claim 1, wherein determination of the best estimate of glucose level in real time uses a recursive process comprising the acts of: computing a current estimate of the state vector using a preceding best estimate of the state vector and the dynamic process model; computing a current error covariance matrix of the state vector using a preceding error covariance matrix of the state vector, uncertainties associated with the dynamic process model, and the dynamic process model linearized about the current estimate of the state vector; computing a Kalman gain using the current covariance matrix, uncertainties associated with the dynamic measurement model, and the dynamic measurement model linearized about the current estimate of the state vector; computing a new error covariance matrix of the state vector using the current error covariance matrix of the state vector, the Kalman gain, and the dynamic measurement model linearized about the current estimate of the state vector; and computing a new best estimate of the state vector using the current estimate of the state vector, the Kalman gain, the measurement, and the dynamic measurement model.
4. The method of Claim 1, wherein development of the dynamic process model and the dynamic measurement model is an iterative process comprising the acts of: defining estimation variables and uncertain parameters which become elements of the state vector; defining nominal time propagation of the estimation variables and the uncertain parameters; defining a nominal sensor measurement model; defining relationships between the estimation variables, the uncertain parameters, and the measurement; defining uncertainties associated with the estimation variables, the uncertain parameters, and the measurement; and verifying a nominal dynamic process model and a nominal measurement model.
5. The method of Claim 4, wherein a database of measurements is used to empirically verify the nominal dynamic process model and the nominal measurement model.
6. The method of Claim 5, wherein the database has a plurality of sensor measurements with corresponding direct measurements.
7. A real-time glucose estimator comprising: a plurality of measurement inputs, wherein at least one of the measurement inputs is configured to receive an input indicative of glucose level; a plurality of control inputs; and an extended Kalman filter algorithm configured to receive the plurality of measurement inputs and the plurality of control inputs to provide an optimal estimate of glucose level in real time.
8. The real-time glucose estimator of Claim 7, wherein at least two independent glucose sensors are coupled to the measurement inputs respectively.
9. The real-time glucose estimator of Claim 8, wherein at least two independent glucose sensors are of different types.
10. The real-time glucose estimator of Claim 7, wherein at least one of the measurement inputs is configured to receive a direct measurement of glucose.
11. The real-time glucose estimator of Claim 10, wherein the direct measurement is derived from a capillary blood glucose measurement.
12. The real-time glucose estimator of Claim 7, wherein at least one of the control inputs is configured to receive configuration data.
13. The real-time glucose estimator of Claim 7, wherein at least one of the control inputs is configured to receive environment data.
14. The real-time glucose estimator of Claim 7, wherein the extended Kalman filter algorithm produces the optimal estimate of glucose level in real time using a recursive loop formed by a time update module and a measurement update module.
15. The real-time glucose estimator of Claim 14, wherein the plurality of control inputs is processed in real time by the time update module.
16. The real-time glucose estimator of Claim 14, wherein the plurality of measurement inputs is processed in real time by the measurement update module.
17. The real-time glucose estimator of Claim 16, wherein measurement inputs of a substantially identical time interval are processed serially in one cycle of the measurement update module.
18. The real-time glucose estimator of Claim 10, wherein the direct measurement is provided during initialization of the real-time glucose estimator and intermittently thereafter.
19. The real-time glucose estimator of Claim 7, wherein at least one of the measurement inputs is provided periodically to the real-time glucose estimator.
20. The real-time glucose estimator of Claim 7, wherein at least one of the measurement inputs is coupled to a non-invasive glucose sensor.
21. The real-time glucose estimator of Claim 7 further comprising a patient health monitor configured to accept user inputs, to derive configuration data based on the user inputs, and to provide the configuration data in real time to the real-time glucose estimator.
22. The real-time glucose estimator of Claim 7 further comprising a patient health monitor configured to display in real time the optimal estimate of glucose level and to provide an audible alarm when the optimal estimate of glucose level is outside a predetermined range.
23. The real-time glucose estimator of Claim 10 further comprising a patient health monitor configured to display in real time a performance status for the real-time glucose estimator, wherein the performance status is based upon results of statistical testing performed on residuals of the extended Kalman filter.
24. The real-time glucose estimator of Claim 23, wherein the patient health monitor outputs a request for the direct measurement of glucose when the performance status indicates the optimal estimate of glucose level is outside predetermined limits.
25. An estimator for monitoring a physiological parameter comprising: a sensor which outputs a measurement as a function of the physiological parameter; an electronic processor coupled to an output of the sensor, wherein the electronic processor executes an algorithm that implements an extended Kalman filter to estimate the physiological parameter in real time; and an interface coupled to an output of the electronic processor to display the estimate of the physiological parameter in real time.
26. The estimator of Claim 25, wherein the interface accepts user inputs and provides a control signal in real time to the electronic processor based on the user inputs.
27. The estimator of Claim 26, wherein the user inputs describe environmental conditions of the sensor in real time.
28. The estimator of Claim 26, wherein the user inputs describe activities of a patient using the estimator in real time.
29. The estimator of Claim 26, wherein the user inputs describe real-time administration of medication to a patient using the estimator.
30. The estimator of Claim 25, wherein the estimator is a portable device contained in a relatively small package and battery operated.
31. An estimator comprising: means for obtaining a measurement related to a physiological parameter; and means for processing the measurement in real time using a linearized Kalman filter algorithm to provide a real-time estimate of the physiological parameter.
32. The estimator of Claim 31 further comprising means for displaying the real-time estimate of the physiological parameter in real time.
33. The estimator of Claim 31 further comprising means for updating the linearized Kalman filter algorithm in real time with user inputs.
34. The estimator of Claim 31 further comprising means for indicating performance status of the estimator in real time.
35. A method for controlling glucose level in real time comprising the acts of: receiving an estimated glucose level; receiving a reference signal indicative of a desired glucose level; providing the estimated glucose level and the reference signal to a Kalman control algorithm in real time; determining a control command in real time using the Kalman control algorithm; and providing the control command to a dispenser which outputs medication in response to the control command.
36. The method of Claim 35, wherein the estimated glucose level is provided by an optimal estimator implemented using an extended Kalman filter.
37. The method of Claim 35, wherein the reference signal is provided by a patient health monitor which accepts inputs from a user.
38. The method of Claim 35, wherein the reference signal varies with time.
39. The method of Claim 35, wherein the Kalman control algorithm has a dynamic process model forced by the control command and a cost function determining a relative level of control.
40. The method of Claim 39, wherein determination of the control command in real time comprises the acts of: computing a Kalman control gain to minimize the cost function; and adjusting the control command based on the Kalman control gain and a difference between the estimated glucose level and the desired glucose level.
41. The method of Claim 35, wherein the dispenser secretes insulin or glucagon in response to the control command to correct a relatively high or a relatively low estimated glucose level.
42. A method for close-loop control of a physiological parameter comprising the acts of: obtaining a measurement of the physiological parameter from a patient; providing the measurement to an optimal estimator in real time, wherein the optimal estimator outputs a best estimate of the physiological parameter in real time based on the measurements; providing the best estimate of the physiological parameter to an optimal controller in real time, wherein the optimal controller outputs a control command in real time based on the best estimate of the physiological parameter and a control reference; and providing the control command to an actuator, wherein the actuator provides an output to adjust the physiological parameter.
43. The method of Claim 42, wherein the measurement is obtained using a sensor.
44. The method of Claim 42, wherein the optimal estimator is implemented using a linearized Kalman algorithm.
45. The method of Claim 44, wherein the optimal controller is implemented using a Kalman control algorithm.
46. The method of Claim 45, wherein the optimal estimator and the optimal controller have substantially identical dynamic process models forced by the control command.
47. The method of Claim 46, wherein the optimal estimator provides a best estimate state vector to the optimal controller, and the best estimate of the physiological parameter is an element of the best estimate state vector.
48. The method of Claim 46, wherein the optimal controller provides the control command to the optimal estimator.
49. The method of Claim 42 wherein the optimal estimator and the optimal controller are implemented as a joint Kalman algorithm.
50. The method of Claim 42, wherein the control reference is provided by a patient health monitor.
51. The method of Claim 42, wherein the optimal estimator outputs best estimates of additional physiological parameters to the optimal controller in real time, and the optimal controller controls the additional physiological parameters by outputting additional control commands.
52. A real-time optimal glucose controller comprising: a first input configured to receive an estimated glucose level in real time; a second input configured to receive a reference glucose level; a Kalman control algorithm configured to determine a control command based on the estimated glucose level and the reference glucose level, wherein the Kalman control algorithm has a dynamic process model forced by the control command and a cost function defining a desired level of control; and an output configured to provide the control command to a pump, wherein the pump provides medication in response to the control command to minimize a difference between the estimated glucose level and the reference glucose level.
53. An artificial pancreas for controlling glucose level in real time comprising: a glucose sensor to provide a measurement of the glucose level; an optimal glucose estimator, wherein the optimal glucose estimator uses a stochastic model to describe a physiological process relating to the glucose level and uses a linearized Kalman filter to estimate the glucose level in real time based on the measurement from the glucose sensor; an optimal glucose controller, wherein the optimal glucose controller uses a substantially identical stochastic model as the optimal glucose estimator and uses a Kalman control algorithm to determine a control command to adjust the glucose level in real time; and a medical dispenser to provide medication to a patient in response to the control command.
54. The artificial pancreas of Claim 53, wherein the optimal glucose estimator and the optimal glucose controller are implemented using a software algorithm.
55. The artificial pancreas of Claim 53, wherein the artificial pancreas is a portable device.
56. The artificial pancreas of Claim 53 further comprising one or more additional sensors of different types which operate independently to provide respective glucose measurements.
57. The artificial pancreas of Claim 53, wherein the medical dispenser secretes insulin and glucagon to control a relatively high glucose level and a relatively low glucose level respectively.
58. The artificial pancreas of Claim 53 further comprising a patient health monitor with an input/output interface to receive inputs from a user and to display status of the artificial pancreas.
PCT/US2001/029649 2000-09-22 2001-09-21 Method and apparatus for real-time estimation and control of pysiological parameters WO2002024065A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AU2001291189A AU2001291189A1 (en) 2000-09-22 2001-09-21 Method and apparatus for real-time estimation and control of pysiological parameters

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US23463200P 2000-09-22 2000-09-22
US60/234,632 2000-09-22

Publications (1)

Publication Number Publication Date
WO2002024065A1 true WO2002024065A1 (en) 2002-03-28

Family

ID=22882159

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2001/029649 WO2002024065A1 (en) 2000-09-22 2001-09-21 Method and apparatus for real-time estimation and control of pysiological parameters

Country Status (3)

Country Link
US (3) US6572545B2 (en)
AU (1) AU2001291189A1 (en)
WO (1) WO2002024065A1 (en)

Cited By (63)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1518495A1 (en) * 2003-09-23 2005-03-30 Roche Diagnostics GmbH Method and device for continuous monitoring of analyte concentration
WO2006024672A1 (en) * 2004-09-03 2006-03-09 Novo Nordisk A/S System and method for estimating the glucose concentration in blood
WO2007138154A1 (en) * 2006-05-29 2007-12-06 Wristop Technologies Oy Apparatus and method for dosing drug and wireless remote control of a drug pump
DE102007046864A1 (en) * 2007-09-28 2009-04-09 Endress + Hauser Conducta Gesellschaft für Mess- und Regeltechnik mbH + Co. KG Measured value plausibility and quality evaluating method for e.g. control circuit, involves forming and evaluating derivative of measured value according to time, where measured value is compared with predefined warning limit value
US7774145B2 (en) 2003-08-01 2010-08-10 Dexcom, Inc. Transcutaneous analyte sensor
US7778680B2 (en) 2003-08-01 2010-08-17 Dexcom, Inc. System and methods for processing analyte sensor data
US7927274B2 (en) 2003-11-19 2011-04-19 Dexcom, Inc. Integrated receiver for continuous analyte sensor
WO2011092133A1 (en) * 2010-01-26 2011-08-04 Roche Diagnostics Gmbh Methods and systems for processing glucose data measured from a person having diabetes
US8010174B2 (en) 2003-08-22 2011-08-30 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US8808228B2 (en) 2004-02-26 2014-08-19 Dexcom, Inc. Integrated medicament delivery device for use with continuous analyte sensor
US8840553B2 (en) 1998-04-30 2014-09-23 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8886273B2 (en) 2003-08-01 2014-11-11 Dexcom, Inc. Analyte sensor
US8917184B2 (en) 2008-03-21 2014-12-23 Lifescan Scotland Limited Analyte testing method and system
US8915850B2 (en) 2005-11-01 2014-12-23 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8920319B2 (en) 2005-11-01 2014-12-30 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8926585B2 (en) 2004-02-26 2015-01-06 Dexcom, Inc. Integrated delivery device for continuous glucose sensor
US8968198B2 (en) 2006-02-22 2015-03-03 Dexcom, Inc. Analyte sensor
US8974386B2 (en) 1998-04-30 2015-03-10 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8977504B2 (en) 2007-10-12 2015-03-10 Cambridge Enterprise Limited Substance monitoring and control in human or animal bodies
US8986209B2 (en) 2003-08-01 2015-03-24 Dexcom, Inc. Transcutaneous analyte sensor
US9011332B2 (en) 2001-01-02 2015-04-21 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
CN104602599A (en) * 2012-11-01 2015-05-06 泰尔茂株式会社 Sensing device and sensing method
US9060742B2 (en) 2004-07-13 2015-06-23 Dexcom, Inc. Transcutaneous analyte sensor
US9066695B2 (en) 1998-04-30 2015-06-30 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
EP2748747A4 (en) * 2011-08-26 2015-07-08 Univ Virginia Patent Found Method, system and computer readable medium for adaptive advisory control of diabetes
US9078607B2 (en) 2005-11-01 2015-07-14 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9078608B2 (en) 2005-03-10 2015-07-14 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US9107623B2 (en) 2003-12-09 2015-08-18 Dexcom, Inc. Signal processing for continuous analyte sensor
US9135402B2 (en) 2007-12-17 2015-09-15 Dexcom, Inc. Systems and methods for processing sensor data
US9149233B2 (en) 2007-12-17 2015-10-06 Dexcom, Inc. Systems and methods for processing sensor data
US9247900B2 (en) 2004-07-13 2016-02-02 Dexcom, Inc. Analyte sensor
US9282925B2 (en) 2002-02-12 2016-03-15 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9420968B2 (en) 2003-08-22 2016-08-23 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9446194B2 (en) 2009-03-27 2016-09-20 Dexcom, Inc. Methods and systems for promoting glucose management
US9451908B2 (en) 2006-10-04 2016-09-27 Dexcom, Inc. Analyte sensor
CN106139318A (en) * 2016-06-28 2016-11-23 成都市亿泰科技有限公司 A kind of intelligent transfusion monitoring method based on integrated kalman filter
US9717449B2 (en) 2007-10-25 2017-08-01 Dexcom, Inc. Systems and methods for processing sensor data
WO2017129634A1 (en) * 2016-01-25 2017-08-03 Prediktor Medical As Estimating a substance concentration using multiple non-invasive sensors
US9741139B2 (en) 2007-06-08 2017-08-22 Dexcom, Inc. Integrated medicament delivery device for use with continuous analyte sensor
CN107203700A (en) * 2017-07-14 2017-09-26 清华-伯克利深圳学院筹备办公室 A kind of method and device monitored based on continuous blood sugar
US9895089B2 (en) 2003-08-01 2018-02-20 Dexcom, Inc. System and methods for processing analyte sensor data
CN108375381A (en) * 2018-02-08 2018-08-07 北方工业大学 A kind of bionical polarization sensor multi-source error calibrating method based on Extended Kalman filter
US10297350B2 (en) 2016-06-01 2019-05-21 Roche Diabetes Care, Inc. Risk-based control-to-range
US10311976B2 (en) 2016-04-28 2019-06-04 Roche Diabetes Care, Inc. Bolus calculator with probabilistic carbohydrate measurements
US10332632B2 (en) 2016-06-01 2019-06-25 Roche Diabetes Care, Inc. Control-to-range failsafes
US10332633B2 (en) 2016-06-01 2019-06-25 Roche Diabetes Care, Inc. Control-to-range aggressiveness
US10478108B2 (en) 1998-04-30 2019-11-19 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US10478556B2 (en) 2016-03-04 2019-11-19 Roche Diabetes Care, Inc. Probability based controller gain
US10575790B2 (en) 2016-03-02 2020-03-03 Roche Diabetes Care, Inc. Patient diabetes monitoring system with clustering of unsupervised daily CGM profiles (or insulin profiles) and method thereof
US10653835B2 (en) 2007-10-09 2020-05-19 Dexcom, Inc. Integrated insulin delivery system with continuous glucose sensor
US10980461B2 (en) 2008-11-07 2021-04-20 Dexcom, Inc. Advanced analyte sensor calibration and error detection
US11000215B1 (en) 2003-12-05 2021-05-11 Dexcom, Inc. Analyte sensor
CN113503879A (en) * 2021-07-09 2021-10-15 北京航空航天大学 Dynamic adaptive Kalman filter method based on ensemble empirical mode decomposition
US11331022B2 (en) 2017-10-24 2022-05-17 Dexcom, Inc. Pre-connected analyte sensors
US11350862B2 (en) 2017-10-24 2022-06-07 Dexcom, Inc. Pre-connected analyte sensors
US11382539B2 (en) 2006-10-04 2022-07-12 Dexcom, Inc. Analyte sensor
US11399745B2 (en) 2006-10-04 2022-08-02 Dexcom, Inc. Dual electrode system for a continuous analyte sensor
US11432772B2 (en) 2006-08-02 2022-09-06 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US11437141B2 (en) 2009-04-30 2022-09-06 Dexcom, Inc. Performance reports associated with continuous sensor data from multiple analysis time periods
US11559260B2 (en) 2003-08-22 2023-01-24 Dexcom, Inc. Systems and methods for processing analyte sensor data
EP3381364B1 (en) * 2015-11-25 2023-03-29 Nippon Telegraph and Telephone Corporation Respiratory estimation method and device
US11633133B2 (en) 2003-12-05 2023-04-25 Dexcom, Inc. Dual electrode system for a continuous analyte sensor
EP4276851A3 (en) * 2007-06-21 2024-01-24 University Of Virginia Patent Foundation Lqg artificial pancreas control system and related method

Families Citing this family (305)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6036924A (en) 1997-12-04 2000-03-14 Hewlett-Packard Company Cassette of lancet cartridges for sampling blood
US6391005B1 (en) 1998-03-30 2002-05-21 Agilent Technologies, Inc. Apparatus and method for penetration with shaft having a sensor for sensing penetration depth
US6949816B2 (en) 2003-04-21 2005-09-27 Motorola, Inc. Semiconductor component having first surface area for electrically coupling to a semiconductor chip and second surface area for electrically coupling to a substrate, and method of manufacturing same
US6430525B1 (en) * 2000-06-05 2002-08-06 Masimo Corporation Variable mode averager
US6572545B2 (en) * 2000-09-22 2003-06-03 Knobbe, Lim & Buckingham Method and apparatus for real-time control of physiological parameters
US8641644B2 (en) 2000-11-21 2014-02-04 Sanofi-Aventis Deutschland Gmbh Blood testing apparatus having a rotatable cartridge with multiple lancing elements and testing means
US7041468B2 (en) 2001-04-02 2006-05-09 Therasense, Inc. Blood glucose tracking apparatus and methods
US7316700B2 (en) 2001-06-12 2008-01-08 Pelikan Technologies, Inc. Self optimizing lancing device with adaptation means to temporal variations in cutaneous properties
US7344507B2 (en) 2002-04-19 2008-03-18 Pelikan Technologies, Inc. Method and apparatus for lancet actuation
CA2448905C (en) 2001-06-12 2010-09-07 Pelikan Technologies, Inc. Blood sampling apparatus and method
ATE497731T1 (en) 2001-06-12 2011-02-15 Pelikan Technologies Inc DEVICE FOR INCREASING THE SUCCESS RATE OF BLOOD YIELD OBTAINED BY A FINGER PICK
US8337419B2 (en) 2002-04-19 2012-12-25 Sanofi-Aventis Deutschland Gmbh Tissue penetration device
US7025774B2 (en) 2001-06-12 2006-04-11 Pelikan Technologies, Inc. Tissue penetration device
US7981056B2 (en) 2002-04-19 2011-07-19 Pelikan Technologies, Inc. Methods and apparatus for lancet actuation
ATE485766T1 (en) 2001-06-12 2010-11-15 Pelikan Technologies Inc ELECTRICAL ACTUATING ELEMENT FOR A LANCET
WO2002100254A2 (en) 2001-06-12 2002-12-19 Pelikan Technologies, Inc. Method and apparatus for lancet launching device integrated onto a blood-sampling cartridge
US9427532B2 (en) 2001-06-12 2016-08-30 Sanofi-Aventis Deutschland Gmbh Tissue penetration device
US9795747B2 (en) 2010-06-02 2017-10-24 Sanofi-Aventis Deutschland Gmbh Methods and apparatus for lancet actuation
US9226699B2 (en) 2002-04-19 2016-01-05 Sanofi-Aventis Deutschland Gmbh Body fluid sampling module with a continuous compression tissue interface surface
US7451126B2 (en) * 2001-10-18 2008-11-11 Omron Corporation State space navigation system, user system and business methods for machine to machine business
US7892183B2 (en) 2002-04-19 2011-02-22 Pelikan Technologies, Inc. Method and apparatus for body fluid sampling and analyte sensing
US7297122B2 (en) 2002-04-19 2007-11-20 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US9795334B2 (en) 2002-04-19 2017-10-24 Sanofi-Aventis Deutschland Gmbh Method and apparatus for penetrating tissue
US7547287B2 (en) 2002-04-19 2009-06-16 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US8579831B2 (en) 2002-04-19 2013-11-12 Sanofi-Aventis Deutschland Gmbh Method and apparatus for penetrating tissue
US8702624B2 (en) 2006-09-29 2014-04-22 Sanofi-Aventis Deutschland Gmbh Analyte measurement device with a single shot actuator
US7901362B2 (en) 2002-04-19 2011-03-08 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7674232B2 (en) 2002-04-19 2010-03-09 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7909778B2 (en) 2002-04-19 2011-03-22 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7226461B2 (en) 2002-04-19 2007-06-05 Pelikan Technologies, Inc. Method and apparatus for a multi-use body fluid sampling device with sterility barrier release
US7291117B2 (en) 2002-04-19 2007-11-06 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7717863B2 (en) 2002-04-19 2010-05-18 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7976476B2 (en) 2002-04-19 2011-07-12 Pelikan Technologies, Inc. Device and method for variable speed lancet
US7491178B2 (en) 2002-04-19 2009-02-17 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7648468B2 (en) 2002-04-19 2010-01-19 Pelikon Technologies, Inc. Method and apparatus for penetrating tissue
US8267870B2 (en) 2002-04-19 2012-09-18 Sanofi-Aventis Deutschland Gmbh Method and apparatus for body fluid sampling with hybrid actuation
US9248267B2 (en) 2002-04-19 2016-02-02 Sanofi-Aventis Deustchland Gmbh Tissue penetration device
US7371247B2 (en) 2002-04-19 2008-05-13 Pelikan Technologies, Inc Method and apparatus for penetrating tissue
US7229458B2 (en) 2002-04-19 2007-06-12 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US9314194B2 (en) 2002-04-19 2016-04-19 Sanofi-Aventis Deutschland Gmbh Tissue penetration device
US8221334B2 (en) 2002-04-19 2012-07-17 Sanofi-Aventis Deutschland Gmbh Method and apparatus for penetrating tissue
US8784335B2 (en) 2002-04-19 2014-07-22 Sanofi-Aventis Deutschland Gmbh Body fluid sampling device with a capacitive sensor
US7331931B2 (en) 2002-04-19 2008-02-19 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
US7232451B2 (en) 2002-04-19 2007-06-19 Pelikan Technologies, Inc. Method and apparatus for penetrating tissue
JP4764626B2 (en) 2002-05-06 2011-09-07 ベクトン・ディキンソン・アンド・カンパニー Method and device for controlling the pharmacokinetics of a drug
US8574895B2 (en) 2002-12-30 2013-11-05 Sanofi-Aventis Deutschland Gmbh Method and apparatus using optical techniques to measure analyte levels
EP1578262A4 (en) 2002-12-31 2007-12-05 Therasense Inc Continuous glucose monitoring system and methods of use
US7016715B2 (en) * 2003-01-13 2006-03-21 Nellcorpuritan Bennett Incorporated Selection of preset filter parameters based on signal quality
US20040253736A1 (en) * 2003-06-06 2004-12-16 Phil Stout Analytical device with prediction module and related methods
WO2004107964A2 (en) 2003-06-06 2004-12-16 Pelikan Technologies, Inc. Blood harvesting device with electronic control
US8066639B2 (en) 2003-06-10 2011-11-29 Abbott Diabetes Care Inc. Glucose measuring device for use in personal area network
WO2006001797A1 (en) 2004-06-14 2006-01-05 Pelikan Technologies, Inc. Low pain penetrating
US8423113B2 (en) 2003-07-25 2013-04-16 Dexcom, Inc. Systems and methods for processing sensor data
US20100168543A1 (en) 2003-08-01 2010-07-01 Dexcom, Inc. System and methods for processing analyte sensor data
US7494465B2 (en) 2004-07-13 2009-02-24 Dexcom, Inc. Transcutaneous analyte sensor
US8845536B2 (en) 2003-08-01 2014-09-30 Dexcom, Inc. Transcutaneous analyte sensor
US7986986B2 (en) 2003-08-01 2011-07-26 Dexcom, Inc. System and methods for processing analyte sensor data
US8369919B2 (en) 2003-08-01 2013-02-05 Dexcom, Inc. Systems and methods for processing sensor data
US8160669B2 (en) 2003-08-01 2012-04-17 Dexcom, Inc. Transcutaneous analyte sensor
EP1671096A4 (en) 2003-09-29 2009-09-16 Pelikan Technologies Inc Method and apparatus for an improved sample capture device
US9123077B2 (en) 2003-10-07 2015-09-01 Hospira, Inc. Medication management system
US8065161B2 (en) 2003-11-13 2011-11-22 Hospira, Inc. System for maintaining drug information and communicating with medication delivery devices
EP1680014A4 (en) 2003-10-14 2009-01-21 Pelikan Technologies Inc Method and apparatus for a variable user interface
US20050119866A1 (en) * 2003-11-14 2005-06-02 Zaleski John R. Medical parameter processing system
US8287453B2 (en) 2003-12-05 2012-10-16 Dexcom, Inc. Analyte sensor
EP2239567B1 (en) 2003-12-05 2015-09-02 DexCom, Inc. Calibration techniques for a continuous analyte sensor
US7822454B1 (en) 2005-01-03 2010-10-26 Pelikan Technologies, Inc. Fluid sampling device with improved analyte detecting member configuration
EP1706026B1 (en) 2003-12-31 2017-03-01 Sanofi-Aventis Deutschland GmbH Method and apparatus for improving fluidic flow and sample capture
US20050182500A1 (en) * 2004-02-17 2005-08-18 Continuous Control Solutions, Inc. Time delay definition
WO2005089103A2 (en) 2004-02-17 2005-09-29 Therasense, Inc. Method and system for providing data communication in continuous glucose monitoring and management system
WO2005081171A2 (en) * 2004-02-19 2005-09-01 Edward Henry Mathews Insulin bolus calculator for mobile communication device
US8792955B2 (en) 2004-05-03 2014-07-29 Dexcom, Inc. Transcutaneous analyte sensor
US8828203B2 (en) 2004-05-20 2014-09-09 Sanofi-Aventis Deutschland Gmbh Printable hydrogels for biosensors
EP1765194A4 (en) 2004-06-03 2010-09-29 Pelikan Technologies Inc Method and apparatus for a fluid sampling device
US7755503B2 (en) * 2004-06-16 2010-07-13 Dell Products L.P. System and method for error messaging of an incorrect information handling system graphics cable connection
US8452368B2 (en) 2004-07-13 2013-05-28 Dexcom, Inc. Transcutaneous analyte sensor
EP3524142B1 (en) * 2004-07-13 2021-04-28 Dexcom, Inc. Transcutaneous analyte sensor
US7857760B2 (en) 2004-07-13 2010-12-28 Dexcom, Inc. Analyte sensor
US20060270922A1 (en) 2004-07-13 2006-11-30 Brauker James H Analyte sensor
US7783333B2 (en) 2004-07-13 2010-08-24 Dexcom, Inc. Transcutaneous medical device with variable stiffness
US8565848B2 (en) 2004-07-13 2013-10-22 Dexcom, Inc. Transcutaneous analyte sensor
US20080027679A1 (en) * 2004-07-21 2008-01-31 Dror Shklarski Wearable Device, System and Method for Measuring Physiological and/or Environmental Parameters
US8396670B2 (en) * 2004-08-16 2013-03-12 Venture Milling, Inc. Process, system and method for improving the determination of digestive effects upon an ingestable substance
FR2874807B1 (en) * 2004-09-03 2006-11-10 Memscap Sa DEVICE FOR ANALYZING THE PHYSICO-CHEMICAL PROPERTIES OF THE SKIN
US20070191716A1 (en) * 2004-09-29 2007-08-16 Daniel Goldberger Blood monitoring system
US7608042B2 (en) * 2004-09-29 2009-10-27 Intellidx, Inc. Blood monitoring system
US20060229531A1 (en) * 2005-02-01 2006-10-12 Daniel Goldberger Blood monitoring system
US20060069436A1 (en) * 2004-09-30 2006-03-30 Depuy Spine, Inc. Trial disk implant
US9636450B2 (en) 2007-02-19 2017-05-02 Udo Hoss Pump system modular components for delivering medication and analyte sensing at seperate insertion sites
US8652831B2 (en) 2004-12-30 2014-02-18 Sanofi-Aventis Deutschland Gmbh Method and apparatus for analyte measurement test time
EP1850734A4 (en) 2005-01-13 2009-08-26 Welch Allyn Inc Vital signs monitor
US7545272B2 (en) 2005-02-08 2009-06-09 Therasense, Inc. RF tag on test strips, test strip vials and boxes
US7509156B2 (en) * 2005-05-18 2009-03-24 Clarian Health Partners, Inc. System for managing glucose levels in patients with diabetes or hyperglycemia
US8880138B2 (en) 2005-09-30 2014-11-04 Abbott Diabetes Care Inc. Device for channeling fluid and methods of use
WO2007047279A1 (en) * 2005-10-18 2007-04-26 Richards Cynthia C Dispenser having a first pump for insulin and a second pump for glucose or glucagon
US7766829B2 (en) 2005-11-04 2010-08-03 Abbott Diabetes Care Inc. Method and system for providing basal profile modification in analyte monitoring and management systems
US20070123801A1 (en) * 2005-11-28 2007-05-31 Daniel Goldberger Wearable, programmable automated blood testing system
US20080200838A1 (en) * 2005-11-28 2008-08-21 Daniel Goldberger Wearable, programmable automated blood testing system
CA2631870C (en) 2005-11-29 2021-07-06 Venture Gain, L.L.C. Residual-based monitoring of human health
US7826879B2 (en) 2006-02-28 2010-11-02 Abbott Diabetes Care Inc. Analyte sensors and methods of use
EP1991110B1 (en) 2006-03-09 2018-11-07 DexCom, Inc. Systems and methods for processing analyte sensor data
US7620438B2 (en) 2006-03-31 2009-11-17 Abbott Diabetes Care Inc. Method and system for powering an electronic device
US8473022B2 (en) 2008-01-31 2013-06-25 Abbott Diabetes Care Inc. Analyte sensor with time lag compensation
US7801582B2 (en) 2006-03-31 2010-09-21 Abbott Diabetes Care Inc. Analyte monitoring and management system and methods therefor
US8226891B2 (en) 2006-03-31 2012-07-24 Abbott Diabetes Care Inc. Analyte monitoring devices and methods therefor
US7618369B2 (en) 2006-10-02 2009-11-17 Abbott Diabetes Care Inc. Method and system for dynamically updating calibration parameters for an analyte sensor
US7653425B2 (en) 2006-08-09 2010-01-26 Abbott Diabetes Care Inc. Method and system for providing calibration of an analyte sensor in an analyte monitoring system
US8374668B1 (en) 2007-10-23 2013-02-12 Abbott Diabetes Care Inc. Analyte sensor with lag compensation
US9392969B2 (en) 2008-08-31 2016-07-19 Abbott Diabetes Care Inc. Closed loop control and signal attenuation detection
US8346335B2 (en) 2008-03-28 2013-01-01 Abbott Diabetes Care Inc. Analyte sensor calibration management
US8140312B2 (en) 2007-05-14 2012-03-20 Abbott Diabetes Care Inc. Method and system for determining analyte levels
US20070243567A1 (en) * 2006-04-13 2007-10-18 Syhhong Chang Beta cell mimicking control algorithm for artificial pancreas
US8092385B2 (en) * 2006-05-23 2012-01-10 Intellidx, Inc. Fluid access interface
US20080071158A1 (en) 2006-06-07 2008-03-20 Abbott Diabetes Care, Inc. Analyte monitoring system and method
ES2382449T3 (en) * 2006-06-13 2012-06-08 Carefusion 303, Inc. System and method for optimizing control of PCA and PCEA systems
EP1875933A1 (en) * 2006-07-06 2008-01-09 Debiotech S.A. Medical device for delivery of a solution
US8388534B2 (en) * 2006-10-11 2013-03-05 Samsung Electronics Co., Ltd. Apparatus for providing skin care information by measuring skin moisture content and method and medium for the same
US20080091466A1 (en) 2006-10-16 2008-04-17 Hospira, Inc. System and method for comparing and utilizing activity information and configuration information from multiple device management systems
US8721584B2 (en) * 2006-11-02 2014-05-13 The University Of North Carolina At Chapel Hill Methods and systems for determining an intravenous insulin infusion rate to correct hyperglycemia of a patient, to maintain euglycemia of a patient, and to prevent hypoglycemia of a patient
US20080306353A1 (en) * 2006-11-03 2008-12-11 Douglas Joel S Calculation device for metabolic control of critically ill and/or diabetic patients
US20080214919A1 (en) * 2006-12-26 2008-09-04 Lifescan, Inc. System and method for implementation of glycemic control protocols
US10154804B2 (en) 2007-01-31 2018-12-18 Medtronic Minimed, Inc. Model predictive method and system for controlling and supervising insulin infusion
US20080199894A1 (en) 2007-02-15 2008-08-21 Abbott Diabetes Care, Inc. Device and method for automatic data acquisition and/or detection
US8930203B2 (en) 2007-02-18 2015-01-06 Abbott Diabetes Care Inc. Multi-function analyte test device and methods therefor
US8732188B2 (en) 2007-02-18 2014-05-20 Abbott Diabetes Care Inc. Method and system for providing contextual based medication dosage determination
US8123686B2 (en) 2007-03-01 2012-02-28 Abbott Diabetes Care Inc. Method and apparatus for providing rolling data in communication systems
EP4108162A1 (en) 2007-04-14 2022-12-28 Abbott Diabetes Care, Inc. Method and apparatus for providing data processing and control in medical communication system
US9615780B2 (en) 2007-04-14 2017-04-11 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in medical communication system
WO2008130895A2 (en) 2007-04-14 2008-10-30 Abbott Diabetes Care, Inc. Method and apparatus for providing dynamic multi-stage signal amplification in a medical device
EP2146625B1 (en) 2007-04-14 2019-08-14 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in medical communication system
WO2008130896A1 (en) 2007-04-14 2008-10-30 Abbott Diabetes Care, Inc. Method and apparatus for providing data processing and control in medical communication system
US8140142B2 (en) * 2007-04-14 2012-03-20 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in medical communication system
US8461985B2 (en) 2007-05-08 2013-06-11 Abbott Diabetes Care Inc. Analyte monitoring system and methods
US7928850B2 (en) 2007-05-08 2011-04-19 Abbott Diabetes Care Inc. Analyte monitoring system and methods
US8456301B2 (en) 2007-05-08 2013-06-04 Abbott Diabetes Care Inc. Analyte monitoring system and methods
US8665091B2 (en) 2007-05-08 2014-03-04 Abbott Diabetes Care Inc. Method and device for determining elapsed sensor life
US8103471B2 (en) 2007-05-14 2012-01-24 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8600681B2 (en) 2007-05-14 2013-12-03 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8560038B2 (en) 2007-05-14 2013-10-15 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8239166B2 (en) 2007-05-14 2012-08-07 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US9125548B2 (en) 2007-05-14 2015-09-08 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8444560B2 (en) 2007-05-14 2013-05-21 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US10002233B2 (en) 2007-05-14 2018-06-19 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8260558B2 (en) 2007-05-14 2012-09-04 Abbott Diabetes Care Inc. Method and apparatus for providing data processing and control in a medical communication system
US8417311B2 (en) 2008-09-12 2013-04-09 Optiscan Biomedical Corporation Fluid component analysis system and method for glucose monitoring and control
DK2006786T3 (en) * 2007-06-18 2018-08-06 Hoffmann La Roche Method and glucose monitoring system to monitor individual metabolic response and to generate nutrient feedback
US8834366B2 (en) 2007-07-31 2014-09-16 Abbott Diabetes Care Inc. Method and apparatus for providing analyte sensor calibration
US20090143725A1 (en) * 2007-08-31 2009-06-04 Abbott Diabetes Care, Inc. Method of Optimizing Efficacy of Therapeutic Agent
MX2010003855A (en) 2007-10-10 2010-04-27 Optiscan Biomedical Corp Fluid component analysis system and method for glucose monitoring and control.
US8377031B2 (en) 2007-10-23 2013-02-19 Abbott Diabetes Care Inc. Closed loop control system with safety parameters and methods
US8409093B2 (en) 2007-10-23 2013-04-02 Abbott Diabetes Care Inc. Assessing measures of glycemic variability
US8517990B2 (en) 2007-12-18 2013-08-27 Hospira, Inc. User interface improvements for medical devices
US20090164239A1 (en) 2007-12-19 2009-06-25 Abbott Diabetes Care, Inc. Dynamic Display Of Glucose Information
US8229535B2 (en) 2008-02-21 2012-07-24 Dexcom, Inc. Systems and methods for blood glucose monitoring and alert delivery
US8396528B2 (en) 2008-03-25 2013-03-12 Dexcom, Inc. Analyte sensor
WO2009121392A1 (en) * 2008-03-31 2009-10-08 Scibase Ab Method and device for non-invasive determination of the concentration of a substance in a body fluid
US8600682B2 (en) 2008-04-04 2013-12-03 Hygieia, Inc. Apparatus for optimizing a patient's insulin dosage regimen
US9220456B2 (en) 2008-04-04 2015-12-29 Hygieia, Inc. Systems, methods and devices for achieving glycemic balance
US10624577B2 (en) 2008-04-04 2020-04-21 Hygieia, Inc. Systems, devices, and methods for alleviating glucotoxicity and restoring pancreatic beta-cell function in advanced diabetes mellitus
EP2265324B1 (en) 2008-04-11 2015-01-28 Sanofi-Aventis Deutschland GmbH Integrated analyte measurement system
US20090281519A1 (en) * 2008-05-12 2009-11-12 Rao R Harsha Automated system and method for diabetes control
US8876755B2 (en) 2008-07-14 2014-11-04 Abbott Diabetes Care Inc. Closed loop control system interface and methods
JPWO2010007720A1 (en) * 2008-07-16 2012-01-05 パナソニック株式会社 Battery pack
WO2010019919A1 (en) * 2008-08-14 2010-02-18 University Of Toledo Multifunctional neural network system and uses thereof for glycemic forecasting
US7959598B2 (en) 2008-08-20 2011-06-14 Asante Solutions, Inc. Infusion pump systems and methods
US8622988B2 (en) 2008-08-31 2014-01-07 Abbott Diabetes Care Inc. Variable rate closed loop control and methods
US9943644B2 (en) 2008-08-31 2018-04-17 Abbott Diabetes Care Inc. Closed loop control with reference measurement and methods thereof
US8734422B2 (en) 2008-08-31 2014-05-27 Abbott Diabetes Care Inc. Closed loop control with improved alarm functions
US20100057040A1 (en) 2008-08-31 2010-03-04 Abbott Diabetes Care, Inc. Robust Closed Loop Control And Methods
US8986208B2 (en) 2008-09-30 2015-03-24 Abbott Diabetes Care Inc. Analyte sensor sensitivity attenuation mitigation
US20100174228A1 (en) * 2008-10-24 2010-07-08 Bruce Buckingham Hypoglycemia prediction and control
WO2010056718A2 (en) 2008-11-11 2010-05-20 Hygieia, Inc. Apparatus and system for diabetes management
CN102405011A (en) * 2008-11-26 2012-04-04 弗吉尼亚大学专利基金会 Method, system, and computer program product for tracking of blood glucose variability in diabetes
US8103456B2 (en) 2009-01-29 2012-01-24 Abbott Diabetes Care Inc. Method and device for early signal attenuation detection using blood glucose measurements
US9375169B2 (en) 2009-01-30 2016-06-28 Sanofi-Aventis Deutschland Gmbh Cam drive for managing disposable penetrating member actions with a single motor and motor and control system
FR2942387A1 (en) * 2009-02-23 2010-08-27 Ecole Superieure Electricite SYSTEM FOR REGULATING BLOOD GLUCOSE WITH LIVING DIABETIC SUBJECT AND CORRESPONDING COMPUTER PROGRAM PRODUCT.
KR20120047841A (en) * 2009-02-26 2012-05-14 몰 리서치 어플리케이션스 엘티디 Method and system for automatic monitoring of diabetes related treatment
US8753290B2 (en) * 2009-03-27 2014-06-17 Intellectual Inspiration, Llc Fluid transfer system and method
US8497777B2 (en) 2009-04-15 2013-07-30 Abbott Diabetes Care Inc. Analyte monitoring system having an alert
US8271106B2 (en) 2009-04-17 2012-09-18 Hospira, Inc. System and method for configuring a rule set for medical event management and responses
DE102009018649A1 (en) 2009-04-23 2010-10-28 Fresenius Medical Care Deutschland Gmbh Dialysis machine, method and apparatus for a dialysis machine
WO2010127050A1 (en) 2009-04-28 2010-11-04 Abbott Diabetes Care Inc. Error detection in critical repeating data in a wireless sensor system
US8483967B2 (en) 2009-04-29 2013-07-09 Abbott Diabetes Care Inc. Method and system for providing real time analyte sensor calibration with retrospective backfill
US20100279418A1 (en) * 2009-05-04 2010-11-04 Loren Robert Larson Glucose meter adaptable for use with handheld devices, and associated communication network
WO2010138856A1 (en) 2009-05-29 2010-12-02 Abbott Diabetes Care Inc. Medical device antenna systems having external antenna configurations
GB2471066A (en) * 2009-06-10 2010-12-22 Dna Electronics Ltd A glucagon pump controller
ES2888427T3 (en) 2009-07-23 2022-01-04 Abbott Diabetes Care Inc Real-time management of data related to the physiological control of glucose levels
EP3689237B1 (en) 2009-07-23 2021-05-19 Abbott Diabetes Care, Inc. Method of manufacturing and system for continuous analyte measurement
WO2011014851A1 (en) 2009-07-31 2011-02-03 Abbott Diabetes Care Inc. Method and apparatus for providing analyte monitoring system calibration accuracy
EP2473098A4 (en) 2009-08-31 2014-04-09 Abbott Diabetes Care Inc Analyte signal processing device and methods
ES2952361T3 (en) 2009-08-31 2023-10-31 Abbott Diabetes Care Inc Displays for a medical device
US8993331B2 (en) 2009-08-31 2015-03-31 Abbott Diabetes Care Inc. Analyte monitoring system and methods for managing power and noise
WO2011041469A1 (en) 2009-09-29 2011-04-07 Abbott Diabetes Care Inc. Method and apparatus for providing notification function in analyte monitoring systems
KR20120102048A (en) 2009-09-30 2012-09-17 모르 리서치 애플리케이션즈 리미티드 Monitoring device for management of insulin delivery
US8538740B2 (en) * 2009-10-14 2013-09-17 International Business Machines Corporation Real-time performance modeling of software systems with multi-class workload
EP2319394A1 (en) * 2009-11-05 2011-05-11 Rus Medical Technology S.A. Non-invasive device and method for monitoring analytes in biological samples
US8771251B2 (en) 2009-12-17 2014-07-08 Hospira, Inc. Systems and methods for managing and delivering patient therapy through electronic drug delivery systems
WO2011091336A1 (en) * 2010-01-22 2011-07-28 Abbott Diabetes Care Inc. Method and apparatus for providing notification in analyte monitoring systems
EP4183326A1 (en) * 2010-03-26 2023-05-24 University Of Virginia Patent Foundation Method, system, and computer program product for improving the accuracy of glucose sensors using insulin delivery observation in diabetes
US9089292B2 (en) * 2010-03-26 2015-07-28 Medtronic Minimed, Inc. Calibration of glucose monitoring sensor and/or insulin delivery system
FR2958096B1 (en) * 2010-03-29 2012-04-13 Commissariat Energie Atomique METHOD AND SYSTEM FOR CALIBRATION, RECORDING MEDIUM FOR THIS METHOD
US8965476B2 (en) 2010-04-16 2015-02-24 Sanofi-Aventis Deutschland Gmbh Tissue penetration device
US8636711B2 (en) * 2010-06-14 2014-01-28 Legacy Emanuel Hospital & Health Center Stabilized glucagon solutions and uses therefor
EP2624745A4 (en) 2010-10-07 2018-05-23 Abbott Diabetes Care, Inc. Analyte monitoring devices and methods
AU2012299169B2 (en) 2011-08-19 2017-08-24 Icu Medical, Inc. Systems and methods for a graphical interface including a graphical representation of medical data
US9594875B2 (en) 2011-10-21 2017-03-14 Hospira, Inc. Medical device update system
US8710993B2 (en) 2011-11-23 2014-04-29 Abbott Diabetes Care Inc. Mitigating single point failure of devices in an analyte monitoring system and methods thereof
US9317656B2 (en) 2011-11-23 2016-04-19 Abbott Diabetes Care Inc. Compatibility mechanisms for devices in a continuous analyte monitoring system and methods thereof
WO2013078426A2 (en) 2011-11-25 2013-05-30 Abbott Diabetes Care Inc. Analyte monitoring system and methods of use
US20130138094A1 (en) * 2011-11-30 2013-05-30 Amo Development, Llc System and method for ophthalmic surface measurements based on sequential estimates
US9510747B2 (en) * 2011-11-30 2016-12-06 Amo Development, Llc System and method for ophthalmic surface measurements based on objective quality estimation
US10022498B2 (en) 2011-12-16 2018-07-17 Icu Medical, Inc. System for monitoring and delivering medication to a patient and method of using the same to minimize the risks associated with automated therapy
US20130231543A1 (en) * 2012-03-05 2013-09-05 Andrea Facchinetti Method to improve safety monitoring in type-1 diabetic patients by detecting in real-time failures of the glucose
US9995611B2 (en) 2012-03-30 2018-06-12 Icu Medical, Inc. Air detection system and method for detecting air in a pump of an infusion system
ES2743160T3 (en) 2012-07-31 2020-02-18 Icu Medical Inc Patient care system for critical medications
US10132793B2 (en) 2012-08-30 2018-11-20 Abbott Diabetes Care Inc. Dropout detection in continuous analyte monitoring data during data excursions
US9968306B2 (en) 2012-09-17 2018-05-15 Abbott Diabetes Care Inc. Methods and apparatuses for providing adverse condition notification with enhanced wireless communication range in analyte monitoring systems
WO2014055718A1 (en) 2012-10-04 2014-04-10 Aptima, Inc. Clinical support systems and methods
US9119529B2 (en) 2012-10-30 2015-09-01 Dexcom, Inc. Systems and methods for dynamically and intelligently monitoring a host's glycemic condition after an alert is triggered
WO2014109898A1 (en) * 2013-01-14 2014-07-17 The Regents Of University Of California Model-based personalization scheme of an artificial pancreas for type i diabetes applications
US9227014B2 (en) * 2013-02-07 2016-01-05 The Board Of Trustees Of The Laland Stanford Junior University Kalman filter based on-off switch for insulin pump
EP2964079B1 (en) 2013-03-06 2022-02-16 ICU Medical, Inc. Medical device communication method
WO2014190264A1 (en) 2013-05-24 2014-11-27 Hospira, Inc. Multi-sensor infusion system for detecting air or an occlusion in the infusion system
CA2913918C (en) 2013-05-29 2022-02-15 Hospira, Inc. Infusion system and method of use which prevents over-saturation of an analog-to-digital converter
AU2014274146B2 (en) 2013-05-29 2019-01-24 Icu Medical, Inc. Infusion system which utilizes one or more sensors and additional information to make an air determination regarding the infusion system
US9561324B2 (en) 2013-07-19 2017-02-07 Bigfoot Biomedical, Inc. Infusion pump system and method
EP3022668A1 (en) 2013-07-19 2016-05-25 Dexcom, Inc. Time averaged basal rate optimizer
US20150066531A1 (en) 2013-08-30 2015-03-05 James D. Jacobson System and method of monitoring and managing a remote infusion regimen
US9662436B2 (en) 2013-09-20 2017-05-30 Icu Medical, Inc. Fail-safe drug infusion therapy system
FI3060276T3 (en) 2013-10-24 2023-06-07 Infusion system for preventing mischanneling of multiple medicaments
US10311972B2 (en) 2013-11-11 2019-06-04 Icu Medical, Inc. Medical device system performance index
JP2016537175A (en) 2013-11-19 2016-12-01 ホスピーラ インコーポレイテッド Infusion pump automation system and method
US10569015B2 (en) 2013-12-02 2020-02-25 Bigfoot Biomedical, Inc. Infusion pump system and method
GB2523989B (en) 2014-01-30 2020-07-29 Insulet Netherlands B V Therapeutic product delivery system and method of pairing
US9898585B2 (en) 2014-01-31 2018-02-20 Aseko, Inc. Method and system for insulin management
US9486580B2 (en) 2014-01-31 2016-11-08 Aseko, Inc. Insulin management
EP3110474B1 (en) 2014-02-28 2019-12-18 ICU Medical, Inc. Infusion system and method which utilizes dual wavelength optical air-in-line detection
US10722650B2 (en) 2014-03-28 2020-07-28 Roche Diabetes Care, Inc. System and method for adjusting therapy based on risk associated with a glucose state
CA2945647C (en) 2014-04-30 2023-08-08 Hospira, Inc. Patient care system with conditional alarm forwarding
AU2015266706B2 (en) 2014-05-29 2020-01-30 Icu Medical, Inc. Infusion system and pump with configurable closed loop delivery rate catch-up
US9724470B2 (en) 2014-06-16 2017-08-08 Icu Medical, Inc. System for monitoring and delivering medication to a patient and method of using the same to minimize the risks associated with automated therapy
EP2979623A1 (en) * 2014-07-28 2016-02-03 Roche Diagnostics GmbH Medical sensor assembly
AU2015301454A1 (en) * 2014-08-14 2017-03-09 University Of Virginia Patent Foundation Improved accuracy continuous glucose monitoring method, system, and device
US9539383B2 (en) 2014-09-15 2017-01-10 Hospira, Inc. System and method that matches delayed infusion auto-programs with manually entered infusion programs and analyzes differences therein
US11081226B2 (en) 2014-10-27 2021-08-03 Aseko, Inc. Method and controller for administering recommended insulin dosages to a patient
US9892234B2 (en) 2014-10-27 2018-02-13 Aseko, Inc. Subcutaneous outpatient management
US20160147976A1 (en) * 2014-11-26 2016-05-26 RxAdvance Corporation Medication Identification, Tracking And Adherence Management
US10016551B2 (en) 2014-12-18 2018-07-10 Gambro Lundia Ab Method of displaying a predicted state, medical apparatus and computer program
US11344668B2 (en) 2014-12-19 2022-05-31 Icu Medical, Inc. Infusion system with concurrent TPN/insulin infusion
EP3258991B1 (en) 2015-02-18 2020-10-21 Insulet Corporation Fluid delivery and infusion devices, and methods of use thereof
US10850024B2 (en) 2015-03-02 2020-12-01 Icu Medical, Inc. Infusion system, device, and method having advanced infusion features
US9878097B2 (en) 2015-04-29 2018-01-30 Bigfoot Biomedical, Inc. Operating an infusion pump system
EP3304370B1 (en) 2015-05-26 2020-12-30 ICU Medical, Inc. Infusion pump system and method with multiple drug library editor source capability
EP3302664B1 (en) * 2015-06-02 2020-08-05 Koninklijke Philips N.V. Non-invasive method for monitoring patient respiratory status via successive parameter estimation
CA2991058A1 (en) 2015-07-08 2017-01-12 Trustees Of Boston University Infusion system and components thereof
US11553883B2 (en) 2015-07-10 2023-01-17 Abbott Diabetes Care Inc. System, device and method of dynamic glucose profile response to physiological parameters
US11464456B2 (en) 2015-08-07 2022-10-11 Aptima, Inc. Systems and methods to support medical therapy decisions
EP3337402A4 (en) 2015-08-20 2019-04-24 Aseko, Inc. Diabetes management therapy advisor
CA3009409A1 (en) 2016-01-05 2017-07-13 Bigfoot Biomedical, Inc. Operating multi-modal medicine delivery systems
US10449294B1 (en) 2016-01-05 2019-10-22 Bigfoot Biomedical, Inc. Operating an infusion pump system
US10275573B2 (en) 2016-01-13 2019-04-30 Bigfoot Biomedical, Inc. User interface for diabetes management system
US10780223B2 (en) 2016-01-14 2020-09-22 Bigfoot Biomedical, Inc. Adjusting insulin delivery rates
AU2017216393A1 (en) * 2016-02-01 2018-06-28 Dexcom, Inc. System and method for decision support using lifestyle factors
CA3211813A1 (en) 2016-05-02 2017-11-09 Dexcom, Inc. System and method for providing alerts optimized for a user
AU2017264784B2 (en) 2016-05-13 2022-04-21 Icu Medical, Inc. Infusion pump system and method with common line auto flush
EP3468635A4 (en) 2016-06-10 2019-11-20 ICU Medical, Inc. Acoustic flow sensor for continuous medication flow measurements and feedback control of infusion
AU2017295722B2 (en) 2016-07-14 2022-08-11 Icu Medical, Inc. Multi-communication path selection and security system for a medical device
EP3515535A1 (en) 2016-09-23 2019-07-31 Insulet Corporation Fluid delivery device with sensor
US10783801B1 (en) 2016-12-21 2020-09-22 Aptima, Inc. Simulation based training system for measurement of team cognitive load to automatically customize simulation content
RU2757312C2 (en) 2017-01-06 2021-10-13 Трастис Оф Бостон Юниверсити Infusion system and its components
WO2018175489A1 (en) 2017-03-21 2018-09-27 Abbott Diabetes Care Inc. Methods, devices and system for providing diabetic condition diagnosis and therapy
CN116531603A (en) 2017-05-05 2023-08-04 伊莱利利公司 Closed loop control of physiological glucose
WO2019090016A1 (en) * 2017-11-02 2019-05-09 Boston Scientific Scimed, Inc. Systems and methods for graded glucose control
WO2019090009A1 (en) 2017-11-02 2019-05-09 Boston Scientific Scimed, Inc. System and method for providing glucose control therapy
US11901060B2 (en) 2017-12-21 2024-02-13 Ypsomed Ag Closed loop control of physiological glucose
US10089055B1 (en) 2017-12-27 2018-10-02 Icu Medical, Inc. Synchronized display of screen content on networked devices
JP7443235B2 (en) 2018-02-09 2024-03-05 デックスコム・インコーポレーテッド Decision support systems and methods
WO2019213493A1 (en) 2018-05-04 2019-11-07 Insulet Corporation Safety constraints for a control algorithm-based drug delivery system
US11443649B2 (en) * 2018-06-29 2022-09-13 Cadwell Laboratories, Inc. Neurophysiological monitoring training simulator
US10861592B2 (en) 2018-07-17 2020-12-08 Icu Medical, Inc. Reducing infusion pump network congestion by staggering updates
US11483403B2 (en) 2018-07-17 2022-10-25 Icu Medical, Inc. Maintaining clinical messaging during network instability
AU2019306490A1 (en) 2018-07-17 2021-02-04 Icu Medical, Inc. Updating infusion pump drug libraries and operational software in a networked environment
WO2020018389A1 (en) 2018-07-17 2020-01-23 Icu Medical, Inc. Systems and methods for facilitating clinical messaging in a network environment
US10692595B2 (en) 2018-07-26 2020-06-23 Icu Medical, Inc. Drug library dynamic version management
CA3107315C (en) 2018-07-26 2023-01-03 Icu Medical, Inc. Drug library management system
CN112789070A (en) 2018-09-28 2021-05-11 英赛罗公司 Mode of activity of the artificial pancreas System
EP3864668A1 (en) 2018-10-11 2021-08-18 Insulet Corporation Event detection for drug delivery system
EP3997712A4 (en) * 2019-07-09 2023-07-12 Eli Lilly And Company Method and system for determining glucose change in a subject
DE112020003400T5 (en) 2019-07-16 2022-03-31 Beta Bionics, Inc. OUTPATIENT INSTITUTION AND ITS COMPONENTS
US20210045640A1 (en) * 2019-08-16 2021-02-18 Poltorak Technologies, LLC Device and method for medical diagnostics
US11801344B2 (en) 2019-09-13 2023-10-31 Insulet Corporation Blood glucose rate of change modulation of meal and correction insulin bolus quantity
US11278671B2 (en) 2019-12-04 2022-03-22 Icu Medical, Inc. Infusion pump with safety sequence keypad
US11833329B2 (en) 2019-12-20 2023-12-05 Insulet Corporation Techniques for improved automatic drug delivery performance using delivery tendencies from past delivery history and use patterns
US11551802B2 (en) 2020-02-11 2023-01-10 Insulet Corporation Early meal detection and calorie intake detection
US11547800B2 (en) 2020-02-12 2023-01-10 Insulet Corporation User parameter dependent cost function for personalized reduction of hypoglycemia and/or hyperglycemia in a closed loop artificial pancreas system
US11324889B2 (en) 2020-02-14 2022-05-10 Insulet Corporation Compensation for missing readings from a glucose monitor in an automated insulin delivery system
US11278661B2 (en) 2020-03-10 2022-03-22 Beta Bionics, Inc. Infusion system and components thereof
US11607493B2 (en) 2020-04-06 2023-03-21 Insulet Corporation Initial total daily insulin setting for user onboarding
WO2022020184A1 (en) 2020-07-21 2022-01-27 Icu Medical, Inc. Fluid transfer devices and methods of use
US11684716B2 (en) 2020-07-31 2023-06-27 Insulet Corporation Techniques to reduce risk of occlusions in drug delivery systems
US20220189603A1 (en) 2020-12-07 2022-06-16 Beta Bionics, Inc. Medicament pumps and control systems for managing glucose control therapy data of a subject
US11135360B1 (en) 2020-12-07 2021-10-05 Icu Medical, Inc. Concurrent infusion with common line auto flush
US11904140B2 (en) 2021-03-10 2024-02-20 Insulet Corporation Adaptable asymmetric medicament cost component in a control system for medicament delivery
US11738144B2 (en) 2021-09-27 2023-08-29 Insulet Corporation Techniques enabling adaptation of parameters in aid systems by user input
US11439754B1 (en) 2021-12-01 2022-09-13 Insulet Corporation Optimizing embedded formulations for drug delivery

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1997028737A1 (en) * 1996-02-12 1997-08-14 Nokia Mobile Phones Limited A method for monitoring the health of a patient by measuring and predicting the glucose level of the patient's blood sample
WO1999056613A1 (en) * 1998-04-30 1999-11-11 Therasense, Inc. Analyte monitoring device and methods of use
WO2000015102A1 (en) * 1998-09-10 2000-03-23 Spectrx, Inc. Attribute compensation for analyte detection and/or continuous monitoring

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4655225A (en) 1985-04-18 1987-04-07 Kurabo Industries Ltd. Spectrophotometric method and apparatus for the non-invasive
US5086229A (en) 1989-01-19 1992-02-04 Futrex, Inc. Non-invasive measurement of blood glucose
US5298022A (en) 1989-05-29 1994-03-29 Amplifon Spa Wearable artificial pancreas
DE3927768A1 (en) * 1989-08-23 1991-02-28 Perkin Elmer Corp METHOD FOR DETERMINING CONCENTRATION BY MEANS OF ATOMEMISSION SPECTROSCOPY
US5070874A (en) 1990-01-30 1991-12-10 Biocontrol Technology, Inc. Non-invasive determination of glucose concentration in body of patients
US5490505A (en) 1991-03-07 1996-02-13 Masimo Corporation Signal processing apparatus
US5632272A (en) 1991-03-07 1997-05-27 Masimo Corporation Signal processing apparatus
EP0930045A3 (en) 1991-03-07 1999-10-27 Masimo Corporation Signal processing apparatus and method for an oximeter
EP0522674B1 (en) 1991-07-12 1998-11-11 Mark R. Robinson Oximeter for reliable clinical determination of blood oxygen saturation in a fetus
DE69218536T2 (en) * 1991-10-04 1997-07-03 Perkin Elmer Corp Method and device for comparing spectra
US5387327A (en) * 1992-10-19 1995-02-07 Duquesne University Of The Holy Ghost Implantable non-enzymatic electrochemical glucose sensor
US5956501A (en) * 1997-01-10 1999-09-21 Health Hero Network, Inc. Disease simulation system and method
US5497772A (en) 1993-11-19 1996-03-12 Alfred E. Mann Foundation For Scientific Research Glucose monitoring system
US5459317A (en) 1994-02-14 1995-10-17 Ohio University Method and apparatus for non-invasive detection of physiological chemicals, particularly glucose
US5995860A (en) 1995-07-06 1999-11-30 Thomas Jefferson University Implantable sensor and system for measurement and control of blood constituent levels
US5853364A (en) 1995-08-07 1998-12-29 Nellcor Puritan Bennett, Inc. Method and apparatus for estimating physiological parameters using model-based adaptive filtering
US6152876A (en) 1997-04-18 2000-11-28 Rio Grande Medical Technologies, Inc. Method for non-invasive blood analyte measurement with improved optical interface
US5687733A (en) * 1995-10-26 1997-11-18 Baxter International Inc. System and method for estimating cardiac output
CA2218623A1 (en) 1996-02-23 1997-08-28 Circe Biomedical, Inc. Novel artificial pancreas
US5792047A (en) * 1997-01-15 1998-08-11 Coggins; George Physiological parameter monitoring and bio-feedback apparatus
EP1006863B1 (en) 1997-03-21 2003-10-15 Nellcor Puritan Bennett Incorporated Method and apparatus for adaptively averaging data signals
US6558351B1 (en) * 1999-06-03 2003-05-06 Medtronic Minimed, Inc. Closed loop system for controlling insulin infusion
US6216690B1 (en) * 1997-10-15 2001-04-17 Datex-Ohmeda, Inc. Method and apparatus for rapid control of set inspired gas concentration in anesthesia delivery systems
EP0910023A2 (en) 1997-10-17 1999-04-21 Siemens Aktiengesellschaft Method and device for the neuronal modelling of a dynamic system with non-linear stochastic behavior
US6368272B1 (en) * 1998-04-10 2002-04-09 Proactive Metabolics Company Equipment and method for contemporaneous decision supporting metabolic control
JP4571317B2 (en) * 1999-06-01 2010-10-27 マサチューセッツ インスティテュート オブ テクノロジー Pressureless continuous blood pressure monitoring device
US6572545B2 (en) * 2000-09-22 2003-06-03 Knobbe, Lim & Buckingham Method and apparatus for real-time control of physiological parameters
US6544212B2 (en) * 2001-07-31 2003-04-08 Roche Diagnostics Corporation Diabetes management system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1997028737A1 (en) * 1996-02-12 1997-08-14 Nokia Mobile Phones Limited A method for monitoring the health of a patient by measuring and predicting the glucose level of the patient's blood sample
WO1999056613A1 (en) * 1998-04-30 1999-11-11 Therasense, Inc. Analyte monitoring device and methods of use
WO2000015102A1 (en) * 1998-09-10 2000-03-23 Spectrx, Inc. Attribute compensation for analyte detection and/or continuous monitoring

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PARKER R S ET AL: "Model predictive control for infusion pump insulin delivery", PROCEEDINGS OF THE 18TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. BRIDGING DISCIPLINES FOR BIOMEDICINE' (CAT. NO.96CH36036), PROCEEDINGS OF 18TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING I, 1997, New York, NY, USA, IEEE, USA, pages 1822 - 1823 vol.5, XP002190918, ISBN: 0-7803-3811-1 *
WOODRUFF E A ET AL: "The closed-loop regulation of blood glucose in diabetics", PROCEEDINGS OF THE FOURTEENTH ANNUAL NORTHEAST BIOENGINEERING CONFERENCE (IEEE CAT. NO.88-CH2666-6), DURHAM, NH, USA, 10-11 MARCH 1988, 1988, New York, NY, USA, IEEE, USA, pages 54 - 57, XP010078511 *
WU X ET AL: "Fast on-line data evaluation of flow-injection analysis signals based on parameter estimation by an extended Kalman filter", JOURNAL OF BIOTECHNOLOGY, ELSEVIER SCIENCE PUBLISHERS, AMSTERDAM, NL, vol. 62, no. 1, 11 June 1998 (1998-06-11), pages 11 - 28, XP004127150, ISSN: 0168-1656 *

Cited By (188)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9066697B2 (en) 1998-04-30 2015-06-30 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9326714B2 (en) 1998-04-30 2016-05-03 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US10478108B2 (en) 1998-04-30 2019-11-19 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9014773B2 (en) 1998-04-30 2015-04-21 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8880137B2 (en) 1998-04-30 2014-11-04 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8974386B2 (en) 1998-04-30 2015-03-10 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8840553B2 (en) 1998-04-30 2014-09-23 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9011331B2 (en) 1998-04-30 2015-04-21 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9042953B2 (en) 1998-04-30 2015-05-26 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9066695B2 (en) 1998-04-30 2015-06-30 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9072477B2 (en) 1998-04-30 2015-07-07 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9066694B2 (en) 1998-04-30 2015-06-30 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9011332B2 (en) 2001-01-02 2015-04-21 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9498159B2 (en) 2001-01-02 2016-11-22 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9610034B2 (en) 2001-01-02 2017-04-04 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9282925B2 (en) 2002-02-12 2016-03-15 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US8442610B2 (en) 2003-08-01 2013-05-14 Dexcom, Inc. System and methods for processing analyte sensor data
US10786185B2 (en) 2003-08-01 2020-09-29 Dexcom, Inc. System and methods for processing analyte sensor data
US7826981B2 (en) 2003-08-01 2010-11-02 Dexcom, Inc. System and methods for processing analyte sensor data
US7797028B2 (en) 2003-08-01 2010-09-14 Dexcom, Inc. System and methods for processing analyte sensor data
US7778680B2 (en) 2003-08-01 2010-08-17 Dexcom, Inc. System and methods for processing analyte sensor data
US7774145B2 (en) 2003-08-01 2010-08-10 Dexcom, Inc. Transcutaneous analyte sensor
US8060173B2 (en) 2003-08-01 2011-11-15 Dexcom, Inc. System and methods for processing analyte sensor data
US8886273B2 (en) 2003-08-01 2014-11-11 Dexcom, Inc. Analyte sensor
US8052601B2 (en) 2003-08-01 2011-11-08 Dexcom, Inc. System and methods for processing analyte sensor data
US9895089B2 (en) 2003-08-01 2018-02-20 Dexcom, Inc. System and methods for processing analyte sensor data
US8915849B2 (en) 2003-08-01 2014-12-23 Dexcom, Inc. Transcutaneous analyte sensor
US8986209B2 (en) 2003-08-01 2015-03-24 Dexcom, Inc. Transcutaneous analyte sensor
US11589823B2 (en) 2003-08-22 2023-02-28 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US11559260B2 (en) 2003-08-22 2023-01-24 Dexcom, Inc. Systems and methods for processing analyte sensor data
US9427183B2 (en) 2003-08-22 2016-08-30 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9149219B2 (en) 2003-08-22 2015-10-06 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9750460B2 (en) 2003-08-22 2017-09-05 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9420968B2 (en) 2003-08-22 2016-08-23 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US8843187B2 (en) 2003-08-22 2014-09-23 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US8010174B2 (en) 2003-08-22 2011-08-30 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9585607B2 (en) 2003-08-22 2017-03-07 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9724045B1 (en) 2003-08-22 2017-08-08 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9649069B2 (en) 2003-08-22 2017-05-16 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US9510782B2 (en) 2003-08-22 2016-12-06 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
EP1518495A1 (en) * 2003-09-23 2005-03-30 Roche Diagnostics GmbH Method and device for continuous monitoring of analyte concentration
US7389133B1 (en) 2003-09-23 2008-06-17 Roche Diagnostics Operations, Inc. Method and device for continuous monitoring of the concentration of an analyte
US11564602B2 (en) 2003-11-19 2023-01-31 Dexcom, Inc. Integrated receiver for continuous analyte sensor
US9538946B2 (en) 2003-11-19 2017-01-10 Dexcom, Inc. Integrated receiver for continuous analyte sensor
US8282550B2 (en) 2003-11-19 2012-10-09 Dexcom, Inc. Integrated receiver for continuous analyte sensor
US7927274B2 (en) 2003-11-19 2011-04-19 Dexcom, Inc. Integrated receiver for continuous analyte sensor
US11633133B2 (en) 2003-12-05 2023-04-25 Dexcom, Inc. Dual electrode system for a continuous analyte sensor
US11000215B1 (en) 2003-12-05 2021-05-11 Dexcom, Inc. Analyte sensor
US11020031B1 (en) 2003-12-05 2021-06-01 Dexcom, Inc. Analyte sensor
US9107623B2 (en) 2003-12-09 2015-08-18 Dexcom, Inc. Signal processing for continuous analyte sensor
US10898113B2 (en) 2003-12-09 2021-01-26 Dexcom, Inc. Signal processing for continuous analyte sensor
US9420965B2 (en) 2003-12-09 2016-08-23 Dexcom, Inc. Signal processing for continuous analyte sensor
US9192328B2 (en) 2003-12-09 2015-11-24 Dexcom, Inc. Signal processing for continuous analyte sensor
US9364173B2 (en) 2003-12-09 2016-06-14 Dexcom, Inc. Signal processing for continuous analyte sensor
US9498155B2 (en) 2003-12-09 2016-11-22 Dexcom, Inc. Signal processing for continuous analyte sensor
US11638541B2 (en) 2003-12-09 2023-05-02 Dexconi, Inc. Signal processing for continuous analyte sensor
US9750441B2 (en) 2003-12-09 2017-09-05 Dexcom, Inc. Signal processing for continuous analyte sensor
US9351668B2 (en) 2003-12-09 2016-05-31 Dexcom, Inc. Signal processing for continuous analyte sensor
US8808228B2 (en) 2004-02-26 2014-08-19 Dexcom, Inc. Integrated medicament delivery device for use with continuous analyte sensor
US10835672B2 (en) 2004-02-26 2020-11-17 Dexcom, Inc. Integrated insulin delivery system with continuous glucose sensor
US10278580B2 (en) 2004-02-26 2019-05-07 Dexcom, Inc. Integrated medicament delivery device for use with continuous analyte sensor
US9937293B2 (en) 2004-02-26 2018-04-10 Dexcom, Inc. Integrated delivery device for continuous glucose sensor
US10966609B2 (en) 2004-02-26 2021-04-06 Dexcom, Inc. Integrated medicament delivery device for use with continuous analyte sensor
US11246990B2 (en) 2004-02-26 2022-02-15 Dexcom, Inc. Integrated delivery device for continuous glucose sensor
US8926585B2 (en) 2004-02-26 2015-01-06 Dexcom, Inc. Integrated delivery device for continuous glucose sensor
US10993642B2 (en) 2004-07-13 2021-05-04 Dexcom, Inc. Analyte sensor
US10813576B2 (en) 2004-07-13 2020-10-27 Dexcom, Inc. Analyte sensor
US10827956B2 (en) 2004-07-13 2020-11-10 Dexcom, Inc. Analyte sensor
US10918314B2 (en) 2004-07-13 2021-02-16 Dexcom, Inc. Analyte sensor
US10799158B2 (en) 2004-07-13 2020-10-13 Dexcom, Inc. Analyte sensor
US9247900B2 (en) 2004-07-13 2016-02-02 Dexcom, Inc. Analyte sensor
US10918315B2 (en) 2004-07-13 2021-02-16 Dexcom, Inc. Analyte sensor
US10918313B2 (en) 2004-07-13 2021-02-16 Dexcom, Inc. Analyte sensor
US10932700B2 (en) 2004-07-13 2021-03-02 Dexcom, Inc. Analyte sensor
US10799159B2 (en) 2004-07-13 2020-10-13 Dexcom, Inc. Analyte sensor
US10980452B2 (en) 2004-07-13 2021-04-20 Dexcom, Inc. Analyte sensor
US9668677B2 (en) 2004-07-13 2017-06-06 Dexcom, Inc. Analyte sensor
US9060742B2 (en) 2004-07-13 2015-06-23 Dexcom, Inc. Transcutaneous analyte sensor
US10722152B2 (en) 2004-07-13 2020-07-28 Dexcom, Inc. Analyte sensor
US10993641B2 (en) 2004-07-13 2021-05-04 Dexcom, Inc. Analyte sensor
US10709363B2 (en) 2004-07-13 2020-07-14 Dexcom, Inc. Analyte sensor
US10709362B2 (en) 2004-07-13 2020-07-14 Dexcom, Inc. Analyte sensor
US11026605B1 (en) 2004-07-13 2021-06-08 Dexcom, Inc. Analyte sensor
US11045120B2 (en) 2004-07-13 2021-06-29 Dexcom, Inc. Analyte sensor
US10524703B2 (en) 2004-07-13 2020-01-07 Dexcom, Inc. Transcutaneous analyte sensor
US9801572B2 (en) 2004-07-13 2017-10-31 Dexcom, Inc. Transcutaneous analyte sensor
US11064917B2 (en) 2004-07-13 2021-07-20 Dexcom, Inc. Analyte sensor
US11883164B2 (en) 2004-07-13 2024-01-30 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10314525B2 (en) 2004-07-13 2019-06-11 Dexcom, Inc. Analyte sensor
WO2006024672A1 (en) * 2004-09-03 2006-03-09 Novo Nordisk A/S System and method for estimating the glucose concentration in blood
US10918318B2 (en) 2005-03-10 2021-02-16 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10918317B2 (en) 2005-03-10 2021-02-16 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US11051726B2 (en) 2005-03-10 2021-07-06 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US11000213B2 (en) 2005-03-10 2021-05-11 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US9078608B2 (en) 2005-03-10 2015-07-14 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10925524B2 (en) 2005-03-10 2021-02-23 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10918316B2 (en) 2005-03-10 2021-02-16 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US9220449B2 (en) 2005-03-10 2015-12-29 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US9918668B2 (en) 2005-03-10 2018-03-20 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US9314196B2 (en) 2005-03-10 2016-04-19 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10898114B2 (en) 2005-03-10 2021-01-26 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10856787B2 (en) 2005-03-10 2020-12-08 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10743801B2 (en) 2005-03-10 2020-08-18 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10716498B2 (en) 2005-03-10 2020-07-21 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10709364B2 (en) 2005-03-10 2020-07-14 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10617336B2 (en) 2005-03-10 2020-04-14 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10610137B2 (en) 2005-03-10 2020-04-07 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10610135B2 (en) 2005-03-10 2020-04-07 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10610136B2 (en) 2005-03-10 2020-04-07 Dexcom, Inc. System and methods for processing analyte sensor data for sensor calibration
US10813577B2 (en) 2005-06-21 2020-10-27 Dexcom, Inc. Analyte sensor
US10952652B2 (en) 2005-11-01 2021-03-23 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US10201301B2 (en) 2005-11-01 2019-02-12 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US11911151B1 (en) 2005-11-01 2024-02-27 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8920319B2 (en) 2005-11-01 2014-12-30 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9078607B2 (en) 2005-11-01 2015-07-14 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US8915850B2 (en) 2005-11-01 2014-12-23 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US11399748B2 (en) 2005-11-01 2022-08-02 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US11363975B2 (en) 2005-11-01 2022-06-21 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US11103165B2 (en) 2005-11-01 2021-08-31 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US10231654B2 (en) 2005-11-01 2019-03-19 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9326716B2 (en) 2005-11-01 2016-05-03 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US11272867B2 (en) 2005-11-01 2022-03-15 Abbott Diabetes Care Inc. Analyte monitoring device and methods of use
US9724028B2 (en) 2006-02-22 2017-08-08 Dexcom, Inc. Analyte sensor
US8968198B2 (en) 2006-02-22 2015-03-03 Dexcom, Inc. Analyte sensor
WO2007138154A1 (en) * 2006-05-29 2007-12-06 Wristop Technologies Oy Apparatus and method for dosing drug and wireless remote control of a drug pump
US11432772B2 (en) 2006-08-02 2022-09-06 Dexcom, Inc. Systems and methods for replacing signal artifacts in a glucose sensor data stream
US10349873B2 (en) 2006-10-04 2019-07-16 Dexcom, Inc. Analyte sensor
US11382539B2 (en) 2006-10-04 2022-07-12 Dexcom, Inc. Analyte sensor
US11399745B2 (en) 2006-10-04 2022-08-02 Dexcom, Inc. Dual electrode system for a continuous analyte sensor
US9451908B2 (en) 2006-10-04 2016-09-27 Dexcom, Inc. Analyte sensor
US11373347B2 (en) 2007-06-08 2022-06-28 Dexcom, Inc. Integrated medicament delivery device for use with continuous analyte sensor
US10403012B2 (en) 2007-06-08 2019-09-03 Dexcom, Inc. Integrated medicament delivery device for use with continuous analyte sensor
US9741139B2 (en) 2007-06-08 2017-08-22 Dexcom, Inc. Integrated medicament delivery device for use with continuous analyte sensor
EP4276851A3 (en) * 2007-06-21 2024-01-24 University Of Virginia Patent Foundation Lqg artificial pancreas control system and related method
DE102007046864A1 (en) * 2007-09-28 2009-04-09 Endress + Hauser Conducta Gesellschaft für Mess- und Regeltechnik mbH + Co. KG Measured value plausibility and quality evaluating method for e.g. control circuit, involves forming and evaluating derivative of measured value according to time, where measured value is compared with predefined warning limit value
US11744943B2 (en) 2007-10-09 2023-09-05 Dexcom, Inc. Integrated insulin delivery system with continuous glucose sensor
US11160926B1 (en) 2007-10-09 2021-11-02 Dexcom, Inc. Pre-connected analyte sensors
US10653835B2 (en) 2007-10-09 2020-05-19 Dexcom, Inc. Integrated insulin delivery system with continuous glucose sensor
US8977504B2 (en) 2007-10-12 2015-03-10 Cambridge Enterprise Limited Substance monitoring and control in human or animal bodies
EP2223251B1 (en) * 2007-10-12 2018-07-18 Cambridge Enterprise Limited Method and apparatus for monitoring a substance in human or animal bodies
US10182751B2 (en) 2007-10-25 2019-01-22 Dexcom, Inc. Systems and methods for processing sensor data
US11272869B2 (en) 2007-10-25 2022-03-15 Dexcom, Inc. Systems and methods for processing sensor data
US9717449B2 (en) 2007-10-25 2017-08-01 Dexcom, Inc. Systems and methods for processing sensor data
US11342058B2 (en) 2007-12-17 2022-05-24 Dexcom, Inc. Systems and methods for processing sensor data
US9135402B2 (en) 2007-12-17 2015-09-15 Dexcom, Inc. Systems and methods for processing sensor data
US9149233B2 (en) 2007-12-17 2015-10-06 Dexcom, Inc. Systems and methods for processing sensor data
US9149234B2 (en) 2007-12-17 2015-10-06 Dexcom, Inc. Systems and methods for processing sensor data
US9901307B2 (en) 2007-12-17 2018-02-27 Dexcom, Inc. Systems and methods for processing sensor data
US9339238B2 (en) 2007-12-17 2016-05-17 Dexcom, Inc. Systems and methods for processing sensor data
US9839395B2 (en) 2007-12-17 2017-12-12 Dexcom, Inc. Systems and methods for processing sensor data
US10827980B2 (en) 2007-12-17 2020-11-10 Dexcom, Inc. Systems and methods for processing sensor data
US10506982B2 (en) 2007-12-17 2019-12-17 Dexcom, Inc. Systems and methods for processing sensor data
US9626480B2 (en) 2008-03-21 2017-04-18 Lifescan Scotland Limited Analyte testing method and system
US8917184B2 (en) 2008-03-21 2014-12-23 Lifescan Scotland Limited Analyte testing method and system
US10980461B2 (en) 2008-11-07 2021-04-20 Dexcom, Inc. Advanced analyte sensor calibration and error detection
US10537678B2 (en) 2009-03-27 2020-01-21 Dexcom, Inc. Methods and systems for promoting glucose management
US10610642B2 (en) 2009-03-27 2020-04-07 Dexcom, Inc. Methods and systems for promoting glucose management
US9446194B2 (en) 2009-03-27 2016-09-20 Dexcom, Inc. Methods and systems for promoting glucose management
US10675405B2 (en) 2009-03-27 2020-06-09 Dexcom, Inc. Methods and systems for simulating glucose response to simulated actions
US11437141B2 (en) 2009-04-30 2022-09-06 Dexcom, Inc. Performance reports associated with continuous sensor data from multiple analysis time periods
EP2762073A1 (en) * 2010-01-26 2014-08-06 Roche Diagnostics GmbH Methods and systems for processing glucose data measured from a person having diabetes
WO2011092133A1 (en) * 2010-01-26 2011-08-04 Roche Diagnostics Gmbh Methods and systems for processing glucose data measured from a person having diabetes
US8843321B2 (en) 2010-01-26 2014-09-23 Roche Diagnostics Operations, Inc. Methods and systems for processing glucose data measured from a person having diabetes
CN104605862B (en) * 2010-01-26 2019-01-18 霍夫曼-拉罗奇有限公司 For handling the method and system of the glucose data measured from tool person with diabetes
CN102740770A (en) * 2010-01-26 2012-10-17 霍夫曼-拉罗奇有限公司 Methods and systems for processing glucose data measured from a person having diabetes
CN104605862A (en) * 2010-01-26 2015-05-13 霍夫曼-拉罗奇有限公司 Methods and systems for processing glucose data measured from a person having diabetes
EP3799073A1 (en) * 2011-08-26 2021-03-31 The University of Virginia Patent Foundation Method, system and computer readable medium for adaptive advisory control of diabetes
EP2748747A4 (en) * 2011-08-26 2015-07-08 Univ Virginia Patent Found Method, system and computer readable medium for adaptive advisory control of diabetes
US11024429B2 (en) 2011-08-26 2021-06-01 University Of Virginia Patent Foundation Method, system and computer readable medium for adaptive and advisory control of diabetes
CN104602599A (en) * 2012-11-01 2015-05-06 泰尔茂株式会社 Sensing device and sensing method
EP3381364B1 (en) * 2015-11-25 2023-03-29 Nippon Telegraph and Telephone Corporation Respiratory estimation method and device
WO2017129634A1 (en) * 2016-01-25 2017-08-03 Prediktor Medical As Estimating a substance concentration using multiple non-invasive sensors
CN107924715A (en) * 2016-01-25 2018-04-17 普雷迪科特医疗有限公司 Estimate material concentration using multiple non-invasive sensors
US10575790B2 (en) 2016-03-02 2020-03-03 Roche Diabetes Care, Inc. Patient diabetes monitoring system with clustering of unsupervised daily CGM profiles (or insulin profiles) and method thereof
US10478556B2 (en) 2016-03-04 2019-11-19 Roche Diabetes Care, Inc. Probability based controller gain
US10311976B2 (en) 2016-04-28 2019-06-04 Roche Diabetes Care, Inc. Bolus calculator with probabilistic carbohydrate measurements
US10332632B2 (en) 2016-06-01 2019-06-25 Roche Diabetes Care, Inc. Control-to-range failsafes
US10332633B2 (en) 2016-06-01 2019-06-25 Roche Diabetes Care, Inc. Control-to-range aggressiveness
US10297350B2 (en) 2016-06-01 2019-05-21 Roche Diabetes Care, Inc. Risk-based control-to-range
CN106139318A (en) * 2016-06-28 2016-11-23 成都市亿泰科技有限公司 A kind of intelligent transfusion monitoring method based on integrated kalman filter
CN107203700A (en) * 2017-07-14 2017-09-26 清华-伯克利深圳学院筹备办公室 A kind of method and device monitored based on continuous blood sugar
CN107203700B (en) * 2017-07-14 2020-05-05 清华-伯克利深圳学院筹备办公室 Method and device based on continuous blood glucose monitoring
US11382540B2 (en) 2017-10-24 2022-07-12 Dexcom, Inc. Pre-connected analyte sensors
US11706876B2 (en) 2017-10-24 2023-07-18 Dexcom, Inc. Pre-connected analyte sensors
US11350862B2 (en) 2017-10-24 2022-06-07 Dexcom, Inc. Pre-connected analyte sensors
US11331022B2 (en) 2017-10-24 2022-05-17 Dexcom, Inc. Pre-connected analyte sensors
CN108375381A (en) * 2018-02-08 2018-08-07 北方工业大学 A kind of bionical polarization sensor multi-source error calibrating method based on Extended Kalman filter
CN113503879A (en) * 2021-07-09 2021-10-15 北京航空航天大学 Dynamic adaptive Kalman filter method based on ensemble empirical mode decomposition

Also Published As

Publication number Publication date
AU2001291189A1 (en) 2002-04-02
US20020099282A1 (en) 2002-07-25
US6572545B2 (en) 2003-06-03
US20030195404A1 (en) 2003-10-16
US20020111547A1 (en) 2002-08-15
US6575905B2 (en) 2003-06-10

Similar Documents

Publication Publication Date Title
US6575905B2 (en) Method and apparatus for real-time estimation of physiological parameters
US20040034295A1 (en) Method and apparatus for real-time estimation and control of physiological parameters
US11189381B2 (en) Moving-horizon state-initializer for control applications
EP2516671B1 (en) Systems for managing drug delivery devices
EP2223251B1 (en) Method and apparatus for monitoring a substance in human or animal bodies
Incremona et al. Model predictive control with integral action for artificial pancreas
EP1881786B1 (en) Fully automated control system for type 1 diabetes
US20100262117A1 (en) Predictive control based system and method for control of insulin delivery in diabetes using glucose sensing
JP7042808B2 (en) An automated system to control a patient&#39;s blood glucose level
US11213249B2 (en) Automated system for controlling the blood glucose level of a patient
Barcelo-Rico et al. Adaptive calibration algorithm for plasma glucose estimation in continuous glucose monitoring
KR102601111B1 (en) Closed-loop blood glucose control system and method
Boiroux et al. Parameter estimation in type 1 diabetes models for model-based control applications
US20220211329A1 (en) Method and system for enhancing glucose prediction
Kirchsteiger et al. LMI-based approaches for the calibration of continuous glucose measurement sensors
US11883162B2 (en) System for predicting a patient&#39;s blood glucose level
Garnica et al. Noise spectral analysis and error estimation of continuous glucose monitors under real-life conditions of diabetes patients
Kirchsteiger et al. A novel online recalibration strategy for continuous glucose measurement sensors employing LMI techniques
Gireesh et al. Blood glucose prediction algorithms for hypoglycemic and/or hyperglycemic alerts
Al Ahdab et al. Stochastic Modeling and Inference for Type 2 Diabetes
Aguirre-Zapata et al. Glycemic monitoring in critical care using nonlinear state estimators
Halvorsen Robust Estimation of Blood Glucose Level Based on Subcutaneous Sensors and Analysis of Blood Samples
Sutradhar et al. On-line Identification of ARX model for Glucose-Insulin Interaction in a Diabetic Patient from Input-Output Data
US20210260285A1 (en) Automated system for monitoring a patient&#39;s blood sugar
JP2021529059A (en) An automated system to control a patient&#39;s blood sugar

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ CZ DE DE DK DK DM DZ EC EE EE ES FI FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ PH PL PT RO RU SD SE SG SI SK SK SL TJ TM TR TT TZ UA UG UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

122 Ep: pct application non-entry in european phase
NENP Non-entry into the national phase

Ref country code: JP