CA2479107A1 - Analytical device with prediction module and related methods - Google Patents
Analytical device with prediction module and related methods Download PDFInfo
- Publication number
- CA2479107A1 CA2479107A1 CA002479107A CA2479107A CA2479107A1 CA 2479107 A1 CA2479107 A1 CA 2479107A1 CA 002479107 A CA002479107 A CA 002479107A CA 2479107 A CA2479107 A CA 2479107A CA 2479107 A1 CA2479107 A1 CA 2479107A1
- Authority
- CA
- Canada
- Prior art keywords
- isf
- subject
- analyte concentration
- rate
- series
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/14—Devices for taking samples of blood ; Measuring characteristics of blood in vivo, e.g. gas concentration within the blood, pH-value of blood
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring 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/14532—Measuring 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/15—Devices for taking samples of blood
- A61B5/150007—Details
- A61B5/150015—Source of blood
- A61B5/150022—Source of blood for capillary blood or interstitial fluid
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/15—Devices for taking samples of blood
- A61B5/150007—Details
- A61B5/150358—Strips for collecting blood, e.g. absorbent
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
- A61B2560/04—Constructional details of apparatus
- A61B2560/0431—Portable apparatus, e.g. comprising a handle or case
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0295—Strip shaped analyte sensors for apparatus classified in A61B5/145 or A61B5/157
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/022—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring 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/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring 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/1486—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using enzyme electrodes, e.g. with immobilised oxidase
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/15—Devices for taking samples of blood
- A61B5/157—Devices characterised by integrated means for measuring characteristics of blood
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6824—Arm or wrist
-
- Y—GENERAL 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
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10T—TECHNICAL SUBJECTS COVERED BY FORMER US CLASSIFICATION
- Y10T436/00—Chemistry: analytical and immunological testing
- Y10T436/11—Automated chemical analysis
Abstract
An analytical device for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid (ISF) analyte concentration includes an ISF sampling module, an analysis module and a prediction module. The ISF sampling module is configured to sequentially extract a plurality of ISF samples from a subject. The analysis module is configured to sequentially determining an ISF analyte concentration (e.g., ISF
glucose concentration) in each of the ISF samples, resulting in a series of ISF
analyte concentrations.
The prediction module is configured for storing the series of ISF analyte concentrations and predicting the subject's whole blood analyte concentration based on the series by performing at least one algorithm. A method for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid analyte concentration includes extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner and sequentially determining an ISF analyte concentration in each of the plurality of ISF
samples to create a series of ISF analyte concentrations. The subject's blood analyte concentration is than predicted based on the series of ISF analyte concentrations by performing at least one algorithm.
glucose concentration) in each of the ISF samples, resulting in a series of ISF
analyte concentrations.
The prediction module is configured for storing the series of ISF analyte concentrations and predicting the subject's whole blood analyte concentration based on the series by performing at least one algorithm. A method for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid analyte concentration includes extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner and sequentially determining an ISF analyte concentration in each of the plurality of ISF
samples to create a series of ISF analyte concentrations. The subject's blood analyte concentration is than predicted based on the series of ISF analyte concentrations by performing at least one algorithm.
Description
ANALYTICAL DEVICE WITH PREDICTION MODULE
AND RELATED METHODS
BACKGROUND OF INVENTION
[QOOI] This application is a continuation-in part of U.S. Application No.
10/653,023, filed August 28, 2003, which claims the benefit of U.S.
Provisional Application No. 60/476,733, filed June 6, 2003.
1. Field of the Invention [0002] The present invention relates, in general, to analytical devices and, in particular, to analytical devices and associated methods for predicting a subject's blood analyte concentration from a subject's interstitial fluid (ISF) analyte concentration.
Description of the Related Art [0003] In the field of analyte (e.g., glucose) monitoring, continuous or semi-continuous analytical devices and methods are advantageous in that they provide enhanced insight into analyte concentration trends, a subject's overall analyte control and the effect of food, exercise and/or medication on an analyte's concentration. In practice, however, such analytical devices can have drawbacks. For example, interstitial fluid (ISF) analytical devices can suffer inaccuracies due to, for instance, physiological lag (i.e., the time-dependent difference between a subject's ISF analyte concentration and a subject's blood analyte concentration) andlor bias effects (i.e., the fluid characteristic-dependent difference between a subject's ISF anaiyte concentration and a subject's blood anaIyte concentration).
AND RELATED METHODS
BACKGROUND OF INVENTION
[QOOI] This application is a continuation-in part of U.S. Application No.
10/653,023, filed August 28, 2003, which claims the benefit of U.S.
Provisional Application No. 60/476,733, filed June 6, 2003.
1. Field of the Invention [0002] The present invention relates, in general, to analytical devices and, in particular, to analytical devices and associated methods for predicting a subject's blood analyte concentration from a subject's interstitial fluid (ISF) analyte concentration.
Description of the Related Art [0003] In the field of analyte (e.g., glucose) monitoring, continuous or semi-continuous analytical devices and methods are advantageous in that they provide enhanced insight into analyte concentration trends, a subject's overall analyte control and the effect of food, exercise and/or medication on an analyte's concentration. In practice, however, such analytical devices can have drawbacks. For example, interstitial fluid (ISF) analytical devices can suffer inaccuracies due to, for instance, physiological lag (i.e., the time-dependent difference between a subject's ISF analyte concentration and a subject's blood analyte concentration) andlor bias effects (i.e., the fluid characteristic-dependent difference between a subject's ISF anaiyte concentration and a subject's blood anaIyte concentration).
[0004] Conventional ISF analytical devices can employ ISF samples obtained from various sites on a subject's body and from variaus penetration depths in a subject's skin. The use of various sites and penetration depths for obtaining an ISF sample cawbe a contxibuting factor in an ISF analytical devices' inaccuracy. In addition, other analytically relevant properties of an ISF
sample can be influenced by the site and/or penetration depth at which the 1SF sample is collected. For example, ISF collected from the subcutaneous region of a subject's skin can be more prone to containing contaminating substances such as triglycerides, which can affect analyte analysis in terms of volume ezxor and sensor fouling.
sample can be influenced by the site and/or penetration depth at which the 1SF sample is collected. For example, ISF collected from the subcutaneous region of a subject's skin can be more prone to containing contaminating substances such as triglycerides, which can affect analyte analysis in terms of volume ezxor and sensor fouling.
[0005] Furthermore, conventional ISF analytical devices can require inconvenient and cumbersome calibration procedures involving samples of capillary blood.
[0006] Still needed in the field, therefore, is an analytical device and associated method with reduced inaccuracy due to physiological lag and bias effects. In addition, the analytical device and associated methods should not require samples of capillary blood for calibration.
SUMMARY OF INVENTION
SUMMARY OF INVENTION
[0007] Embodiments of the present invention include analytical devices and methods that accurately account for physiological lag and bias effects. In addition, the analytical device and associated methods do not require samples of capillary blood for calibration.
[0008] An analytical device for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid (ISF) analyte concentration according to an exemplary embodiment of the present invention includes an ISF
sampling module, an analysis module and a prediction module.
sampling module, an analysis module and a prediction module.
[0009] The ISF sampling module is configured to extract a plurality of ISF
samples from a subject in a sequential manner. The analysis module is configured to sequentially determining an ISF analyte concentration (e.g., ISF glucose concentration) in each of the plurality of ISF samples. The result of this sequential determination is a series of ISF analyte concentrations. The prediction module is configured for storing the series of ISF analyte concentrations and predicting the subject's whole blood analyte concentration based on the series of ISF analyte concentrations by performing at least one algorithm.
samples from a subject in a sequential manner. The analysis module is configured to sequentially determining an ISF analyte concentration (e.g., ISF glucose concentration) in each of the plurality of ISF samples. The result of this sequential determination is a series of ISF analyte concentrations. The prediction module is configured for storing the series of ISF analyte concentrations and predicting the subject's whole blood analyte concentration based on the series of ISF analyte concentrations by performing at least one algorithm.
[00010] An exemplary embodiment of a method for predicting a subject's whole blood analyte concentration based on the subject's interstitial analyte concentration according to the present invention includes extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner and determining an ISF
analyte concentration in each of the plurality of ISF samples in a sequential manner to create a series of ISF analyte concentrations. The subject's blood analyte concentration is then predicted based on the series of ISF analyze concentrations by performing at least one algorithm.
analyte concentration in each of the plurality of ISF samples in a sequential manner to create a series of ISF analyte concentrations. The subject's blood analyte concentration is then predicted based on the series of ISF analyze concentrations by performing at least one algorithm.
[00011] Embodiments of analytical devices and methods according to the present invention predict a subject's blood analyte concentration based solely on a series of ISF analyte concentrations derived from ISF samples extracted in a continuous or semi-continuous manner. The analytical devices and methods do so using an algorithm that predicts the subject's blood analyte concentration based on the series of ISF analyte concentrations. The algorithm accounts for physiological lag and bias effects. In addition, the analytical device does not require calibration using capillary blood.
BRIEF DESCRIPTION OF DRAWINGS
BRIEF DESCRIPTION OF DRAWINGS
[00012] A better understanding of the features and advantages of the present invention will be obtained by reference to the following detailed description that sets forth illustrative embodiments, in which the principles of the invention are utilized, and the accompanying drawings of which:
FIG. 1 is a block diagram of an analytical device for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid (ISF) analyte concentration according to an exemplary embodiment of the present invention;
FIG. 2 is a Clarke Error Grid Plot for interpolate finger blood glucose (reference) versus ISF glucose concentration (ISFo);
FIG. 3 is a Clarke Error Grid Plot for interpolated finger blood glucose versus predicted finger blood glucose for an algorithm (i.e., Eqn 1) that can be' employed in analytical devices and methods according to the present invention;
FIG. 4 is a Clarke Error Grid Plot for interpolated finger blood glucose versus predicted finger blood glucose for another algorithm (i.e., Eqn 2) that can be employed in analytical devices and methods according to the present invention;
and FIG. 5 is a flow chart illustrating a sequence of steps in a process according to an exemplary embodiment of the present invention.
DETAILED DESCRIPTION OF THE INVENTION
FIG. 1 is a block diagram of an analytical device for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid (ISF) analyte concentration according to an exemplary embodiment of the present invention;
FIG. 2 is a Clarke Error Grid Plot for interpolate finger blood glucose (reference) versus ISF glucose concentration (ISFo);
FIG. 3 is a Clarke Error Grid Plot for interpolated finger blood glucose versus predicted finger blood glucose for an algorithm (i.e., Eqn 1) that can be' employed in analytical devices and methods according to the present invention;
FIG. 4 is a Clarke Error Grid Plot for interpolated finger blood glucose versus predicted finger blood glucose for another algorithm (i.e., Eqn 2) that can be employed in analytical devices and methods according to the present invention;
and FIG. 5 is a flow chart illustrating a sequence of steps in a process according to an exemplary embodiment of the present invention.
DETAILED DESCRIPTION OF THE INVENTION
[00013] FIG. 1 is a block diagram of an analytical device 104 (within the dashed lines) for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid (ISF) analyte concentration according to an exemplary embodiment of the present invention. Analytical device 100 includes an interstitial fluid (ISF) sampling module 1,10, an analysis module 120 and a prediction module 130. The arrows of FIG. 1 indicate operative communication between the ISF sampling, analysis and prediction modules.
[00014] ISF sampling module 1 I0 is configured to extract a plurality of ISF
samples from the subject in a sequential manner and to deliver the ISF samples to analysis module 120. ISF sampling module 110 can, for example, extract ISF samples in a sequential manner with a time interval between samples in the range of 5'minutes to 15 minutes.
[00015] Analysis module 120 is adapted for sequentially determining an ISF
analyte concentration in each of the plurality of ISF samples extracted by ISF
sampling module 1 I0. The result of such a sequential determination is a series of ISF
analyte concentrations. , Such a series of ISF analyte concentrations will typically include a plurality of ISF analyte concentrations (also refewed to as "values"), each associated with a time that corresponds to the time at which the ISF sample was extracted by ISF
sampling module 110. Analysis module 120 can also be configured to transfer the series of ISF analyte concentrations to prediction module 130, either individually or as a group. ' , [00016] Interstitial sampling module 110 can take any suitable form known to one skilled in the art including, but not limited to, sampling modules and ISF
extraction devices described in co-pending U.S. Provisional Patent Application No.60/476,733 (f led June 6, 2003}, and sampling modules described in International Application PCT/GB01/05634 (as published as International Publication Na. WO 02/49507 A1 on June 27, 2002), both of which are hereby fully incorporated herein by reference.
samples from the subject in a sequential manner and to deliver the ISF samples to analysis module 120. ISF sampling module 110 can, for example, extract ISF samples in a sequential manner with a time interval between samples in the range of 5'minutes to 15 minutes.
[00015] Analysis module 120 is adapted for sequentially determining an ISF
analyte concentration in each of the plurality of ISF samples extracted by ISF
sampling module 1 I0. The result of such a sequential determination is a series of ISF
analyte concentrations. , Such a series of ISF analyte concentrations will typically include a plurality of ISF analyte concentrations (also refewed to as "values"), each associated with a time that corresponds to the time at which the ISF sample was extracted by ISF
sampling module 110. Analysis module 120 can also be configured to transfer the series of ISF analyte concentrations to prediction module 130, either individually or as a group. ' , [00016] Interstitial sampling module 110 can take any suitable form known to one skilled in the art including, but not limited to, sampling modules and ISF
extraction devices described in co-pending U.S. Provisional Patent Application No.60/476,733 (f led June 6, 2003}, and sampling modules described in International Application PCT/GB01/05634 (as published as International Publication Na. WO 02/49507 A1 on June 27, 2002), both of which are hereby fully incorporated herein by reference.
[00017] Furthermore, analysis module 120 can also take any suitable form known to one skilled in the art including those described in co-pending U.S.
Provisional Patent Application No. 601476,733 (filed June 6, 2003) and in International Application PCT/GBOl/05634 (published as International Publication No. WO 02/49507 A1 on June 27, 2002). In addition, the analysis module can include any suitable analyte sensor known to those of skill in the art including, but not limited to, glucose analyte sensors based on photometric or electrochemical analytical techniques.
[oOOlsJ Prediction module I30 is canfigured for storing the series of ISF
analyte concentrations created by the analysis module and predicting the subject's whole blood analyte concentration by performing at least one algorithm of the following general form {referred to as Eqn 1):
PC = f (ISF;k, ratel, significant interaction terms) (Eqn 1) where:
PC is a predicted subject's whole blood analyte concentration;
i is an integer with predetermined values selected from the values of, for example, 0,1, 2, 3, 4 and 5;
j is an integer with predetermined values selected from the values of, for example, 1, 2, 3, 4 and 5;
k is an integers) with predetermined values selected from the values of, for example, l and 2;
ISF; is a measured ISF analyte concentration in the series of ISF analyze concentrations, with the subscript (i) indicating which ISF value is being refen~ed to, i.e., i = 0 indicates the most recently measured ISF analyte, i =1 indicates one value back in the series of ISF analyte concentrations, 2 = two values back in the series of ISF
analyte concentrations, etc.;
rated is the rate of change between immediately adjacent ISF~analyte concentrations in the series of ISF analyte concentrations (calculated as the difference between the immediately adjacent ISF concentrations divided by the time difference between when the immediately adjacent ISF concentrations were measured by the analysis module) with the subscript (j) referring to which immediately adjacent ISF
values are used to calculate the rate, i.e., when j =1 indicates the rate between current ISF value and the previous ISF value, j = 2 is indicative of the rate between the ISF
values one previous and two previous relative to the current ISF value, etc.;
and significant interaction terms = interaction ternns involving at least two of ISF;k, rate.
[00419] The mathematical form of the function (~ employed in Eqn 1 can be any suitable mathematical form that accounts for physiological lag between ISF
analyte concentration and blood analyte concentration as well as any bias effect between ISF
and blood analyte concentrations. However, it has been determined that such a relationship is suitably accurate when measured ISF analyte concentrations (i.e., ISF;~'), rates {ratei) and interaction terms are included. The use of rates is particularly beneficial in providing an accurate algorithm, and hence an accurate analytical device, since the time interval between the ISF analyte concentrations can be non-uniform (e.g., the time interval could vary between 5 minutes and 15 minutes).
[00020] The form of the function (~ can determined by, for example, a least squares regression analysis of a statistically relevant number of ISF analyte concentrations and associated blood anaIyte concentrations. Those skilled in the art will appreciate that any number of mathematical methods (e.g., mathematical modeling methods) can be used to analyze such data and arrive at a suitable function ~. For example, linear and polynomial regression analysis, time series analysis, or neural networks can be used.
In the circumstance that the analyte is glucose, ISF glucose concentrations and blood glucose concentrations can be determined from ISF and blood samples extracted from diabetic subjects that have ingested glucose.
[00021] If desired, a suitable algorithm can be obtained using a mathematical modeling method that includes weighting factors to provide for greater accuracy at lower analyte concentrations values, to account for curvature in the response, and/or to account for noise in the modeled data. Weighting of input observations can also be similarly beneficial in such mathematical modeling methods.
[00022] The determination of suitable weighting factors can be, for example, an iterative process in which a weighting factors) is applied in a model, the weighting factor's effect on model results observed, and the weighting factors) adjusted based on model error reduction. The choice of weighting factors in the mathematical modeling method can also be determined, for example, by the relative importance of data ranges andlor trending direction. For example, when glucose is the analyte of interest, greater accuracy for the low end of the physiological glucose concentration range may be deemed important, and thus a weighting factor that enhances the importance of lower glucose concentrations can be employed. Such an enhancement can be accomplished, far example, by multiplying observed glucose concentrations by the inverse of the observed value raised to a predetermined power. Similarly, weighting factors can be determined which will enhance the importance of certain events or trends in observed values, such as a magnitude of the gradient of an observed rate and/or a change in direction. Furthermore, prospective weighting factors can also be arbitrarily chosen with suitable weighting factors chosen from the prospective weighting factors based on their effect on model error reduction.
[00023] Prediction module 130 can take any suitable forn~ known to one skilled in the art including, but not limited to, the remote control modules described in co-pending U.S. Provisional Patent Application No. 60/476,733 (filed June 6, 2003), which is hereby fully incorporated herein by reference.
[00024] As an alternative to the use of Eqn 1 above, prediction module 130 can be configured for storing the series of ISF analyte concentrations created by the analysis module and predicting the subject's whole blood analyte concentration by performing at least one algorithm of the following general form (referred to as Eqn 2):
PC = f (ISF;k, ratei, ma~,rate~,~, significant interaction terms) (Eqn 2) where:
n and m are integers with predetermined values selected from the values of, for example, 1, 2 and 3;
p is an integers) with predetermined values selected from the values of, for example, l and 2; and ma"ratem is the moving average rate between immediately adjacent averages of groupings of ISF values, with the subscript (n) referring to the number of ISF
values included in the moving average and the subscript (m) referring to the time position of the adjacent average values relative to the current values as follows (with a further explanation in the next paragraph):
n+1 = the number of points used in the moving average rate;
m-1 = first point included in the moving average. If m-1= 0, then the current ISF value is used as the first point in the moving average calculation.
n + m always adds up to the number of points back (or removed from the current ISF value) that will be needed for calculating the moving average calculation.
n and m are integers with predetermined values selected from the values of, for example, 1, 2 and 3; and significant interaction terms = interaction terms involving at least two of ISF;k, ratej, and ma"ratemp. , (00025] The following example illustrates the concept of the moving average rate (ma"ratero) employed in Eqn 2 above. For the exemplary moving average rate malratel, the moving average rate is a 2-point moving average rate {since m+n =1+1 = 2) with the moving average rate calculated iietween the average of the grouping that includes the most recent iSF concentration and the ISF concentration one point back and the average of the grouping that includes the ISF concentrations one and two points prior to the most recent ISF concentration.
(00026) Eqn 2 includes moving average rates (i.e., ma"rate"~ to smooth the data (i:e., the series of ISF analyte concentrations and/or rates) with respect to both rate and the trending direction of an analyte concentration, thereby removing noise from the data and increasing fhe analytical device's accuracy. Although significant (i.e., major) changes in adjacent ISF values can be regarded as important in terms of algorithm accuracy, significant changes can also be due to noise that can adversely effect an algorithm's accuracy. The moving average rate, which is the rate of change between the means of adjacent (and overlapping) groupings of ISF values, dampens noise caused by outlier values that do not represent a true trend in the data.
(00027] Examples of suitable algorithms and the techniques used to derive the algorithms are includes in the examples below.
Example 1: Predictive Algorithm for a Glucose Analytical Device Utilizing ISF
k, ratel, and Significant Interaction Terms A data set (i.e., a series of ISF glucose concentrations) was generated using an experimental ISF sampling and analysis modules. The ISF sampling module and analysis module employed to generate the data set were configured to extract an ISF
sample from a subject's dermal layer of skin (i.e., dermis), for example from a subject's forearm, and to measure the glucose concentration in the ISF sample.
The ISF sampling module and analysis module were an integrated unit comprising a one-piece sampling module and a modified OneTouch~ Ultra glucose meter with test strip.
The sampling utilized a 30-gauge cannula and a penetraeion depth of about 1 to 2 mm.
It should be noted that an ISF sample collected from the dermis is considered to have a beneficially reduced physiological lag in comparison to an ISF sample collected from a subcutaneous layer due to the dermis being closer to vascular capillary beds than the subcutaneous layer. ' j00028] The ISF sampling module extracted an approximately 1 p.I. ISF sample from a subject's dermis via the cannula and deposited the ISF sample automatically and directly into a measurement zone of the test strip. After a brief electrochemistry development period, the meter displayed the ISF glucose concentration.
[00029] Prior to the ISF samples being extracted, 2 to 4 pounds of pressure was applied to a subject's dermis for 34 seconds, followed by a 5 minute waiting period to allow blood to perfuse (flow into) the sampling area from which the ISF would be extracted.
This elevated blood-flow in the sampling area has the desirable effect of mitigating the physiological lag between blood glucose concentration and ISF glucose concentration, simply because the sampling area is better perfused with flowing blood.
[00030] Finger stick blood glucose measurements in mg/dL (i.e., blood glucose concentrations) were taken from 20 subjects, followed by measurements of glucose in forearm interstitial fluid (i.e., ISF glucose concentrations) as described above. The finger stick blood measurements were taken approximately 15 minutes apart and each was followed approximately 5 minutes later by an ISF sample extraction and ISF
glucose concentration measurement.
[00031] Approximately thirty (3Q) pairs of observations (i.e., pairs of blood glucose concentration measurements and ISF glucose concentration measurements) were obtained fox each of the 20 subjects. The observations were collected over the course of one day for each~subject, in whom a change in glucose concentration was induced through the ingestion of 75g of glucose. The blood glucose concentration for each to observation represents a finger stick draw occurring approximately five minutes prior to the ISF draw.. Blood glucose concentration at the time of the ISF sample extraction was, therefore, linearly interpolated, with the linearly interpolated value used as a response variable in developing the algorithm below. The final ISF glucose concentration for each subject was excluded during the development of the algorithm due to the inability to accurately interpolate a blood glucose concentration.
[40032] An algorithm of the form identified above as Eqn 2 was developed from the data set using multiple linear regression. The algorithm thus developed weighted lower ISF analyte concentrations more heavily, primarily due to the relative importance of accurately predicting glucose at lower concentrations. The weight used was ISF-0. In the absence of such weighting, higher ISF glucose concentration values produced undesirable variability in the residuals.
[00033] The parameters, estimates, errors, t-values and Pr values for the model were as follows:
Parameter Estimate Error t ValuePr > Its ISFO 0.964114574 0.00900642107.05 <.0001 rate2 3.564454310 0.791431254.50 <.0001 rate2*ratel 1.032526146 0.276843433.73 0.0002 rate3 2.115098820 0.572521873.69 0.0002 ratel*rate3 0.728563905 0.325075672.24 0.0255 rate2*rate3 0.993732089 0.362936362.74 0.0064 rate4 2.620714810 0.480338955.46 <.0001 rate2*rate4 1.149236162 0.380759903.02 0.0027 rate2*rate3*rate4 0.140279472.99 0.0029 0.419884620 rates 1.704279771 0.409084594.17 <.0001 R-squared=.98 [00034] The algorithm, following estimators therefore, has form are the when the employed with significant two decimal places:
. PC = 0.96ISFo +3.56rate2 + 1.03 (rate2*rate~ ) + 2.I lrate3 +0.72 (rates *rate3) .
+0.99rate 4 +1.14(rate2*rate4) + 0.42(rate2*rate2*rate4) + 1.70rate5.
[00035] One skilled in the art will recognize,that the above equation is of the form of Eqn. 1 above with:
~..._ ., .., .<nxe~.._c r. -~av ~"-.~rs_a."-.ymrao~~c>°a ar;a-R~.~,,"".
s,.~~~E. ~~,:~,~h =n.~.n°.~xamP=. ~-n m~ .~.~e.....r......~~,~...._.__.
..,....._..._.-......___.,....,.....w....._...N~"._..",.,._._.-__",..._.~,..,.._..._,.....__.~.
I=0 k=1 j=2,3,4and5 and interaction terms = rate2*ratel,ratet*rate3,rate2*rateq, arid rate2*rate2*rate4.
[00036] A Clarke Error Grid analysis can be employed to determuine the accuracy and suitability of an algorithm for the prediction of a subject's blood glucose concentration. The error grid of such an analysis categorizes an analytical'device's response against a reference value into one of five (5) clinical accuracy zones (i.e., zones A-E). Zone A indicates clinically accurate results, zone B indicates results that are not clinically accurate but pose minimal risk to patient health, and zones C through E indicate clinically inaccurate results that pose increasing potential risk to patient health (see Clarke, William L. et al., Evaluating Clinical Accuracy of Systems for Sedf Monitoring of Blood Glucose, Diabetes Care, Vol: 10 hTo. 5, 622-628 [1987]).
An effective~and accurate blood glucose monitoring device should have greater than approximately 85-90% of the data in the A and B zones of the Clark Error Grid analysis, with a majority of the data in the A zone (Claxk et al., supra).
[00037] A Clarke Error Grid Analysis for the prediction of a subject's blood glucose concentration based solely on a single measurement of the subject's ISF
glucose concentration is depicted in FIG. 2. FIG. 3 is a Clarke Error Grid Analysis for the prediction of a subject's blood glucose concentration based on a series of ISF
glucose concentrations and the algorithm immediately above. Both FIG. 3 and FIG. 4 were obtained using the data set described above.
[00038] Referring to FIGS. 2 and 3, it is evident that use of a series of ISF
glucose concentrations and the algorithm above beneficially increased the percentage of predicted blood glucose concentrations in zone A to 88.2%a compared to 79.5%
when a sole ISF glucose concentration was employed to predict blood glucose concentration.
"..,... a ,n: :.,.. , ;,~;.~",, n~,~.~,,>x....;ammxr.«-.~mm"~.-.~..,:~:vmurar.s,~~a~nn~a~.,~sn:z~s..$:w"~..,rmr~..~:.,~ma~... ~-... n ,m..N
~>«n..rw.,..r.. -,...,.. .",."~.........,...,-<-.....___-,~._~....-...,.-....~.__~"....,. _.._.._...._..."..~..__,...,.,...~.
Example 2: Predictive Algorithm for a Glucose Analytical Device ITtilizing ISF
~', rated, ma"rate,~p; Significant Interaction Terms [00039] Employing the same data set as in Example 1 above, algorithms employing ISF;k, rate, ma"ratemP, significant interaction terms were developed as described below. The algorithms employed smoothing variables of the general form ma~rate~, (discussed above) using two to four point moving averages. Weighting variables were also included to improve the algorithms' ability to accurately predict blood glucose concentration from the series of ISF glucose concentrations. The weighting algorithm used was as follows (in SAS~ code):
weight4=ISF* *-4;
tiewweight=200;
if malratel < 0 and ma3ratel <= 0 then do;
if ratel <= 0 then newweight=weight4*(-1*ratel+1)**2;
if ratel > 0 then newweight=(weight4*(abs(malratel}+1)**2)/(1+ratel);
end;
if malratel > 0 and ma3ratel >= 0 then do;
if ratel >= 0 then newweight=weight4*{1*ratel+1)**2;
if ratel < 0 then newweight={weight4*(abs(malratel}+1)**2)/(1+ abs(ratel));
end;
if malratel <= 0 and ma3ratel > 0 then do;
if ratel >= 0 then newweight={weight4*(1*ratel+1)**2)!4;
if ratel < 0 then newweight=(weight4*(-1.*ratel+1)**2}/2;
end;
if ma 1 rate 1 >= 0 and ma3rate 1 < 0 then do;
if ratel > 0 then newweight=(weight4*(1*ratel+1)**2)I2;
if ratel <= 0 then newweight=(weight4*(-1*ratel+1)**2)!4;
end;
if newweight ne 200 then do;
newweight=10000000000*newweight;
end;
[00040] Separate equations were developed for increasing (rising) and decreasing (falling) ISF glucose concentration trends in order to provide analytical devices and methods of superior accuracy. For data series that indicate a decreasing (falling) ISF
.4.... . e.~. ~ _. ._~ ~ r ,~M-~~ ..~b._ .4.. ~.k~.-~~~.~ A~~~ ~~ya.~..~~ ,~~, .~x~a~... ~a~.. ~.__~ ~.~....~~,.~~__..~~..a~.~
glucose concentrations, the following model was obtained by least squares regression analysis using SAS~ version 8.02 and N=278 data points:
PC = 8.23malratel + 0.88ISF3 + 12.04malrate2 + 10.S4ratel + 1.71ratel*rate2 0.056ISF*ratel + 0.71(ratel)2 + 0.68(rate2)2 + 0.0014(ISF)2 - 0.001 I (TSF3)2 Fox data series that indicate an increasing (rising) ISF glucose concentration, the following model was obtained by least squares regression analysis with SAS~
version 8.02 and N=I80 data points:
PC = 4.I3ISF - l .S IISFI - I.69ISF3 - 37.06malrate2 + I3.67ma3ratel -28.3Sratel -3.S6ratel *rate2 + 0. lOISF*ratel + 0.1 SISF*rate2 + 0.47ratel *rate2*rate3 -1.13(rate3)2 - 0.0061 (ISF)2 + 0.0060(ISF~Z
[00041] FIG. 4 is a CIarke Error Grid Analysis for the precliction of a subject's blood glucose concentration based on a series of ISF glucose concentrations and the algorithms immediately above. FIG. 4 was obtained using the data set described above with respect to Example 1.
[fl0042] Another measure of device accuracy is the mean absolute % error (MPE(%)) which is determined from the mean of individual % error (PE) given by the following function:
PE = (PG~-BG~)BGc where:
BG~ = the reference glucose measurement at time t, and PGt = the predicted glucose measurement at time t.
The MPE(%) results for the use of no algorithm (i.e., simply predicting that subject's blood glucose concentration is equal to a subject's ISF glucose concentration) and the two algorithms described inunediately above are depicted in Table 1 along with selected results from FIG. 4.
Yet another measure of device accuracy is average percent bias (Avg Bias(%)).
Bias (%) is determined by the following equation:
Bias(%) = [(PG~-BG~)/BG~]* 100 Avg Bias(°k) = [sum of all Bias(%)]/total number of measurements Effective measurements should have an Avg Bias(%) of about 10% or less. Table shows that the Avg Bias (%1 criterion is benef5cially decreased by use of the predictive algorithm.
[00043) The correlation between calculated and measured blood glucose values was also assessed. The correlation coefficient values (R) also presented in Table 1. below.
Effective measurements should have R values of greater than about 0.85. As can be seen, the predictive algorithm of the present invention provides for improved correlation between actual and predicted values.
Algorithms N MPE(%) A (%) B (%) Other (%) R Avg Bias (°k) None 458 14 79.9 17.7 2.4 0.94 6.76 Example 2 458 10 88.4 9.6 2.0 0.96 0.46 [00044] FIG. S is a flow chart illustrating a sequence of steps in a process 500 for predicting a subject's blood analyte concentration based on the subjects' ISF
analyte concentration according to an exemplary embodiment of the present invention.
Process 500 includes extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner, as set forth in step 510, and sequentially determining an ISF analyte concentration in each of the plurality of ISF samples, as set forth in step 520. The result of step 520 is the creation of a series of ISF analyte concentrations.
[00045] Steps S 10 and 520 of process 500 can be accomglished using any suitable techniques including those described above with respect to sampling modules and analysis modules of analytical devices according to the present invention.
[00046] Next, the subject's blood analyte concentration is predicted based on the series of ISF analyte concentrations by performing at least one algorithm of the forms) described above with respect to analytical devices according to the present invention, as set forth in step 530.
[00047] While preferred embodiments of the present invention have been shown and described herein, it will be obvious to those skilled in the art that such embodiments are provided by way of example only. Numerous variations, changes, and substitutions will now occur to hose skilled in the art without departing from the invention.
[00048] It should be understood that various alternatives to the embodiments of the invention described herein may be employed in practicing the invention. It is intended that the following claims define the scope of the invention and that methods and structures within the scope of these claims and their equivalents be covered thereby.
Provisional Patent Application No. 601476,733 (filed June 6, 2003) and in International Application PCT/GBOl/05634 (published as International Publication No. WO 02/49507 A1 on June 27, 2002). In addition, the analysis module can include any suitable analyte sensor known to those of skill in the art including, but not limited to, glucose analyte sensors based on photometric or electrochemical analytical techniques.
[oOOlsJ Prediction module I30 is canfigured for storing the series of ISF
analyte concentrations created by the analysis module and predicting the subject's whole blood analyte concentration by performing at least one algorithm of the following general form {referred to as Eqn 1):
PC = f (ISF;k, ratel, significant interaction terms) (Eqn 1) where:
PC is a predicted subject's whole blood analyte concentration;
i is an integer with predetermined values selected from the values of, for example, 0,1, 2, 3, 4 and 5;
j is an integer with predetermined values selected from the values of, for example, 1, 2, 3, 4 and 5;
k is an integers) with predetermined values selected from the values of, for example, l and 2;
ISF; is a measured ISF analyte concentration in the series of ISF analyze concentrations, with the subscript (i) indicating which ISF value is being refen~ed to, i.e., i = 0 indicates the most recently measured ISF analyte, i =1 indicates one value back in the series of ISF analyte concentrations, 2 = two values back in the series of ISF
analyte concentrations, etc.;
rated is the rate of change between immediately adjacent ISF~analyte concentrations in the series of ISF analyte concentrations (calculated as the difference between the immediately adjacent ISF concentrations divided by the time difference between when the immediately adjacent ISF concentrations were measured by the analysis module) with the subscript (j) referring to which immediately adjacent ISF
values are used to calculate the rate, i.e., when j =1 indicates the rate between current ISF value and the previous ISF value, j = 2 is indicative of the rate between the ISF
values one previous and two previous relative to the current ISF value, etc.;
and significant interaction terms = interaction ternns involving at least two of ISF;k, rate.
[00419] The mathematical form of the function (~ employed in Eqn 1 can be any suitable mathematical form that accounts for physiological lag between ISF
analyte concentration and blood analyte concentration as well as any bias effect between ISF
and blood analyte concentrations. However, it has been determined that such a relationship is suitably accurate when measured ISF analyte concentrations (i.e., ISF;~'), rates {ratei) and interaction terms are included. The use of rates is particularly beneficial in providing an accurate algorithm, and hence an accurate analytical device, since the time interval between the ISF analyte concentrations can be non-uniform (e.g., the time interval could vary between 5 minutes and 15 minutes).
[00020] The form of the function (~ can determined by, for example, a least squares regression analysis of a statistically relevant number of ISF analyte concentrations and associated blood anaIyte concentrations. Those skilled in the art will appreciate that any number of mathematical methods (e.g., mathematical modeling methods) can be used to analyze such data and arrive at a suitable function ~. For example, linear and polynomial regression analysis, time series analysis, or neural networks can be used.
In the circumstance that the analyte is glucose, ISF glucose concentrations and blood glucose concentrations can be determined from ISF and blood samples extracted from diabetic subjects that have ingested glucose.
[00021] If desired, a suitable algorithm can be obtained using a mathematical modeling method that includes weighting factors to provide for greater accuracy at lower analyte concentrations values, to account for curvature in the response, and/or to account for noise in the modeled data. Weighting of input observations can also be similarly beneficial in such mathematical modeling methods.
[00022] The determination of suitable weighting factors can be, for example, an iterative process in which a weighting factors) is applied in a model, the weighting factor's effect on model results observed, and the weighting factors) adjusted based on model error reduction. The choice of weighting factors in the mathematical modeling method can also be determined, for example, by the relative importance of data ranges andlor trending direction. For example, when glucose is the analyte of interest, greater accuracy for the low end of the physiological glucose concentration range may be deemed important, and thus a weighting factor that enhances the importance of lower glucose concentrations can be employed. Such an enhancement can be accomplished, far example, by multiplying observed glucose concentrations by the inverse of the observed value raised to a predetermined power. Similarly, weighting factors can be determined which will enhance the importance of certain events or trends in observed values, such as a magnitude of the gradient of an observed rate and/or a change in direction. Furthermore, prospective weighting factors can also be arbitrarily chosen with suitable weighting factors chosen from the prospective weighting factors based on their effect on model error reduction.
[00023] Prediction module 130 can take any suitable forn~ known to one skilled in the art including, but not limited to, the remote control modules described in co-pending U.S. Provisional Patent Application No. 60/476,733 (filed June 6, 2003), which is hereby fully incorporated herein by reference.
[00024] As an alternative to the use of Eqn 1 above, prediction module 130 can be configured for storing the series of ISF analyte concentrations created by the analysis module and predicting the subject's whole blood analyte concentration by performing at least one algorithm of the following general form (referred to as Eqn 2):
PC = f (ISF;k, ratei, ma~,rate~,~, significant interaction terms) (Eqn 2) where:
n and m are integers with predetermined values selected from the values of, for example, 1, 2 and 3;
p is an integers) with predetermined values selected from the values of, for example, l and 2; and ma"ratem is the moving average rate between immediately adjacent averages of groupings of ISF values, with the subscript (n) referring to the number of ISF
values included in the moving average and the subscript (m) referring to the time position of the adjacent average values relative to the current values as follows (with a further explanation in the next paragraph):
n+1 = the number of points used in the moving average rate;
m-1 = first point included in the moving average. If m-1= 0, then the current ISF value is used as the first point in the moving average calculation.
n + m always adds up to the number of points back (or removed from the current ISF value) that will be needed for calculating the moving average calculation.
n and m are integers with predetermined values selected from the values of, for example, 1, 2 and 3; and significant interaction terms = interaction terms involving at least two of ISF;k, ratej, and ma"ratemp. , (00025] The following example illustrates the concept of the moving average rate (ma"ratero) employed in Eqn 2 above. For the exemplary moving average rate malratel, the moving average rate is a 2-point moving average rate {since m+n =1+1 = 2) with the moving average rate calculated iietween the average of the grouping that includes the most recent iSF concentration and the ISF concentration one point back and the average of the grouping that includes the ISF concentrations one and two points prior to the most recent ISF concentration.
(00026) Eqn 2 includes moving average rates (i.e., ma"rate"~ to smooth the data (i:e., the series of ISF analyte concentrations and/or rates) with respect to both rate and the trending direction of an analyte concentration, thereby removing noise from the data and increasing fhe analytical device's accuracy. Although significant (i.e., major) changes in adjacent ISF values can be regarded as important in terms of algorithm accuracy, significant changes can also be due to noise that can adversely effect an algorithm's accuracy. The moving average rate, which is the rate of change between the means of adjacent (and overlapping) groupings of ISF values, dampens noise caused by outlier values that do not represent a true trend in the data.
(00027] Examples of suitable algorithms and the techniques used to derive the algorithms are includes in the examples below.
Example 1: Predictive Algorithm for a Glucose Analytical Device Utilizing ISF
k, ratel, and Significant Interaction Terms A data set (i.e., a series of ISF glucose concentrations) was generated using an experimental ISF sampling and analysis modules. The ISF sampling module and analysis module employed to generate the data set were configured to extract an ISF
sample from a subject's dermal layer of skin (i.e., dermis), for example from a subject's forearm, and to measure the glucose concentration in the ISF sample.
The ISF sampling module and analysis module were an integrated unit comprising a one-piece sampling module and a modified OneTouch~ Ultra glucose meter with test strip.
The sampling utilized a 30-gauge cannula and a penetraeion depth of about 1 to 2 mm.
It should be noted that an ISF sample collected from the dermis is considered to have a beneficially reduced physiological lag in comparison to an ISF sample collected from a subcutaneous layer due to the dermis being closer to vascular capillary beds than the subcutaneous layer. ' j00028] The ISF sampling module extracted an approximately 1 p.I. ISF sample from a subject's dermis via the cannula and deposited the ISF sample automatically and directly into a measurement zone of the test strip. After a brief electrochemistry development period, the meter displayed the ISF glucose concentration.
[00029] Prior to the ISF samples being extracted, 2 to 4 pounds of pressure was applied to a subject's dermis for 34 seconds, followed by a 5 minute waiting period to allow blood to perfuse (flow into) the sampling area from which the ISF would be extracted.
This elevated blood-flow in the sampling area has the desirable effect of mitigating the physiological lag between blood glucose concentration and ISF glucose concentration, simply because the sampling area is better perfused with flowing blood.
[00030] Finger stick blood glucose measurements in mg/dL (i.e., blood glucose concentrations) were taken from 20 subjects, followed by measurements of glucose in forearm interstitial fluid (i.e., ISF glucose concentrations) as described above. The finger stick blood measurements were taken approximately 15 minutes apart and each was followed approximately 5 minutes later by an ISF sample extraction and ISF
glucose concentration measurement.
[00031] Approximately thirty (3Q) pairs of observations (i.e., pairs of blood glucose concentration measurements and ISF glucose concentration measurements) were obtained fox each of the 20 subjects. The observations were collected over the course of one day for each~subject, in whom a change in glucose concentration was induced through the ingestion of 75g of glucose. The blood glucose concentration for each to observation represents a finger stick draw occurring approximately five minutes prior to the ISF draw.. Blood glucose concentration at the time of the ISF sample extraction was, therefore, linearly interpolated, with the linearly interpolated value used as a response variable in developing the algorithm below. The final ISF glucose concentration for each subject was excluded during the development of the algorithm due to the inability to accurately interpolate a blood glucose concentration.
[40032] An algorithm of the form identified above as Eqn 2 was developed from the data set using multiple linear regression. The algorithm thus developed weighted lower ISF analyte concentrations more heavily, primarily due to the relative importance of accurately predicting glucose at lower concentrations. The weight used was ISF-0. In the absence of such weighting, higher ISF glucose concentration values produced undesirable variability in the residuals.
[00033] The parameters, estimates, errors, t-values and Pr values for the model were as follows:
Parameter Estimate Error t ValuePr > Its ISFO 0.964114574 0.00900642107.05 <.0001 rate2 3.564454310 0.791431254.50 <.0001 rate2*ratel 1.032526146 0.276843433.73 0.0002 rate3 2.115098820 0.572521873.69 0.0002 ratel*rate3 0.728563905 0.325075672.24 0.0255 rate2*rate3 0.993732089 0.362936362.74 0.0064 rate4 2.620714810 0.480338955.46 <.0001 rate2*rate4 1.149236162 0.380759903.02 0.0027 rate2*rate3*rate4 0.140279472.99 0.0029 0.419884620 rates 1.704279771 0.409084594.17 <.0001 R-squared=.98 [00034] The algorithm, following estimators therefore, has form are the when the employed with significant two decimal places:
. PC = 0.96ISFo +3.56rate2 + 1.03 (rate2*rate~ ) + 2.I lrate3 +0.72 (rates *rate3) .
+0.99rate 4 +1.14(rate2*rate4) + 0.42(rate2*rate2*rate4) + 1.70rate5.
[00035] One skilled in the art will recognize,that the above equation is of the form of Eqn. 1 above with:
~..._ ., .., .<nxe~.._c r. -~av ~"-.~rs_a."-.ymrao~~c>°a ar;a-R~.~,,"".
s,.~~~E. ~~,:~,~h =n.~.n°.~xamP=. ~-n m~ .~.~e.....r......~~,~...._.__.
..,....._..._.-......___.,....,.....w....._...N~"._..",.,._._.-__",..._.~,..,.._..._,.....__.~.
I=0 k=1 j=2,3,4and5 and interaction terms = rate2*ratel,ratet*rate3,rate2*rateq, arid rate2*rate2*rate4.
[00036] A Clarke Error Grid analysis can be employed to determuine the accuracy and suitability of an algorithm for the prediction of a subject's blood glucose concentration. The error grid of such an analysis categorizes an analytical'device's response against a reference value into one of five (5) clinical accuracy zones (i.e., zones A-E). Zone A indicates clinically accurate results, zone B indicates results that are not clinically accurate but pose minimal risk to patient health, and zones C through E indicate clinically inaccurate results that pose increasing potential risk to patient health (see Clarke, William L. et al., Evaluating Clinical Accuracy of Systems for Sedf Monitoring of Blood Glucose, Diabetes Care, Vol: 10 hTo. 5, 622-628 [1987]).
An effective~and accurate blood glucose monitoring device should have greater than approximately 85-90% of the data in the A and B zones of the Clark Error Grid analysis, with a majority of the data in the A zone (Claxk et al., supra).
[00037] A Clarke Error Grid Analysis for the prediction of a subject's blood glucose concentration based solely on a single measurement of the subject's ISF
glucose concentration is depicted in FIG. 2. FIG. 3 is a Clarke Error Grid Analysis for the prediction of a subject's blood glucose concentration based on a series of ISF
glucose concentrations and the algorithm immediately above. Both FIG. 3 and FIG. 4 were obtained using the data set described above.
[00038] Referring to FIGS. 2 and 3, it is evident that use of a series of ISF
glucose concentrations and the algorithm above beneficially increased the percentage of predicted blood glucose concentrations in zone A to 88.2%a compared to 79.5%
when a sole ISF glucose concentration was employed to predict blood glucose concentration.
"..,... a ,n: :.,.. , ;,~;.~",, n~,~.~,,>x....;ammxr.«-.~mm"~.-.~..,:~:vmurar.s,~~a~nn~a~.,~sn:z~s..$:w"~..,rmr~..~:.,~ma~... ~-... n ,m..N
~>«n..rw.,..r.. -,...,.. .",."~.........,...,-<-.....___-,~._~....-...,.-....~.__~"....,. _.._.._...._..."..~..__,...,.,...~.
Example 2: Predictive Algorithm for a Glucose Analytical Device ITtilizing ISF
~', rated, ma"rate,~p; Significant Interaction Terms [00039] Employing the same data set as in Example 1 above, algorithms employing ISF;k, rate, ma"ratemP, significant interaction terms were developed as described below. The algorithms employed smoothing variables of the general form ma~rate~, (discussed above) using two to four point moving averages. Weighting variables were also included to improve the algorithms' ability to accurately predict blood glucose concentration from the series of ISF glucose concentrations. The weighting algorithm used was as follows (in SAS~ code):
weight4=ISF* *-4;
tiewweight=200;
if malratel < 0 and ma3ratel <= 0 then do;
if ratel <= 0 then newweight=weight4*(-1*ratel+1)**2;
if ratel > 0 then newweight=(weight4*(abs(malratel}+1)**2)/(1+ratel);
end;
if malratel > 0 and ma3ratel >= 0 then do;
if ratel >= 0 then newweight=weight4*{1*ratel+1)**2;
if ratel < 0 then newweight={weight4*(abs(malratel}+1)**2)/(1+ abs(ratel));
end;
if malratel <= 0 and ma3ratel > 0 then do;
if ratel >= 0 then newweight={weight4*(1*ratel+1)**2)!4;
if ratel < 0 then newweight=(weight4*(-1.*ratel+1)**2}/2;
end;
if ma 1 rate 1 >= 0 and ma3rate 1 < 0 then do;
if ratel > 0 then newweight=(weight4*(1*ratel+1)**2)I2;
if ratel <= 0 then newweight=(weight4*(-1*ratel+1)**2)!4;
end;
if newweight ne 200 then do;
newweight=10000000000*newweight;
end;
[00040] Separate equations were developed for increasing (rising) and decreasing (falling) ISF glucose concentration trends in order to provide analytical devices and methods of superior accuracy. For data series that indicate a decreasing (falling) ISF
.4.... . e.~. ~ _. ._~ ~ r ,~M-~~ ..~b._ .4.. ~.k~.-~~~.~ A~~~ ~~ya.~..~~ ,~~, .~x~a~... ~a~.. ~.__~ ~.~....~~,.~~__..~~..a~.~
glucose concentrations, the following model was obtained by least squares regression analysis using SAS~ version 8.02 and N=278 data points:
PC = 8.23malratel + 0.88ISF3 + 12.04malrate2 + 10.S4ratel + 1.71ratel*rate2 0.056ISF*ratel + 0.71(ratel)2 + 0.68(rate2)2 + 0.0014(ISF)2 - 0.001 I (TSF3)2 Fox data series that indicate an increasing (rising) ISF glucose concentration, the following model was obtained by least squares regression analysis with SAS~
version 8.02 and N=I80 data points:
PC = 4.I3ISF - l .S IISFI - I.69ISF3 - 37.06malrate2 + I3.67ma3ratel -28.3Sratel -3.S6ratel *rate2 + 0. lOISF*ratel + 0.1 SISF*rate2 + 0.47ratel *rate2*rate3 -1.13(rate3)2 - 0.0061 (ISF)2 + 0.0060(ISF~Z
[00041] FIG. 4 is a CIarke Error Grid Analysis for the precliction of a subject's blood glucose concentration based on a series of ISF glucose concentrations and the algorithms immediately above. FIG. 4 was obtained using the data set described above with respect to Example 1.
[fl0042] Another measure of device accuracy is the mean absolute % error (MPE(%)) which is determined from the mean of individual % error (PE) given by the following function:
PE = (PG~-BG~)BGc where:
BG~ = the reference glucose measurement at time t, and PGt = the predicted glucose measurement at time t.
The MPE(%) results for the use of no algorithm (i.e., simply predicting that subject's blood glucose concentration is equal to a subject's ISF glucose concentration) and the two algorithms described inunediately above are depicted in Table 1 along with selected results from FIG. 4.
Yet another measure of device accuracy is average percent bias (Avg Bias(%)).
Bias (%) is determined by the following equation:
Bias(%) = [(PG~-BG~)/BG~]* 100 Avg Bias(°k) = [sum of all Bias(%)]/total number of measurements Effective measurements should have an Avg Bias(%) of about 10% or less. Table shows that the Avg Bias (%1 criterion is benef5cially decreased by use of the predictive algorithm.
[00043) The correlation between calculated and measured blood glucose values was also assessed. The correlation coefficient values (R) also presented in Table 1. below.
Effective measurements should have R values of greater than about 0.85. As can be seen, the predictive algorithm of the present invention provides for improved correlation between actual and predicted values.
Algorithms N MPE(%) A (%) B (%) Other (%) R Avg Bias (°k) None 458 14 79.9 17.7 2.4 0.94 6.76 Example 2 458 10 88.4 9.6 2.0 0.96 0.46 [00044] FIG. S is a flow chart illustrating a sequence of steps in a process 500 for predicting a subject's blood analyte concentration based on the subjects' ISF
analyte concentration according to an exemplary embodiment of the present invention.
Process 500 includes extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner, as set forth in step 510, and sequentially determining an ISF analyte concentration in each of the plurality of ISF samples, as set forth in step 520. The result of step 520 is the creation of a series of ISF analyte concentrations.
[00045] Steps S 10 and 520 of process 500 can be accomglished using any suitable techniques including those described above with respect to sampling modules and analysis modules of analytical devices according to the present invention.
[00046] Next, the subject's blood analyte concentration is predicted based on the series of ISF analyte concentrations by performing at least one algorithm of the forms) described above with respect to analytical devices according to the present invention, as set forth in step 530.
[00047] While preferred embodiments of the present invention have been shown and described herein, it will be obvious to those skilled in the art that such embodiments are provided by way of example only. Numerous variations, changes, and substitutions will now occur to hose skilled in the art without departing from the invention.
[00048] It should be understood that various alternatives to the embodiments of the invention described herein may be employed in practicing the invention. It is intended that the following claims define the scope of the invention and that methods and structures within the scope of these claims and their equivalents be covered thereby.
Claims (16)
1. An analytical device for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid analyte concentration, the analytical device comprising:
an interstitial fluid sampling module for extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner;
an analysis module for sequentially determining an ISF analyte concentration in each of the plurality of ISF samples, thereby creating a series of ISF
analyte concentrations; and a prediction module for storing the series of ISF analyte concentrations and predicting the subject's whole blood analyte concentration based on the series of ISF
analyte concentrations by performing at least one algorithm of the following general form:
PC = f(ISF i k, rate j, significant interaction terms) where:
PC = the predicted subject's whole blood analyte concentration;
i is an integer with predetermined values selected from the values of 0, 1, 2, 3, 4 and 5;
j is an integer with predetermined values selected from the values of 1, 2, 3, 4 and 5;
k is an integer with predetermined values selected from the values of 1 and 2;
ISF i is a measured ISF analyte concentration in the series of ISF analyte concentrations;
rated is a rate of change between immediately adjacent ISF analyte concentrations in the series of ISF analyte concentrations; and significant interaction terms = statistically significant interaction terms involving terms selected from the group consisting of ISF i k and rate j.
an interstitial fluid sampling module for extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner;
an analysis module for sequentially determining an ISF analyte concentration in each of the plurality of ISF samples, thereby creating a series of ISF
analyte concentrations; and a prediction module for storing the series of ISF analyte concentrations and predicting the subject's whole blood analyte concentration based on the series of ISF
analyte concentrations by performing at least one algorithm of the following general form:
PC = f(ISF i k, rate j, significant interaction terms) where:
PC = the predicted subject's whole blood analyte concentration;
i is an integer with predetermined values selected from the values of 0, 1, 2, 3, 4 and 5;
j is an integer with predetermined values selected from the values of 1, 2, 3, 4 and 5;
k is an integer with predetermined values selected from the values of 1 and 2;
ISF i is a measured ISF analyte concentration in the series of ISF analyte concentrations;
rated is a rate of change between immediately adjacent ISF analyte concentrations in the series of ISF analyte concentrations; and significant interaction terms = statistically significant interaction terms involving terms selected from the group consisting of ISF i k and rate j.
2. The analytical device of claim 1, wherein i = 0, k = 1, j = 2, 3, 4 and 5 and interaction terms = rate2*rate1, rate1*rate3, rate2*rate4, and rate2*rate2*rate4.
3. The analytical device of claim 1, wherein the analyte is glucose.
4. The analytical device of claim 1, wherein the predicting the subject's whole blood analyte concentration is based on the series of ISF analyte concentrations by performing at least one algorithm of the following general form:
PC = f(ISF i k, rate j, ma n rate m p, significant interaction terms) where:
p is an integer with predetermined values selected from the values of 1 and 2;
n and m are integers with predetermined values selected from the values of 1, 2 and 3;
ma n rate m is the moving average rate between adjacent averages of groupings of ISF values; and significant interaction terms = statistically significant interaction terms involving terms selected from the group consisting of ISF i k, rate j, and ma n rate m p.
PC = f(ISF i k, rate j, ma n rate m p, significant interaction terms) where:
p is an integer with predetermined values selected from the values of 1 and 2;
n and m are integers with predetermined values selected from the values of 1, 2 and 3;
ma n rate m is the moving average rate between adjacent averages of groupings of ISF values; and significant interaction terms = statistically significant interaction terms involving terms selected from the group consisting of ISF i k, rate j, and ma n rate m p.
5. The analytical device of claim 4, wherein the prediction module predicts the subject's whole blood analyte concentration by determining whether the series of ISF analyte concentrations is indicative of a rising ISF analyte concentration or a falling ISF analyte concentration, selecting an algorithm based on the determination and performing the selected algorithm.
6. The analytical device of claim 5, wherein the prediction module predicts the subject's whole blood analyte concentration by determining whether the series of ISF analyte concentrations is indicative of a rising ISF analyte concentration or a falling ISF analyte concentration based on an ma n rate m.
7. The analytical device of claim 6, wherein the algorithm employed for a falling ISF analyte concentration is:
PC = 8.23ma1rate1 + 0.88ISF3 + 12.04ma1rate2 + 10.54rate1 +
1.71rate1*rate2 - 0.056ISF*rate1 + 0.71(rate1)2 + 0.68(rate2)2 +
0.0014(ISF)2_sq -0.0011 (ISF3)2.
PC = 8.23ma1rate1 + 0.88ISF3 + 12.04ma1rate2 + 10.54rate1 +
1.71rate1*rate2 - 0.056ISF*rate1 + 0.71(rate1)2 + 0.68(rate2)2 +
0.0014(ISF)2_sq -0.0011 (ISF3)2.
8. The analytical device of claim 6, wherein the algorithm employed for a rising ISF analyte concentration is:
PC = 4.13ISF - 1.51ISF1 - 1.69ISF3 - 37.06ma1rate2 + 13.67ma3rate1 -28.35rate1 -3.56rate1*rate2 + 0.10ISF*rate1 + 0.15ISF*rate2 + 0.47rate1*rate2*rate3 -1.13(rate3)2 - 0.0061(ISF)2 + 0.0060(ISF2)2.
PC = 4.13ISF - 1.51ISF1 - 1.69ISF3 - 37.06ma1rate2 + 13.67ma3rate1 -28.35rate1 -3.56rate1*rate2 + 0.10ISF*rate1 + 0.15ISF*rate2 + 0.47rate1*rate2*rate3 -1.13(rate3)2 - 0.0061(ISF)2 + 0.0060(ISF2)2.
9. The analytical device of claim 4, wherein the analyte is glucose.
10. The analytical device of claim 1, wherein the series of ISF analyte concentrations includes five ISF analyte concentrations.
11. The analytical device of claim 1, wherein the sampling module extracts the plurality of ISF samples at a time interval in the range of five to fifteen minutes.
12. A method for predicting a subject's whole blood analyte concentration based on the subject's interstitial fluid analyte concentration, the method comprising:
extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner;
sequentially determining an ISF analyte concentration in each of the plurality of ISF samples, thereby creating a series of ISF analyte concentrations; and predicting the subject's blood analyze concentration based on the series of ISF analyte concentrations by performing at least one algorithm of the following form:
PC = f(ISF i k, rate j, significant interaction terms) where:
PC = the predicted subject's whole blood analyte concentration;
i is an integer with predetermined values selected from the values of 0, 1, 2, 3, 4 and 5;
j is an integer with predetermined values selected from the values of 1, 2, 3, 4 and 5;
k is an integer with predetermined values selected from the values of 1 and 2;
ISF i is a measured ISF analyte concentration in the series of ISF analyte concentrations;
rated is the rate of change between adjacent ISF analyte concentrations in the series of ISF analyte concentrations; and significant interaction terms = statistically significant interaction terms involving terms selected from the group consisting of ISF i k and rate j.
extracting a plurality of interstitial fluid (ISF) samples from a subject in a sequential manner;
sequentially determining an ISF analyte concentration in each of the plurality of ISF samples, thereby creating a series of ISF analyte concentrations; and predicting the subject's blood analyze concentration based on the series of ISF analyte concentrations by performing at least one algorithm of the following form:
PC = f(ISF i k, rate j, significant interaction terms) where:
PC = the predicted subject's whole blood analyte concentration;
i is an integer with predetermined values selected from the values of 0, 1, 2, 3, 4 and 5;
j is an integer with predetermined values selected from the values of 1, 2, 3, 4 and 5;
k is an integer with predetermined values selected from the values of 1 and 2;
ISF i is a measured ISF analyte concentration in the series of ISF analyte concentrations;
rated is the rate of change between adjacent ISF analyte concentrations in the series of ISF analyte concentrations; and significant interaction terms = statistically significant interaction terms involving terms selected from the group consisting of ISF i k and rate j.
13. The method of claim 12, wherein the predicting set employs an algorithm of the form:
PC = f(ISF i k, rate j, ma n rate m p, significant interaction terms) where:
p is an integer with predetermined values selected from the values of 1 and 2;
m and n are integers with predetermined values selected from the values of 1, 2 and 3;
ma n rate m is the moving average rate between adjacent averages of groupings of ISF values; and significant interaction terms = statistically significant interaction terms involving terms selected from the group consisting of ISF i k, rate j, and ma n rate m p.
PC = f(ISF i k, rate j, ma n rate m p, significant interaction terms) where:
p is an integer with predetermined values selected from the values of 1 and 2;
m and n are integers with predetermined values selected from the values of 1, 2 and 3;
ma n rate m is the moving average rate between adjacent averages of groupings of ISF values; and significant interaction terms = statistically significant interaction terms involving terms selected from the group consisting of ISF i k, rate j, and ma n rate m p.
14. The method of claim 13, wherein the predicting step predicts the subject's whole blood analyte concentration by determining whether the series of ISF
analyte concentrations is indicative of a rising ISF analyte concentration or a falling ISF analyte concentration, selecting the algorithm based on the determination and performing the selected algorithm.
analyte concentrations is indicative of a rising ISF analyte concentration or a falling ISF analyte concentration, selecting the algorithm based on the determination and performing the selected algorithm.
15. The method of claim 13, wherein the predicting step predicts the subject's whole blood analyte concentration by determining whether the series of ISF
analyte concentrations is indicative of a rising ISF analyte concentration or a falling ISF analyte concentration based on an ma n rate m.
analyte concentrations is indicative of a rising ISF analyte concentration or a falling ISF analyte concentration based on an ma n rate m.
16. The method of claim 12, wherein the extracting step extracts a plurality of ISF samples from a subject's dermis.
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/653,023 US7258673B2 (en) | 2003-06-06 | 2003-08-28 | Devices, systems and methods for extracting bodily fluid and monitoring an analyte therein |
US10/653,023 | 2003-08-28 | ||
US10/652,464 US20040253736A1 (en) | 2003-06-06 | 2003-08-28 | Analytical device with prediction module and related methods |
US10/652,464 | 2003-08-28 |
Publications (1)
Publication Number | Publication Date |
---|---|
CA2479107A1 true CA2479107A1 (en) | 2005-02-28 |
Family
ID=34108168
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002479107A Abandoned CA2479107A1 (en) | 2003-08-28 | 2004-08-26 | Analytical device with prediction module and related methods |
Country Status (12)
Country | Link |
---|---|
US (1) | US20040253736A1 (en) |
EP (1) | EP1510171A1 (en) |
JP (1) | JP2005077412A (en) |
KR (1) | KR20050021955A (en) |
CN (1) | CN1644163A (en) |
AU (1) | AU2004205264A1 (en) |
CA (1) | CA2479107A1 (en) |
IL (1) | IL163729A0 (en) |
NO (1) | NO20043583L (en) |
RU (1) | RU2004126182A (en) |
SG (2) | SG131928A1 (en) |
TW (1) | TW200513647A (en) |
Families Citing this family (63)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8465425B2 (en) | 1998-04-30 | 2013-06-18 | Abbott Diabetes Care Inc. | Analyte monitoring device and methods of use |
US8688188B2 (en) | 1998-04-30 | 2014-04-01 | 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 |
US8480580B2 (en) | 1998-04-30 | 2013-07-09 | Abbott Diabetes Care Inc. | Analyte monitoring device and methods of use |
US8346337B2 (en) | 1998-04-30 | 2013-01-01 | Abbott Diabetes Care Inc. | Analyte monitoring device and methods of use |
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 |
US6175752B1 (en) | 1998-04-30 | 2001-01-16 | Therasense, 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 |
US6560471B1 (en) | 2001-01-02 | 2003-05-06 | Therasense, Inc. | Analyte monitoring device and methods of use |
EP1397068A2 (en) | 2001-04-02 | 2004-03-17 | Therasense, Inc. | Blood glucose tracking apparatus and methods |
US7993108B2 (en) | 2002-10-09 | 2011-08-09 | Abbott Diabetes Care Inc. | Variable volume, shape memory actuated insulin dispensing pump |
DE60336834D1 (en) | 2002-10-09 | 2011-06-01 | Abbott Diabetes Care Inc | FUEL FEEDING DEVICE, SYSTEM AND METHOD |
US7727181B2 (en) | 2002-10-09 | 2010-06-01 | Abbott Diabetes Care Inc. | Fluid delivery device with autocalibration |
AU2003303597A1 (en) | 2002-12-31 | 2004-07-29 | Therasense, Inc. | Continuous glucose monitoring system and methods of use |
US7679407B2 (en) | 2003-04-28 | 2010-03-16 | Abbott Diabetes Care Inc. | Method and apparatus for providing peak detection circuitry for data communication systems |
US8066639B2 (en) | 2003-06-10 | 2011-11-29 | Abbott Diabetes Care Inc. | Glucose measuring device for use in personal area network |
EP1718198A4 (en) | 2004-02-17 | 2008-06-04 | Therasense Inc | Method and system for providing data communication in continuous glucose monitoring and management system |
US20060036187A1 (en) | 2004-06-30 | 2006-02-16 | Hester Vos | Devices, systems and methods for extracting bodily fluid and monitoring an analyte therein |
CA2601441A1 (en) | 2005-03-21 | 2006-09-28 | Abbott Diabetes Care Inc. | Method and system for providing integrated medication infusion and analyte monitoring system |
US8112240B2 (en) | 2005-04-29 | 2012-02-07 | Abbott Diabetes Care Inc. | Method and apparatus for providing leak detection in data monitoring and management systems |
US7768408B2 (en) | 2005-05-17 | 2010-08-03 | Abbott Diabetes Care Inc. | Method and system for providing data management in data monitoring system |
US7620437B2 (en) | 2005-06-03 | 2009-11-17 | Abbott Diabetes Care Inc. | Method and apparatus for providing rechargeable power in data monitoring and management systems |
US20070038147A1 (en) | 2005-08-11 | 2007-02-15 | Joel Mechelke | Method for extracting interstitial fluid |
US20070060844A1 (en) * | 2005-08-29 | 2007-03-15 | Manuel Alvarez-Icaza | Applied pressure sensing cap for a lancing device |
US7756561B2 (en) | 2005-09-30 | 2010-07-13 | Abbott Diabetes Care Inc. | Method and apparatus for providing rechargeable power in data monitoring and management systems |
US7583190B2 (en) | 2005-10-31 | 2009-09-01 | Abbott Diabetes Care Inc. | Method and apparatus for providing data communication in data monitoring and management systems |
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 |
US8344966B2 (en) | 2006-01-31 | 2013-01-01 | Abbott Diabetes Care Inc. | Method and system for providing a fault tolerant display unit in an electronic device |
US7620438B2 (en) | 2006-03-31 | 2009-11-17 | Abbott Diabetes Care Inc. | Method and system for powering an electronic device |
US8226891B2 (en) | 2006-03-31 | 2012-07-24 | Abbott Diabetes Care Inc. | Analyte monitoring devices and methods therefor |
US7920907B2 (en) | 2006-06-07 | 2011-04-05 | Abbott Diabetes Care Inc. | Analyte monitoring system and method |
US8579853B2 (en) | 2006-10-31 | 2013-11-12 | Abbott Diabetes Care Inc. | Infusion devices and methods |
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 |
US8461985B2 (en) | 2007-05-08 | 2013-06-11 | 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 |
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 |
GB0800981D0 (en) * | 2008-01-18 | 2008-02-27 | Plaque Attack Ltd | Catheter |
US20110144463A1 (en) * | 2008-02-27 | 2011-06-16 | Benny Pesach | Device, system and method for modular analyte monitoring |
US20110160555A1 (en) * | 2008-07-31 | 2011-06-30 | Jacques Reifman | Universal Models for Predicting Glucose Concentration in Humans |
US8103456B2 (en) | 2009-01-29 | 2012-01-24 | Abbott Diabetes Care Inc. | Method and device for early signal attenuation detection using blood glucose measurements |
US8560082B2 (en) | 2009-01-30 | 2013-10-15 | Abbott Diabetes Care Inc. | Computerized determination of insulin pump therapy parameters using real time and retrospective data processing |
US9226701B2 (en) | 2009-04-28 | 2016-01-05 | Abbott Diabetes Care Inc. | Error detection in critical repeating data in a wireless sensor system |
WO2010129375A1 (en) | 2009-04-28 | 2010-11-11 | Abbott Diabetes Care Inc. | Closed loop blood glucose control algorithm analysis |
WO2010138856A1 (en) | 2009-05-29 | 2010-12-02 | Abbott Diabetes Care Inc. | Medical device antenna systems having external antenna configurations |
EP2456351B1 (en) | 2009-07-23 | 2016-10-12 | Abbott Diabetes Care, Inc. | Real time management of data relating to physiological control of glucose levels |
WO2011026147A1 (en) | 2009-08-31 | 2011-03-03 | Abbott Diabetes Care Inc. | Analyte signal processing device and methods |
WO2011026148A1 (en) | 2009-08-31 | 2011-03-03 | Abbott Diabetes Care Inc. | Analyte monitoring system and methods for managing power and noise |
US9320461B2 (en) | 2009-09-29 | 2016-04-26 | Abbott Diabetes Care Inc. | Method and apparatus for providing notification function in analyte monitoring systems |
WO2011105178A1 (en) * | 2010-02-26 | 2011-09-01 | アークレイ株式会社 | Analysis device, analysis method, and analysis system |
US9089292B2 (en) * | 2010-03-26 | 2015-07-28 | Medtronic Minimed, Inc. | Calibration of glucose monitoring sensor and/or insulin delivery system |
JP6443802B2 (en) | 2011-11-07 | 2018-12-26 | アボット ダイアベティス ケア インコーポレイテッドAbbott Diabetes Care Inc. | Analyte monitoring apparatus and method |
GB2502287B (en) * | 2012-05-21 | 2016-11-23 | Dermal Diagnostics Ltd | Cumulative measurement of an analyte |
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 |
AU2015211258A1 (en) | 2014-01-31 | 2016-08-11 | Trustees Of Boston University | Offline glucose control based on preceding periods |
EP3957233A3 (en) | 2015-08-07 | 2022-05-11 | Trustees of Boston University | Glucose control system with automatic adaptation of glucose target |
CN109690318B (en) * | 2016-06-30 | 2023-08-29 | 拜克门寇尔特公司 | Core calibration of an analyzer |
CN107137093B (en) * | 2017-04-20 | 2019-06-07 | 浙江大学 | A kind of continuous blood sugar monitoring device comprising abnormal plasma glucose probability alarm |
EP4000075A4 (en) | 2019-07-16 | 2023-10-04 | Beta Bionics, Inc. | Blood glucose control system |
CA3146965A1 (en) | 2019-07-16 | 2021-02-21 | Beta Bionics, Inc. | Blood glucose control system |
CN114081484B (en) * | 2021-11-24 | 2024-02-27 | 上海微创生命科技有限公司 | Continuous blood glucose detection method, system and readable storage medium |
Family Cites Families (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5002054A (en) * | 1987-02-25 | 1991-03-26 | Ash Medical Systems, Inc. | Interstitial filtration and collection device and method for long-term monitoring of physiological constituents of the body |
NL8702370A (en) * | 1987-10-05 | 1989-05-01 | Groningen Science Park | METHOD AND SYSTEM FOR GLUCOSE DETERMINATION AND USEABLE MEASURING CELL ASSEMBLY. |
US5139023A (en) * | 1989-06-02 | 1992-08-18 | Theratech Inc. | Apparatus and method for noninvasive blood glucose monitoring |
US6040194A (en) * | 1989-12-14 | 2000-03-21 | Sensor Technologies, Inc. | Methods and device for detecting and quantifying substances in body fluids |
US5231975A (en) * | 1990-02-23 | 1993-08-03 | Cygnus Therapeutic Systems | Ultrasound-enhanced delivery of materials into and through the skin |
US5593852A (en) * | 1993-12-02 | 1997-01-14 | Heller; Adam | Subcutaneous glucose electrode |
US5795325A (en) * | 1991-07-16 | 1998-08-18 | Heartport, Inc. | Methods and apparatus for anchoring an occluding member |
US5165418B1 (en) * | 1992-03-02 | 1999-12-14 | Nikola I Tankovich | Blood sampling device and method using a laser |
US5956501A (en) * | 1997-01-10 | 1999-09-21 | Health Hero Network, Inc. | Disease simulation system and method |
US5582184A (en) * | 1993-10-13 | 1996-12-10 | Integ Incorporated | Interstitial fluid collection and constituent measurement |
US5458140A (en) * | 1993-11-15 | 1995-10-17 | Non-Invasive Monitoring Company (Nimco) | Enhancement of transdermal monitoring applications with ultrasound and chemical enhancers |
US6240306B1 (en) * | 1995-08-09 | 2001-05-29 | Rio Grande Medical Technologies, Inc. | Method and apparatus for non-invasive blood analyte measurement with fluid compartment equilibration |
US5879367A (en) * | 1995-09-08 | 1999-03-09 | Integ, Inc. | Enhanced interstitial fluid collection |
US5706806A (en) * | 1996-04-26 | 1998-01-13 | Bioanalytical Systems, Inc. | Linear microdialysis probe with support fiber |
DE19618597B4 (en) * | 1996-05-09 | 2005-07-21 | Institut für Diabetestechnologie Gemeinnützige Forschungs- und Entwicklungsgesellschaft mbH an der Universität Ulm | Method for determining the concentration of tissue glucose |
US5951493A (en) * | 1997-05-16 | 1999-09-14 | Mercury Diagnostics, Inc. | Methods and apparatus for expressing body fluid from an incision |
US6332871B1 (en) * | 1996-05-17 | 2001-12-25 | Amira Medical | Blood and interstitial fluid sampling device |
JP2002515786A (en) * | 1996-06-28 | 2002-05-28 | ソントラ メディカル,エル.ピー. | Ultrasound enhancement of transdermal delivery |
US6558321B1 (en) * | 1997-03-04 | 2003-05-06 | Dexcom, Inc. | Systems and methods for remote monitoring and modulation of medical devices |
ATE386934T1 (en) * | 1997-06-04 | 2008-03-15 | Sensor Technologies Inc | METHOD AND DEVICE FOR DETECTING OR QUANTIFYING COMPOUNDS CONTAINING CARBOHYDRATES |
AU8031898A (en) * | 1997-06-16 | 1999-01-04 | Elan Medical Technologies Limited | Methods of calibrating and testing a sensor for (in vivo) measurement of an analyte and devices for use in such methods |
US6706000B2 (en) * | 1997-11-21 | 2004-03-16 | Amira Medical | Methods and apparatus for expressing body fluid from an incision |
US5964718A (en) * | 1997-11-21 | 1999-10-12 | Mercury Diagnostics, Inc. | Body fluid sampling device |
US6155992A (en) * | 1997-12-02 | 2000-12-05 | Abbott Laboratories | Method and apparatus for obtaining interstitial fluid for diagnostic tests |
US6022316A (en) * | 1998-03-06 | 2000-02-08 | Spectrx, Inc. | Apparatus and method for electroporation of microporated tissue for enhancing flux rates for monitoring and delivery applications |
DE69910007T2 (en) * | 1998-05-13 | 2004-04-22 | Cygnus, Inc., Redwood City | DEVICE FOR PREDICTING PHYSIOLOGICAL MEASUREMENTS |
US6233471B1 (en) * | 1998-05-13 | 2001-05-15 | Cygnus, Inc. | Signal processing for measurement of physiological analysis |
WO2000018289A1 (en) * | 1998-09-30 | 2000-04-06 | Cygnus, Inc. | Method and device for predicting physiological values |
US6468229B1 (en) * | 1998-10-20 | 2002-10-22 | Abbott Laboratories | Apparatus and method for the collection of interstitial fluids |
CA2361062A1 (en) * | 1999-02-04 | 2000-08-10 | Integ, Inc. | Needle for body fluid tester |
US6424847B1 (en) * | 1999-02-25 | 2002-07-23 | Medtronic Minimed, Inc. | Glucose monitor calibration methods |
EP1048265A1 (en) * | 1999-04-30 | 2000-11-02 | V.Lilienfeld-Toal, Hermann, Prof. Dr. med. | Apparatus and method for detecting a substance |
DE19935165A1 (en) * | 1999-07-28 | 2001-02-01 | Roche Diagnostics Gmbh | Method and arrangement for determining the concentration of glucose in a body fluid |
KR20020035136A (en) * | 1999-09-13 | 2002-05-09 | 세야 히로미치 | Tetrafluoroethylene/ethylene copolymer and film thereof |
AU2606901A (en) * | 1999-12-28 | 2001-07-09 | Pindi Products, Inc. | Method and apparatus for non-invasive analysis of blood glucose |
US6477392B1 (en) * | 2000-07-14 | 2002-11-05 | Futrex Inc. | Calibration of near infrared quantitative measurement device using optical measurement cross-products |
WO2002024065A1 (en) * | 2000-09-22 | 2002-03-28 | Knobbe, Lim & Buckingham | Method and apparatus for real-time estimation and control of pysiological parameters |
GB0030929D0 (en) * | 2000-12-19 | 2001-01-31 | Inverness Medical Ltd | Analyte measurement |
US6549796B2 (en) * | 2001-05-25 | 2003-04-15 | Lifescan, Inc. | Monitoring analyte concentration using minimally invasive devices |
US6501976B1 (en) * | 2001-06-12 | 2002-12-31 | Lifescan, Inc. | Percutaneous biological fluid sampling and analyte measurement devices and methods |
US6702857B2 (en) * | 2001-07-27 | 2004-03-09 | Dexcom, Inc. | Membrane for use with implantable devices |
US20030108976A1 (en) * | 2001-10-09 | 2003-06-12 | Braig James R. | Method and apparatus for improving clinical accuracy of analyte measurements |
US6966880B2 (en) * | 2001-10-16 | 2005-11-22 | Agilent Technologies, Inc. | Universal diagnostic platform |
US20030212344A1 (en) * | 2002-05-09 | 2003-11-13 | Vadim Yuzhakov | Physiological sample collection devices and methods of using the same |
DK2260759T3 (en) * | 2005-06-17 | 2015-08-17 | Hoffmann La Roche | Feel and conditioning device and method for monitoring a connection, in particular glucose, in body tissue |
-
2003
- 2003-08-28 US US10/652,464 patent/US20040253736A1/en not_active Abandoned
-
2004
- 2004-08-24 SG SG200702279-1A patent/SG131928A1/en unknown
- 2004-08-24 SG SG200404739A patent/SG109580A1/en unknown
- 2004-08-25 IL IL16372904A patent/IL163729A0/en unknown
- 2004-08-26 CA CA002479107A patent/CA2479107A1/en not_active Abandoned
- 2004-08-26 AU AU2004205264A patent/AU2004205264A1/en not_active Abandoned
- 2004-08-27 TW TW093125822A patent/TW200513647A/en unknown
- 2004-08-27 RU RU2004126182/15A patent/RU2004126182A/en not_active Application Discontinuation
- 2004-08-27 NO NO20043583A patent/NO20043583L/en not_active Application Discontinuation
- 2004-08-27 CN CNA2004100851804A patent/CN1644163A/en active Pending
- 2004-08-27 KR KR1020040068181A patent/KR20050021955A/en not_active Application Discontinuation
- 2004-08-27 EP EP04255207A patent/EP1510171A1/en not_active Withdrawn
- 2004-08-27 JP JP2004249038A patent/JP2005077412A/en active Pending
Also Published As
Publication number | Publication date |
---|---|
US20040253736A1 (en) | 2004-12-16 |
SG131928A1 (en) | 2007-05-28 |
RU2004126182A (en) | 2006-02-10 |
CN1644163A (en) | 2005-07-27 |
TW200513647A (en) | 2005-04-16 |
EP1510171A1 (en) | 2005-03-02 |
IL163729A0 (en) | 2005-12-18 |
NO20043583L (en) | 2005-02-28 |
KR20050021955A (en) | 2005-03-07 |
JP2005077412A (en) | 2005-03-24 |
SG109580A1 (en) | 2005-03-30 |
AU2004205264A1 (en) | 2005-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2479107A1 (en) | Analytical device with prediction module and related methods | |
AU2022200642B2 (en) | Analysis of glucose median, variability, and hypoglycemia risk for therapy guidance | |
US8718943B2 (en) | Method and device for utilizing analyte levels to assist in the treatment of diabetes | |
PL198154B1 (en) | Hba1c check method basing on patient glucose concentration in blood, hypoglycemia evaluation method and diabetes mellitus patient check instrument | |
US20020133063A1 (en) | Use of targeted glycemic profiles in the calibration of a noninvasive blood glucose monitor | |
EP2051620A2 (en) | Method and device for monitoring a physiological parameter | |
WO2000015102A1 (en) | Attribute compensation for analyte detection and/or continuous monitoring | |
EP3236851A1 (en) | Method and system for determining a blood glucose level for a patient | |
CN101594823A (en) | The non-invasive methods that uses metabolic measurements estimation blood sugar level to change | |
Zhao et al. | An automatic glucose monitoring signal denoising method with noise level estimation and responsive filter updating | |
MXPA04008357A (en) | Analytical device with prediction module and related methods. | |
CN115064261A (en) | Blood glucose statistical system and method based on artificial intelligence | |
US20190246964A1 (en) | Combined Non Invasive Blood Glucose Monitor Device | |
US20220039698A1 (en) | Method and System for Enhancing Monitoring Glucose | |
US20170363565A1 (en) | Analyte Measurement | |
US20230218204A1 (en) | Personalized calibration for glucose sensing | |
Biswas et al. | Learning the Sensing Delay for Personalized Continuous Health Monitoring | |
Gioe et al. | How well do self-monitored capillary glucose measurements predict the mean blood glucose from 24-hour continuous monitoring in endocrine practice? | |
US20170363566A1 (en) | Method of Determining Parameters of a Test Fluid | |
US20160100779A1 (en) | Combined non invasive blood glucose monitor device | |
陳榮興 et al. | Accuracy of the FreeStyle Method: Measuring Blood-Glucose by Skin-Prick Blood Extraction From Forearm | |
Thomaseth | Optimal experiment design for assessing plasma clearance of/sup 125/I-Iothalamate in peritoneal dialysis | |
Shi et al. | Modeling of relationship between glucose concentration in blood and glucose concentration in interstitial fluid | |
Khoury | Analysis of transient failures in continuous glucose monitoring sensors: data modeling and simulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FZDE | Discontinued |