CA2504661C - Method and device for monitoring an analyte concentration in the living body of a human or animal - Google Patents
Method and device for monitoring an analyte concentration in the living body of a human or animal Download PDFInfo
- Publication number
- CA2504661C CA2504661C CA002504661A CA2504661A CA2504661C CA 2504661 C CA2504661 C CA 2504661C CA 002504661 A CA002504661 A CA 002504661A CA 2504661 A CA2504661 A CA 2504661A CA 2504661 C CA2504661 C CA 2504661C
- Authority
- CA
- Canada
- Prior art keywords
- delta
- analyte
- values
- determined
- function
- 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.)
- Expired - Fee Related
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/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
Abstract
The invention relates to a method for continuous monitoring of an analyte concentration by determining its progression in the living body of a human or animal, in which analyte values y(t n) correlating with the concentration to be determined are determined for consecutive points in time t n. For predicting a prediction value for an analyte value y(t n0+.DELTA.t) over a prediction period .DELTA.t a function F (t k, t k-.DELTA.n, t k-2.DELTA.n, ..., t k-(m- 2).DELTA.n, t k-(m-1).DELTA.n) , which depends on analyte values y(t k) , y(t k- .DELTA.n) , y(t k-2.DELTA.n) , ..., y(t k-(m-2).DELTA.n, y(t k-(m-1).DELTA.n) is determined, which can be used to approximate the progression of the analyte values y(t n) at time t n0 in the vicinity U of an analyte value y(t n0) with a pre-determined accuracy .sigma., such that (see formula I) whereby .DELTA.n is an integer, and-subsequently the function F is used to calculate the prediction value y(t n0 + .DELTA.t) =F(t n0 , t n0 - .DELTA.n, ... ..., t(n0 - (m - 2).DELTA.n),t(n0 - (m - 1).DELTA.n))
Description
Applicant: Roche Diagnostics GmbH, Mannheim, Germany Method and Device for Monitoring an Analyte Concentration in the Living Body of a Human or Animal The invention concerns a method and a device for monitoring an analyte concentration in the living body of a human or animal, in which analyte values y(tn) correlating with the concentration to be determined are determined for consecutive points in time tn.
A method of this type and a suitable device are described in US 5,507,288, for example. This concerns in particular the continuous monitoring of the glucose concentration in the body of a patient, which is of crucial medical significance. Study results have demonstrated that extremely serious long-term corollaries of diabetes mellitus (for example blinding of the patient due to retinopathy) can be prevented by continuous in vivo monitoring of the time profile of the glucose concentration. This allows the dosage of the medication required (insulin) to be adjusted exactly at any point in time and thus the blood sugar levels to be kept within narrow limits at all times like in a healthy person.
In order for the medication process to be as efficient as possible it is desired to be able to determine not only a current value of the concentration of an analyte, but also predict future concentration values. The prediction of future values of the concentration of an analyte is of particular significance in the case of diabetes, not only in order to be able to counteract an increase of the blood glucose concentration beyond a critical value by administering insulin, but also to alert the patient to an imminent decrease of the blood glucose concentration below a critical value and the associated danger of fainting such that this can be prevented by timely uptake io of carbohydrates.
The prediction of a future concentration of an analyte value in the living body of a human or animal is made significantly more complicated by the fact that concentration values depend not only on historic values, but also on other parameters, such as food and medication intake or physical exercise.
A method for monitoring the blood glucose concentration is known from US 6,272,480 B1, in which an artificial neuronal network is used to derive a prediction of future blood glucose values from a time series of glucose concentrations that were measured. However, a significant disadvantage of artificial neuronal networks is that they need to be trained prior to their use and yield unreliable results when confronted with time series data outside of their training range.
From EP 1102194 A2, a method is known in which blood glucose concentration values measured at consecutive points in time are extrapolated taking into consideration patient-specific data, such as, in particular, any carbohydrate intake and insulin administration. For practical applications, the need to account for such patient-specific data is a major disadvantage, since accounting for such data is difficult and often error-prone.
The invention seeks to provide means for continuous monitoring of the concentration of an analyte in the living body of a human or animal and predicting future concentration values. In particular, this method shall, aside from data of historic concentration values, not require any patient-specific data, such as food or insulin intake, for the prediction of future concentration values.
In one aspect, the invention provides a method for continuous monitoring of an analyte concentration by determining its progression in the living body of a human or animal, in which analyte values y(tn) correlating with the concentration to be determined are determined for consecutive points in time tn, wherein for predicting an analyte value y(tno+Ot) over a prediction period At, a function F(tk, tk-An, tk-20n, ... , tk-(m-2)4ni tk-(m-1)An) 1S
determined, which depends on analyte values y(tk) , y(tk-on) Y( tk-2on) ,===,Y( tk- (m-2)on, Y( tk- (m-1) on) and which can be used to approximate the progression of the analyte values y(tn) at time tno in a vicinity U of an analyte value y(tno) with a predetermined accuracy a, such that . . .
62 [Y(tk) - F(tk-At, tk-On - At, tk-20n - At, tk e U
. . . , t(k - (m - 2)On) - At, t(k - (m -1)On) - Ot ) ]2 (eq. 1) , whereby An, no and m are integers, and calculation of the prediction value using the following equation:
Y(tnO + Ot) =F(tnO, tn0 - An, . . .
= = = I t(n0 - (m - 2)On),t(n0 - (m - l)On)) (eq. 2).
In another aspect, the invention provides a device for continuous monitoring of a concentration of an analyte by determining its progression in the living body of a human or animal, comprising a measuring unit, by means of which measured values of a measuring parameter that is correlated with the concentration to be determined are measured as measuring signals at consecutive points in time tn, and an analytical unit, by means of which analyte values y(tn) that are correlated with the concentration are determined from the measured values, the analytical unit being adapted for predicting an analyte value y(tno+At) over a prediction period Ot by determining a function F(tk, tk-An, tk-20n, tk-(m-2)4ni tk-(m-1)On) , which depends on analyte values y(tk) , y ( tk-on) , y( tk-2An) i ... . Y( tk- (m-2) Ant Y( tk- (m-1) On) , and can be used to approximate the progression of the analyte values y(tn) at time tno in a vicinity U of an analyte value y(tno) with a predetermined accuracy a, such that 6 2 1 tk) - ~ tk-At, tk-dn - At, tk-2An - At, ...
L( tk E U
. . . , t(k-(m-2)4n) - At, t(k-(m-1)On) - Ot) ]2 whereby An, no and m are integers, and subsequently the function F is used by the analytical unit to calculate a prediction value Y(tnO + Ot) =F(tnO, tn0 - On, . . .
A method of this type and a suitable device are described in US 5,507,288, for example. This concerns in particular the continuous monitoring of the glucose concentration in the body of a patient, which is of crucial medical significance. Study results have demonstrated that extremely serious long-term corollaries of diabetes mellitus (for example blinding of the patient due to retinopathy) can be prevented by continuous in vivo monitoring of the time profile of the glucose concentration. This allows the dosage of the medication required (insulin) to be adjusted exactly at any point in time and thus the blood sugar levels to be kept within narrow limits at all times like in a healthy person.
In order for the medication process to be as efficient as possible it is desired to be able to determine not only a current value of the concentration of an analyte, but also predict future concentration values. The prediction of future values of the concentration of an analyte is of particular significance in the case of diabetes, not only in order to be able to counteract an increase of the blood glucose concentration beyond a critical value by administering insulin, but also to alert the patient to an imminent decrease of the blood glucose concentration below a critical value and the associated danger of fainting such that this can be prevented by timely uptake io of carbohydrates.
The prediction of a future concentration of an analyte value in the living body of a human or animal is made significantly more complicated by the fact that concentration values depend not only on historic values, but also on other parameters, such as food and medication intake or physical exercise.
A method for monitoring the blood glucose concentration is known from US 6,272,480 B1, in which an artificial neuronal network is used to derive a prediction of future blood glucose values from a time series of glucose concentrations that were measured. However, a significant disadvantage of artificial neuronal networks is that they need to be trained prior to their use and yield unreliable results when confronted with time series data outside of their training range.
From EP 1102194 A2, a method is known in which blood glucose concentration values measured at consecutive points in time are extrapolated taking into consideration patient-specific data, such as, in particular, any carbohydrate intake and insulin administration. For practical applications, the need to account for such patient-specific data is a major disadvantage, since accounting for such data is difficult and often error-prone.
The invention seeks to provide means for continuous monitoring of the concentration of an analyte in the living body of a human or animal and predicting future concentration values. In particular, this method shall, aside from data of historic concentration values, not require any patient-specific data, such as food or insulin intake, for the prediction of future concentration values.
In one aspect, the invention provides a method for continuous monitoring of an analyte concentration by determining its progression in the living body of a human or animal, in which analyte values y(tn) correlating with the concentration to be determined are determined for consecutive points in time tn, wherein for predicting an analyte value y(tno+Ot) over a prediction period At, a function F(tk, tk-An, tk-20n, ... , tk-(m-2)4ni tk-(m-1)An) 1S
determined, which depends on analyte values y(tk) , y(tk-on) Y( tk-2on) ,===,Y( tk- (m-2)on, Y( tk- (m-1) on) and which can be used to approximate the progression of the analyte values y(tn) at time tno in a vicinity U of an analyte value y(tno) with a predetermined accuracy a, such that . . .
62 [Y(tk) - F(tk-At, tk-On - At, tk-20n - At, tk e U
. . . , t(k - (m - 2)On) - At, t(k - (m -1)On) - Ot ) ]2 (eq. 1) , whereby An, no and m are integers, and calculation of the prediction value using the following equation:
Y(tnO + Ot) =F(tnO, tn0 - An, . . .
= = = I t(n0 - (m - 2)On),t(n0 - (m - l)On)) (eq. 2).
In another aspect, the invention provides a device for continuous monitoring of a concentration of an analyte by determining its progression in the living body of a human or animal, comprising a measuring unit, by means of which measured values of a measuring parameter that is correlated with the concentration to be determined are measured as measuring signals at consecutive points in time tn, and an analytical unit, by means of which analyte values y(tn) that are correlated with the concentration are determined from the measured values, the analytical unit being adapted for predicting an analyte value y(tno+At) over a prediction period Ot by determining a function F(tk, tk-An, tk-20n, tk-(m-2)4ni tk-(m-1)On) , which depends on analyte values y(tk) , y ( tk-on) , y( tk-2An) i ... . Y( tk- (m-2) Ant Y( tk- (m-1) On) , and can be used to approximate the progression of the analyte values y(tn) at time tno in a vicinity U of an analyte value y(tno) with a predetermined accuracy a, such that 6 2 1 tk) - ~ tk-At, tk-dn - At, tk-2An - At, ...
L( tk E U
. . . , t(k-(m-2)4n) - At, t(k-(m-1)On) - Ot) ]2 whereby An, no and m are integers, and subsequently the function F is used by the analytical unit to calculate a prediction value Y(tnO + Ot) =F(tnO, tn0 - On, . . .
5 = = = I t(n0 - (m - 2)On),t(n0 - (m - l)On)).
According to the method of the invention, the living body of a human or animal and the concentration of an analyte that is present therein are mathematically treated as a deterministic dynamic system with some possible interference by noise. The development of this system over time is then described as a phase trajectory in a phase space whose dimensionality depends on the number of independent variables influencing the concentration of the analyte to be determined. In the case of diabetes, for example, the number of ingested bread units and the quantity of insulin administered are variables of this type, on which future values of the blood glucose concentration depend.
For practical purposes, the number of variables, on which the analyte concentration to be determined depends, is usually unknown and consequently the dimensionality of the phase space, in which the development of the system over time can be described by means of the phase trajectory, is also unknown.
The invention utilizes the fact that there is no need to know the entire phase trajectory of the system in a possibly high-dimensional phase space in order to be able to predict an analyte value. For extrapolating a point on the phase trajectory it is only necessary to know the phase trajectory in a vicinity of that point. In particular, there is no need to consider the entire phase space of the system for extrapolation, but rather it is sufficient to study the dynamics in an artificial, reconstructed phase space of significantly lower dimensionality.
It is sufficient if the dimensionality is high enough to Zo account for all dominant degrees of freedom driving the dynamics of the system. Degrees of freedom with a negligible effect on the dynamics of the system can be left unaccounted for such that the dimensionality of the reconstructed phase space can be selected accordingly lower.
Negligible degrees of freedom are found especially for dynamic systems, in which a sub-system is weakly coupled to the remaining system. Considering diabetes as an example, a sub-system, that depends on the administered quantity of insulin, ingested bread units and physical exercise, is weakly coupled to hormonal regulation mechanisms, which also exert an influence on the blood sugar levels.
A time series of analyte values y(tn) determined for consecutive points in time tn contains in principle all of the information regarding the dynamic behavior of the system. Therefore, from a finite number of analyte values y(tn) determined for consecutive points in time tn - often as few as approx. 10 analyte values in the vicinity of the point to be extrapolated are sufficient - a local approximation of the phase trajectory in the vicinity of an analyte value y(tõo) at point in time tno can be determined and extrapolated into the future by a prediction period At such that a prediction value of an analyte value y(tno+Ot) can be calculated.
It is important that the analyte values used to determine the function F and thus a local approximation of the phase trajectory were not determined for points in time tn that are too close to each other. Otherwise there is a risk that an observed progression is based mainly on measuring errors and noise effects and does not reflect io the dynamics over time of a concentration change of the analyte.
In order to be able to reduce the influence of inevitable measuring errors to the extent possible by statistical i.s methods it is desired to determine analyte values y(tn) for consecutive points in time tn that are as close to each other as possible. Therefore, it may be necessary not to use all analyte values y(tn) determined for directly consecutive points in time tn for determining a 20 local approximation of the phase trajectory, but rather only every second, third or n-th such analyte value that has been determined.
For this reason, equation 1 includes a parameter An to 25 ensure that the analyte values y(tn) that are used to determine the function F were measured for consecutive points in time tn that are separated by sufficiently large time intervals.
30 Whether or not a function F as determined can be used to approximate the time profile of the analyte values y(tn) with a pre-determined accuracy, can be tested by means of equation 1. In principle, it is necessary to sum up over all points in time tk from the vicinity U of the analyte 35 value y(tno) at time tno. However, this does not exclude that individual analyte values y(tn) recognized as erroneous, for example by means of applying suitable filters, are ignored. The function F can be a non-linear function and depend not only on the analyte values, but in addition on their first and second time derivative.
However, it is preferable for the function F to be a linear function. It has proven to be particularly favorable to represent the function F using coefficients ao to am as follows:
m F = ao+ E y(t(k_(m_j) An))aj j=l In practice, the coefficients, ao to am, can be determined by minimizing the following sum:
m 2 Y(tk) - ao - Y(t(k-(m-i) An) - Ot) ai i=1 tkEU
In doing so, the vicinity U must be selected such that it contains a sufficient number of points in time tn for the determination of the coefficients, ao to am.
The method according to the invention can be applied in several steps in order to determine from a prediction value y( tno+At ) another prediction value y( to+20t ) even further in the future. A prediction value y(to+nAt) can thus be determined by multiple application of the method.
In the process, it is important that the coefficients, ao to am, used to approximate the progression of the analyte values in the vicinity U are re-determined in each step for a vicinity U that is shifted by the prediction period At.
According to the method of the invention, the living body of a human or animal and the concentration of an analyte that is present therein are mathematically treated as a deterministic dynamic system with some possible interference by noise. The development of this system over time is then described as a phase trajectory in a phase space whose dimensionality depends on the number of independent variables influencing the concentration of the analyte to be determined. In the case of diabetes, for example, the number of ingested bread units and the quantity of insulin administered are variables of this type, on which future values of the blood glucose concentration depend.
For practical purposes, the number of variables, on which the analyte concentration to be determined depends, is usually unknown and consequently the dimensionality of the phase space, in which the development of the system over time can be described by means of the phase trajectory, is also unknown.
The invention utilizes the fact that there is no need to know the entire phase trajectory of the system in a possibly high-dimensional phase space in order to be able to predict an analyte value. For extrapolating a point on the phase trajectory it is only necessary to know the phase trajectory in a vicinity of that point. In particular, there is no need to consider the entire phase space of the system for extrapolation, but rather it is sufficient to study the dynamics in an artificial, reconstructed phase space of significantly lower dimensionality.
It is sufficient if the dimensionality is high enough to Zo account for all dominant degrees of freedom driving the dynamics of the system. Degrees of freedom with a negligible effect on the dynamics of the system can be left unaccounted for such that the dimensionality of the reconstructed phase space can be selected accordingly lower.
Negligible degrees of freedom are found especially for dynamic systems, in which a sub-system is weakly coupled to the remaining system. Considering diabetes as an example, a sub-system, that depends on the administered quantity of insulin, ingested bread units and physical exercise, is weakly coupled to hormonal regulation mechanisms, which also exert an influence on the blood sugar levels.
A time series of analyte values y(tn) determined for consecutive points in time tn contains in principle all of the information regarding the dynamic behavior of the system. Therefore, from a finite number of analyte values y(tn) determined for consecutive points in time tn - often as few as approx. 10 analyte values in the vicinity of the point to be extrapolated are sufficient - a local approximation of the phase trajectory in the vicinity of an analyte value y(tõo) at point in time tno can be determined and extrapolated into the future by a prediction period At such that a prediction value of an analyte value y(tno+Ot) can be calculated.
It is important that the analyte values used to determine the function F and thus a local approximation of the phase trajectory were not determined for points in time tn that are too close to each other. Otherwise there is a risk that an observed progression is based mainly on measuring errors and noise effects and does not reflect io the dynamics over time of a concentration change of the analyte.
In order to be able to reduce the influence of inevitable measuring errors to the extent possible by statistical i.s methods it is desired to determine analyte values y(tn) for consecutive points in time tn that are as close to each other as possible. Therefore, it may be necessary not to use all analyte values y(tn) determined for directly consecutive points in time tn for determining a 20 local approximation of the phase trajectory, but rather only every second, third or n-th such analyte value that has been determined.
For this reason, equation 1 includes a parameter An to 25 ensure that the analyte values y(tn) that are used to determine the function F were measured for consecutive points in time tn that are separated by sufficiently large time intervals.
30 Whether or not a function F as determined can be used to approximate the time profile of the analyte values y(tn) with a pre-determined accuracy, can be tested by means of equation 1. In principle, it is necessary to sum up over all points in time tk from the vicinity U of the analyte 35 value y(tno) at time tno. However, this does not exclude that individual analyte values y(tn) recognized as erroneous, for example by means of applying suitable filters, are ignored. The function F can be a non-linear function and depend not only on the analyte values, but in addition on their first and second time derivative.
However, it is preferable for the function F to be a linear function. It has proven to be particularly favorable to represent the function F using coefficients ao to am as follows:
m F = ao+ E y(t(k_(m_j) An))aj j=l In practice, the coefficients, ao to am, can be determined by minimizing the following sum:
m 2 Y(tk) - ao - Y(t(k-(m-i) An) - Ot) ai i=1 tkEU
In doing so, the vicinity U must be selected such that it contains a sufficient number of points in time tn for the determination of the coefficients, ao to am.
The method according to the invention can be applied in several steps in order to determine from a prediction value y( tno+At ) another prediction value y( to+20t ) even further in the future. A prediction value y(to+nAt) can thus be determined by multiple application of the method.
In the process, it is important that the coefficients, ao to am, used to approximate the progression of the analyte values in the vicinity U are re-determined in each step for a vicinity U that is shifted by the prediction period At.
Instead of predicting an analyte value y(tn) for a point in time tn in the future by means of multiple application of the method in multiple steps, it is also feasible to use a suitably larger prediction period Ot and coefficients, ao to am, determined accordingly to predict in a single step an analyte value y(tn+Ot) in the distant future. However, it has been found that stepwise application of the method provides for a more accurate prediction over extended periods of time.
The finding that the dynamic behavior of a deterministic system can be described by a finite number of analyte values y(tn) that are determined for consecutive points in time~tn is based on the "Takens time delay embedding"
method described in the publication by T. Schreiber in Phys. Rep. 85, 1-64 (1999). According to the terminology used in the publication of T. Schreiber, the interval between consecutive points in time, for which analyte values y(tn) are determined, is called delay time.
Reference is made to the publication of T. Schreiber in a supplementing fashion.
The invention and preferred embodiments are illustrated in more detail in the following with reference to the figures. The features illustrated and described in the following can be used individually or in combination in order to create preferred further refinements of the invention. In the figures:
Fig. 1 shows a block diagram of a device according to the invention;
Fig. 2 shows a schematic sketch of a sensor suited for the invention;
Fig. 3 shows a measuring signal of a sensor according 5 to Fig. 2;
Fig. 4 shows measured blood glucose concentration values and blood glucose concentration values predicted for a period of 30 minutes; and Fig. 5 shows a schematic example of a phase space with data points.
The essential components of a device according to the invention are shown in Fig. 1. A sensor 1 measures values at consecutive points in time tn. The resulting measurement signal is then transmitted - wireless in the case shown - to a receiver 2 which in turn transmits the measurement signal to an analytical unit 3 containing a microprocessor 4 and a data storage unit 5.
Alternatively, both data and commands can be transmitted to the analytical unit 3 by an input unit 6. The results are output by an output unit 7 which may include a display and other common output means. It is self-evident that data processing in the analytical unit 3 is performed in a digital fashion, and corresponding transformers for transformating analogous signals into digital signals are provided. More detailed explanations are not required in this context, since the basic design of devices of this type is known (from US 5,507,288, for example) and the invention is suitable for a broad range of measuring techniques, in which - as described above -different measuring signals that are correlated with the desired working signal are measured.
Fig. 2 shows in the form of a schematic sketch a sensor 1, in which an implantable catheter 10 is used to aspirate interstitial fluid from the subcutaneous fatty tissue by means of a pump 11 and to aspirate the fluid through a photometric measuring unit 12 into a disposal container 13. Inside the photometric measuring unit 12, the tubing 14, through which the fluid is being transported, contains a transparent measuring cell 15, into which primary light emitted by a light sensor 16 is guided. The secondary light resulting from passage of the measuring cell 15 is then measured by means of a photo-detector 17 and the resulting signal processed by means of measuring electronics (not shown) to generate a raw signal, which - as is shown in Fig. 1 in an exemplary fashion - is transmitted to an analytical unit 3. The analytical unit 3 is used to determine from the raw signal analyte values y(tn) that are correlated with the concentration to be determined for consecutive points in time tn. In the scope of the exemplary embodiment, these analyte values y(tn) represent the blood glucose level of a subject at a time tn.
It is self-evident that the analyte values y(tn) can also be determined by other methods, for example methods based on an enzymatic conversion of glucose and subsequent electrochemical detection of the hydrogen peroxide formed in the process.
Fig. 3 shows a curve representing the progression of the blood glucose level of a subject. The two crosses drawn in Fig. 3 denote, in an exemplary fashion, one prediction value each that was determined for a prediction period of 30 minutes by the method to be described in the following. As shown by the good agreement of the left cross and the curve represented as a continuous line, the accuracy of the predicted values is good. It is self-evident that this applies only if the blood glucose level is not affected by an administration of insulin or uptake of carbohydrates during the prediction interval. The right cross in Fig. 3 shows a prediction value that is below a critical threshold value. In the case of a prediction value being this low, a warning signal is generated in order to alert a patient to a dangerous decrease in the blood glucose level such that the io decrease can be counteracted by ingestion of carbohydrates.
The analyte values y(tn) shown in Fig. 3 in the form of a curve were determined from measured values that were recorded at intervals of 1 second each. To use a measuring rate this high is advantageous in that sufficient data is available for statistical processing of the measured values, but it is not obligatory. For many applications, it is sufficient to determine measured values at intervals of a few minutes, e.g. every 2 to 5 minutes.
In the embodiment shown, the measured values were processed by means of a filter algorithm that includes an operation, in which the influence of a current measured value on the working signal is weighted by means of a weighting factor ("controllable filter algorithm"), whereby a signal variation parameter (related in each case to the current point in time, i.e. changing over time) is determined on the basis of signal variations detected in a close temporal relationship to the measurement of the current measured value, and the weighting factor is adapted dynamically as a function of the signal variation parameter determined for the current measured value for the point in time. Filter algorithms of this type are known as Kalmann filter algorithms and have been described, for example, in DE 10343863.7.
For the purposes of the method described in the following, extensive statistical processing of the measured values is not obligatory. Alternatively, it is also feasible, for example, to average over several measured values to obtain an analyte value y(tn) or to obtain exactly one analyte value y(tn) from each measured value.
The analyte values y(tn) determined are then used to determine first the function, F(tk,tk-en,tk-2Qn,===rtk-(m-2)On.tk-(m-1)An), which depends on analyte values y(tk) , Y(tk-On) , Y(tk-2An) .====- Y(tk-(m-2)An) , Y(tk-(m-1)An) and approximates the progression of the analyte values in a vicinity U of a point in time tno at a pre-determined accuracy 62. Although the progression of the analyte values y(tn) is usually not linear, it has been found that a linear function of the type m F- a0+ Y_ y [t(k(f))Jai j=1 with coefficients ao to am suffices. The value of the parameter, m, of function F(tk,tk-An,tk-2An,--=. tk- (m-2)On.tk-(m-1)On) (e.g. in am and tk-(m-1)An) describes the dimensionality of the phase space, in which the dynamical behavior of the system is approximated.
The coefficients, ao to am, can be determined by -minimizing the equation m 2 62 Y(tk) - ao - Y(t(k-(m-i) On) - Ot) ai tk EU i=1 for various analyte values y(tk) from the vicinity U. If measured values are detected at a high measuring rate and analyte values y(tn) are determined for accordingly tightly spaced points in time tn, it is advantageous not to include each and every analyte value that was determined in the equation above for the determination of the coefficients ao bis am, but instead, for example, only every second or third analyte value such that An is not necessarily equal to 1, but rather can be a different integer. Although, in principle, the sum for all points in time tk in the vicinity U of the time tno should be obtained in the equation stated above, this does not exclude that analyte values y(tk) recognized to be outliers, for example by means of suitable filters, are not included in the analysis.
Ideally, the times tn, for which the analyte values y(tn) are determined, should be at equidistant time intervals from each other. However, analyte values in a living body of a human or animal possess such high auto-correlation when viewed over extended periods of time, that minor deviations from a constant measuring rate affect the accuracy of the prediction values only to a small degree.
For example, the concentrations of lactate - under aerobic conditions - or creatinine display high auto-correlation over a period of several hours. Even relatively rapidly varying analytes, such as glucose, show an auto-correlation coefficient in excess of 0.8 over a period of 30 minutes.
The value of the parameter m in the function F(tk,tk-On, tk-2An, =.., tk- (m-2)Qn, tk- (m-1)pn) , i. e. the dimensionality of the relevant phase space, cannot be stated generally for all possible cases, in which the concentration of an analyte in the living body of a human or animal must be monitored. This is due mainly to the fact that the 5 concentration of the various analytes of medical interest depend on different numbers of independent variables. In the case of the blood glucose concentration of a diabetes patient, two such variables are known, i.e. the quantity of bread units ingested and the quantity of insulin 10 administered. In general, it can be said that the parameter does not need to be larger than twice the number of independent variables plus 1. Therefore, m = 5 in the present exemplary embodiment illustrated by means of Fig. 3.
As has been mentioned above, the value of the parameter m corresponds to the dimensionality of a phase space reconstructed according to the "Takens time delay embedding" method, in which the dynamic behavior of the system can be described. For practical purposes, m usually has at most a value between 3 and 10, preferably between 5 and 8, and can be determined empirically. One option for determining the dimensionality of the phase space and thus the value of this parameter is provided by the method of so-called "false nearest neighbour"
analysis as described in a publication by M.B. Kennel et al., Phys. Rev. A 45 3403 (1995). Reference is made to this publication in a supplementing fashion.
Once the coefficients, ao to am, have been determined they can be used to calculate a prediction value m y (tn0 + At) = a0 + I Y (t(nO_(m_j) An)l aj j=1 Subsequently, a prediction value determined accordingly is compared to a predetermined lower and/or upper threshold value and a warning signal is generated if that value is below or above the threshold value, respectively.
When a patient responds to a warning signal of this type by ingesting carbohydrates or administering a dose of insulin, it takes a while before the blood glucose concentration is affected by such ingestion of carbohydrates or administration of insulin. For this reason, it is usually necessary to predict analyte value 10 to 90 minutes, preferably 20 to 40 min, into the future such that there is sufficient time to counteract a dangerous increase or decrease of the analyte concentration. It has been found that analyte values are better predicted over extended periods of time if, rather than selecting an accordingly large prediction period (e.g. At = 30 min) in a single step, the procedure is applied in multiple smaller steps of 5 minutes, for example. In each additional step coefficients ao to am are determined for an already calculated prediction value y(tno+At), which coefficients are used to approximate the progression of analyte values y(tn) in a vicinity U of the analyte value at point in time tno+Ot. The results are used to calculate another prediction value y(tno+20t). In a stepwise procedure of this type, the prediction period At of the individual steps should be selected to be equal to the interval between times tn, for which the analyte values y(tn) used for calculating the prediction values were determined. If the time interval between consecutive points in time tn equals dt, this means that the prediction period is given by At = dt An.
The method described above cannot only be used to predict future analyte values, but also to test analyte values y(tn) that were determined for the past and identify false analyte values y(tn) which are due to outliers or measuring errors. The generally very good agreement between the measured analyte values y(tn) and the prediction values is shown in Figure 4 in an exemplary fashion. The smooth continuous curve represents the time profile of the prediction values, whereas the curve showing strong short-term variations represents the time profile of the measured values. The prediction period was 2 min, and a total of 15 steps were performed consecutively such that analyte values were predicted for a period of 30 minutes. Therefore, the y axis does not represent the blood glucose concentration in Figure 4, but rather the intensity of the measuring signal on which the blood glucose determination is based as detected electrochemically (in units of nanoampere).
The analyte values y (tk) , y (tk-An) , y (tk-2,&n) , = = = = , Y (tk- (m-2)On) ,Y (tk- (m-1),&n) , on which the function F depends, are considered in the "Takens time delay embedding" method to represent the coordinates of a reconstructed phase space, in which the dynamic behavior of the system is approximated and extrapolated. These coordinates are also called delay coordinates. For the calculation, it can be advantageous to transform the coordinates y(tk),y(tk-On) , Y (tk-2Qn) , = = = = , Y (tk- (m-2)On) , Y (tk- (m-1)On) and determine F
as a function of a set of transformed coordinates Ty(tk), Ty (tk-Qn) , = = = , Ty (tk- (m-1)On) =
It is preferable to perform a linear transformation on the set of coordinates. Hereby, it is useful for the transformed coordinates to be so-called latent coordinates aligned along the main axes of a so-called ellipsoid of variance, as is illustrated in Figure S. An ellipsoid of variance of this type is drawn in a two-dimensional phase space in Figure 5, which phase space is represented by its delay coordinates, y(t) and y(t-k).
io The ellipsoid of variance surrounds data points that are indicated by black circles, which, when connected according to their sequence in time, form a phase trajectory.
The ellipsoid of variance can be calculated from the scattering of the data points in the phase space, whereby the main axes of the ellipsoid indicate the variance of the data points in the corresponding direction of the phase space. Figure 5 also shows the new phase space obtained by transformation along with its coordinates, T(1) and T(2). The first latent coordinate T(1) extends along the largest axis of inertia of the ellipsoid of variance, the second latent coordinate T(2) along the second-to-largest axis of inertia. Numerically, a transformation to latent coordinates can be performed for example by so-called principal component analysis (M.
Otto, Chemometrics, Wiley, Weinheim (Germany), 1999, pages 124 to 133) or singular value decomposition (W.H.
Press et al., Numerical Recipes, Cambridge University Press, Cambridge (UK), 1989, pages 52 to 61).
By performing a transformation to latent coordinates of this type, it is possible to achieve that one or several of the transformed coordinates have a negligible effect on the function F at the predetermined accuracy 6 such that only a correspondingly smaller number of transformed coordinates needs to be included in the analysis in the calculation of function F.
These and other features and advantages of the present invention should be more fully understood from the above detailed description of the invention taken together with the accompanying claims. It is noted that the scope of the claims is defined by the recitations therein and not by the specific discussion of the features and advantages set forth in the present description.
It is noted that terms like "preferably", "commonly" and "typically" are not utilized therein to limit the scope of the claimed invention or to imply that certain features are critical, essential or even important to applications of the claimed invention. Rather, these terms are merely intended to highlight alternative or additional features that may or may not be utilized in a particular embodiment of the present invention.
The finding that the dynamic behavior of a deterministic system can be described by a finite number of analyte values y(tn) that are determined for consecutive points in time~tn is based on the "Takens time delay embedding"
method described in the publication by T. Schreiber in Phys. Rep. 85, 1-64 (1999). According to the terminology used in the publication of T. Schreiber, the interval between consecutive points in time, for which analyte values y(tn) are determined, is called delay time.
Reference is made to the publication of T. Schreiber in a supplementing fashion.
The invention and preferred embodiments are illustrated in more detail in the following with reference to the figures. The features illustrated and described in the following can be used individually or in combination in order to create preferred further refinements of the invention. In the figures:
Fig. 1 shows a block diagram of a device according to the invention;
Fig. 2 shows a schematic sketch of a sensor suited for the invention;
Fig. 3 shows a measuring signal of a sensor according 5 to Fig. 2;
Fig. 4 shows measured blood glucose concentration values and blood glucose concentration values predicted for a period of 30 minutes; and Fig. 5 shows a schematic example of a phase space with data points.
The essential components of a device according to the invention are shown in Fig. 1. A sensor 1 measures values at consecutive points in time tn. The resulting measurement signal is then transmitted - wireless in the case shown - to a receiver 2 which in turn transmits the measurement signal to an analytical unit 3 containing a microprocessor 4 and a data storage unit 5.
Alternatively, both data and commands can be transmitted to the analytical unit 3 by an input unit 6. The results are output by an output unit 7 which may include a display and other common output means. It is self-evident that data processing in the analytical unit 3 is performed in a digital fashion, and corresponding transformers for transformating analogous signals into digital signals are provided. More detailed explanations are not required in this context, since the basic design of devices of this type is known (from US 5,507,288, for example) and the invention is suitable for a broad range of measuring techniques, in which - as described above -different measuring signals that are correlated with the desired working signal are measured.
Fig. 2 shows in the form of a schematic sketch a sensor 1, in which an implantable catheter 10 is used to aspirate interstitial fluid from the subcutaneous fatty tissue by means of a pump 11 and to aspirate the fluid through a photometric measuring unit 12 into a disposal container 13. Inside the photometric measuring unit 12, the tubing 14, through which the fluid is being transported, contains a transparent measuring cell 15, into which primary light emitted by a light sensor 16 is guided. The secondary light resulting from passage of the measuring cell 15 is then measured by means of a photo-detector 17 and the resulting signal processed by means of measuring electronics (not shown) to generate a raw signal, which - as is shown in Fig. 1 in an exemplary fashion - is transmitted to an analytical unit 3. The analytical unit 3 is used to determine from the raw signal analyte values y(tn) that are correlated with the concentration to be determined for consecutive points in time tn. In the scope of the exemplary embodiment, these analyte values y(tn) represent the blood glucose level of a subject at a time tn.
It is self-evident that the analyte values y(tn) can also be determined by other methods, for example methods based on an enzymatic conversion of glucose and subsequent electrochemical detection of the hydrogen peroxide formed in the process.
Fig. 3 shows a curve representing the progression of the blood glucose level of a subject. The two crosses drawn in Fig. 3 denote, in an exemplary fashion, one prediction value each that was determined for a prediction period of 30 minutes by the method to be described in the following. As shown by the good agreement of the left cross and the curve represented as a continuous line, the accuracy of the predicted values is good. It is self-evident that this applies only if the blood glucose level is not affected by an administration of insulin or uptake of carbohydrates during the prediction interval. The right cross in Fig. 3 shows a prediction value that is below a critical threshold value. In the case of a prediction value being this low, a warning signal is generated in order to alert a patient to a dangerous decrease in the blood glucose level such that the io decrease can be counteracted by ingestion of carbohydrates.
The analyte values y(tn) shown in Fig. 3 in the form of a curve were determined from measured values that were recorded at intervals of 1 second each. To use a measuring rate this high is advantageous in that sufficient data is available for statistical processing of the measured values, but it is not obligatory. For many applications, it is sufficient to determine measured values at intervals of a few minutes, e.g. every 2 to 5 minutes.
In the embodiment shown, the measured values were processed by means of a filter algorithm that includes an operation, in which the influence of a current measured value on the working signal is weighted by means of a weighting factor ("controllable filter algorithm"), whereby a signal variation parameter (related in each case to the current point in time, i.e. changing over time) is determined on the basis of signal variations detected in a close temporal relationship to the measurement of the current measured value, and the weighting factor is adapted dynamically as a function of the signal variation parameter determined for the current measured value for the point in time. Filter algorithms of this type are known as Kalmann filter algorithms and have been described, for example, in DE 10343863.7.
For the purposes of the method described in the following, extensive statistical processing of the measured values is not obligatory. Alternatively, it is also feasible, for example, to average over several measured values to obtain an analyte value y(tn) or to obtain exactly one analyte value y(tn) from each measured value.
The analyte values y(tn) determined are then used to determine first the function, F(tk,tk-en,tk-2Qn,===rtk-(m-2)On.tk-(m-1)An), which depends on analyte values y(tk) , Y(tk-On) , Y(tk-2An) .====- Y(tk-(m-2)An) , Y(tk-(m-1)An) and approximates the progression of the analyte values in a vicinity U of a point in time tno at a pre-determined accuracy 62. Although the progression of the analyte values y(tn) is usually not linear, it has been found that a linear function of the type m F- a0+ Y_ y [t(k(f))Jai j=1 with coefficients ao to am suffices. The value of the parameter, m, of function F(tk,tk-An,tk-2An,--=. tk- (m-2)On.tk-(m-1)On) (e.g. in am and tk-(m-1)An) describes the dimensionality of the phase space, in which the dynamical behavior of the system is approximated.
The coefficients, ao to am, can be determined by -minimizing the equation m 2 62 Y(tk) - ao - Y(t(k-(m-i) On) - Ot) ai tk EU i=1 for various analyte values y(tk) from the vicinity U. If measured values are detected at a high measuring rate and analyte values y(tn) are determined for accordingly tightly spaced points in time tn, it is advantageous not to include each and every analyte value that was determined in the equation above for the determination of the coefficients ao bis am, but instead, for example, only every second or third analyte value such that An is not necessarily equal to 1, but rather can be a different integer. Although, in principle, the sum for all points in time tk in the vicinity U of the time tno should be obtained in the equation stated above, this does not exclude that analyte values y(tk) recognized to be outliers, for example by means of suitable filters, are not included in the analysis.
Ideally, the times tn, for which the analyte values y(tn) are determined, should be at equidistant time intervals from each other. However, analyte values in a living body of a human or animal possess such high auto-correlation when viewed over extended periods of time, that minor deviations from a constant measuring rate affect the accuracy of the prediction values only to a small degree.
For example, the concentrations of lactate - under aerobic conditions - or creatinine display high auto-correlation over a period of several hours. Even relatively rapidly varying analytes, such as glucose, show an auto-correlation coefficient in excess of 0.8 over a period of 30 minutes.
The value of the parameter m in the function F(tk,tk-On, tk-2An, =.., tk- (m-2)Qn, tk- (m-1)pn) , i. e. the dimensionality of the relevant phase space, cannot be stated generally for all possible cases, in which the concentration of an analyte in the living body of a human or animal must be monitored. This is due mainly to the fact that the 5 concentration of the various analytes of medical interest depend on different numbers of independent variables. In the case of the blood glucose concentration of a diabetes patient, two such variables are known, i.e. the quantity of bread units ingested and the quantity of insulin 10 administered. In general, it can be said that the parameter does not need to be larger than twice the number of independent variables plus 1. Therefore, m = 5 in the present exemplary embodiment illustrated by means of Fig. 3.
As has been mentioned above, the value of the parameter m corresponds to the dimensionality of a phase space reconstructed according to the "Takens time delay embedding" method, in which the dynamic behavior of the system can be described. For practical purposes, m usually has at most a value between 3 and 10, preferably between 5 and 8, and can be determined empirically. One option for determining the dimensionality of the phase space and thus the value of this parameter is provided by the method of so-called "false nearest neighbour"
analysis as described in a publication by M.B. Kennel et al., Phys. Rev. A 45 3403 (1995). Reference is made to this publication in a supplementing fashion.
Once the coefficients, ao to am, have been determined they can be used to calculate a prediction value m y (tn0 + At) = a0 + I Y (t(nO_(m_j) An)l aj j=1 Subsequently, a prediction value determined accordingly is compared to a predetermined lower and/or upper threshold value and a warning signal is generated if that value is below or above the threshold value, respectively.
When a patient responds to a warning signal of this type by ingesting carbohydrates or administering a dose of insulin, it takes a while before the blood glucose concentration is affected by such ingestion of carbohydrates or administration of insulin. For this reason, it is usually necessary to predict analyte value 10 to 90 minutes, preferably 20 to 40 min, into the future such that there is sufficient time to counteract a dangerous increase or decrease of the analyte concentration. It has been found that analyte values are better predicted over extended periods of time if, rather than selecting an accordingly large prediction period (e.g. At = 30 min) in a single step, the procedure is applied in multiple smaller steps of 5 minutes, for example. In each additional step coefficients ao to am are determined for an already calculated prediction value y(tno+At), which coefficients are used to approximate the progression of analyte values y(tn) in a vicinity U of the analyte value at point in time tno+Ot. The results are used to calculate another prediction value y(tno+20t). In a stepwise procedure of this type, the prediction period At of the individual steps should be selected to be equal to the interval between times tn, for which the analyte values y(tn) used for calculating the prediction values were determined. If the time interval between consecutive points in time tn equals dt, this means that the prediction period is given by At = dt An.
The method described above cannot only be used to predict future analyte values, but also to test analyte values y(tn) that were determined for the past and identify false analyte values y(tn) which are due to outliers or measuring errors. The generally very good agreement between the measured analyte values y(tn) and the prediction values is shown in Figure 4 in an exemplary fashion. The smooth continuous curve represents the time profile of the prediction values, whereas the curve showing strong short-term variations represents the time profile of the measured values. The prediction period was 2 min, and a total of 15 steps were performed consecutively such that analyte values were predicted for a period of 30 minutes. Therefore, the y axis does not represent the blood glucose concentration in Figure 4, but rather the intensity of the measuring signal on which the blood glucose determination is based as detected electrochemically (in units of nanoampere).
The analyte values y (tk) , y (tk-An) , y (tk-2,&n) , = = = = , Y (tk- (m-2)On) ,Y (tk- (m-1),&n) , on which the function F depends, are considered in the "Takens time delay embedding" method to represent the coordinates of a reconstructed phase space, in which the dynamic behavior of the system is approximated and extrapolated. These coordinates are also called delay coordinates. For the calculation, it can be advantageous to transform the coordinates y(tk),y(tk-On) , Y (tk-2Qn) , = = = = , Y (tk- (m-2)On) , Y (tk- (m-1)On) and determine F
as a function of a set of transformed coordinates Ty(tk), Ty (tk-Qn) , = = = , Ty (tk- (m-1)On) =
It is preferable to perform a linear transformation on the set of coordinates. Hereby, it is useful for the transformed coordinates to be so-called latent coordinates aligned along the main axes of a so-called ellipsoid of variance, as is illustrated in Figure S. An ellipsoid of variance of this type is drawn in a two-dimensional phase space in Figure 5, which phase space is represented by its delay coordinates, y(t) and y(t-k).
io The ellipsoid of variance surrounds data points that are indicated by black circles, which, when connected according to their sequence in time, form a phase trajectory.
The ellipsoid of variance can be calculated from the scattering of the data points in the phase space, whereby the main axes of the ellipsoid indicate the variance of the data points in the corresponding direction of the phase space. Figure 5 also shows the new phase space obtained by transformation along with its coordinates, T(1) and T(2). The first latent coordinate T(1) extends along the largest axis of inertia of the ellipsoid of variance, the second latent coordinate T(2) along the second-to-largest axis of inertia. Numerically, a transformation to latent coordinates can be performed for example by so-called principal component analysis (M.
Otto, Chemometrics, Wiley, Weinheim (Germany), 1999, pages 124 to 133) or singular value decomposition (W.H.
Press et al., Numerical Recipes, Cambridge University Press, Cambridge (UK), 1989, pages 52 to 61).
By performing a transformation to latent coordinates of this type, it is possible to achieve that one or several of the transformed coordinates have a negligible effect on the function F at the predetermined accuracy 6 such that only a correspondingly smaller number of transformed coordinates needs to be included in the analysis in the calculation of function F.
These and other features and advantages of the present invention should be more fully understood from the above detailed description of the invention taken together with the accompanying claims. It is noted that the scope of the claims is defined by the recitations therein and not by the specific discussion of the features and advantages set forth in the present description.
It is noted that terms like "preferably", "commonly" and "typically" are not utilized therein to limit the scope of the claimed invention or to imply that certain features are critical, essential or even important to applications of the claimed invention. Rather, these terms are merely intended to highlight alternative or additional features that may or may not be utilized in a particular embodiment of the present invention.
Claims (16)
1. Method for continuous monitoring of an analyte concentration by determining its progression in the living body of a human or animal, in which analyte values y(t n) correlating with the concentration to be determined are determined for consecutive points in time t n, wherein for predicting an analyte value y(t n0+.DELTA.t) over a prediction period .DELTA.t a function F(t k, t k-.DELTA.n, t k-2.DELTA.n, ..., t k-(m-2).DELTA.n, t k-(m-1).DELTA.n) is determined, which depends on analyte values y(t k) , y(t k-.DELTA.n), y(t k-2.DELTA.n), ..., y(t k-(m-2).DELTA.n, y(t k-(m-1).DELTA.n) and which can be used to approximate the progression of the analyte values y(t n) at time t n0 in a vicinity U of an analyte value y(t n0) with a pre-determined accuracy .sigma., such that whereby .DELTA.n, n0 and m are integers, and calculating a prediction for the analyte value y(t n0+.DELTA.t) by using the following equation y(t n0 + .DELTA.t) = F(t n0, t n0 - .DELTA.n, ...
..., t(n0-(m-2).DELTA.n),t(n0-(m-1).DELTA.n)).
..., t(n0-(m-2).DELTA.n),t(n0-(m-1).DELTA.n)).
2. Method according to claim 1, wherein the function F
depends on analyte values y(t k), y(t k-.DELTA.n), y(t k-2.DELTA.n),...., y(t k-(m-2).DELTA.n), y(t k-(m-1).DELTA.n) and on their first time derivatives.
depends on analyte values y(t k), y(t k-.DELTA.n), y(t k-2.DELTA.n),...., y(t k-(m-2).DELTA.n), y(t k-(m-1).DELTA.n) and on their first time derivatives.
3. Method according to claim 1 or 2, wherein the function F depends on analyte values y(t k), y(t k-.DELTA.n), y(t k-2.DELTA.n),...., y(t k-(m-2).DELTA.n), y(t k-(m-1).DELTA.n) and on their first and second time derivatives.
4. Method according to any one of claims 1 to 3, wherein the function F contains linear or square terms.
5. Method according to any one of claims 1 to 4, wherein the function F is a linear function.
6. Method according to any one of claims 1 to 5, wherein the function F can be represented by coefficients a0 to a m as follows
7. Method according to claim 6, wherein the coefficients, a0 to a m, are determined by minimizing the sum
8. Method according to claim 6 or 7, wherein the coefficients, a0 to a m, are determined by solving a system of linear equations which contains one equation each of the type for at least m+1 different analyte values y(t k) from the vicinity U of the point in time t n0.
9. Method according to claim 8, wherein the system of equations contains more than m+1 equations and is solved numerically by approximation to determine the coefficients, a0 to a m.
10. Method according to any one of claims 1 to 9, wherein the analyte values y(t n) used to determine the function F are obtained by numerical processing, in particular by filtering, of measured values.
11. Method according to any one of claims 1 to 10, wherein in an additional step the function F(t k, t k-.DELTA.n, t k-2.DELTA.n,..., t k-(m-2).DELTA.n, t k-(m-1).DELTA.n) is adapted in a vicinity U of an analyte value y(t n0+.DELTA.t) at a point in time t n =t n0+.DELTA.t for a calculated prediction value y(t n0+.DELTA.t), such that the adapted function F
approximates the progression of analyte values y(t n) in the vicinity U with a predetermined accuracy and an additional prediction value y(t n0+2.DELTA.t) is calculated therefrom.
approximates the progression of analyte values y(t n) in the vicinity U with a predetermined accuracy and an additional prediction value y(t n0+2.DELTA.t) is calculated therefrom.
12. Method according to any one of claims 1 to 11, wherein the analyte values y(t n) used to calculate the prediction value are determined for consecutive points in time t n which are separated by intervals of 30 seconds to 5 minutes, preferably by 1 to 3 minutes.
13. Method according to any one of claim 1 to 12, wherein the analyte values y(t n) used to calculate the prediction value or one of the prediction values are determined for times t n which are separated by an interval corresponding to the prediction period .DELTA.t.
14. Method according to any one of claims 1 to 13, wherein the number of the coefficients, a0 to a m, is 4 to 11, preferably 6 to 9.
15. Method according to any one of claims 1 to 14, wherein the function F is determined as a function of values Ty (t k), T y (t k-.DELTA.n), Ty(t k-2.DELTA.n),...., Ty(t k-(m-2).DELTA.n), Ty(t k-(m-1).DELTA.n), which were determined by means of a transformation, preferably a linear transformation, from analyte values y(t k), y(t k-.DELTA.n), y(t k-2.DELTA.n),...., y(t k-(m-2).DELTA.n), y(t k-(m-1).DELTA.n), whereby the transformation is selected such that at least one of the values Ty(t k), T t k-.DELTA.n), Ty(t k-2.DELTA.n),...., Ty(t k-(m-2).DELTA.n), Ty(t k-(m-1).DELTA.n) has a negligible influence on the function F at the predetermined accuracy .sigma..
16. Device for continuous monitoring of an analyte concentration by determining its progression in the living body of a human or animal, comprising a measuring unit (12), by means of which values of a measuring parameter that is correlated with the concentration to be determined are measured as measuring signals at consecutive points in time t n, and an analytical unit (3), by means of which analyte values y(t n) that are correlated with the concentration are determined from the measured values, the analytical unit (3) being adapted for predicting an analyte value y(t n0+.DELTA.t) over a prediction period .DELTA.t by determining a function F(t k, t k-.DELTA.n, t k-2.DELTA.n, ..., t k-(m-2).DELTA.n, t k-(m-1).DELTA.n), which depends on analyte values y(t k), y(t k-.DELTA.n), y(t k-2.DELTA.n), ..., y(t k-(m-2).DELTA.n, y(t k-(m-1).DELTA.n) and can be used to approximate the progression of the analyte values y(t n) at time t n0 in a vicinity U of an analyte value y(t n0) with a pre-determined accuracy .sigma., such that whereby .DELTA.n, n0 and m are integers, and subsequently the function F is used by the analytical unit (3) to calculate a prediction for the analyte value y(t n0+.DELTA.t) by using the following equation y(t n0 + .DELTA.t) = F(t n0, t n0 - .DELTA.n, ...
..., t(n0-(m-2).DELTA.n),t(n0-(m-1).DELTA.n)).
..., t(n0-(m-2).DELTA.n),t(n0-(m-1).DELTA.n)).
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102004020160A DE102004020160A1 (en) | 2004-04-24 | 2004-04-24 | Method and device for monitoring a concentration of an analyte in the living body of a human or animal |
DE102004020160.9 | 2004-04-24 |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2504661A1 CA2504661A1 (en) | 2005-10-24 |
CA2504661C true CA2504661C (en) | 2009-08-18 |
Family
ID=34934888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002504661A Expired - Fee Related CA2504661C (en) | 2004-04-24 | 2005-04-21 | Method and device for monitoring an analyte concentration in the living body of a human or animal |
Country Status (5)
Country | Link |
---|---|
US (2) | US7188034B2 (en) |
EP (1) | EP1589438B1 (en) |
JP (1) | JP4693474B2 (en) |
CA (1) | CA2504661C (en) |
DE (1) | DE102004020160A1 (en) |
Families Citing this family (63)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6391005B1 (en) | 1998-03-30 | 2002-05-21 | Agilent Technologies, Inc. | Apparatus and method for penetration with shaft having a sensor for sensing penetration depth |
US8641644B2 (en) | 2000-11-21 | 2014-02-04 | Sanofi-Aventis Deutschland Gmbh | Blood testing apparatus having a rotatable cartridge with multiple lancing elements and testing means |
US9795747B2 (en) | 2010-06-02 | 2017-10-24 | Sanofi-Aventis Deutschland Gmbh | Methods and apparatus for lancet actuation |
US7749174B2 (en) | 2001-06-12 | 2010-07-06 | Pelikan Technologies, Inc. | Method and apparatus for lancet launching device intergrated onto a blood-sampling cartridge |
US9226699B2 (en) | 2002-04-19 | 2016-01-05 | Sanofi-Aventis Deutschland Gmbh | Body fluid sampling module with a continuous compression tissue interface surface |
US7981056B2 (en) | 2002-04-19 | 2011-07-19 | Pelikan Technologies, Inc. | Methods and apparatus for lancet actuation |
US7025774B2 (en) | 2001-06-12 | 2006-04-11 | Pelikan Technologies, Inc. | Tissue penetration device |
US8337419B2 (en) | 2002-04-19 | 2012-12-25 | Sanofi-Aventis Deutschland Gmbh | Tissue penetration device |
ATE485766T1 (en) | 2001-06-12 | 2010-11-15 | Pelikan Technologies Inc | ELECTRICAL ACTUATING ELEMENT FOR A LANCET |
DE60234598D1 (en) | 2001-06-12 | 2010-01-14 | Pelikan Technologies Inc | SELF-OPTIMIZING LANZET DEVICE WITH ADAPTANT FOR TEMPORAL FLUCTUATIONS OF SKIN PROPERTIES |
US9427532B2 (en) | 2001-06-12 | 2016-08-30 | Sanofi-Aventis Deutschland Gmbh | Tissue penetration device |
US9314194B2 (en) | 2002-04-19 | 2016-04-19 | Sanofi-Aventis Deutschland Gmbh | Tissue penetration device |
US9795334B2 (en) | 2002-04-19 | 2017-10-24 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus for penetrating tissue |
US7229458B2 (en) | 2002-04-19 | 2007-06-12 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US7175642B2 (en) | 2002-04-19 | 2007-02-13 | Pelikan Technologies, Inc. | Methods and apparatus for lancet actuation |
US8267870B2 (en) | 2002-04-19 | 2012-09-18 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus for body fluid sampling with hybrid actuation |
US7976476B2 (en) | 2002-04-19 | 2011-07-12 | Pelikan Technologies, Inc. | Device and method for variable speed lancet |
US7297122B2 (en) | 2002-04-19 | 2007-11-20 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US9248267B2 (en) | 2002-04-19 | 2016-02-02 | Sanofi-Aventis Deustchland Gmbh | Tissue penetration device |
US7232451B2 (en) | 2002-04-19 | 2007-06-19 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US7331931B2 (en) | 2002-04-19 | 2008-02-19 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US7909778B2 (en) | 2002-04-19 | 2011-03-22 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US7491178B2 (en) | 2002-04-19 | 2009-02-17 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US7674232B2 (en) | 2002-04-19 | 2010-03-09 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US7226461B2 (en) | 2002-04-19 | 2007-06-05 | Pelikan Technologies, Inc. | Method and apparatus for a multi-use body fluid sampling device with sterility barrier release |
US8784335B2 (en) | 2002-04-19 | 2014-07-22 | Sanofi-Aventis Deutschland Gmbh | Body fluid sampling device with a capacitive sensor |
US7892183B2 (en) | 2002-04-19 | 2011-02-22 | Pelikan Technologies, Inc. | Method and apparatus for body fluid sampling and analyte sensing |
US8702624B2 (en) | 2006-09-29 | 2014-04-22 | Sanofi-Aventis Deutschland Gmbh | Analyte measurement device with a single shot actuator |
US8579831B2 (en) | 2002-04-19 | 2013-11-12 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus for penetrating tissue |
US8360992B2 (en) | 2002-04-19 | 2013-01-29 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus for penetrating tissue |
US7901362B2 (en) | 2002-04-19 | 2011-03-08 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US8221334B2 (en) | 2002-04-19 | 2012-07-17 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus for penetrating tissue |
US7892185B2 (en) | 2002-04-19 | 2011-02-22 | Pelikan Technologies, Inc. | Method and apparatus for body fluid sampling and analyte sensing |
US7547287B2 (en) | 2002-04-19 | 2009-06-16 | Pelikan Technologies, Inc. | Method and apparatus for penetrating tissue |
US8574895B2 (en) | 2002-12-30 | 2013-11-05 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus using optical techniques to measure analyte levels |
EP1628567B1 (en) | 2003-05-30 | 2010-08-04 | Pelikan Technologies Inc. | Method and apparatus for fluid injection |
DK1633235T3 (en) | 2003-06-06 | 2014-08-18 | Sanofi Aventis Deutschland | Apparatus for sampling body fluid and detecting analyte |
WO2006001797A1 (en) | 2004-06-14 | 2006-01-05 | Pelikan Technologies, Inc. | Low pain penetrating |
US8282576B2 (en) | 2003-09-29 | 2012-10-09 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus for an improved sample capture device |
EP1680014A4 (en) | 2003-10-14 | 2009-01-21 | Pelikan Technologies Inc | Method and apparatus for a variable user interface |
US7822454B1 (en) | 2005-01-03 | 2010-10-26 | Pelikan Technologies, Inc. | Fluid sampling device with improved analyte detecting member configuration |
EP1706026B1 (en) | 2003-12-31 | 2017-03-01 | Sanofi-Aventis Deutschland GmbH | Method and apparatus for improving fluidic flow and sample capture |
US8828203B2 (en) | 2004-05-20 | 2014-09-09 | Sanofi-Aventis Deutschland Gmbh | Printable hydrogels for biosensors |
US9775553B2 (en) | 2004-06-03 | 2017-10-03 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus for a fluid sampling device |
EP1765194A4 (en) | 2004-06-03 | 2010-09-29 | Pelikan Technologies Inc | Method and apparatus for a fluid sampling device |
US8652831B2 (en) | 2004-12-30 | 2014-02-18 | Sanofi-Aventis Deutschland Gmbh | Method and apparatus for analyte measurement test time |
EP1728468A1 (en) * | 2005-06-04 | 2006-12-06 | Roche Diagnostics GmbH | Evaluation of blood glucose concentration values for adaptation of insulin dosage |
ATE465674T1 (en) * | 2006-08-08 | 2010-05-15 | Koninkl Philips Electronics Nv | METHOD AND DEVICE FOR MONITORING A PHYSIOLOGICAL PARAMETER |
US8439837B2 (en) * | 2006-10-31 | 2013-05-14 | Lifescan, Inc. | Systems and methods for detecting hypoglycemic events having a reduced incidence of false alarms |
ES2683852T3 (en) * | 2007-01-23 | 2018-09-28 | Ascensia Diabetes Care Holdings Ag | Analyte test device |
EP1967132B8 (en) | 2007-03-09 | 2015-06-03 | Roche Diagnostics GmbH | Medical data transmission system |
US20090160656A1 (en) * | 2007-10-11 | 2009-06-25 | Mahesh Seetharaman | Analyte monitoring system alarms |
DE102007046864A1 (en) * | 2007-09-28 | 2009-04-09 | Endress + Hauser Conducta Gesellschaft für Mess- und Regeltechnik mbH + Co. KG | Measured value plausibility and quality evaluating method for e.g. control circuit, involves forming and evaluating derivative of measured value according to time, where measured value is compared with predefined warning limit value |
US7731659B2 (en) * | 2007-10-18 | 2010-06-08 | Lifescan Scotland Limited | Method for predicting a user's future glycemic state |
US7695434B2 (en) * | 2007-10-19 | 2010-04-13 | Lifescan Scotland, Ltd. | Medical device for predicting a user's future glycemic state |
WO2009126900A1 (en) | 2008-04-11 | 2009-10-15 | Pelikan Technologies, Inc. | Method and apparatus for analyte detecting device |
US20110160555A1 (en) * | 2008-07-31 | 2011-06-30 | Jacques Reifman | Universal Models for Predicting Glucose Concentration in Humans |
US9375169B2 (en) | 2009-01-30 | 2016-06-28 | Sanofi-Aventis Deutschland Gmbh | Cam drive for managing disposable penetrating member actions with a single motor and motor and control system |
CN102458234B (en) | 2009-04-17 | 2014-02-26 | 爱科来株式会社 | Data provision system, data provision method, program, and server device |
US8843321B2 (en) * | 2010-01-26 | 2014-09-23 | Roche Diagnostics Operations, Inc. | Methods and systems for processing glucose data measured from a person having diabetes |
US8965476B2 (en) | 2010-04-16 | 2015-02-24 | Sanofi-Aventis Deutschland Gmbh | Tissue penetration device |
EP2700031A2 (en) * | 2011-04-20 | 2014-02-26 | Novo Nordisk A/S | Glucose predictor based on regularization networks with adaptively chosen kernels and regularization parameters |
KR101962315B1 (en) * | 2017-03-29 | 2019-07-17 | 이화여자대학교 산학협력단 | Method of prediction mercury concentration in blood |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FI931234A (en) * | 1993-03-19 | 1994-09-20 | Instrumentarium Oy | Procedure for monitoring a patient's condition |
KR0142483B1 (en) * | 1994-02-28 | 1998-08-17 | 고지마 게이지 | Non-linear time sequential data predicting device |
DE4415896A1 (en) | 1994-05-05 | 1995-11-09 | Boehringer Mannheim Gmbh | Analysis system for monitoring the concentration of an analyte in the blood of a patient |
EP0910023A2 (en) * | 1997-10-17 | 1999-04-21 | Siemens Aktiengesellschaft | Method and device for the neuronal modelling of a dynamic system with non-linear stochastic behavior |
JPH11296598A (en) * | 1998-04-07 | 1999-10-29 | Seizaburo Arita | System and method for predicting blood-sugar level and record medium where same method is recorded |
ATE246356T1 (en) * | 1998-05-13 | 2003-08-15 | Cygnus Therapeutic Systems | DEVICE FOR PREDICTING PHYSIOLOGICAL MEASUREMENTS |
US6925393B1 (en) | 1999-11-18 | 2005-08-02 | Roche Diagnostics Gmbh | System for the extrapolation of glucose concentration |
JP3563333B2 (en) | 2000-08-04 | 2004-09-08 | Necエレクトロニクス株式会社 | High-side switch for inrush current and overcurrent control method |
WO2002015777A1 (en) | 2000-08-18 | 2002-02-28 | Cygnus, Inc. | Methods and devices for prediction of hypoglycemic events |
DE60334365D1 (en) | 2002-03-22 | 2010-11-11 | Animas Technologies Llc | INCREASED PERFORMANCE OF AN ANALYSIS MONITORING DEVICE |
US7278983B2 (en) | 2002-07-24 | 2007-10-09 | Medtronic Minimed, Inc. | Physiological monitoring device for controlling a medication infusion device |
WO2004015539A2 (en) | 2002-08-13 | 2004-02-19 | University Of Virginia Patent Foundation | Managing and processing self-monitoring blood glucose |
DE10343863A1 (en) * | 2003-09-23 | 2005-04-14 | Roche Diagnostics Gmbh | Method and device for continuously monitoring the concentration of an analyte |
US20070060869A1 (en) | 2005-08-16 | 2007-03-15 | Tolle Mike C V | Controller device for an infusion pump |
US20070060870A1 (en) | 2005-08-16 | 2007-03-15 | Tolle Mike Charles V | Controller device for an infusion pump |
-
2004
- 2004-04-24 DE DE102004020160A patent/DE102004020160A1/en not_active Withdrawn
-
2005
- 2005-04-08 EP EP05007701.5A patent/EP1589438B1/en active Active
- 2005-04-14 JP JP2005117542A patent/JP4693474B2/en not_active Expired - Fee Related
- 2005-04-21 CA CA002504661A patent/CA2504661C/en not_active Expired - Fee Related
- 2005-04-25 US US11/113,606 patent/US7188034B2/en active Active
-
2007
- 2007-02-15 US US11/675,580 patent/US7650244B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
EP1589438A3 (en) | 2006-03-22 |
JP4693474B2 (en) | 2011-06-01 |
JP2005308742A (en) | 2005-11-04 |
US20050240356A1 (en) | 2005-10-27 |
CA2504661A1 (en) | 2005-10-24 |
US7650244B2 (en) | 2010-01-19 |
EP1589438B1 (en) | 2015-11-04 |
EP1589438A2 (en) | 2005-10-26 |
DE102004020160A1 (en) | 2005-11-10 |
US7188034B2 (en) | 2007-03-06 |
US20070148778A1 (en) | 2007-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2504661C (en) | Method and device for monitoring an analyte concentration in the living body of a human or animal | |
Joshi et al. | iGLU 2.0: A new wearable for accurate non-invasive continuous serum glucose measurement in IoMT framework | |
US6968221B2 (en) | Low-cost method and apparatus for non-invasively measuring blood glucose levels | |
JP3242346U (en) | Devices and systems for non-invasive monitoring of physiological measurements | |
US10825569B2 (en) | Universal non-invasive blood glucose estimation method based on time series analysis | |
CA2481627C (en) | Method and device for continuous monitoring of the concentration of an analyte | |
EP1658578B1 (en) | A system and method for detecting signal artifacts | |
EP2408354B1 (en) | Systems and methods for monitoring heart rate and blood pressure correlation | |
US20140309507A1 (en) | Reducing nuisance alarms | |
US8571623B2 (en) | System and method for detection of venous pulsation | |
US6064896A (en) | Non-invasive measurement of blood glucose using instruments that have less precise detection capability | |
KR20060061838A (en) | A method of monitoring glucose level | |
WO2005096922A1 (en) | Physiological assessment system | |
US20110021929A1 (en) | Systems and methods for continuous non-invasive blood pressure monitoring | |
Champseix et al. | A python package for heart rate variability analysis and signal preprocessing | |
Vettoretti et al. | Modeling the error of factory-calibrated continuous glucose monitoring sensors: application to Dexcom G6 sensor data | |
CA2597004A1 (en) | Method and apparatus for monitoring a sedated patient | |
DK1793321T3 (en) | Evaluation method and analysis system of an analyte in the body fluid of a human or animal | |
KR20110094405A (en) | Method for non-invasive glucose measurement and non-invasive glucose measuring apparatus using the same method | |
RU2230485C2 (en) | Method for determination of blood glucose concentration in humans | |
US20190104975A1 (en) | Device And Method For Measurement of Glycated Hemoglobin (A1c) | |
Kossowski et al. | Robust IR attenuation measurement for non-invasive glucose level analysis | |
US20220142520A1 (en) | A non-invasive glucometer | |
CN116327186A (en) | Noninvasive blood glucose detection system, noninvasive blood glucose detection method, noninvasive blood glucose detection device and noninvasive blood glucose detection storage medium | |
CN116965807A (en) | System for realizing dynamic continuous monitoring of blood sugar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
MKLA | Lapsed |
Effective date: 20170421 |
|
MKLA | Lapsed |
Effective date: 20170421 |